Title :
Hybrid Poisson/polynomial objective functions for tomographic image reconstruction from transmission scans
Author :
Fessler, Jeffrey A.
Author_Institution :
Dept. of Electr. Eng. & Comput. Sci., Michigan Univ., Ann Arbor, MI, USA
fDate :
10/1/1995 12:00:00 AM
Abstract :
This paper describes rapidly converging algorithms for computing attenuation maps from Poisson transmission measurements using penalized-likelihood objective functions. We demonstrate that an under-relaxed cyclic coordinate-ascent algorithm converges faster than the convex algorithm of Lange (see ibid., vol.4, no.10, p.1430-1438, 1995), which in turn converges faster than the expectation-maximization (EM) algorithm for transmission tomography. To further reduce computation, one could replace the log-likelihood objective with a quadratic approximation. However, we show with simulations and analysis that the quadratic objective function leads to biased estimates for low-count measurements. Therefore we introduce hybrid Poisson/polynomial objective functions that use the exact Poisson log-likelihood for detector measurements with low counts, but use computationally efficient quadratic or cubic approximations for the high-count detector measurements. We demonstrate that the hybrid objective functions reduce computation time without increasing estimation bias
Keywords :
approximation theory; convergence of numerical methods; image reconstruction; maximum likelihood estimation; medical image processing; polynomials; positron emission tomography; signal detection; stochastic processes; Poisson transmission measurements; attenuation maps; computation time reduction; convex algorithm; cubic approximations; detector measurements; estimation bias; exact Poisson log-likelihood; expectation-maximization algorithm; hybrid Poisson/polynomial objective functions; log-likelihood objective; low-count measurements; penalized-likelihood objective functions; positron emission tomography; quadratic approximation; quadratic approximations; quadratic objective function; rapidly converging algorithms; tomographic image reconstruction; transmission scans; transmission tomography; under-relaxed cyclic coordinate-ascent algorithm; Attenuation measurement; Biomedical imaging; Computational modeling; Detectors; Image converters; Image reconstruction; Lungs; Polynomials; Positron emission tomography; Thorax;
Journal_Title :
Image Processing, IEEE Transactions on