Maximum likelihood, least squares, and penalized least squares for PET

The EM algorithm is the basic approach used to maximize the log likelihood objective function for the reconstruction problem in positron emission tomography (PET). The EM algorithm is a scaled steepest ascent algorithm that elegantly handles the nonnegativity constraints of the problem. It is shown that the same scaled steepest descent algorithm can be applied to the least squares merit function, and that it can be accelerated using the conjugate gradient approach. The experiments suggest that one can cut the computation by about a factor of 3 by using this technique. The results are applied to various penalized least squares functions which might be used to produce a smoother image.

[1]  O. Nalcioglu,et al.  Constrained Iterative Reconstruction by the Conjugate Gradient Method , 1985, IEEE Transactions on Medical Imaging.

[2]  M. Bierlaire,et al.  On iterative algorithms for linear least squares problems with bound constraints , 1991 .

[3]  G T Herman,et al.  Performance evaluation of an iterative image reconstruction algorithm for positron emission tomography. , 1991, IEEE transactions on medical imaging.

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

[5]  Gabor T. Herman,et al.  Evaluation And Optimization Of Iterative Reconstruction Techniques , 1990, [1990] Proceedings of the Twelfth Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[6]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[7]  O. Axelsson.  Conjugate gradient type methods for unsymmetric and inconsistent systems of linear equations , 1980 .

[8]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[9]  M I Miller,et al.  Maximum likelihood SPECT in clinical computation times using mesh-connected parallel computers. , 1991, IEEE transactions on medical imaging.

[10]  C. M. Reeves,et al.  Function minimization by conjugate gradients , 1964, Comput. J..

[11]  Eric C. Frey,et al.  Comparison between ML-EM and WLS-CG algorithms for SPECT image reconstruction , 1991 .

[12]  E. Levitan,et al.  A Maximum a Posteriori Probability Expectation Maximization Algorithm for Image Reconstruction in Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

[13]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[14]  A. Dhawan,et al.  A multigrid expectation maximization reconstruction algorithm for positron emission tomography. , 1988, IEEE transactions on medical imaging.

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

[16]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[17]  T. Pan,et al.  Numerical study of multigrid implementations of some iterative image reconstruction algorithms , 1991, Conference Record of the 1991 IEEE Nuclear Science Symposium and Medical Imaging Conference.

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

[19]  Linda Kaufman,et al.  Solving emission tomography problems on vector machines , 1990 .

[20]  E. Polak.  Introduction to linear and nonlinear programming , 1973 .

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

[22]  Linda Kaufman,et al.  Implementing and Accelerating the EM Algorithm for Positron Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

[23]  Åke Bjürck.  METHODS FOR SPARSE LINEAR LEAST SQUARES PROBLEMS , 1976 .

[24]  L. Shepp,et al.  A Statistical Model for Positron Emission Tomography , 1985 .

[25]  Zang-Hee Cho,et al.  Parallelization of the EM algorithm for 3-D PET image reconstruction , 1991 .

[26]  T. Hebert,et al.  A generalized EM algorithm for 3-D Bayesian reconstruction from Poisson data using Gibbs priors. , 1989, IEEE transactions on medical imaging.