Modified MLEM Algorithm for Artifact Suppression in CT

Iterative methods for CT image reconstruction are widely disregarded in clinical routine due to the inherent massive computational effort. On the other hand, the continuously growing computational power of today's standard computers has led to a rediscovery of these methods. Currently, filtered backprojection (FBP) is the reconstruction method of choice, because it is very fast, especially on dedicated hardware. However, a disadvantage of FBP is that it ignores some of the physical effects like beam hardening, scatter and noise, leading to artifacts in the reconstructed images. Artifacts of this type are particularly dominant if metal objects are inside the patient, because FBP interprets the corresponding projection data as inconsistent. For the method proposed here, the inconsistency of the metal-influenced projection data is reduced with an interpolation scheme working in the Radon-space domain. As the interpolation yields a vague estimate of the real soft-tissue projection values, residual inconsistencies remain. Those residuals can be treated by a modified maximum likelihood expectation maximization (MLEM) method presented in this paper. Within the novel approach, different projection lines through the object can be associated with appropriate weightings that decrease the influence of the residual inconsistencies - resulting in an overall image-quality enhancement. The approach is demonstrated on real CT measurements of a torso phantom equipped with metal markers that allows an evaluation with the corresponding ground truth data.

[1]  P. Joseph,et al.  A Method for Correcting Bone Induced Artifacts in Computed Tomography Scanners , 1978, Journal of computer assisted tomography.

[2]  H. Skinner,et al.  CT image correction for beam hardening using simulated projection data , 1990 .

[3]  P M Joseph,et al.  The exponential edge-gradient effect in x-ray computed tomography. , 1981, Physics in medicine and biology.

[4]  Donald Geman,et al.  Bayesian Image Analysis , 1986 .

[5]  Bruno De Man,et al.  Iterative reconstruction for reduction of metal artifacts in computed tomography , 2001 .

[6]  P M Joseph,et al.  The effects of sampling on CT images. , 1980, Computerized tomography.

[7]  A. Macovski,et al.  Nonlinear Polychromatic and Noise Artifacts in X‐Ray Computed Tomography Images , 1979, Journal of computer assisted tomography.

[8]  K. Lange,et al.  A Theoretical Study of Some Maximum Likelihood Algorithms for Emission and Transmission Tomography , 1987, IEEE Transactions on Medical Imaging.

[9]  Patrick Dupont,et al.  An iterative maximum-likelihood polychromatic algorithm for CT , 2001, IEEE Transactions on Medical Imaging.

[10]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[11]  S. Zhao,et al.  X-ray CT metal artifact reduction using wavelets: an application for imaging total hip prostheses , 2000, IEEE Transactions on Medical Imaging.

[12]  P. Rüegsegger,et al.  CT BONE DENSITOMETRY OF THE ANCHORAGE OF ARTIFICIAL KNEE JOINTS , 1985 .

[13]  N J Pelc,et al.  Nonlinear partial volume artifacts in x-ray computed tomography. , 1980, Medical physics.

[14]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[15]  J. Nuyts,et al.  Reduction of attenuation correction artifacts in PET-CT , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[16]  P M Joseph,et al.  A method for simultaneous correction of spectrum hardening artifacts in CT images containing both bone and iodine. , 1997, Medical physics.

[17]  A. Macovski,et al.  Polychromatic streak artifacts in computed tomography images. , 1978, Journal of computer assisted tomography.

[18]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[19]  Jeffrey A. Fessler,et al.  Iterative tomographic image reconstruction using Fourier-based forward and back-projectors , 2004, IEEE Transactions on Medical Imaging.

[20]  Gary H. Glover,et al.  Compton scatter effects in CT reconstructions , 1982 .

[21]  Joseph A. O'Sullivan,et al.  Iterative deblurring for CT metal artifact reduction , 1996, IEEE Trans. Medical Imaging.

[22]  Patrick Dupont,et al.  Reduction of metal streak artifacts in X-ray computed tomography using a transmission maximum a posteriori algorithm , 1999 .

[23]  D. Robertson,et al.  Total hip prosthesis metal-artifact suppression using iterative deblurring reconstruction. , 1997, Journal of computer assisted tomography.

[24]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.

[25]  E. Fishman,et al.  Evaluation of CT techniques for reducing artifacts in the presence of metallic orthopedic implants. , 1988, Journal of computer assisted tomography.

[26]  P. Rüegsegger,et al.  Anchorage of femoral implants visualized by modified computed tomography , 2004, Archives of orthopaedic and traumatic surgery.

[27]  Peter M. Joseph,et al.  View sampling requirements in fan beam computed tomography. , 1980 .

[28]  P. Rüegsegger,et al.  Computed Tomography Reconstruction from Hollow Projections: An Application to In Vivo Evaluation of Artificial Hip Joints , 1979, Journal of computer assisted tomography.

[29]  Victor J. Sank,et al.  IMAGE RECONSTRUCTION FROM PROJECTIONS: ***I , 1978 .

[30]  J A Veiga-Pires,et al.  Artefacts in CT scanning. , 1979, The British journal of radiology.

[31]  Predrag Sukovic,et al.  Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography , 2000, IEEE Transactions on Medical Imaging.

[32]  Chye Hwang Yan,et al.  Reconstruction algorithm for polychromatic CT imaging: application to beam hardening correction , 2000, IEEE Transactions on Medical Imaging.

[33]  William R. Brody,et al.  Iterative convolution backprojection algorithms for image reconstruction from limited data , 1983 .

[34]  P M Joseph,et al.  The effects of scatter in x-ray computed tomography. , 1982, Medical physics.

[35]  C. Adey,et al.  Acute Coronary Intervention , 1988 .

[36]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[37]  K. Kletter,et al.  Analyse des Konvergenzverhaltens von Rekonstruktions-algorithmen anhand lokaler und globaler Bildparameter , 2001 .

[38]  Jeffrey A. Fessler,et al.  Globally Convergent Algorithms for Maximum a , 1995 .

[39]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[40]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[41]  Detlef Zerfowski Kompensation von Metallartefakten in der Computertomographie , 1998, Bildverarbeitung für die Medizin.

[42]  Gary H. Glover,et al.  An algorithm for the reduction of metal clip artifacts in CT reconstructions. , 1981 .

[43]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[44]  Jonas August Decoupling the equations of regularized tomography , 2002, Proceedings IEEE International Symposium on Biomedical Imaging.

[45]  Willi A. Kalender,et al.  Algorithms for the reduction of CT artifacts caused by metallic implants , 1990, Medical Imaging.

[46]  Jeffrey A. Fessler Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[47]  P. Suetens,et al.  Metal streak artifacts in X-ray computed tomography: a simulation study , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[48]  K Freeman,et al.  CT scans through metal scanning technique versus hardware composition. , 1994, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[49]  W. Kalender,et al.  Reduction of CT artifacts caused by metallic implants. , 1987 .

[50]  K. Lange Convergence of EM image reconstruction algorithms with Gibbs smoothing. , 1990, IEEE transactions on medical imaging.

[51]  D D Robertson,et al.  Quantitative bone measurements using x-ray computed tomography with second-order correction. , 1986, Medical physics.

[52]  Ryan J. Murphy,et al.  Deblurring subject to nonnegativity constraints when known functions are present with application to object-constrained computerized tomography , 2001, IEEE Transactions on Medical Imaging.

[53]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[54]  B. E. Oppenheim,et al.  Reconstruction tomography from incomplete projections , 1975 .

[55]  Jeffrey A. Fessler,et al.  Globally convergent ordered subsets algorithms: application to tomography , 2001, 2001 IEEE Nuclear Science Symposium Conference Record (Cat. No.01CH37310).