PET-MR Motion Correction Using a Surrogate-Driven Respiratory Motion Model
(Funded by EU HYPERIMAGE project)
PET imaging is commonly used for detecting cancer and monitoring its progress, but its quantitative accuracy is limited by the problem of motion, such as that caused by respiration. The state-of-the-art in motion correction of PET data is limited to coarse level modelling of the motion and often requires significant extra scanning resources to acquire the necessary imaging data. As a result the translation of motion correction techniques into clinical use has been limited. The recent introduction to the market of simultaneous PET-MR scanners raises the possibility of motion correcting PET data based on motion estimates derived from MR imaging.
Initial work focused on the use of a surrogate-driven motion model. The model was built from dynamic MR images, and captured the relationship between a simple 1-D MR navigator signal and the complex 3-D respiratory motion of the thorax. Once the model was formed, motion estimates could be made by acquiring only the MR navigator and using it as input to the motion model. The resulting motion fields were used to motion correct reconstructed PET gates in simulations.
- A. P. King, C. Tsoumpas, C. Buerger, V. Schulz, P. Marsden, T. Schaeffter, "Real-time respiratory motion correction for simultaneous PET-MR using an MR-derived motion model", Proceedings IEEE NSS-MIC, 2011. (paper)
This work was then developed to use a more sophisticated motion model and imaging data as the surrogate. Specifically, dimensionality reduction using PCA was applied to a set of motion fields derived from dynamic MR imaging. The resulting eigenvectors represented the respiratory motion model. To apply the model, a fast 2D MR slice was acquired, and the modes of variation of the model were optimised to maximise the match between the acquired slice (known as the surrogate) and an exhale MR volume transformed using the model transformation. Thus, a 2D slice was used to estimate a 3D transformation, constrained by the motions seen during the model formation phase. This technique was also evaluated on simulated PET gates.
- A. P. King, C. Buerger, C. Tsoumpas, P. Marsden, T. Schaeffter, "Thoracic Respiratory Motion Estimation from MRI Using a Statistical Model and a 2-D Image Navigator", Medical Image Analysis 16:252-264, 2012. (paper)
PET-MR Motion Correction Constrained by a Respiratory Motion Model
Subsequent work on the use of motion models for PET-MR motion correction demonstrated that it was not actually necessary to acquire the surrogate data at all in order to apply the model. Given a simple surrogate-driven model such as those used in the work described above, reconstructed PET gates can be transformed to a reference motion state using the optimal motion field produced by the model, i.e. the one that maximises the similarity between the transformed PET gate and the reference PET gate. Effectively, this technique used a PET-PET registration that was constrained by the MR-derived respiratory motion model.
- D. R. Balfour, P. K. Marsden, I. Polycarpou, C. Kolbitsch, A. P. King, "Respiratory motion correction of PET using MR-constrained PET-PET registration", Biomedical Engineering Online 14:85, 2015. (open access paper)
Motion Correction for Cardiac Catheterisations
(Funded by EPSRC grants EP/D061471/1 and EP/H046410/1)
Guidance of minimally invasive cardiac catheterisations can be enhanced by using 'roadmaps' derived from preprocedure imaging, such as MR or CT. However, motion, such as that caused by respiration, can cause a misalignment between the static roadmap and the moving anatomy that it represents. This work aimed to estimate respiratory motion transformations for updating the roadmaps in such procedures. An affine surrogate-driven respiratory motion model was formed from dynamic MR images, using an MR pencil-beam navigator positioned on the diaphraghm as the surrogate. During the procedure, the diaphragm was tracked from X-ray images, and this trace was used as input to the motion model. The model then produced an affine transform that was used to update the roadmap. The technique was evaluated on a number of clinical cases, with accuracy assessed by overlaying a roadmap representing a structure with an X-ray opaque catheter positioned inside it. The alignment between roadmap and anatomy was assessed by comparing the catheter position with the roadmap (see movies below).
- A. P. King, R. Boubertakh, K. S. Rhode, Y. L. Ma, P. Chinchapatnam, G. Gao, T. Tangcharoen, M. Ginks, M. Cooklin, J. S. Gill, D. J. Hawkes, R. S. Razavi, T. Schaeffter. "A Subject-Specific Technique for Respiratory Motion Correction in Image-Guided Cardiac Catheterisation Procedures", Medical Image Analysis, 13(3): 419-431, 2009. (paper)
- A. P. King, R. Boubertakh, K. L. Ng, Y. L. Ma, P. Chinchapatnam, G. Gang, T. Schaeffter, D. J. Hawkes, R. Razavi, K. S. Rhode, "A technique for respiratory motion correction in image-guided cardiac catheterization procedures", Proceedings SPIE Medical Imaging, 2008. (paper)
In subsequent work, the motion model was refined to enable it to adapt to changes in breathing pattern, that can occur during some procedures. For example, patients may breathe deeply due to discomfort as the heart is ablated, or breathe more shallowly when they are relaxed due to the effects of sedation. Such breathing pattern changes affect the accuracy of the model, and the new adaptive and predictive model was able to adapt to these changes.
- A. P. King, K. S. Rhode, R. Razavi, T. Schaeffter, "An Adaptive and Predictive Respiratory Motion Model for Image-Guided Interventions: Theory and First Clinical Application", IEEE Transactions on Medical Imaging, 28(12):2020-2032, 2009. (paper)
- A. P. King, K. S. Rhode, R. Razavi, T. Schaeffter, "An Adaptive and Predictive Respiratory Motion Model for Image-Guided Interventions", Proceedings MICCAI Workshop on Image-Guidance and Computer Assistance for Soft-Tissue Interventions, 2008.
Later, it was shown how intraprocedure ultrasound data could be used as input to the motion model. This innovation enabled the motion information contained in the noisy, low quality but real-time ultrasound data to be exploited, whilst being constrained by a robust preprocedure motion model.
- A. P. King, C. Jansen, K. S. Rhode, D. Caulfield, R. S. Razavi, G. P. Penney, "Respiratory Motion Correction for Image-Guided Cardiac Interventions Using 3-D Echocardiography", Medical Image Analysis, 14(1):21-29, 2010. (paper)
- A. P. King, C. Jansen, R. Boubertakh, K. S. Rhode, R. Razavi, G. P. Penney, "Model-based respiratory motion correction using 3-D echocardiography", Proceedings ISBI, pp1465-1468, 2008. (paper)
The ultrasound based approach was then extended to make use of a more flexible model, which permitted some deviation from the motions observed during the model formation. This was made possible using a Bayesian approach, in which the motion model defined a prior distribution of motion states, and the intraprocedure ultrasound data defined the likelihood term. By maximising the posterior probability, the most probably motion state was estimated, resulting in improvements in motion estimation accuracy.
- D. Peressutti, G. P. Penney, R. J. Housden, C. Kolbitsch, A. Gomez, E.-J. Rijkhorst, D. C. Barratt, K. S. Rhode, A. P. King, "A Novel Bayesian Respiratory Motion Model to Estimate and Resolve Uncertainty in Image-Guided Cardiac Interventions", Medical Image Analysis 17(4):488-502, 2013. (paper)
- D. Peressutti, E.-J. Rijkhoorst, D. C. Barratt, G. P. Penney, A. P. King, "Estimating and resolving uncertainty in cardiac respiratory motion modelling", Proceedings ISBI, pp262-265, 2012. (paper)
Population-based motion models are attractive as they obviate the need for a subject-specific scan for model formation. However, to date they have been less accurate and robust than subject-specific models. We aimed to address this weakness by developing an improved personalisation scheme for population-based respiratory motion models. The approach used neighbourhood approximation and learnt anatomical features to enable an accurate motion model to be formed for a new subject using a static anatomical image only. We demonstrated that accuracy was improved compared to the state-of-the-art in population-based models.
- D. Peressutti, G. P. Penney, C. Kolbitsch, A. P. King, "Personalising Population-Based Respiratory Motion Models Of The Heart Using Neighbourhood Approximation Based On Learnt Anatomical Features", Medical Image Analysis 18(7):1015-1025, 2014. (paper)
- D. Peressutti, G. P. Penney, C. Kolbitsch, A. P. King, "Personalising Cross-Population Respiratory Motion Models Using Anatomical Features", Proceedings MIUA, 2013. (proceedings) - AWARDED STUDENT PAPER PRIZE
Autoadaptive Motion Modelling
(Funded by EPSRC grant EP/H046410/1)
In an extension of the simultaneous groupwise manifold alignment (SGA) scheme, we demonstrated how using SGA with 2D motion fields, rather than 2D images, resulted in a highly novel form of surrogate-driven motion model. In this autoadaptive motion model, the data used to form the model and the surrogate data used to apply the model were both of the same type (i.e. 2D motion fields). This meant that the surrogate data could be automatically added to the model each time the model was applied, enabling it to automatically adapt to changes in breathing pattern without any need for interrupting the application of the model.
- C. F. Baumgartner, C. Kolbitsch, J. R. McClelland, D. Rueckert, A. P. King, "Autoadaptive Motion Modelling for MR-Based Respiratory Motion Estimation", Medical Image Analysis 35:83-110. (open access paper)
- C. F. Baumgartner, C. Kolbitsch, J. R. McClelland, D. Rueckert, A. P. King, "Autoadaptive Motion Modelling", Proceedings ISBI, pp457-460, 2014. (paper)