Enhanced Gibbs sampling for an application to X-ray imaging
Document Type
Presentation Abstract
Presentation Date
4-25-2016
Abstract
Image deblurring techniques derived from convolution require, a priori, an estimate for the convolution kernel or point spread function (PSF). Standard techniques for estimating the PSF involve imaging a bright point source, but this is not always feasible (e.g. high energy radiography). This work takes a novel non-parametric approach to modeling a radially symmetric PSF, in which an estimate can be obtained from the calibration image of a vertical edge. Moreover, we employ a hierarchical Bayesian model that in addition to providing a method for estimation, also gives a quantification of uncertainty in the estimate. We will present a recently developed improvement to the Gibbs algorithm for simulating samples from the Bayesian posterior of the hierarchical model, referred to as partial collapse. The improved algorithm has been independently derived in several other works, however, it has been shown that partial collapse may be improperly implemented resulting in a sampling algorithm that that no longer converges to the desired posterior. The algorithm we present is proven to satisfy invariance with respect to the target density.
Recommended Citation
Joyce, Kevin, "Enhanced Gibbs sampling for an application to X-ray imaging" (2016). Colloquia of the Department of Mathematical Sciences. 495.
https://scholarworks.umt.edu/mathcolloquia/495
Additional Details
Monday, April 25, 2016 at 3:10 p.m. in Math 103
Refreshments at 4:00 p.m. in Math Lounge 109