Maximum-likelihood reconstruction of transmission images in emission computed tomography via the EM algorithm
- 1 March 1994
- journal article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Medical Imaging
- Vol. 13 (1) , 89-101
- https://doi.org/10.1109/42.276147
Abstract
The expectation-maximization (EM) algorithm for computing maximum-likelihood estimates of transmission images in positron-emission tomography (PET) (see K. Lange and R. Carson, J. Comput. Assist. Tomogr., vol.8, no.2, p.306-16, 1984) is extended to include measurement error, accidental coincidences and Compton scatter. A method for accomplishing the maximization step using one step of Newton's method is proposed. The algorithm is regularized with the method of sieves. Evaluations using both Monte Carlo simulations and phantom studies on the Siemens 953B scanner suggest that the algorithm yields unbiased images with significantly lower variances than filtered-backprojection when the images are reconstructed to the intrinsic resolution. Large features in the images converge in under 200 iterations while the smallest features required up to 2,000 iterations. All but the smallest features in typical transmission scans converge in approximately 250 iterations. The initial implementation of the algorithm requires 50 sec per iteration on a DECStation 5000.Keywords
This publication has 21 references indexed in Scilit:
- A statistical approach to transmission image reconstruction from ring source calibration measurements in PETPublished by Institute of Electrical and Electronics Engineers (IEEE) ,2003
- Maximum a posteriori estimation for SPECT using regularization techniques on massively parallel computersIEEE Transactions on Medical Imaging, 1993
- A local update strategy for iterative reconstruction from projectionsIEEE Transactions on Signal Processing, 1993
- Bayesian estimation of transmission tomograms using segmentation based optimizationIEEE Transactions on Nuclear Science, 1992
- Practical tradeoffs between noise, quantitation, and number of iterations for maximum likelihood-based reconstructionsIEEE Transactions on Medical Imaging, 1991
- The use of constraints to eliminate artifacts in maximum-likelihood image estimation for emission tomographyIEEE Transactions on Nuclear Science, 1988
- Noise and Edge Artifacts in Maximum-Likelihood Reconstructions for Emission TomographyIEEE Transactions on Medical Imaging, 1987
- Transmission Profile Filtering for Positron Emission TomographyIEEE Transactions on Nuclear Science, 1986
- On the Convergence Properties of the EM AlgorithmThe Annals of Statistics, 1983
- The Fourier reconstruction of a head sectionIEEE Transactions on Nuclear Science, 1974