Applications of MCMC method based on Bayesian hierarchical model in flash radiography reconstruction
-
摘要: 针对闪光图像得到的光程数据,采用贝叶斯分层模型建立了后验概率模型,运用Gibbs抽样动态构造马尔可夫链;进而获得了关于线吸收系数的统计结果及其不确定度,并与约束共轭梯度(CCG)方法进行对比分析。数值实验结果表明,马尔可夫链蒙特卡罗(MCMC)方法对理想光程图像的重建结果与真值近似完全一致;在含模糊和噪声时,重建结果与CCG方法相当;当含模糊且噪声干扰较大时,MCMC方法的重建结果要略优于CCG;更重要的是MCMC方法能够给出重建结果的不确定度。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.
-
Key words:
- MCMC method /
- flash radiography /
- image reconstruction /
- linear inversion /
- uncertainty
-
表 1 重建数据均方根误差比较
Table 1. Root mean square errors by different method
RSN σ/% CCG1 CCG2 MCMC 9.50 11.86 11.65 11.20 12.6 10.73 10.48 10.12 32.3 8.22 7.32 7.69 46.7 7.93 7.00 7.39 ∞ 0 0.79 -
[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.htmJing 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.htmJing 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.