核辐射剂量场分布进行实时成像测量的新方法―阵列式吸收发光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)保持所有投影积分的相等性，这是雷当变换所要求的.

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

　设计的前端探测系统是由多根闪烁光纤构成的，其产生的闪烁光通过光导光纤耦合至CCD进行数据读取.因此，由于各根光纤性能的不一致，端面处理及反光特性的差异，传输效率和光耦合效率的不同，以及可能受到的损伤而引起性能的改变等等，必将会致使相同的输入，有不相同的输出响应.同时，还有光学成像系统的光损失及CCD光敏元的不均匀性等.为此，构造的阵列式前端探测器在实际应用中一般还须进行坪场修正.所谓坪场修正，就是对敏根光纤在CCD上获取的输出响应数据乘以一修正因子，使它们各自的综合性能保持一致，即具有相同的场强与电信号的转换特性.
　　2.机械旋转扫描系统
　　我们提出的用阵列式吸收发光CT法探测剂量场强度分布的构想.为此需要设计运动机架以带动探测器在180°范围内进行等角度旋转扫描.根据机械系统的要求，可选择步进电机作为驱动部件.因为步进电机特点是定位精度高，无累积误差，因此被广泛应用于开环数控系统.设计的步进电机控制电路采用集成模块结构，与微机直接相联，能同时控制两组步进电机，其功能强、响应速度快，可靠性高.原理框图见图3.

图3　步进电机控制电路的原理框图

　　3.定点数据采集系统
　　因在实际应用中仅需对感兴趣相对应传送投影数据的数量较少的均匀光斑进行采集，故可借鉴通常静态图像慢速数据采集系统所采用的方法，提出了定点数据采集方法.定点采集系统通过对视频同步信号的计数控制，产生A/D变换器的启动信号，采集相应时刻的视频数据信号后送入计算机进行处理.系统的硬件主要分为视频信号定点控制和数据变换采集两大部分，其结构框图见图4.

图4　定点采集系统结构框图

　　定点控制是根据监视器屏幕二维空间上某点的位置，确定与其相对应的一维视频信号中该点的时刻.定点控制电路的原理框图如图5所示.

影响测量系统精度的主要因素有：前端探测系统的随机误差；不同滤波函数对重建图像质量的影响；闪烁光纤芯直径大小对重建精度的影响；数据采样速率所产生的影响；探测器旋转中心偏移产生的影响.
　　对于本文所设计的阵列式闪烁光纤探测器(有效探测区域100mm×100mm)，将其放置所构造的剂量场中进行实时成像测量，根据理论推导和实验结果的数据分析，可估算可能导致的各种误差［9］，以综合评估系统的性能.若输入的投影数据，其相对误差不超过±3%，则模拟实验和计算表明，重建误差可控制在3%左右.对于直径为1mm的光纤，获取的投影数据平均相对误差经折算约为0.5%，重建平均相对误差约为0.4%.投影方向数一般小于7个时，则完全不能重建.当方向数增加，则重建图像误差逐渐减小.采样频率的选取同样应满足Niquist定理，否则，会影响图像重建精度及空间分辨.对于一般的剂量场分布，若不考虑各种其他因素的影响，当投影方向数为60，采样间隔等于1mm时，图像重建精度是非常高的.通过模拟运算，其重建场平均相对误差非常小，约百分之零点几.中心偏移对重建图像质量的影响十分大.在制作阵列式闪烁光纤探测器时，一定要切实注意精确地确定其旋转中心位置，否则，会产生很大误差，甚至导致变形.若中心偏差控制在不超过0.1个象素点，则产生的重建误差可控制在2.0%以内.
　　综合各种因素，本系统的成像测量的总体平均相对误差可控制在5%以内，空间分辨率不低于1mm.

图1　系统总体框图

　　1.前端探测系统
　　由纤芯是闪烁材料构成的光纤能适应E＞5kev的X射线、γ射线及其它射线的辐照探测.但直到现在，获得的主要研究成果是涉及在高能粒子物理中的应用［5，6］.闪烁光纤对带电粒子比x射线和γ射线灵敏，这是由于光纤纤芯的直径较x射线或γ射线与其作用产生次级电子的有效射程相比太小，一般仅很少部分能量沉积在光纤纤芯中，以产生闪烁光.而在所涉及测量的能量范围内，主要作用机制是康普顿效应，这是由于构成光纤的材料是低Z所决定的，因此光电效应和电子对效应都相对较弱［7］.
　　根据γ刀及其它剂量场和其与闪烁体相互作用的特点，我们提出了阵列式吸收发光CT测量方法，并据此构造前端探测器.探测器设计为：由若干个一定长度的特种闪烁光纤水平紧密放置构成一平面光纤阵列，其一端端面覆盖反射层，以提高其输出光响应，另一端可耦合至光接收器(CCD).所设计的阵列式闪烁探测器与光接收器CCD，通过光导光纤连一成体，并将其加固，构成前端探测系统.探测器在剂量场中，将所吸收的辐射能转换成光能，经线积分后，再通过光导光纤引出.并在光接收器的光敏区形成了按一定间隔排列的光束，从而将剂量场强度信号转换成视频电信号.
　　由塑料闪烁光纤阵列构成的探测器，具有如下特点［8］
　　(1)较短的衰减时间(即无长余辉)，约2～3ns.
　　(2)性能稳定.探测器是有辐射损伤的，但经实际测量在102GY辐照量以下探测器受到的损伤不甚明显.
　　(3)光传输性能好.光衰减长度可达500cm.
　　(4)结构简单、使用寿命长等.