Fully automatic nonrigid registration-based local motion estimation for motion-corrected iterative cardiac CT reconstruction.

PURPOSE Cardiac computed tomography is a rapidly emerging technique for noninvasive diagnosis of cardiovascular diseases. Nevertheless, the cardiac motion continues to be a limiting factor. Electrocardiogram-gated cardiac computed tomography reconstruction methods yield excellent results, but these are limited in their temporal resolution due to the mechanical movement of the gantry, and lead to residual motion blurring artifacts. If the motion of the cardiac region of interest is determined, motion compensated gated reconstructions can be applied to reduce motion artifacts. In this paper it is shown that elastic image registration methods can be an accurate solution to determine the cardiac motion. A method, which combines elastic registration and iterative computed tomography reconstruction methods delivering motion-corrected images of a chosen cardiac region of interest, is introduced. METHODS Using a gated four-dimensional region of interest image data set, a fully automatic elastic image registration is applied to recover a cardiac displacement field from a reference phase to a number of phases within the RR interval. Here, a stochastic optimizer and multiresolution approach are adopted to speed up the registration process. Subsequently, motion-compensated iterative reconstruction using the determined motion field is carried out. For the image representation volume-adapted spherical basis functions (blobs) are used to take the volume change caused by a divergent motion vector field into account. RESULTS The method is evaluated on phantom data and on four clinical data sets at a strong cardiac motion phase. Comparing the method to standard gated iterative reconstruction results shows that motion compensation strongly improves the image quality in these phases. A qualitative and quantitative accuracy study is presented for the estimated cardiac motion field. For the first time a blob-volume adaptation is applied on clinical data, and in the case of divergent motion it yields improved image quality. CONCLUSIONS A fully automatic local cardiac motion compensated gated iterative method with volume-adapted blobs is proposed. The method leads to excellent motion-corrected images which outperform nonmotion corrected results in phases of strong cardiac motion. In clinical cases, a volume-dependent blob-footprint adaptation proves to be a good solution to take care of the change in the blob volume caused by a divergent motion field.

[1]  Michael Unser,et al.  Fast parametric elastic image registration , 2003, IEEE Trans. Image Process..

[2]  Akram Aldroubi,et al.  B-SPLINE SIGNAL PROCESSING: PART I-THEORY , 1993 .

[3]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[4]  J. Spall Implementation of the simultaneous perturbation algorithm for stochastic optimization , 1998 .

[5]  Marc Kachelrieß,et al.  Advanced single-slice rebinning in cone-beam spiral CT: theoretical considerations and medical applications , 2000, Image Processing.

[6]  T Nielsen,et al.  Helical cardiac cone beam reconstruction using retrospective ECG gating. , 2003, Physics in medicine and biology.

[7]  M Grass,et al.  Motion-compensated iterative cone-beam CT image reconstruction with adapted blobs as basis functions , 2008, Physics in medicine and biology.

[8]  M Grass,et al.  Aperture weighted cardiac reconstruction for cone-beam CT , 2006, Physics in medicine and biology.

[9]  Eldon Hansen,et al.  Numerical Methods for Unconstrained Optimization and Nonlinear Equations (J. E. Dennis, Jr., and Robert B. Schnabel) , 1986 .

[10]  D. Hill,et al.  Non-rigid image registration: theory and practice. , 2004, The British journal of radiology.

[11]  D. Shepard A two-dimensional interpolation function for irregularly-spaced data , 1968, ACM National Conference.

[12]  T Nielsen,et al.  Adaptive temporal resolution optimization in helical cardiac cone beam CT reconstruction. , 2003, Medical physics.

[13]  Ge Wang,et al.  A knowledge-based cone-beam x-ray CT algorithm for dynamic volumetric cardiac imaging. , 2002, Medical physics.

[14]  J. Kiefer,et al.  Stochastic Estimation of the Maximum of a Regression Function , 1952 .

[15]  Josien P. W. Pluim,et al.  Evaluation of Optimization Methods for Nonrigid Medical Image Registration Using Mutual Information and B-Splines , 2007, IEEE Transactions on Image Processing.

[16]  R. Lewitt Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

[17]  Willi A Kalender,et al.  Extended parallel backprojection for standard three-dimensional and phase-correlated four-dimensional axial and spiral cone-beam CT with arbitrary pitch, arbitrary cone-angle, and 100% dose usage. , 2004, Medical physics.

[18]  Jeffrey A. Fessler,et al.  A Simple Regularizer for B-spline Nonrigid Image Registration That Encourages Local Invertibility , 2009, IEEE Journal of Selected Topics in Signal Processing.

[19]  Michael Grass,et al.  Automatic generation of time resolved motion vector fields of coronary arteries and 4D surface extraction using rotational x-ray angiography , 2009, Physics in medicine and biology.

[20]  Michael Grass,et al.  Improvement of cardiac CT reconstruction using local motion vector fields , 2009, Comput. Medical Imaging Graph..

[21]  Nicholas Ayache,et al.  3D tomographic reconstruction of coronary arteries using a precomputed 4D motion field. , 2004, Physics in medicine and biology.

[22]  Günter Lauritsch,et al.  Cardiac C-Arm CT: A Unified Framework for Motion Estimation and Dynamic CT , 2009, IEEE Transactions on Medical Imaging.

[23]  M. Bandstra,et al.  Measurements of Fukushima fallout by the Berkeley Radiological Air and Water Monitoring project , 2011, 2011 IEEE Nuclear Science Symposium Conference Record.

[24]  L Xing,et al.  Motion correction for improved target localization with on-board cone-beam computed tomography , 2006, Physics in medicine and biology.

[25]  Joachim Hornegger,et al.  Multi-modal 2D-3D non-rigid registration , 2006, SPIE Medical Imaging.

[26]  Weiguo Lu,et al.  A simple fixed-point approach to invert a deformation field. , 2007, Medical physics.

[27]  H. Robbins A Stochastic Approximation Method , 1951 .

[28]  J. Halton On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals , 1960 .

[29]  Michael Grass,et al.  Motion compensated iterative reconstruction of a region of interest in cardiac cone-beam CT , 2010, Comput. Medical Imaging Graph..

[30]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[31]  Akram Aldroubi,et al.  B-SPLINE SIGNAL PROCESSING: PART II-EFFICIENT DESIGN AND APPLICATIONS , 1993 .

[32]  Michael Grass,et al.  Motion-compensated and gated cone beam filtered back-projection for 3-D rotational X-ray angiography , 2006, IEEE Transactions on Medical Imaging.

[33]  Lucas C. Parra,et al.  List-Mode Likelihood EM Algorithm and Noise Estimation Demonstrated on 2D-PET , 1998, IEEE Trans. Medical Imaging.

[34]  Katsuyuki Taguchi,et al.  Direct cone-beam cardiac reconstruction algorithm with cardiac banding artifact correction. , 2006, Medical physics.

[35]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Globally Convergent Algorithms for Maximum a Posteriori Transmission Tomography , 2022 .

[36]  Gabor T. Herman,et al.  Algebraic reconstruction techniques can be made computationally efficient [positron emission tomography application] , 1993, IEEE Trans. Medical Imaging.

[37]  Daniel Rueckert,et al.  Nonrigid registration using free-form deformations: application to breast MR images , 1999, IEEE Transactions on Medical Imaging.

[38]  Thomas Köhler,et al.  Efficient projection and backprojection scheme for spherically symmetric basis functions in divergent beam geometry. , 2006, Medical physics.

[39]  Michael Unser,et al.  B-spline signal processing. I. Theory , 1993, IEEE Trans. Signal Process..

[40]  M Grass,et al.  A motion-compensated scheme for helical cone-beam reconstruction in cardiac CT angiography. , 2008, Medical physics.

[41]  J. Spall Multivariate stochastic approximation using a simultaneous perturbation gradient approximation , 1992 .

[42]  Rémy Prost,et al.  Motion Correction for Coronary Stent Reconstruction From Rotational X-ray Projection Sequences , 2007, IEEE Transactions on Medical Imaging.

[43]  Fred L. Bookstein,et al.  Principal Warps: Thin-Plate Splines and the Decomposition of Deformations , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[44]  Michael A. King,et al.  An interior point iterative maximum-likelihood reconstruction algorithm incorporating upper and lower bounds with application to SPECT transmission imaging , 2001, IEEE Transactions on Medical Imaging.

[45]  Michael Unser,et al.  Multiresolution approximation using shifted splines , 1998, IEEE Trans. Signal Process..

[46]  Michael Unser,et al.  Splines: a perfect fit for signal and image processing , 1999, IEEE Signal Process. Mag..

[47]  Torsten Rohlfing,et al.  Volume-preserving nonrigid registration of MR breast images using free-form deformation with an incompressibility constraint , 2003, IEEE Transactions on Medical Imaging.

[48]  Michael Grass,et al.  Iterative reconstruction of a region of interest for transmission tomography. , 2008 .

[49]  Akram Aldroubi,et al.  B-spline signal processing. II. Efficiency design and applications , 1993, IEEE Trans. Signal Process..

[50]  R M Lewitt,et al.  Multidimensional digital image representations using generalized Kaiser-Bessel window functions. , 1990, Journal of the Optical Society of America. A, Optics and image science.

[51]  A. Boyer,et al.  Radiation dose reduction in four-dimensional computed tomography. , 2005, Medical physics.

[52]  Robert M. Lewitt,et al.  Practical considerations for 3-D image reconstruction using spherically symmetric volume elements , 1996, IEEE Trans. Medical Imaging.

[53]  Michael Unser,et al.  The L2-Polynomial Spline Pyramid , 1993, IEEE Trans. Pattern Anal. Mach. Intell..

[54]  M Desco,et al.  Cardiac motion analysis from ultrasound sequences using nonrigid registration: validation against Doppler tissue velocity. , 2006, Ultrasound in medicine & biology.

[55]  E. Fishman,et al.  Toward Time Resolved Cardiac CT Images with Patient Dose Reduction: Image-based Motion Estimation , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[56]  Günter Lauritsch,et al.  Towards cardiac C-arm computed tomography , 2006, IEEE Transactions on Medical Imaging.

[57]  Jeongtae Kim Intensity based image registration using robust similarity measure and constrained optimization: Applications for radiation therapy. , 2004 .

[58]  Olivier Ecabert,et al.  Local Cardiac Wall Motion Estimation from Retrospectively Gated CT Images , 2009, FIMH.

[59]  Mark Holden,et al.  A Review of Geometric Transformations for Nonrigid Body Registration , 2008, IEEE Transactions on Medical Imaging.

[60]  T Nielsen,et al.  Cardiac cone-beam CT volume reconstruction using ART. , 2005, Medical physics.