Доступ предоставлен для: Guest
Портал Begell Электронная Бибилиотека e-Книги Журналы Справочники и Сборники статей Коллекции
International Journal for Uncertainty Quantification
Импакт фактор: 3.259 5-летний Импакт фактор: 2.547 SJR: 0.531 SNIP: 0.8 CiteScore™: 1.52

ISSN Печать: 2152-5080
ISSN Онлайн: 2152-5099

Свободный доступ

International Journal for Uncertainty Quantification

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

A METROPOLIS-HASTINGS METHOD FOR LINEAR INVERSE PROBLEMS WITH POISSON LIKELIHOOD AND GAUSSIAN PRIOR

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.


Articles with similar content:

A MULTI-INDEX MARKOV CHAIN MONTE CARLO METHOD
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 1
Ajay Jasra, Yan Zhou, Kengo Kamatani, Kody J. H. Law
DATA ASSIMILATION FOR NAVIER-STOKES USING THE LEAST-SQUARES FINITE-ELEMENT METHOD
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 5
Richard P. Dwight, Alexander Schwarz
VARIABLE-SEPARATION BASED ITERATIVE ENSEMBLE SMOOTHER FOR BAYESIAN INVERSE PROBLEMS IN ANOMALOUS DIFFUSION REACTION MODELS
International Journal for Uncertainty Quantification, Vol.9, 2019, issue 3
Yuming Ba, Na Ou, Lijian Jiang
A NOVEL HYBRID APPROACH FOR SIMPLIFIED NEUTROSOPHIC DECISION-MAKING WITH COMPLETELY UNKNOWN WEIGHT INFORMATION
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 2
R?dvan Şahin, Gökçe Dilek Küçük
ALGORITHM FOR PICTURE FUZZY MULTIPLE ATTRIBUTE DECISION-MAKING BASED ON NEW DISTANCE MEASURE
International Journal for Uncertainty Quantification, Vol.7, 2017, issue 2
Jingguo Dai, Xindong Peng