Motion-compensated iterative cone-beam CT image reconstruction with adapted blobs as basis functions

This paper presents a three-dimensional method to reconstruct moving objects from cone-beam X-ray projections using an iterative reconstruction algorithm and a given motion vector field. For the image representation, adapted blobs are used, which can be implemented efficiently as basis functions. Iterative reconstruction requires the calculation of line integrals (forward projections) through the image volume, which are compared with the actual measurements to update the image volume. In the existence of a divergent motion vector field, a change in the volumes of the blobs has to be taken into account in the forward and backprojections. An efficient method to calculate the line integral through the adapted blobs is proposed. It solves the problem, how to compensate for the divergence in the motion vector field on a grid of basis functions. The method is evaluated on two phantoms, which are subject to three different known motions. Moreover, a motion-compensated filtered back-projection reconstruction method is used, and the reconstructed images are compared. Using the correct motion vector field with the iterative motion-compensated reconstruction, sharp images are obtained, with a quality that is significantly better than gated reconstructions.

[1]  Wojciech Zbijewski,et al.  Characterization and suppression of edge and aliasing artefacts in iterative x-ray CT reconstruction. , 2004, Physics in medicine and biology.

[2]  J. Hornegger,et al.  Cardiac C-arm CT: Efficient Motion Correction for 4D-FBP , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

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

[4]  T. Kohler,et al.  SNR-weighted ART applied to transmission tomography , 2003, 2003 IEEE Nuclear Science Symposium. Conference Record (IEEE Cat. No.03CH37515).

[5]  Alvaro R. De Pierro,et al.  Fast EM-like methods for maximum "a posteriori" estimates in emission tomography , 2001, IEEE Transactions on Medical Imaging.

[6]  R. Bracewell.  The Fourier transform. , 1989, Scientific American.

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

[8]  M. Defrise,et al.  Cone-beam filtered-backprojection algorithm for truncated helical data. , 1998, Physics in medicine and biology.

[9]  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.

[10]  W A Kalender,et al.  Electrocardiogram-correlated image reconstruction from subsecond spiral computed tomography scans of the heart. , 1998, Medical physics.

[11]  Pierre Grangeat,et al.  Exact reconstruction in 2D dynamic CT: compensation of time-dependent affine deformations. , 2004, Physics in medicine and biology.

[12]  Robert M. Lewitt,et al.  Application of the row action maximum likelihood algorithm with spherical basis functions to clinical PET imaging , 2001 .

[13]  R Proksa,et al.  Noise and resolution in images reconstructed with FBP and OSC algorithms for CT. , 2007, Medical physics.

[14]  B. De Man,et al.  A study of noise and spatial resolution for 2D and 3D filtered backprojection reconstruction , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

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

[16]  Georges Voronoi.  Nouvelles applications des paramètres continus à la théorie des formes quadratiques. Premier mémoire. Sur quelques propriétés des formes quadratiques positives parfaites. , 1908 .

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

[18]  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.

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

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

[21]  K. Mueller,et al.  Anti-aliased three-dimensional cone-beam reconstruction of low-contrast objects with algebraic methods , 1999, IEEE Transactions on Medical Imaging.

[22]  Yongmin Kim,et al.  Correction of computed tomography motion artifacts using pixel-specific back-projection , 1996, IEEE Trans. Medical Imaging.

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

[24]  Á. R. De Pierro,et al.  Fast EM-like methods for maximum "a posteriori" estimates in emission tomography. , 2001, IEEE transactions on medical imaging.

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

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

[27]  Jed D. Pack,et al.  Dynamic computed tomography with known motion field , 2004, SPIE Medical Imaging.

[28]  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.

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

[30]  L. J. Thomas,et al.  Noise and Edge Artifacts in Maximum-Likelihood Reconstructions for Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

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

[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]  Michael Grass,et al.  ECG gated circular cone-beam multi-cycle short-scan reconstruction algorithm , 2007, SPIE Medical Imaging.

[34]  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.

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

[36]  Alexander Katsevich,et al.  Theoretically exact FBP-type inversion algorithm for spiral CT , 2001 .

[37]  Ronald N. Bracewell,et al.  The Fourier Transform and Its Applications , 1966 .