Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

The ML-EM algorithm in continuum: sparse measure solutions

Abstract : Linear inverse problems $A \mu = y$ with Poisson noise and non-negative unknown $\mu \geq 0$ are ubiquitous in applications where light is involved, such as Positron Emission Tomography (PET) in medical imaging. The associated maximum likelihood problem is most often solved by the celebrated iterative ML-EM algorithm, which is known to yield results which look spiky even when stopped early. This work provides an explanation for this phenomenon by going into continuum in image space when the operator $A$ is integral, such as in PET. This is done by considering the image $\mu$ as a measure rather than an element of a Lebesgue space. We prove that if the measurements $y$ are not in the cone $\{A \mu, \mu \geq 0\}$, which is typical of high noise or short exposure times, likelihood maximisers as well as ML-EM cluster points must be sparse, i.e., typically a sum of Dirac masses. In the better case where $y$ is in the cone, we prove that cluster points of ML-EM will be measures without singular part, as long as the initial guess is smooth. Finally, we provide concentration bounds for the probability to be in the undesirable sparse case, proving that it is exponentially small with exposure time, with a rate controlled by the distance of the noiseless data to the boundary of the cone.
Document type :
Preprints, Working Papers, ...
Complete list of metadata

Cited literature [32 references]  Display  Hide  Download
Contributor : Camille Pouchol Connect in order to contact the contributor
Submitted on : Monday, August 12, 2019 - 3:55:05 PM
Last modification on : Thursday, March 5, 2020 - 3:31:33 PM


Files produced by the author(s)


  • HAL Id : hal-02192480, version 2



Camille Pouchol, Olivier Verdier. The ML-EM algorithm in continuum: sparse measure solutions. 2019. ⟨hal-02192480v2⟩



Record views


Files downloads