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

　　一、引　　言

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

　　二、基本原理

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

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