Reversible jump Markov chain Monte Carlo for Bayesian deconvolution of point sources

Guillaume Stawinski; Arnaud Doucet; Patrick Duvaut

doi:10.1117/12.323798

22 September 1998 Reversible jump Markov chain Monte Carlo for Bayesian deconvolution of point sources

Guillaume Stawinski, Arnaud Doucet, Patrick Duvaut

Proceedings Volume 3459, Bayesian Inference for Inverse Problems; (1998) https://doi.org/10.1117/12.323798
Event: SPIE's International Symposium on Optical Science, Engineering, and Instrumentation, 1998, San Diego, CA, United States

Abstract

In this article, we address the problem of Bayesian deconvolution of point sources in nuclear imaging under the assumption of Poissonian statistics. The observed image is the result of the convolution by a known point spread function of an unknown number of point sources with unknown parameters. To detect the number of sources and estimate their parameters we follow a Bayesian approach. However, instead of using a classical low level prior model based on Markov random fields, we prose a high-level model which describes the picture as a list of its constituent objects, rather than as a list of pixels on which the data are recorded. More precisely, each source is assumed to have a circular Gaussian shape and we set a prior distribution on the number of sources, on their locations and on the amplitude and width deviation of the Gaussian shape. This high-level model has far less parameters than a Markov random field model as only s small number of sources are usually present. The Bayesian model being defined, all inference is based on the resulting posterior distribution. This distribution does not admit any closed-form analytical expression. We present here a Reversible Jump MCMC algorithm for its estimation. This algorithm is tested on both synthetic and real data.

Citation Download Citation

Guillaume Stawinski, Arnaud Doucet, and Patrick Duvaut "Reversible jump Markov chain Monte Carlo for Bayesian deconvolution of point sources", Proc. SPIE 3459, Bayesian Inference for Inverse Problems, (22 September 1998); https://doi.org/10.1117/12.323798

ACCESS THE FULL ARTICLE

INSTITUTIONAL
Select your institution to access the SPIE Digital Library.

SELECT YOUR INSTITUTION

PERSONAL
Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.

PERSONAL SIGN IN

No SPIE Account? Create one

PURCHASE THIS CONTENT

SUBSCRIBE TO DIGITAL LIBRARY

50 downloads per 1-year subscription

Members: $195

Non-members: $335 ADD TO CART

25 downloads per 1 - year subscription

Members: $145

Non-members: $250 ADD TO CART

PURCHASE SINGLE ARTICLE

Includes PDF, HTML & Video, when available