Interior-point methodology for 3-D PET reconstruction

Interior-point methods have been successfully applied to a wide variety of linear and nonlinear programming applications. This paper presents a class of algorithms, based on path-following interior-point methodology, for performing regularized maximum-likelihood (ML) reconstructions on three-dimensional (3-D) emission tomography data. The algorithms solve a sequence of subproblems that converge to the regularized maximum likelihood solution from the interior of the feasible region (the nonnegative orthant). The authors propose two methods, a primal method which updates only the primal image variables and a primal-dual method which simultaneously updates the primal variables and the Lagrange multipliers. A parallel implementation permits the interior-point methods to scale to very large reconstruction problems. Termination is based on well-defined convergence measures, namely, the Karush-Kuhn-Tucker first-order necessary conditions for optimality. The authors demonstrate the rapid convergence of the path-following interior-point methods using both data from a small animal scanner and Monte Carlo simulated data. The proposed methods can readily be applied to solve the regularized, weighted least squares reconstruction problem.

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

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

[3]  B. Tsui,et al.  Noise properties of the EM algorithm: II. Monte Carlo simulations. , 1994, Physics in medicine and biology.

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

[5]  D S Lalush,et al.  A fast and stable maximum a posteriori conjugate gradient reconstruction algorithm. , 1995, Medical physics.

[6]  Michael V. Green,et al.  Characteristics of a pair of small field-of-view LSO scintillation cameras , 1996 .

[7]  Alvaro R. De Pierro,et al.  On the convergence of an EM-type algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[8]  P. Green On Use of the EM Algorithm for Penalized Likelihood Estimation , 1990 .

[9]  Jeffrey A. Fessler,et al.  Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction , 1999, IEEE Trans. Image Process..

[10]  P M Bloomfield,et al.  PET scanners for small animals. , 1995, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[11]  S. Strother,et al.  Practical tradeoffs between noise, quantitation, and number of iterations for maximum likelihood-based reconstructions. , 1991, IEEE transactions on medical imaging.

[12]  Chung-Ming Chen,et al.  On Parallelizing the EM Algorithm for PET Image Reconstruction , 1994, IEEE Trans. Parallel Distributed Syst..

[13]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[14]  S. Nash Preconditioning of Truncated-Newton Methods , 1985 .

[15]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[16]  D. Lalush,et al.  The importance of preconditioners in fast Poisson-based iterative reconstruction algorithms for SPECT , 1995, 1995 IEEE Nuclear Science Symposium and Medical Imaging Conference Record.

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

[18]  Garry Chinn,et al.  A general class of preconditioners for statistical iterative reconstruction of emission computed tomography , 1997, IEEE Transactions on Medical Imaging.

[19]  S. Nash,et al.  Assessing a search direction within a truncated-newton method , 1990 .

[20]  B. Tsui,et al.  Noise properties of filtered-backprojection and ML-EM reconstructed emission tomographic images , 1992, IEEE Conference on Nuclear Science Symposium and Medical Imaging.

[21]  Alfred O. Hero,et al.  Ieee Transactions on Image Processing: to Appear Penalized Maximum-likelihood Image Reconstruction Using Space-alternating Generalized Em Algorithms , 2022 .

[22]  Ariela Sofer,et al.  Evaluation of 3D reconstruction algorithms for a small animal PET camera , 1996 .

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

[24]  Gabor T. Herman,et al.  Evaluation of reconstruction algorithms , 1991 .

[25]  P. Toint,et al.  A primal-dual algorithm for minimizing a non-convex function subject to bound and linear equality constraints , 2000 .

[26]  Gabor T. Herman,et al.  Image reconstruction from projections : the fundamentals of computerized tomography , 1980 .

[27]  Trond Steihaug,et al.  Truncated-newtono algorithms for large-scale unconstrained optimization , 1983, Math. Program..

[28]  A R De Pierro,et al.  On the convergence of an EM-type algorithm for penalized likelihood estimation in emission tomography. , 1995, IEEE transactions on medical imaging.

[29]  M. Kojima,et al.  A primal-dual interior point algorithm for linear programming , 1988 .

[30]  Stephen G. Nash,et al.  A Barrier Method for Large-Scale Constrained Optimization , 1993, INFORMS J. Comput..

[31]  W. Murray,et al.  Analytical expressions for the eigenvalues and eigenvectors of the Hessian matrices of barrier and penalty functions , 1971 .

[32]  Linda Kaufman,et al.  Maximum likelihood, least squares, and penalized least squares for PET , 1993, IEEE Trans. Medical Imaging.

[33]  T. K. Narayan,et al.  Evaluation of task-oriented performance of several fully 3D PET reconstruction algorithms. , 1994, Physics in medicine and biology.

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

[35]  Donald W. Wilson,et al.  Noise properties of the EM algorithm. I. Theory , 1994 .

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

[37]  Simon R. Cherry,et al.  Fast gradient-based methods for Bayesian reconstruction of transmission and emission PET images , 1994, IEEE Trans. Medical Imaging.

[38]  E U Mumcuoğlu,et al.  Bayesian reconstruction of PET images: methodology and performance analysis. , 1996, Physics in medicine and biology.

[39]  Anders Forsgren,et al.  Primal-Dual Interior Methods for Nonconvex Nonlinear Programming , 1998, SIAM J. Optim..

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

[41]  M. J. D. Powell,et al.  Nonlinear Programming—Sequential Unconstrained Minimization Techniques , 1969 .

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

[43]  Tinsu Pan,et al.  Preconditioning methods for improved convergence rates in iterative reconstructions , 1993, IEEE Trans. Medical Imaging.

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

[45]  Stephen G. Nash,et al.  Block truncated-Newton methods for parallel optimization , 1989, Math. Program..

[46]  Yuchen Yan,et al.  A system for the 3D reconstruction of retracted-septa PET data using the EM algorithm , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[47]  Roy E. Marsten,et al.  Feature Article - Interior Point Methods for Linear Programming: Computational State of the Art , 1994, INFORMS J. Comput..

[48]  J. M. Ollinger,et al.  Positron Emission Tomography , 2018, Handbook of Small Animal Imaging.

[49]  Paul Kinahan,et al.  Analytic 3D image reconstruction using all detected events , 1989 .

[50]  Stephen G. Nash,et al.  A General-Purpose Parallel Algorithm for Unconstrained Optimization , 1991, SIAM J. Optim..

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