The convergence of object dependent resolution in maximum likelihood based tomographic image reconstruction.

Study of the maximum likelihood by EM algorithm (ML) with a reconstruction kernel equal to the intrinsic detector resolution and sieve regularization has demonstrated that any image improvements over filtered backprojection (FBP) are a function of image resolution. Comparing different reconstruction algorithms potentially requires measuring and matching the image resolution. Since there are no standard methods for describing the resolution of images from a nonlinear algorithm such as ML, we have defined measures of effective local Gaussian resolution (ELGR) and effective global Gaussian resolution (EGGR) and examined their behaviour in FBP images and in ML images using two different measurement techniques. For FBP these two resolution measures are equal and exhibit the standard convolution behaviour of linear systems. For ML, the FWHM of the ELGR monotonically increased with decreasing Gaussian object size due to slower convergence rates for smaller objects. For the simple simulated phantom used, this resolution dependence is independent of object position. With increasing object size, number of iterations and sieve size the object size dependence of the ELGR decreased. The FWHM of the EGGR converged after approximately 200 iterations, masking the fact that the ELGR for small objects was far from convergence. When FBP is compared to a nonlinear algorithm such as ML, it is recommended that at least the EGGR be matched; for ML this requires more than the number of iterations (e.g., < 100) that are typically run to minimize the mean square error or to satisfy a feasibility or similar stopping criterion. For many tasks, matching the EGGR of ML to FBP images may be insufficient and >> 200 iterations may be needed, particularly for small objects in the ML image because their ELGR has not yet converged.

[1]  E. Veklerov,et al.  Stopping Rule for the MLE Algorithm Based on Statistical Hypothesis Testing , 1987, IEEE Transactions on Medical Imaging.

[2]  Eugene Veklerov,et al.  The feasibility of images reconstructed with the methods of sieves , 1990 .

[3]  Eiichi Tanaka,et al.  Improved iterative image reconstruction with automatic noise artifact suppression , 1992, IEEE Trans. Medical Imaging.

[4]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[5]  T. Pan,et al.  Numerical study of multigrid implementations of some iterative image reconstruction algorithms , 1991, Conference Record of the 1991 IEEE Nuclear Science Symposium and Medical Imaging Conference.

[6]  S. Holte,et al.  Iterative image reconstruction for positron emission tomography: a study of convergence and quantitation problems , 1990 .

[7]  R. F. Wagner,et al.  Toward a unified view of radiological imaging systems. Part II: Noisy images. , 1977, Medical physics.

[8]  G T Herman,et al.  Performance evaluation of an iterative image reconstruction algorithm for positron emission tomography. , 1991, IEEE transactions on medical imaging.

[9]  L. J. Thomas,et al.  Noise and Edge Artifacts in Maximum-Likelihood Reconstructions for Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

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

[11]  R. Carson A Maximum Likelihood Method for Region-of-Interest Evaluation in Emission Tomography , 1986, Journal of computer assisted tomography.

[12]  David G. Politte Image improvements in positron-emission tomography due to measuring differential time-of-flight and using maximum-likelihood estimation , 1990 .

[13]  S. Strother,et al.  Practical tradeoffs between noise, quantitation, and number of iterations for maximum likelihood-based reconstructions. , 1991, IEEE transactions on medical imaging.

[14]  R. Kessler,et al.  Analysis of emission tomographic scan data: limitations imposed by resolution and background. , 1984, Journal of computer assisted tomography.

[15]  Kenneth F. Koral,et al.  Object-dependent performance comparison of two iterative reconstruction algorithms , 1988 .