A unified approach to statistical tomography using coordinate descent optimization
- 1 March 1996
- journal article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Image Processing
- Vol. 5 (3), 480-492
- https://doi.org/10.1109/83.491321
Abstract
Over the past years there has been considerable interest in statistically optimal reconstruction of cross-sectional images from tomographic data. In particular, a variety of such algorithms have been proposed for maximum a posteriori (MAP) reconstruction from emission tomographic data. While MAP estimation requires the solution of an optimization problem, most existing reconstruction algorithms take an indirect approach based on the expectation maximization (EM) algorithm. We propose a new approach to statistically optimal image reconstruction based on direct optimization of the MAP criterion. The key to this direct optimization approach is greedy pixel-wise computations known as iterative coordinate decent (ICD). We propose a novel method for computing the ICD updates, which we call ICD/Newton-Raphson. We show that ICD/Newton-Raphson requires approximately the same amount of computation per iteration as EM-based approaches, but the new method converges much more rapidly (in our experiments, typically five to ten iterations). Other advantages of the ICD/Newton-Raphson method are that it is easily applied to MAP estimation of transmission tomograms, and typical convex constraints, such as positivity, are easily incorporated.Keywords
This publication has 36 references indexed in Scilit:
- Provably convergent coordinate descent in statistical tomographic reconstructionPublished by Institute of Electrical and Electronics Engineers (IEEE) ,2002
- A modified expectation maximization algorithm for penalized likelihood estimation in emission tomographyIEEE Transactions on Medical Imaging, 1995
- Acceleration and filtering in the generalized Landweber iteration using a variable shaping matrixIEEE Transactions on Medical Imaging, 1993
- Preconditioning methods for improved convergence rates in iterative reconstructionsIEEE Transactions on Medical Imaging, 1993
- Bayesian estimation of transmission tomograms using segmentation based optimizationIEEE Transactions on Nuclear Science, 1992
- The GEM MAP algorithm with 3-D SPECT system responseIEEE Transactions on Medical Imaging, 1992
- Bayesian reconstructions from emission tomography data using a modified EM algorithmIEEE Transactions on Medical Imaging, 1990
- Convergence of EM image reconstruction algorithms with Gibbs smoothingIEEE Transactions on Medical Imaging, 1990
- Fast MLE for SPECT using an intermediate polar representation and a stopping criterionIEEE Transactions on Nuclear Science, 1988
- Stopping Rule for the MLE Algorithm Based on Statistical Hypothesis TestingIEEE Transactions on Medical Imaging, 1987