Abstract
Respiratory motion poses significant challenges in image-guided interventions. In emerging treatments such as MR-guided HIFU or MR-guided radiotherapy, it may cause significant misalignments between interventional road maps obtained pre-procedure and the anatomy during the treatment, and may affect intra-procedural imaging such as MR-thermometry. Patient specific respiratory motion models provide a solution to this problem. They establish a correspondence between the patient motion and simpler surrogate data which can be acquired easily during the treatment. Patient motion can then be estimated during the treatment by acquiring only the simpler surrogate data. In the majority of classical motion modelling approaches once the correspondence between the surrogate data and the patient motion is established it cannot be changed unless the model is recalibrated. However, breathing patterns are known to significantly change in the time frame of MR-guided interventions. Thus, the classical motion modelling approach may yield inaccurate motion estimations when the relation between the motion and the surrogate data changes over the duration of the treatment and frequent recalibration may not be feasible. We propose a novel methodology for motion modelling which has the ability to automatically adapt to new breathing patterns. This is achieved by choosing the surrogate data in such a way that it can be used to estimate the current motion in 3D as well as to update the motion model. In particular, in this work, we use 2D MR slices from different slice positions to build as well as to apply the motion model. We implemented such an autoadaptive motion model by extending our previous work on manifold alignment. We demonstrate a proof-of-principle of the proposed technique on cardiac gated data of the thorax and evaluate its adaptive behaviour on realistic synthetic data containing two breathing types generated from 6 volunteers, and real data from 4 volunteers. On synthetic data the autoadaptive motion model yielded 21.45% more accurate motion estimations compared to a non-adaptive motion model 10min after a change in breathing pattern. On real data we demonstrated the method's ability to maintain motion estimation accuracy despite a drift in the respiratory baseline. Due to the cardiac gating of the imaging data, the method is currently limited to one update per heart beat and the calibration requires approximately 12min of scanning. Furthermore, the method has a prediction latency of 800ms. These limitations may be overcome in future work by altering the acquisition protocol.
Highlights
Recent advances in magnetic resonance (MR) compatible materials and the development of fast parallel computational tech-C.F
Both the calibration as well as the surrogate data are 2D MR slices acquired from variable imaging planes
We extend our previously proposed simultaneous groupwise manifold alignment (SGA) technique to build an autoadaptive motion model from multiple 2D motion fields derived from sagittal 2D MR slices acquired at different anatomical positions
Summary
Recent advances in magnetic resonance (MR) compatible materials and the development of fast parallel computational tech-C.F. In MRg-RT, magnetic resonance imaging is used to accurately identify and track the target and to prevent the irradiation of healthy tissue in organs at risk (Crijns et al, 2012). For treatments targeting organs affected by breathing motion such as the lungs, the liver, the kidneys or the heart, accurate knowledge of the respiratory motion is essential. Apart from ensuring the irradiation or ablation of the intended target and sparing of the organs at risk, knowledge of respiratory motion is crucial to correct for motion-induced image-artefacts and for adjusting accumulated dose calculations such as temperature maps in MRgHIFU or dose simulations in MRg-RT (De Senneville et al, 2015). We will discuss how correspondences in the low-dimensional embedded space can be established for data obtained from two slice positions (Section 2.2). For an overview of the mathematical notation used in the remainder of this paper refer to Table 1
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.