Bayesian inference and uncertainty quantification for medical image reconstruction with Poisson data

Qingping Zhou, Tengchao Yu, Xiaoqun Zhang, Jinglai Li

Research output: Contribution to journalArticlepeer-review

3 Citations (Scopus)
297 Downloads (Pure)

Abstract

We provide a complete framework for performing infinite dimensional Bayesian inference and uncertainty quantification for image reconstruction with Poisson data. In particular, we address the following issues to make the Bayesian framework applicable in practice. We first introduce a positivity-preserving reparametrization, and we prove that under the reparametrization and a hybrid prior, the posterior distribution is well-posed in the infinite dimensional setting. Second, we provide a dimension-independent Markov chain Monte Carlo algorithm, based on the preconditioned Crank--Nicolson Langevin method, in which we use a primal-dual scheme to compute the offset direction. Third, we give a method combining the model discrepancy method and maximum likelihood estimation to determine the regularization parameter in the hybrid prior. Finally we propose to use the obtained posterior distribution to detect artifacts in a recovered image. We provide an example to demonstrate the effectiveness of the proposed method.
Original languageEnglish
Pages (from-to)29–52
Number of pages24
JournalSIAM Journal on Imaging Sciences
Volume13
Issue number1
DOIs
Publication statusPublished - 7 Jan 2020

Keywords

  • Poisson distribution
  • Bayesian inference
  • image reconstruction
  • uncertainty quantification
  • Markov chain Monte Carlo
  • positron emission tomography

Fingerprint

Dive into the research topics of 'Bayesian inference and uncertainty quantification for medical image reconstruction with Poisson data'. Together they form a unique fingerprint.

Cite this