Alternating Minimization Algorithms for Transmission Tomography

A family of alternating minimization algorithms for finding maximum-likelihood estimates of attenuation functions in transmission X-ray tomography is described. The model from which the algorithms are derived includes polyenergetic photon spectra, background events, and nonideal point spread functions. The maximum-likelihood image reconstruction problem is reformulated as a double minimization of the I-divergence. A novel application of the convex decomposition lemma results in an alternating minimization algorithm that monotonically decreases the objective function. Each step of the minimization is in closed form. The family of algorithms includes variations that use ordered subset techniques for increasing the speed of convergence. Simulations demonstrate the ability to correct the cupping artifact due to beam hardening and the ability to reduce streaking artifacts that arise from beam hardening and background events

[1]  Joseph A O'Sullivan,et al.  Image reconstruction for transmission tomography when projection data are incomplete , 2006, Physics in medicine and biology.

[2]  Alvaro R. De Pierro,et al.  A row-action alternative to the EM algorithm for maximizing likelihood in emission tomography , 1996, IEEE Trans. Medical Imaging.

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

[4]  Lasse Riis Østergaard,et al.  Early detection of AD using cortical thickness measurements , 2007, SPIE Medical Imaging.

[5]  I. Csiszár.  Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems , 1991 .

[6]  Joseph A. O'Sullivan,et al.  Alternating minimization algorithm for dual energy X-ray CT , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[7]  H. V. Trees.  Detection, Estimation, And Modulation Theory , 2001 .

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

[9]  J. O’Sullivan.  Alternating Minimization Algorithms: From Blahut-Arimoto to Expectation-Maximization , 1998 .

[10]  Joseph A. O'Sullivan,et al.  Alternating minimization algorithms for x-ray computed tomography: multigrid acceleration and dual energy application , 2005 .

[11]  Jeffrey A. Fessler,et al.  Statistical x-ray-computed tomography image reconstruction with beam- hardening correction , 2001, SPIE Medical Imaging.

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

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

[14]  John M. Ollinger,et al.  Maximum-likelihood reconstruction of transmission images in emission computed tomography via the EM algorithm , 1994, IEEE Trans. Medical Imaging.

[15]  Joseph A. O'Sullivan,et al.  Roughness penalties on finite domains , 1995, IEEE Trans. Image Process..

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

[17]  Yair Censor,et al.  Proximity Function Minimization Using Multiple Bregman Projections, with Applications to Split Feasibility and Kullback–Leibler Distance Minimization , 2001, Ann. Oper. Res..

[18]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 2002, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[19]  Jeffrey A. Fessler,et al.  Globally convergent image reconstruction for emission tomography using relaxed ordered subsets algorithms , 2003, IEEE Transactions on Medical Imaging.

[20]  Charles L. Byrne,et al.  Iterative image reconstruction algorithms based on cross-entropy minimization , 1993, IEEE Trans. Image Process..

[21]  J. P. Stonestrom,et al.  A Framework for Spectral Artifact Corrections in X-Ray CT , 1981, IEEE Transactions on Biomedical Engineering.

[22]  Joseph A. O'Sullivan,et al.  Pose estimation of known objects during transmission tomographic image reconstruction , 2006, IEEE Transactions on Medical Imaging.

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

[24]  Joseph A. O'Sullivan.  Iterative Algorithms for Maximum Likelihood Sequence Detection , 2002 .

[25]  Donald L. Snyder,et al.  Transmission tomographic image reconstruction using alternating minimization algorithms , 2004 .

[26]  J. Darroch,et al.  Generalized Iterative Scaling for Log-Linear Models , 1972 .

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

[28]  M.I. Miller,et al.  The role of likelihood and entropy in incomplete-data problems: Applications to estimating point-process intensities and toeplitz constrained covariances , 1987, Proceedings of the IEEE.

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

[30]  J. O’Sullivan,et al.  Properties of preprocessed sinogram data in x-ray computed tomography. , 2006, Medical physics.

[31]  M. Defrise,et al.  Iterative reconstruction for helical CT: a simulation study. , 1998, Physics in medicine and biology.

[32]  Joseph A. O'Sullivan,et al.  Information-Theoretic Image Formation , 1998, IEEE Trans. Inf. Theory.

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

[34]  Freek J. Beekman,et al.  Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm , 1998, IEEE Transactions on Medical Imaging.

[35]  A. R. De Pierro,et al.  On the relation between the ISRA and the EM algorithm for positron emission tomography , 1993, IEEE Trans. Medical Imaging.

[36]  J. Fessler.  Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[37]  J. Fitzpatrick,et al.  Medical image processing and analysis , 2000 .

[38]  Joseph A O'Sullivan,et al.  Prospects for quantitative computed tomography imaging in the presence of foreign metal bodies using statistical image reconstruction. , 2002, Medical physics.

[39]  Bruce R. Whiting,et al.  Signal statistics in x-ray computed tomography , 2002, SPIE Medical Imaging.

[40]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[41]  Joseph A. O'Sullivan,et al.  Deblurring subject to nonnegativity constraints , 1992, IEEE Trans. Signal Process..

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

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

[44]  B. Whiting,et al.  On two-parameter models of photon cross sections: application to dual-energy CT imaging. , 2006, Medical physics.