Regularization parameter selection and an efficient algorithm for total variation-regularized positron emission tomography

Research output: Contribution to journalArticlepeer-review

Abstract

In positron emission tomography, image data corresponds to measurements of emitted photons from a radioactive tracer in the subject. Such count data is typically modeled using a Poisson random variable, leading to the use of the negative-log Poisson likelihood fit-to-data function. Regularization is needed, however, in order to guarantee reconstructions with minimal artifacts. Given that tracer densities are primarily smoothly varying, but also contain sharp jumps (or edges), total variation regularization is a natural choice. However, the resulting computational problem is quite challenging. In this paper, we present an efficient computational method for this problem. Convergence of the method has been shown for quadratic regularization functions and here convergence is shown for total variation regularization. We also present three regularization parameter choice methods for use on total variation-regularized negative-log Poisson likelihood problems. We test the computational and regularization parameter selection methods on two synthetic data sets.

Original languageEnglish
Pages (from-to)255-271
Number of pages17
JournalNumerical Algorithms
Volume57
Issue number2
DOIs
StatePublished - Jun 2011

Keywords

  • Inverse problems
  • Positron emission tomography
  • Regularization parameter selection
  • Statistical methods
  • Total variation

Fingerprint

Dive into the research topics of 'Regularization parameter selection and an efficient algorithm for total variation-regularized positron emission tomography'. Together they form a unique fingerprint.

Cite this