A Maximum a Posteriori Probability Expectation Maximization Algorithm for Image Reconstruction in Emission Tomography

The expectation maximization method for maximum likelihood image reconstruction in emission tomography, based on the Poisson distribution of the statistically independent components of the image and measurement vectors, is extended to a maximum aposteriori image reconstruction using a multivariate Gaussian a priori probability distribution of the image vector. The approach is equivalent to a penalized maximum likelihood estimation with a special choice of the penalty function. The expectation maximization method is applied to find the a posteriori probability maximizer. A simple iterative formula is derived for a penalty function that is a weighted sum of the squared deviations of image vector components from their a priori mean values. The method is demonstrated to be superior to pure likelihood maximization, in that the penalty function prevents the occurrence of irregular high amplitude patterns in the image with a large number of iterations (the so-called "checkerboard effect" or "noise artifact").

[1]  Abderrahmane Haddad,et al.  Estimation theory with applications to communications and control , 1972 .

[2]  H. Hurwitz,et al.  Entropy reduction in Bayesian analysis of measurements , 1975 .

[3]  Gabor T. Herman,et al.  A Computer Implementation of a Bayesian Analysis of Image Reconstruction , 1976, Inf. Control..

[4]  Azriel Rosenfeld,et al.  Digital Picture Processing , 1976 .

[5]  B. R. Hunt,et al.  Digital Image Restoration , 1977 .

[6]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[7]  R. Tapia,et al.  Nonparametric Probability Density Estimation , 1978 .

[8]  S. Rowland,et al.  Computer implementation of image reconstruction formulas , 1979 .

[9]  Gabor T. Herman,et al.  On the Bayesian Approach to Image Reconstruction , 1979, Inf. Control..

[10]  Gabor T. Herman,et al.  Image reconstruction from projections : the fundamentals of computerized tomography , 1980 .

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

[12]  G. W. Wecksung,et al.  Bayesian approach to limited-angle reconstruction in computed tomography , 1983 .

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

[14]  Michael I. Miller,et al.  The Use of Sieves to Stabilize Images Produced with the EM Algorithm for Emission Tomography , 1985, IEEE Transactions on Nuclear Science.

[15]  L. Shepp,et al.  A Statistical Model for Positron Emission Tomography , 1985 .

[16]  Gabor T. Herman.  Application of Maximum Entropy and Bayesian Optimization Methods to Image Reconstruction from Projections , 1985 .

[17]  J. Skilling,et al.  Maximum-entropy and Bayesian methods in inverse problems , 1985 .

[18]  Stephen M. Moore,et al.  An Evaluation of the Use of Sieves for Producing Estimates. Of Radioactivity Distributions with the EM Algorithm for PET , 1986, IEEE Transactions on Nuclear Science.

[19]  B. Yandell,et al.  Automatic Smoothing of Regression Functions in Generalized Linear Models , 1986 .

[20]  Jorge Llacer,et al.  On The Applicability Of The Maximum Likelihood Estimator Algorithm For Image Recovery In Accelerated Positron Emitter Beam Injection , 1986, Other Conferences.

[21]  Z. Liang,et al.  Bayesian image processing of data from constrained source distributions—II. valued, uncorrelated and correlated constraints , 1987 .

[22]  Z. Liang,et al.  Bayesian image processing of data from constrained source distributions—I. non-valued, uncorrelated and correlated constraints , 1987 .

[23]  M. Melamed.  Detection , 2021, SETI: Astronomy as a Contact Sport.