Noise properties of the EM algorithm. II. Monte Carlo simulations

• Published in: Physics in medicine & biology Vol. 39; no. 5; pp. 847 - 871
• Main Authors: Wilson, D W, Tsui, B M W, Barrett, H H
• Format: Journal Article
• Language: English
• Published: England IOP Publishing 01.05.1994
• Subjects: Algorithms; Brain - diagnostic imaging; Computer Simulation; Image Enhancement - methods; Image Interpretation, Computer-Assisted - methods; Likelihood Functions; Models, Biological; Models, Statistical; Monte Carlo Method; Phantoms, Imaging; Positron-Emission Tomography - instrumentation; Positron-Emission Tomography - methods; Reproducibility of Results; Sensitivity and Specificity; Stochastic Processes
• ISSN: 0031-9155; 1361-6560
• DOI: 10.1088/0031-9155/39/5/005

In an earlier paper we derived a theoretical formulation for estimating the statistical properties of images reconstructed using the iterative ML-EM algorithm. To gain insight into this complex problem, two levels of approximation were considered in the theory. These techniques revealed the dependence of the variance and covariance of the reconstructed image noise on the source distribution, imaging system transfer function, and iteration number. In this paper a Monte Carlo approach was taken to study the noise properties of the ML-EM algorithm and to test the predictions of the theory. The study also served to evaluate the approximations used in the theory. Simulated data from phantoms were used in the Monte Carlo experiments. The ML-EM statistical properties were calculated from sample averages of a large number of images with different noise realizations. The agreement between the more exact form of the theoretical formulation and the Monte Carlo formulation was better than 10% in most cases examined, and for many situations the agreement was within the expected error of the Monte Carlo experiments. Results from the studies provide valuable information about the noise characteristics of ML-EM reconstructed images. Furthermore, the studies demonstrate the power of the theoretical and Monte Carlo approaches for investigating noise properties of statistical reconstruction algorithms.