Scalable Block Gibbs Sampling for Image Deblurring in X-Ray Radiography
AuthorAdams, Jesse John William
Luttman, Aaron B.
MetadataShow full item record
PublisherThe University of Arizona.
RightsCopyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction, presentation (such as public display or performance) of protected items is prohibited except with permission of the author.
AbstractQuantitative image analysis in the security sciences formulates an image deblurring problem as a Bayesian inverse problem to reduce and quantify noise and blur. We consider images of size 16 megapixels and, since each pixel represents an unknown, the dimension of the Bayesian inverse problem is on the order of 107. The large dimension poses numerical and computational difficulties for two reasons. First, Markov chain Monte Carlo (MCMC), typically used to solve a Bayesian inverse problem, is generally slow to converge in high dimensions. Second, even generating one step in a Markov chain is challenging at this size. We present a Gibbs sampler that is scalable to the large dimension required in the security sciences and its scalability is achieved in two steps. We (i) accelerate MCMC convergence by exploring banded structure in the posterior precision matrix; and (ii) use a matrix-free implementation, because constructing and storing even sparse matrices is infeasible in our target application.
Degree ProgramGraduate College