ML-EM in continuum: sparse measure solutions - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Pré-Publication, Document De Travail Année : 2019

ML-EM in continuum: sparse measure solutions

Résumé

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.
Fichier principal
Vignette du fichier
MLEM_Arxiv.pdf (819.38 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Loading...

Dates et versions

hal-02192480 , version 1 (23-07-2019)
hal-02192480 , version 2 (12-08-2019)

Identifiants

  • HAL Id : hal-02192480 , version 1

Citer

Camille Pouchol, Olivier Verdier. ML-EM in continuum: sparse measure solutions. 2019. ⟨hal-02192480v1⟩
80 Consultations
230 Téléchargements

Partager

Gmail Facebook X LinkedIn More