核辐射剂量场分布进行实时成像测量的新方法—阵列式吸收发光CT法

本文利用闪烁体的吸收发光特点，并结合计算机断层扫描技术，提出了对核辐射剂量场分布进行实时成像测量的新方法——阵列式吸收发光CT法，研制出闪烁光纤阵列构成的核探测器及其伺服控系统.采用高灵敏度的电荷耦合器件(CCD)拾取探测器产生的微弱闪烁光信号，并用定点采集的方法对视频信号进行数据的快速采集.在图像重建方法上，提出了迭代滤波反投影重建方法和利用非完全投影进行数据修复，对获取的投影数据进行变换和处理.通过选择滤波函数及其适当的参数，获得最佳滤波效果，以重建剂量场的二维场分布，从而建立了核辐照剂量场的实时、高精度的成像测量系统.实验及模拟实验结果均很好地证实了，该测量方法的可行性和测量系统工作的可靠性.
　　The Study on Real-Time Radiography System in Dose Field

SHENG Jin-hua
(China Academy of Telecommunications Technology,Beijing 100083,China)
YIN Ze-jie
(Department of Modern Physics,University of Science and Technology of China,Hefei 230027,China)

Abstract:By the comprehensive analysis of existing mcasuring methods,a new method is proposed in the paper to measure the nuclear radiation dose field with array absorption-emission computed tomography.The detector based on scintillation optical fiber array and its servo-control system are designed.The weak optical signal is measured by Charge Coupled Devices (CCD) and the formed visual signal is quickly sampled with the method of testing fixed signed points.The reconstruction method of iterative filter back projection and the data renovation based on noncomplete projection are proposed to reconstruct dose field much more accurately.The filter function is also an important factor in the image reconstruction.The Butterworth filter function is chosen.By adjusting its parameters,we have obtained optimum result.Finally,all the experimentation show that the methods proposed in paper are proper and whole measuring system is advanced and feasible.
　　Key words:nuclear radiation dose field;scintillation;computed tomography;iterative filter back projection;image reconstruction

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

二、基本原理
　　对各种核辐剂量射场的探测，其原理跟它们与物质的相互作用是密切相关的.我们所提出的阵列式吸收发光CT法是一种将闪烁体发光特性与图像重建技术相结合的核探测方法［1］.该方法依据射线与物质相互作用的机制，及闪烁体自身的发光特点进行能量变换，将剂量场上各点的强度转换成与其成线性关系的闪烁光强度.借鉴计算机断层扫描技术的基本思想，设计出阵列式探测器.对探测平面各象素点上形成的闪烁光，分别沿轴线方向线积分后接收，采用特定的测量方式，可获得不同方向、不同位置的完备的投影数据.最后利用图像重建技术，可有效地实现对剂量场进行实时成像测量.
　　图像重建的主要方法一般有：直接反投影法、傅立叶变换重建法、卷积反投影重建法、代数迭代法等.目前所采用的各种重建方法，都还存在着一些不足的地方，而算法对重建图像的质量与速度起着关键的作用.对不同的目标应用不同的算法与之相适应.在卷积反投影算法中，选择不同的卷积函数，对重建图像的质量影响是很大的，需根据不同的情况，作相应的调整［2］.本文对卷积反投影重建法进行延伸，使其不仅能在空间域进行卷积处理，而且能方便地选择适当的滤波函数和参数在频域进行频谱修正，达到最佳的处理效果，从而使重建精度和空间分辨率都得到进一步的提高.综合卷积反投影法和代数迭代法的长处，我们提出了一种迭代滤波反投影法，可更好地实现图像重建［1，3］.迭代滤波反投影重建法是一种迭代优化的过程：在每次迭代运算中，首先根据上次的重建结果，依次在每个投影方向上计算重建图像的投影，再同实测的投影数据相比较，将差值再滤波反投影在图像上，以修正重建结果，即完成一次迭代运算，并将该次的运算结果作为下一次迭代的初值.重复上述过程，直到投影误差总和小于给定的阈值或设定的迭代数，从而结束重建过程.
　　另外，在实际应用中可根据需要考虑利用非完全投影重建法来进行适当的数据处理.假设在整个(s,θ)平面上，投影函数P(s,θ)是解析的，即使有部分数据丢失，可根据其解析特性，将所需要的数据有效地估算出.由于闪烁光纤直径的限制，影响了空间分辨率的进一步提高.采样的投影数据对于s变量是离散的，但也应注意到投影数据是投影方向上各点数值的积分，其隐含着该方向上各点数据的连续性.故采用非完全投影重建法对一些点的投影数据进行估计和补齐，补齐缺少的数据一般必须满足三个条件［1，4］：
　　(1)在数据缺少区域，根据已测数据的变化规律，使欲补齐的数据与其保持连续、光滑，由被测场的解析特性决定.
　　(2)对于各方向的投影数据，保持其积分的相等.
　　(3)保持所有投影积分的相等性，这是雷当变换所要求的.