核辐射剂量场实时成像测量系统的研究

一、引言

　　辐射剂量学是试图探讨射线能量的传递及生物组织对其能量的吸收，并用实验的方法测定辐射量值.从早期使用X射线起，人们就开始采用感光胶片进行剂量测定，以后发展了量热剂量学、化学剂量学以及利用热释光现象的剂量测定技术.近年来，又研制出电离室探测器和微型半导体探测器.但从目前国内外所采用的剂量测量方法和仪器来看，各有许多不足之处，如都不能进行剂量场的强度分布及形态的实时测量；一些方法还需要另外的数据读出设备，致使无法现场获取结果；大多也仅能进行点测量，且材料的一致性差；除感光胶片法外，其他空间分辨都较差.因此，现在尚无能够比较好地适合用于剂量场的强度分布及形态的实时测量系统. 对于剂量场的测量，包括束流的总强度及其横截面上的强度分布及形态的实时成像测量，目前在国内外沿属空白(未见报道).而对ts32-1.gif (83 bytes)刀、x刀、医学射线加速器等这些对射线有一定汇聚要求的辐照源的射线汇聚束斑空间形状的强度分布的实时成像测量，具有特别重要和迫切的实用价值.

　　二、基本原理

　　对各种核辐剂量射场的探测，其原理跟它们与物质的相互作用是密切相关的.我们所提出的阵列式吸收发光CT法是一种将闪烁体发光特性与图像重建技术相结合的核探测方法［1］.该方法依据射线与物质相互作用的机制，及闪烁体自身的发光特点进行能量变换，将剂量场上各点的强度转换成与其成线性关系的闪烁光强度.借鉴计算机断层扫描技术的基本思想，设计出阵列式探测器.对探测平面各象素点上形成的闪烁光，分别沿轴线方向线积分后接收，采用特定的测量方式，可获得不同方向、不同位置的完备的投影数据.最后利用图像重建技术，可有效地实现对剂量场进行实时成像测量.

　　图像重建的主要方法一般有：直接反投影法、傅立叶变换重建法、卷积反投影重建法、代数迭代法等.目前所采用的各种重建方法，都还存在着一些不足的地方，而算法对重建图像的质量与速度起着关键的作用.对不同的目标应用不同的算法与之相适应.在卷积反投影算法中，选择不同的卷积函数，对重建图像的质量影响是很大的，需根据不同的情况，作相应的调整［2］.本文对卷积反投影重建法进行延伸，使其不仅能在空间域进行卷积处理，而且能方便地选择适当的滤波函数和参数在频域进行频谱修正，达到最佳的处理效果，从而使重建精度和空间分辨率都得到进一步的提高.综合卷积反投影法和代数迭代法的长处，我们提出了一种迭代滤波反投影法，可更好地实现图像重建［1，3］.迭代滤波反投影重建法是一种迭代优化的过程：在每次迭代运算中，首先根据上次的重建结果，依次在每个投影方向上计算重建图像的投影，再同实测的投影数据相比较，将差值再滤波反投影在图像上，以修正重建结果，即完成一次迭代运算，并将该次的运算结果作为下一次迭代的初值.重复上述过程，直到投影误差总和小于给定的阈值或设定的迭代数，从而结束重建过程.

　　另外，在实际应用中可根据需要考虑利用非完全投影重建法来进行适当的数据处理.假设在整个(s,θ)平面上，投影函数P(s,θ)是解析的，即使有部分数据丢失，可根据其解析特性，将所需要的数据有效地估算出.由于闪烁光纤直径的限制，影响了空间分辨率的进一步提高.采样的投影数据对于s变量是离散的，但也应注意到投影数据是投影方向上各点数值的积分，其隐含着该方向上各点数据的连续性.故采用非完全投影重建法对一些点的投影数据进行估计和补齐，补齐缺少的数据一般必须满足三个条件［1，4］：

　　(1)在数据缺少区域，根据已测数据的变化规律，使欲补齐的数据与其保持连续、光滑，由被测场的解析特性决定.

　　(2)对于各方向的投影数据，保持其积分的相等.

　　(3)保持所有投影积分的相等性，这是雷当变换所要求的.

　　三、系统构成

　　整个系统分为四个部分：前端探测系统、机械旋转扫描系统、定点数据采集系统、数据处理及图像重建系统，见图1.

图1　系统总体框图

　　1.前端探测系统

　　由纤芯是闪烁材料构成的光纤能适应E＞5kev的X射线、γ射线及其它射线的辐照探测.但直到现在，获得的主要研究成果是涉及在高能粒子物理中的应用［5，6］.闪烁光纤对带电粒子比x射线和γ射线灵敏，这是由于光纤纤芯的直径较x射线或γ射线与其作用产生次级电子的有效射程相比太小，一般仅很少部分能量沉积在光纤纤芯中，以产生闪烁光.而在所涉及测量的能量范围内，主要作用机制是康普顿效应，这是由于构成光纤的材料是低Z所决定的，因此光电效应和电子对效应都相对较弱［7］.