The ML-EM algorithm in continuum: sparse measure solutions
Abstract
Linear inverse problems with Poisson noise and non-negative unknown are ubiquitous in applications, for instance in Positron Emission Tomography (PET) in medical imaging. The associated maximum likelihood problem is routinely solved using an expectation-maximisation algorithm (ML-EM). This typically results in images which look spiky, even with early stopping. We give an explanation for this phenomenon. We first regard the image as a measure. We prove that if the measurements are not in the cone , which is typical of short exposure times, likelihood maximisers as well as ML-EM cluster points must be sparse, i.e., typically a sum of point masses. On the other hand, in the long exposure regime, we prove that cluster points of ML-EM will be measures without singular part. Finally, we provide concentration bounds for the probability to be in the sparse case.
Cite
@article{arxiv.1909.01966,
title = {The ML-EM algorithm in continuum: sparse measure solutions},
author = {Camille Pouchol and Olivier Verdier},
journal= {arXiv preprint arXiv:1909.01966},
year = {2020}
}