Library Subscription: Guest
Begell Digital Portal Begell Digital Library eBooks Journals References & Proceedings Research Collections
International Journal for Uncertainty Quantification

Impact factor: 1.000

ISSN Print: 2152-5080
ISSN Online: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2016013678
pages 35-55


Johnathan M. Bardsley
Department of Mathematical Sciences, The University of Montana, Missoula, Montana 59812-0864, USA
Aaron Luttman
Signal Processing and Data Analysis, National Security Technologies, LLC, P.O. Box 98521, M/S NLV078, Las Vegas, Nevada, 89193-8521, USA


Poisson noise models arise in a wide range of linear inverse problems in imaging. In the Bayesian setting, the Poisson likelihood function together with a Gaussian prior yields a posterior density that is not of a well-known form and is thus difficult to sample from, especially for large-scale problems. In this work, we present a method for computing samples from posterior density functions with Poisson likelihood and Gaussian prior, using a Gaussian approximation of the posterior as an independence proposal within a Metropolis−Hastings framework. We consider a class of Gaussian priors, some of which are edge-preserving, and which we motivate using Markov random fields. We present two sampling algorithms: one which samples the unknown image alone, leaving the prior scaling (or regularization) parameter alone, and another which samples both the unknown image and the prior scaling parameter. For this paper, we make the assumption that our unknown image is sufficiently positive that proposed samples are always positive, allowing us to ignore the nonnegativity constraint. Results are demonstrated on synthetic data−including a synthetic X-ray radiograph generated from a radiation transport code−and on real images used to calibrate a pulsed power high-energy X-ray source at the U.S. Department of Energy's Nevada National Security Site.