Point spread function estimation in X-ray imaging with partially collapsed Gibbs sampling
From MaRDI portal
Abstract: The point spread function (PSF) of a translation invariant imaging system is its impulse response, which cannot always be measured directly. This is the case in high energy X-ray radiography, and it must be estimated from images of calibration objects indirectly related to the impulse response. When the PSF is assumed to have radial symmetry, it can be estimated from an image of an opaque straight edge. We use a non-parametric Bayesian approach, where the prior probability density for the PSF is modeled as a Gaussian Markov random field and radial symmetry is incorporated in a novel way. Markov Chain Monte Carlo posterior estimation is carried out by adapting a recently developed improvement to the Gibbs sampling algorithm, referred to as partially collapsed Gibbs sampling. Moreover, the algorithm we present is proven to satisfy invariance with respect to the target density. Finally, we demonstrate the efficacy of these methods on radiographic data obtained from a high-energy X-ray diagnostic system at the U.S. Department of Energy's Nevada National Security Site.
Recommendations
- Sampling-based uncertainty quantification in deconvolution of X-ray radiographs
- SURE-type functionals as criteria for parametric PSF estimation
- Point-spread function identification in an image restoration problem
- The non-parametric sub-pixel local point spread function estimation is a well posed problem
- Subpixel point spread function estimation from two photographs at different distances
Cites work
- A METROPOLIS-HASTINGS METHOD FOR LINEAR INVERSE PROBLEMS WITH POISSON LIKELIHOOD AND GAUSSIAN PRIOR
- A primer on space-time modeling from a Bayesian perspective
- A simple method for generating gamma variables
- Analysis of the Gibbs Sampler for Hierarchical Inverse Problems
- Bayesian Abel inversion in quantitative X-ray radiography
- Bayesian data analysis.
- Computational Methods for Inverse Problems
- Deblurring Images
- Discrete inverse problems. Insight and algorithms.
- Fast sampling in a linear-Gaussian inverse problem
- Gaussian Markov Random Fields
- scientific article; zbMATH DE number 49594 (Why is no real title available?)
- scientific article; zbMATH DE number 1350773 (Why is no real title available?)
- scientific article; zbMATH DE number 1124118 (Why is no real title available?)
- Introduction to Bayesian Scientific Computing
- Inverse problems: a Bayesian perspective
- Markov chain Monte Carlo. Stochastic simulation for Bayesian inference.
- MCMC-based image reconstruction with uncertainty quantification
- Numerical Solution of Partial Differential Equations
- Partially Collapsed Gibbs Samplers
- Sampling-based uncertainty quantification in deconvolution of X-ray radiographs
- Statistical and computational inverse problems.
Cited in
(7)- Numerical study of point spread function of a fast neutron radiography system based on scintillating-fiber array
- Low-rank independence samplers in hierarchical Bayesian inverse problems
- Sampling-based uncertainty quantification in deconvolution of X-ray radiographs
- A blocking scheme for dimension-robust Gibbs sampling in large-scale image deblurring
- Efficient Marginalization-Based MCMC Methods for Hierarchical Bayesian Inverse Problems
- Improving PSF calibration in confocal microscopic imaging -- estimating and exploiting bilateral symmetry
- An approach to characterizing spatial aspects of image system blur
This page was built for publication: Point spread function estimation in X-ray imaging with partially collapsed Gibbs sampling
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4641603)
