Accelerated image reconstruction using ordered subsets of projection data

The authors define ordered subset processing for standard algorithms (such as expectation maximization, EM) for image restoration from projections. Ordered subsets methods group projection data into an ordered sequence of subsets (or blocks). An iteration of ordered subsets EM is defined as a single pass through all the subsets, in each subset using the current estimate to initialize application of EM with that data subset. This approach is similar in concept to block-Kaczmarz methods introduced by Eggermont et al. (1981) for iterative reconstruction. Simultaneous iterative reconstruction (SIRT) and multiplicative algebraic reconstruction (MART) techniques are well known special cases. Ordered subsets EM (OS-EM) provides a restoration imposing a natural positivity condition and with close links to the EM algorithm. OS-EM is applicable in both single photon (SPECT) and positron emission tomography (PET). In simulation studies in SPECT, the OS-EM algorithm provides an order-of-magnitude acceleration over EM, with restoration quality maintained.

[1]  M. Miller,et al.  An evaluation of maximum likelihood reconstruction for SPECT. , 1990, IEEE transactions on medical imaging.

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

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

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

[5]  P. Gilbert.  Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

[6]  A. Lent,et al.  Iterative algorithms for large partitioned linear systems, with applications to image reconstruction , 1981 .

[7]  B F Hutton,et al.  Use of 3D reconstruction to correct for patient motion in SPECT. , 1994, Physics in medicine and biology.

[8]  Gerald Q. Maguire,et al.  Fusion of Radiolabeled Monoclonal Antibody SPECT Images with CT/MRI Images , 1992 .

[9]  Marc Teboulle,et al.  A primal-dual iterative algorithm for a maximum likelihood estimation problem , 1992 .

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

[11]  Y. Censor,et al.  On Block-Iterative Entropy Maximization , 1987 .

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

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

[14]  D. Titterington.  Recursive Parameter Estimation Using Incomplete Data , 1984 .

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

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

[17]  Y. Censor.  Finite series-expansion reconstruction methods , 1983, Proceedings of the IEEE.

[18]  Khairil Anwar Notodiputro,et al.  The use of small scale prototypes in image reconstruction from projections , 1989 .