Inscrição na biblioteca: Guest
Portal Digital Begell Biblioteca digital da Begell eBooks Diários Referências e Anais Coleções de pesquisa
International Journal for Uncertainty Quantification
Fator do impacto: 3.259 FI de cinco anos: 2.547 SJR: 0.417 SNIP: 0.8 CiteScore™: 1.52

ISSN Imprimir: 2152-5080
ISSN On-line: 2152-5099

Open Access

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

RESUMO

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:

Stochastic Model of Time-Base Signal Errors in a Decametric Communication Channel, Allowing for their Fractal Properties
Telecommunications and Radio Engineering, Vol.68, 2009, issue 1
S. A. Volkov
MODEL OF ENERGY CONSUMPTION IN THE ENTERPRISE
ENERGETIKA, Vol.2, 2005, issue 1
Piotras Bachanovas
Increasing the Efficiency of Error Correction by the Reed-Solomon Interference-Resistant Codes in Digital Telecommunication Channels
Telecommunications and Radio Engineering, Vol.66, 2007, issue 13
S. I. Egorov, G. Markaryan
FORWARD AND INVERSE UNCERTAINTY QUANTIFICATION USING MULTILEVEL MONTE CARLO ALGORITHMS FOR AN ELLIPTIC NONLOCAL EQUATION
International Journal for Uncertainty Quantification, Vol.6, 2016, issue 6
Ajay Jasra, Yan Zhou, Kody J. H. Law
OPTIMAL SENSOR PLACEMENT FOR THE ESTIMATION OF TURBULENCE MODEL PARAMETERS IN CFD
International Journal for Uncertainty Quantification, Vol.5, 2015, issue 6
Costas Papadimitriou, Dimitrios I. Papadimitriou