蒙特卡罗方法在高能闪光照相图像边界检测中应用.pdfVIP

第十一届全国蒙特卡罗方法及其应用学术交流会论文 蒙特卡罗方法在高能闪光照相图像边界检测中的 应用 张绚 ,刘军，管永红，刘进 （中国工程物理研究院 流体物理研究所，四川 绵阳 621900） 摘摘 要：要：边界退化是高能闪光照相诊断技术中影响边界信息精确检测的关键因素。本文探讨用蒙特卡罗方法来模拟高 摘摘 要要：： 能闪光照相物理过程，提供模拟图像用于边缘检测，依据模拟图像的边界退化来修正实验图像的边界，通过迭代计算 获得较精确的客体边界。利用数值计算来验证和评估方法的精度与可靠性，实验图像处理结果表明：修正后的边界值 相比于修正前精度得到了显著提升，并且该法具有较好收敛性。 关键词：关键词：闪光照相；边界退化；蒙特卡罗方法 关键词关键词：： 1 引言 边界信息是武器初级内爆过程诊断的主要内容之一[1,2]。然而，对于高能闪光照相图像，由于成像 过程的非理想特性，比如光源尺寸、源强分布、散射、探测器模糊和噪声等物理因素的影响，图像质 [3] 量发生退化 。图像质量的退化将会导致采用经典的边缘检测方法（依据图像灰度极值或灰度梯度极 值）[4]确定的边界量与理想的边界位置存在偏移。因此，为了实现边界的精确测量，边界偏移量是一 个需要充分考虑的问题。 图像复原是解决上述问题的一条途径。但是，对于闪光图像而言，引起图像降级的因素带有一定 随机性，很难对这些因素有充分认识。因此，对闪光图像彻底复原很困难。从实验角度可以采用根据 大量静态闪光照相实验的边界检测结果获得不同实验条件下的边界修正量。这种方法为了获得比较准 确的边界修正量需要以大量的实验为基础，提高了实验成本，不同的照相客体、不同的成像参数以及 边界检测方法均对边界修正量有一定的影响。 [5] 本文考虑采用蒙特卡罗方法模拟实际实验 ，通过对模拟图像的边界检测迭代修正边界偏移量，称 之为基于蒙特卡罗模拟照相过程的边界修正方法（以下简称 MC边界修正方法）。该方法可以针对特定 图像设定相应的实验参数，通过蒙特卡罗随机模拟试验精确地将焦斑尺寸、源强分布和散射等因素的 影响考虑进来。 2 照相物理模型 下面介绍高能闪光照相的物理模型。图 1 是高能闪光照相的示意图。其主要由光源、客体、防护 器件和探测器组成。高能闪光照相的基本过程是：经加速器加速的高能电子撞击高 Z 材料（即金属靶）， 电子与高Z 原子中的原子核发生非弹性碰撞而产生轫致辐射，即 X 射线源。从转换靶发射出来的源光 子会按一定的几率直接穿透照相器件和照相客体而到达探测平面，这种 X 射线光子称为直穿光子。同 时，从转换靶发射出来的部分源光子不可避免地会与防护器件和照相客体发生相互作用而被吸收或散 射，其中产生的散射光子同样会到达探测平面。转换屏将接收到的 X 射线转换成可见光，而后利用 CCD 将光信号转换成易于处理的电信号。直穿光子包含客体信息，而散射光子与客体信息没有多大关 系，这样探测器接收的信号中由于包含了散射光子不能真实反映客体的几何特性和物理特性，从而影 响了信息反演的精度。 3 MC边界处理方法原理 高能闪光照相过程涉及的物理因素较多，可采用以下退化模型来描述： ( ,gx )y = f ([X, x )y (⊗,hx )y +( X, x )]y (⊗ ,hx )y ( +, xn)y (1) d s s I 其中，xd 表示理想图像；g 表示退化后的图像；f 为探测器的响应函数，即屏片系统的 HD 曲线或 CCD 系统的G-X曲线；h 为光源的点扩展函数；x 为成像平面上的散射照射量；h 为探测器的点扩展 s s I 函数；n 为噪声。 MC 边界处理方法流程图如图 2 所示。利用图像处理方法提取图像的初始边界信息，利用蒙特卡罗 收稿日期：