The use of a stopping rule in iterative image reconstruction

During a study of the characteristics of the Maximum Likelihood Estimator (MLE) method of image reconstruction from Positron Emission Tomography (PET) data, we have found that the requirement that the reconstructed image, if it were a source, could have generated the original data imposes the stopping of the iterative procedure at some point. This requirement appears to be in contradiction with the concept of a maximum likelihood estimator until one realizes that maximizing the likelihood results in reaching for the top of the measurement probability distribution P(g\f) in Bayes' formula: P(f\g) = P(g\f)P(f)/P(g) where the results measurement vector g can be described by g = Hf+ n, with / being a source, H the response matrix of the measurement system and ? the noise vector in the measurement. From a Bayesian point of view, we should maximize the Maximum a Posteriori Probability (MAP) p(f\g), i.e., the probability that the source be / given a measurement g, and for that purpose one needs to have the a priori source distribution -P(/), which is not usually available with any degree of reliability. The use of the stopping rule does not reach for the MAP solution but it uses some important physical prior information: it requires that the iterations stop as soon as the image obtained has a characteristic that we call "feasibility". We define a feasible image as one that could have given the original data by the physical process that governs the measurement. We show that MLE reconstructions started from a uniform image field pass through a region of feasibility in which the images represent a good compromise between sharpness and smooth regions of high activity, without the "noise artifact" having yet set in. The shape and characteristics of the region of feasibility are described and implications for future work are described.

[1]  E. L. Alpen,et al.  Imaging by Injection of Accelerated Radioactive Particle Beams , 1984, IEEE Transactions on Medical Imaging.

[2]  S. Gull,et al.  Image reconstruction from incomplete and noisy data , 1978, Nature.

[3]  E. Veklerov,et al.  The high sensitivity of the maximum likelihood estimator method of tomographic image reconstruction , 1987 .

[4]  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.

[5]  J. Skilling,et al.  Maximum entropy image reconstruction: general algorithm , 1984 .

[6]  J. P. Burg,et al.  Maximum entropy spectral analysis. , 1967 .

[7]  R. Narayan,et al.  Maximum Entropy Image Restoration in Astronomy , 1986 .

[8]  M.I. Miller,et al.  The role of likelihood and entropy in incomplete-data problems: Applications to estimating point-process intensities and toeplitz constrained covariances , 1987, Proceedings of the IEEE.

[9]  H H Barrett,et al.  Reconstruction of objects from coded images by simulated annealing. , 1983, Optics letters.

[10]  Edward J. Hoffman,et al.  On The Convergence Of The Maximum Likelihood Estimator Method Of Tomographic Image Reconstruction , 1987, Medical Imaging.

[11]  Stuart Geman,et al.  Statistical methods for tomographic image reconstruction , 1987 .

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

[13]  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.

[14]  E. Levitan,et al.  A Maximum a Posteriori Probability Expectation Maximization Algorithm for Image Reconstruction in Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

[15]  J. Pfleiderer,et al.  Improvement of MEM-deconvolution by an additional constraint. , 1986 .