大孔径厚针孔成像数值模拟研究
Numerical Simulation of Large Thick Aperture Imaging
-
摘要: 为解决cm级尺寸射线源的空间分布诊断,本文提出大孔径厚针孔成像方法。采用Geant4软件编写蒙特卡罗程序模拟5、10和15 mm大孔径厚针孔的成像过程,采用Matlab软件实现逆滤波、Wiener滤波和Lucy-Richardson复原算法,以均方根误差σ作为复原效果评判标准,优化了Wiener滤波信噪比参数和Lucy-Richardson算法迭代次数参数。将复原算法应用于不同孔径厚针孔的图像复原,发现理想逆滤波算法会引起噪声放大;Wiener滤波算法图像复原解存在较严重的伪影;Lucy-Richardson算法可获得清晰的复原图像。结合Lucy-Richardson算法的5、10、15 mm大孔径厚针孔成像获取的图像的σ是0.5 mm孔径传统针孔成像结果的1.15、1.21、1.26倍。结果表明,结合图像复原算法的大孔径厚针孔成像技术能获得与传统针孔成像接近的源图像测量结果,初步验证了该方法具有可行性。Abstract: Large thick aperture imaging method was proposed to measure the radiation intensity distribution of radiation source whose size is several centimeters. The new method contains two steps, which are coded imaging and image reconstruction. Geant4 was used to simulate the coded imaging process. Matlab was used to simulate the image reconstruction process using inverse filtering, Wiener filtering and Lucy-Richardson method. The input parameters of Wiener filtering and Lucy-Richardson method were optimized according to the root mean square error. They were applied to reconstruct the coded imaging of 5, 10, and 15 mm aperture respectively. Inverse filtering causes noise amplification, and Wiener filtering produces artifacts. Lucy-Richardson method can give out a clear restoration image. The root mean square error of the reconstruction image using Lucy-Richardson method is 15% larger than that using traditional pin-hole imaging method. It implies that large thick aperture imaging method can get a similar result with traditional pin-hole imaging. The new method is proved to be effective.