The ability to theoretically model the propagation of photon noise through tomographic reconstruction algorithms is crucial in evaluating reconstructed image quality as a function of parameters of the algorithm. Here, the authors show the theoretical expressions for the propagation of Poisson noise through tomographic SPECT reconstructions using regularized EM algorithms with independent Gamma and multivariate Gaussian priors. The authors' analysis extends the work in H.H. Barrett et al., Phys. Med. Biol., vol. 39, p. 833-46 (1994), in which judicious linearizations were used to enable the propagation of a mean image and covariance matrix from one iteration to the next for the (unregularized) EM algorithm. To validate their theoretical analyses, the authors use a methodology in D.W. Wilson et al., Phys. Med. Biol., vol. 39, p. 847-71 (1994) to compare the results of theoretical calculations to Monte Carlo simulations. The authors also demonstrate an application of the theory to the calculation of an optimal smoothing parameter for a regularized reconstruction. The smoothing parameter is optimal in the context of a quantitation task, defined as the minimization of the expected mean-square error of an estimated number of counts in a hot lesion region. The authors' results thus demonstrate how the theory can be applied to a problem of potential practical use.
[1]
H H Barrett,et al.
Objective assessment of image quality: effects of quantum noise and object variability.
,
1990,
Journal of the Optical Society of America. A, Optics and image science.
[2]
F. Ashcroft,et al.
VIII. References
,
1955
.
[3]
Donald W. Wilson,et al.
Noise properties of the EM algorithm. I. Theory
,
1994
.
[4]
P. Green.
Bayesian reconstructions from emission tomography data using a modified EM algorithm.
,
1990,
IEEE transactions on medical imaging.
[5]
K. Lange,et al.
A Theoretical Study of Some Maximum Likelihood Algorithms for Emission and Transmission Tomography
,
1987,
IEEE Transactions on Medical Imaging.
[6]
P. Green.
On Use of the EM Algorithm for Penalized Likelihood Estimation
,
1990
.
[7]
B. Tsui,et al.
Noise properties of the EM algorithm: II. Monte Carlo simulations.
,
1994,
Physics in medicine and biology.