Maximum-likelihood expectation-maximization reconstruction of sinograms with arbitrary noise distribution using NEC-transformations

The maximum-likelihood (ML) expectation-maximization (EM) [ML-EM] algorithm is being widely used for image reconstruction in positron emission tomography. The algorithm is strictly valid if the data are Poisson distributed. However, it is also often applied to processed sinograms that do not meet this requirement. This may sometimes lead to suboptimal results: streak artifacts appear and the algorithm converges toward a lower likelihood value. As a remedy, the authors propose two simple pixel-by-pixel methods [noise equivalent counts (NEC)-scaling and NEC-shifting] in order to transform arbitrary sinogram noise into noise which is approximately Poisson distributed (the first and second moments of the distribution match those of the Poisson distribution). The convergence speed associated with both transformation methods is compared, and the NEC-scaling method is validated with both simulations and clinical data. These new methods extend the ML-EM algorithm to a general purpose nonnegative reconstruction algorithm.