ResearchTrend.AI
  • Papers
  • Communities
  • Events
  • Blog
  • Pricing
Papers
Communities
Social Events
Terms and Conditions
Pricing
Parameter LabParameter LabTwitterGitHubLinkedInBlueskyYoutube

© 2025 ResearchTrend.AI, All rights reserved.

  1. Home
  2. Papers
  3. 1909.01966
19
7
v1v2 (latest)

The ML-EM algorithm in continuum: sparse measure solutions

4 September 2019
Camille Pouchol
Olivier Verdier
ArXiv (abs)PDFHTML
Abstract

Linear inverse problems Aμ=δA \mu = \deltaAμ=δ with Poisson noise and non-negative unknown μ≥0\mu \geq 0μ≥0 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 μ\muμ as a measure. We prove that if the measurements δ\deltaδ are not in the cone {Aμ,μ≥0}\{A \mu, \mu \geq 0\}{Aμ,μ≥0}, 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.

View on arXiv
Comments on this paper