Sampling-based Uncertainty Quantification in Deconvolution of X-ray Radiographs
- NSTec
In imaging applications that focus on quantitative analysis{such as X-ray radiography in the security sciences--it is necessary to be able to reliably estimate the uncertainties in the processing algorithms applied to the image data, and deconvolving the system blur out of the image is usually an essential step. In this work we solve the deconvolution problem within a Bayesian framework for edge-enhancing reconstruction with uncertainty quantification. The likelihood is a normal approximation to the Poisson likelihood, and the prior is generated from a classical total variation regularized Poisson deconvolution. Samples from the corresponding posterior distribution are computed using a Markov chain Monte Carlo approach, giving a pointwise measure of uncertainty in the final, deconvolved signal. We demonstrate the results on real data used to calibrate a high-energy X-ray source and show that this approach gives reconstructions as good as classical regularization methods, while mitigating many of their drawbacks.
- Research Organization:
- Nevada Test Site (NTS), Mercury, NV (United States)
- Sponsoring Organization:
- USDOE National Nuclear Security Administration (NNSA)
- DOE Contract Number:
- DE-AC52-06NA25946
- OSTI ID:
- 1134211
- Report Number(s):
- DOE/NV/25946-1904
- Journal Information:
- Journal of Computational and Applied Mathematics, Vol. 270; Conference: Fourth International Conference on Finite Element Methods in Engineering and Sciences (FEMTEC 2013)
- Country of Publication:
- United States
- Language:
- English
Similar Records
Bayesian Abel inversion in quantitative X-ray radiography
A stochastic approach to quantifying the blur with uncertainty estimation for high-energy X-ray imaging systems