Convergence of the maximum likelihood reconstruction algorithm for emission computed tomography.

Convergence properties of the maximum likelihood estimator (MLE) for emission computed tomographic (ECT) image reconstruction are evaluated as a function of Poisson noise, precision of the assumed system resolution model and iteration number up to 10,000 iterations. In the ECT reconstruction problem, the photon-emitting source distribution is to be estimated from measurements of projections of the emitted photon flux. The MLE algorithm seeks a source distribution which will maximise the maximum likelihood function relating the estimated and the measured projections. A Monte Carlo model of the system transfer function of a single photon emission computed tomographic (SPECT) system allowed realistic projection data to be simulated from a known source distribution. Poisson noise was added to the Monte Carlo simulations. By using projection data from a known source distribution generated through a known system transfer function, we were able to simultaneously evaluate the convergence of both the projection estimations as well as the source distribution estimations. As predicted by theory, the estimates of the projections did continue to improve (or remain the same) for all combinations of Poisson noise (up to 10% RMS) and system resolution (+/- 10% of true value) tested. Convergence of source distribution estimates to the true value was found for up to 10,000 iterations only for low noise (0.1% RMS) with the correct resolution function. For all other combinations, there was some optimum iteration (between 30 and 400) after which the source estimate was degraded even though the estimate of the projections was improved.