A unified noise analysis for iterative image estimation.

Iterative image estimation methods have been widely used in emission tomography. Accurate estimation of the uncertainty of the reconstructed images is essential for quantitative applications. While both iteration-based noise analysis and fixed-point noise analysis have been developed, current iteration-based results are limited to only a few algorithms that have an explicit multiplicative update equation and some may not converge to the fixed-point result. This paper presents a theoretical noise analysis that is applicable to a wide range of preconditioned gradient-type algorithms. Under a certain condition, the proposed method does not require an explicit expression of the preconditioner. By deriving the fixed-point expression from the iteration-based result, we show that the proposed iteration-based noise analysis is consistent with fixed-point analysis. Examples in emission tomography and transmission tomography are shown. The results are validated using Monte Carlo simulations.

[1]  Alvaro R. De Pierro,et al.  A row-action alternative to the EM algorithm for maximizing likelihood in emission tomography , 1996, IEEE Trans. Medical Imaging.

[2]  J. D. Wilson,et al.  A smoothed EM approach to indirect estimation problems, with particular reference to stereology and emission tomography , 1990 .

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

[4]  Jeffrey A. Fessler,et al.  Channelized hotelling observer performance for penalized-likelihood image reconstruction , 2002, 2002 IEEE Nuclear Science Symposium Conference Record.

[5]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[6]  Charles L. Byrne,et al.  Convergent block-iterative algorithms for image reconstruction from inconsistent data , 1997, IEEE Trans. Image Process..

[7]  Richard M. Leahy,et al.  Resolution and noise properties of MAP reconstruction for fully 3-D PET , 2000, IEEE Transactions on Medical Imaging.

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

[9]  R. Huesman.  A new fast algorithm for the evaluation of regions of interest and statistical uncertainty in computed tomography. , 1984, Physics in medicine and biology.

[10]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[11]  K Thielemans,et al.  Enhanced 3D PET OSEM reconstruction using inter-update Metz filtering. , 2000, Physics in medicine and biology.

[12]  J. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[13]  G. Gindi,et al.  Noise propagation from attenuation correction into PET reconstructions , 2000 .

[14]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[15]  Yuxiang Xing,et al.  RAPID CALCULATION OF DETECTABILITY IN BAYESIAN SPECT , 2004 .

[16]  R H Huesman,et al.  Dynamic PET Data Analysis , 1986, Journal of computer assisted tomography.

[17]  G. Gindi,et al.  Noise analysis of MAP-EM algorithms for emission tomography. , 1997, Physics in medicine and biology.

[18]  Peter Herscovitch,et al.  An approximation formula for the variance of PET region-of-interest values , 1993, IEEE Trans. Medical Imaging.

[19]  P.E. Kinahan,et al.  Modeling and incorporation of system response functions in 3D whole body PET , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[20]  Ronald H. Huesman,et al.  Theoretical study of lesion detectability of MAP reconstruction using computer observers , 2001, IEEE Transactions on Medical Imaging.

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

[22]  Charles L. Byrne,et al.  Noise characterization of block-iterative reconstruction algorithms. I. Theory , 2000, IEEE Transactions on Medical Imaging.

[23]  Jorge Llacer,et al.  Optimization of Maximum Likelihood Estimator Images for PET: II. Statistical Analysis of Human Brain FDG Studies , 1991 .

[24]  Charles L. Byrne,et al.  Accelerating the EMML algorithm and related iterative algorithms by rescaled block-iterative methods , 1998, IEEE Trans. Image Process..

[25]  F. O’Sullivan,et al.  Variability Assessment in Positron Emission Tomography and Related Generalized Deconvolution Models , 1998 .

[26]  Eric C. Frey,et al.  Fast maximum entropy approximation in SPECT using the RBI-MAP algorithm , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

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

[28]  D. Lalush,et al.  Block-iterative techniques for fast 4D reconstruction using a priori motion models in gated cardiac SPECT. , 1998, Physics in medicine and biology.

[29]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 2002, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

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

[31]  Alvaro R. De Pierro,et al.  Fast EM-like methods for maximum "a posteriori" estimates in emission tomography , 2001, IEEE Transactions on Medical Imaging.

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

[33]  Richard M. Leahy,et al.  A theoretical study of the contrast recovery and variance of MAP reconstructions from PET data , 1999, IEEE Transactions on Medical Imaging.

[34]  Jeffrey A. Fessler.  Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): applications to tomography , 1996, IEEE Trans. Image Process..

[35]  Jeffrey A. Fessler,et al.  Regularization for uniform spatial resolution properties in penalized-likelihood image reconstruction , 2000, IEEE Transactions on Medical Imaging.