用于宇宙射线μ子成像的图像重建研究进展.pdfVIP

中国体视学学会cT理论与应用分会、中国工程物理研究院 2008年7月 应用电子学研究所、《CT理论与应用研究》编辑部论文集 用于宇宙射线∥子成像的图像重建研究进展 刘圆圆““，赵自然“6，陈志强。一 (清华大学a．工程物理系：b．教育部粒子技术与辐射成像重点实验室，北京100084) 摘要：宇宙射线F子检测技术具有穿透力强、对高z材料敏感等特点，特别适合检测特殊核材 料，是监控核材料走私的有效方式。图像重建问题的研究是实现该技术的一个关键难点。本文 总结了清华大学在Ⅳ子图像重建和图像评估研究工作的最新进展，包括对已有MLS和MLSD算法 的加速研究，以及对归一化平均绝对距离测量法进行改进，使其适合对Ⅳ子重建图像进行评估。 数值实验结果表明上述方法均达到预期效果，相信这些工作将对宇宙射线Ⅳ子成像技术早日投 入实际应用产生积极意义。 关键词：宇宙射线∥子成像；MLS—OSEM；MLSD--OSEM；图像评估；二值化 引言 目前已公开的可用于宇宙射线Ⅳ子成像的重建算法主要有三种，径迹重建算法(The PointofClosest Likelihood Approach，PoCA)，最大似然散射角重建算法(Maximum Likelihood and Scattering，MLS)和最大似然散射角和位移重建算法(MaximumScattering 分别有各自的特点。POCA，算法简单，容易实现，成像速度快，但是，成像质量较差；而 多；另外，由于MLSD相比MLS增加了位移的信息，成像质量要好于MLS，但是代价却是再 次增加重建耗时。所以这就促使我们去寻找一种重建用时少，成像质量高的重建算法。为 此，我们相继提出了MLS-OSEM¨。和MLSD-OSEMul重建算法。 图像评估是图像重建研究中一个十分关键的环节，通过有效的对重建图像质量进行评 估，可以评价重建算法的性能，进而指导算法的改进。但由于我们的检测目的是在较低误 判率条件下判定核材料的存在，而不是得到核材料分布的纹理特征。所以评价图像质量中 MeanAbsoluteOistanceMeasure， 最常使用的归一化平均绝对距离测量(Normalized NMADM)H。等方法不再适合评估用∥子成像技术获得的重建图像。为了解决这个问题，我们 提出了一种改进的图像评估方法，即图像评估二值化方法(Picture Comparison Binarization Method，PCBM)¨1，有效解决了这个问题。 1 MLS．OSEM和MLSD．OSEM重建算法 为了解决重建耗时和图像质量这对矛盾，我们在2007年首先提出MLS-OSEM重建算法， 并且在此基础上进一步扩展到MLSD重建算法中，得到了重建图像质量更好，重建耗时却基 基金项目：青年自然科学基金；国家自然科学基金。 全国射线数字成像与cT新技术研讨会 本相同的MLSD-OSEM重建算法。 1．1方法概述 1．1．1MLS．OSEM重建算法 假设把探测区域分成NxN个栅格，一共有M个∥子，其中的第i个穿过被探测区域的 每个栅格的路径长度为L(1≤f≤M，1≤．『≤N2)，厶代表第i个∥子通过物体的每一个栅 格的路径长度；s，表示穿过物体以后的角度偏转量；A，是第．『个栅格的散射密度(其中刀是 第，个栅格的散射密度的第一次估计值，而A?+1是第．，个栅格的散射密度的第(n+1)次估计 值)，A为散射角度的标准差。那么，似然函数可用公式(1)表示： ^‘ ^， ^ ．2 “ A=argmin∑【砌(厶五)+立i】，其中乃≥红 (1) i