Applications of MCMC method based on Bayesian hierarchical model in flash radiography reconstruction

Wang Zhongmiao; Liu Jun; Jing Yuefeng; Liu Jin; Guan Yonghong

doi:10.11884/HPLPB201830.180123

Volume 30 Issue 11

Nov. 2018

Turn off MathJax

Article Contents

Article Navigation > High Power Laser and Particle Beams > 2018 > 30(11): 114004

Wang Zhongmiao, Liu Jun, Jing Yuefeng, et al. Applications of MCMC method based on Bayesian hierarchical model in flash radiography reconstruction[J]. High Power Laser and Particle Beams, 2018, 30: 114004. doi: 10.11884/HPLPB201830.180123

Citation:

Wang Zhongmiao, Liu Jun, Jing Yuefeng, et al. Applications of MCMC method based on Bayesian hierarchical model in flash radiography reconstruction[J]. High Power Laser and Particle Beams, 2018, 30: 114004. doi: 10.11884/HPLPB201830.180123

Citation:

PDF( 2042 KB)

Applications of MCMC method based on Bayesian hierarchical model in flash radiography reconstruction

doi: 10.11884/HPLPB201830.180123

1.
Institute of Fluid Physics, CAEP, P. O. Box 919-105, Mianyang 621900, China
2.
Graduate School of CAEP, Mianyang 621900, China

Received Date: 2018-04-25
Rev Recd Date: 2018-07-25
Publish Date: 2018-11-15

Abstract

Abstract

The Markov chain Monte Carlo(MCMC) method combined with Bayesian theory can not only use prior information flexibly, but also give the uncertainty of solution. There is a bright application prospect in quantitative diagnosis of flash radiography. For the optical path length data of the flash radiographic images, a posterior probability model is built by Bayesian hierarchical model, and the Markov chain is dynamically constructed by Gibbs sampling. Then the statistical results of linear attenuation coefficients and their uncertainty are obtained and compared with the constrained conjugate gradient (CCG) method. The results of numerical experiments show that the reconstruction result of MCMC method is approximately the same as the true data for the ideal FTO optical path length image. In the case of blurring and noise, the reconstructed result is equivalent to that by the CCG method. Even when the blurred optical path length data is interfered by high noise, the result of MCMC is slightly better than that of CCG. More importantly, the uncertainty of the reconstruction can be provided by MCMC method. The related work in this paper verifies the feasibility of MCMC reconstruction of flash radiographic images and lays a good foundation for MCMC reconstruction with blurred and noised transmissivity images.
- MCMC method,
- flash radiography,
- image reconstruction,
- linear inversion,
- uncertainty

FullText(HTML)

References(14)

References

[1]	Wong A, Mishra A, Zhang W, et al. Stochastic image denoising based on Markov-chain Monte Carlo sampling[J]. Signal Processing, 2011, 91(8): 2112-2120. doi: 10.1016/j.sigpro.2011.03.021
[2]	Fadili M J, Starck J L. Sparse representation-based image deconvolution by iterative thresholding[C]//Proceedings of Astronomical Data Analysis Ada. 2006.
[3]	Cho M, Shin Y M, Lee K M. Co-recognition of image pairs by data-driven monte carlo image exploration[C]//European Conference on Computer Vision, 2009: 144-157.
[4]	Tian J, Ma K K. A MCMC approach for Bayesian super-resolution image reconstruction[C]//IEEE International Conference on Image Processing. 2005: I-45-8.
[5]	Watson S A, Balzer S, Gossein C. Density measurement errors at DARHT quantifying a decade of progress[R]. Defense Science Quarterly, 2009.
[6]	Gelman A, Carlin J B, Stern H S, et al. Bayesian data analysis[M]. 3rd Ed. Boca Raton: C Press, 2012.
[7]	Bardsley J M. MCMC-based image reconstruction with uncertainty quantification[J]. Siam Journal on Scientific Computing, 2012, 34(3): A1316-A1332. doi: 10.1137/11085760X
[8]	Howard M, Luttman A, Fowler M. Sampling-based uncertainty quantification in deconvolution of X-ray radiographs[J]. Journal of Computational & Applied Mathematics, 2014, 270(270): 43-51.
[9]	Bishop C M, Nasrabadi N M. Pattern recognition and machine learning[M]. Springer, 2006: 461-462.
[10]	Rue H, Held L. Gaussian Markov random fields: Theory and applications[J]. Technometrics, 2005, 48(1): 146-147.
[11]	Howard M, Fowler M, Luttman A, et al. Bayesian Abel inversion in quantitative X-ray radiography[J]. Siam Journal on Scientific Computing, 2016, 38(3): B396-B413. doi: 10.1137/15M1018721
[12]	景越峰, 管永红, 张小琳. 基于约束优化的闪光照相图像重建算法[J]. 强激光与粒子束, 2016, 28: 094002. https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY201609021.htm Jing Yuefeng, Guan Yonghong, Zhang Xiaolin. Constrained optimization reconstruction for flash radiographic image. High Power Laser and Particle Beams, 2016, 28: 094002 https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY201609021.htm
[13]	景越峰, 刘军, 管永红. 改进的约束共轭梯度闪光照相图像重建算法[J]. 强激光与粒子束, 2011, 23(8): 2201-2204. https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY201108048.htm Jing Yuefeng, Liu Jun, Guan Yonghong. Improved constrained conjugate gradient reconstruction algorithm for flash radiographic image. High Power Laser and Particle Beams, 2011, 23(8): 2201-2204 https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY201108048.htm
[14]	Odyniec M, Luttman A, Howardb M. Maximum likelihood estimation and uncertainty quantification for signals with Poisson-Gaussian mixed noise[R]. Site-Directed Research and Development, 2015: 211.