Statistical Modeling of 4D Respiratory Lung Motion Using Diffeomorphic Image Registration

Modeling of respiratory motion has become increasingly important in various applications of medical imaging (e.g., radiation therapy of lung cancer). Current modeling approaches are usually confined to intra-patient registration of 3D image data representing the individual patient's anatomy at different breathing phases. We propose an approach to generate a mean motion model of the lung based on thoracic 4D computed tomography (CT) data of different patients to extend the motion modeling capabilities. Our modeling process consists of three steps: an intra-subject registration to generate subject-specific motion models, the generation of an average shape and intensity atlas of the lung as anatomical reference frame, and the registration of the subject-specific motion models to the atlas in order to build a statistical 4D mean motion model (4D-MMM). Furthermore, we present methods to adapt the 4D mean motion model to a patient-specific lung geometry. In all steps, a symmetric diffeomorphic nonlinear intensity-based registration method was employed. The Log-Euclidean framework was used to compute statistics on the diffeomorphic transformations. The presented methods are then used to build a mean motion model of respiratory lung motion using thoracic 4D CT data sets of 17 patients. We evaluate the model by applying it for estimating respiratory motion of ten lung cancer patients. The prediction is evaluated with respect to landmark and tumor motion, and the quantitative analysis results in a mean target registration error (TRE) of 3.3 ±1.6 mm if lung dynamics are not impaired by large lung tumors or other lung disorders (e.g., emphysema). With regard to lung tumor motion, we show that prediction accuracy is independent of tumor size and tumor motion amplitude in the considered data set. However, tumors adhering to non-lung structures degrade local lung dynamics significantly and the model-based prediction accuracy is lower in these cases. The statistical respiratory motion model is capable of providing valuable prior knowledge in many fields of applications. We present two examples of possible applications in radiation therapy and image guided diagnosis.

[1]  Hans-Ulrich Kauczor,et al.  Measurement of tumor diameter-dependent mobility of lung tumors by dynamic MRI. , 2004, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[2]  Alain Trouvé,et al.  Computing Large Deformation Metric Mappings via Geodesic Flows of Diffeomorphisms , 2005, International Journal of Computer Vision.

[3]  J. Dempsey,et al.  Novel breathing motion model for radiotherapy. , 2005, International journal of radiation oncology, biology, physics.

[4]  V. Boldea,et al.  Simulation of four-dimensional CT images from deformable registration between inhale and exhale breath-hold CT scans. , 2006, Medical physics.

[5]  Guido Gerig,et al.  Unbiased diffeomorphic atlas construction for computational anatomy , 2004, NeuroImage.

[6]  Alejandro F Frangi,et al.  Automatic construction of 3-D statistical deformation models of the brain using nonrigid registration , 2003, IEEE Transactions on Medical Imaging.

[7]  Heinz Handels,et al.  A statistical shape and motion model for the prediction of respiratory lung motion , 2010, Medical Imaging.

[8]  J. Ehrhardt,et al.  An optical flow based method for improved reconstruction of 4D CT data sets acquired during free breathing. , 2007, Medical physics.

[9]  Nicholas Ayache,et al.  A Log-Euclidean Framework for Statistics on Diffeomorphisms , 2006, MICCAI.

[10]  Steve B Jiang,et al.  Synchronized moving aperture radiation therapy (SMART): average tumour trajectory for lung patients. , 2003, Physics in medicine and biology.

[11]  Paul J. Besl,et al.  A Method for Registration of 3-D Shapes , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[12]  R. Castillo,et al.  A framework for evaluation of deformable image registration spatial accuracy using large landmark point sets , 2009, Physics in medicine and biology.

[13]  John Ashburner,et al.  A fast diffeomorphic image registration algorithm , 2007, NeuroImage.

[14]  C. Twining,et al.  Constructing an atlas for the diffeomorphism group of a compact manifold with boundary, with application to the analysis of image registrations , 2008 .

[15]  Nicholas Ayache,et al.  Geometric Means in a Novel Vector Space Structure on Symmetric Positive-Definite Matrices , 2007, SIAM J. Matrix Anal. Appl..

[16]  Steve Webb,et al.  A margin model to account for respiration-induced tumour motion and its variability , 2008, Physics in medicine and biology.

[17]  Eric A. Hoffman,et al.  Registration-based estimates of local lung tissue expansion compared to xenon CT measures of specific ventilation , 2008, Medical Image Anal..

[18]  Jean Meunier,et al.  Average Brain Models: A Convergence Study , 2000, Comput. Vis. Image Underst..

[19]  J. McClelland,et al.  A continuous 4D motion model from multiple respiratory cycles for use in lung radiotherapy. , 2006, Medical physics.

[20]  Monica Hernandez,et al.  Registration of Anatomical Images Using Paths of Diffeomorphisms Parameterized with Stationary Vector Field Flows , 2009, International Journal of Computer Vision.

[21]  Heinz Handels,et al.  Patient-specific finite element modeling of respiratory lung motion using 4D CT image data. , 2009, Medical physics.

[22]  J. Gee,et al.  Geodesic estimation for large deformation anatomical shape averaging and interpolation , 2004, NeuroImage.

[23]  E. Hoffman,et al.  Mass preserving nonrigid registration of CT lung images using cubic B-spline. , 2009, Medical physics.

[24]  Tom Vercauteren,et al.  Diffeomorphic demons: Efficient non-parametric image registration , 2009, NeuroImage.

[25]  Max A. Viergever,et al.  Efficient and reliable schemes for nonlinear diffusion filtering , 1998, IEEE Trans. Image Process..

[26]  Hervé Delingette,et al.  Registration of 4D Time-Series of Cardiac Images with Multichannel Diffeomorphic Demons , 2008, MICCAI.

[27]  G. Christensen,et al.  A method for the reconstruction of four-dimensional synchronized CT scans acquired during free breathing. , 2003, Medical physics.

[28]  Salvador Olmos,et al.  A new algorithm for the computation of the group logarithm of diffeomorphisms , 2008 .

[29]  L. Younes,et al.  Statistics on diffeomorphisms via tangent space representations , 2004, NeuroImage.

[30]  Jean-Philippe Thirion,et al.  Image matching as a diffusion process: an analogy with Maxwell's demons , 1998, Medical Image Anal..

[31]  Eric A. Hoffman,et al.  Automatic lung segmentation for accurate quantitation of volumetric X-ray CT images , 2001, IEEE Transactions on Medical Imaging.

[32]  Nicholas Ayache,et al.  Symmetric Log-Domain Diffeomorphic Registration: A Demons-Based Approach , 2008, MICCAI.

[33]  Jan Modersitzki,et al.  Numerical Methods for Image Registration , 2004 .

[34]  C. Ling,et al.  Respiration-correlated spiral CT: a method of measuring respiratory-induced anatomic motion for radiation treatment planning. , 2002, Medical physics.

[35]  Radhe Mohan,et al.  Four-dimensional radiotherapy planning for DMLC-based respiratory motion tracking. , 2005, Medical physics.

[36]  Xavier Pennec,et al.  Intrinsic Statistics on Riemannian Manifolds: Basic Tools for Geometric Measurements , 2006, Journal of Mathematical Imaging and Vision.

[37]  Tinsu Pan,et al.  Dynamic ventilation imaging from four-dimensional computed tomography , 2006, Physics in medicine and biology.

[38]  R Mohan,et al.  Predicting respiratory motion for four-dimensional radiotherapy. , 2004, Medical physics.

[39]  A. Pevsner,et al.  Quantitation of respiratory motion during 4D-PET/CT acquisition. , 2004, Medical physics.

[40]  Nicholas Ayache,et al.  Iconic feature based nonrigid registration: the PASHA algorithm , 2003, Comput. Vis. Image Underst..

[41]  Daniel Rueckert,et al.  Construction of a 4D Statistical Atlas of the Cardiac Anatomy and Its Use in Classification , 2005, MICCAI.

[42]  H Handels,et al.  Motion Artifact Reducing Reconstruction of 4D CT Image Data for the Analysis of Respiratory Dynamics , 2007, Methods of Information in Medicine.

[43]  Vincent Arsigny,et al.  Processing Data in Lie Groups : An Algebraic Approach. Application to Non-Linear Registration and Diffusion Tensor MRI. (Traitement de données dans les groupes de Lie : une approche algébrique. Application au recalage non-linéaire et à l'imagerie du tenseur de diffusion) , 2006 .

[44]  Michael I. Miller,et al.  Group Actions, Homeomorphisms, and Matching: A General Framework , 2004, International Journal of Computer Vision.

[45]  Sasa Mutic,et al.  Quantitation of the reconstruction quality of a four-dimensional computed tomography process for lung cancer patients. , 2005, Medical physics.

[46]  Alain Trouvé,et al.  Diffeomorphisms Groups and Pattern Matching in Image Analysis , 1998, International Journal of Computer Vision.

[47]  Steve B. Jiang,et al.  The management of respiratory motion in radiation oncology report of AAPM Task Group 76. , 2006, Medical physics.

[48]  Daniel Rueckert,et al.  Spatial transformation of motion and deformation fields using nonrigid registration , 2004, IEEE Transactions on Medical Imaging.

[49]  J. McClelland,et al.  MRI-based measurements of respiratory motion variability and assessment of imaging strategies for radiotherapy planning , 2006, Physics in medicine and biology.

[50]  Brian B. Avants,et al.  A Dynamic Model of Average Lung Deformation Using Capacity-Based Reparameterization and Shape Averaging of Lung MR Images , 2004, MICCAI.

[51]  H Handels,et al.  Structure-preserving Interpolation of Temporal and Spatial Image Sequences Using an Optical Flow-based Method , 2007, Methods of Information in Medicine.

[52]  Daniel Rueckert,et al.  Construction of a Statistical Model for Cardiac Motion Analysis Using Nonrigid Image Registration , 2003, IPMI.

[53]  Quan Chen,et al.  Objective assessment of deformable image registration in radiotherapy: a multi-institution study. , 2007, Medical physics.