电子发烧友App
核辐射剂量场实时成像测量系统的研究
本文利用闪烁体的吸收发光特点，并结合计算机断层扫描技术，提出了对核辐射剂量场分布进行实时成像测量的新方法——阵列式吸收发光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.
|
图1 系统总体框图 1.前端探测系统 |
|
图2 闪烁光纤阵列构成的核探测器及其与光接收器连接 设计的前端探测系统是由多根闪烁光纤构成的，其产生的闪烁光通过光导光纤耦合至CCD进行数据读取.因此，由于各根光纤性能的不一致，端面处理及反光特性的差异，传输效率和光耦合效率的不同，以及可能受到的损伤而引起性能的改变等等，必将会致使相同的输入，有不相同的输出响应.同时，还有光学成像系统的光损失及CCD光敏元的不均匀性等.为此，构造的阵列式前端探测器在实际应用中一般还须进行坪场修正.所谓坪场修正，就是对敏根光纤在CCD上获取的输出响应数据乘以一修正因子，使它们各自的综合性能保持一致，即具有相同的场强与电信号的转换特性. |
|
图3 步进电机控制电路的原理框图 3.定点数据采集系统 |
|
图4 定点采集系统结构框图 定点控制是根据监视器屏幕二维空间上某点的位置，确定与其相对应的一维视频信号中该点的时刻.定点控制电路的原理框图如图5所示. |
|
图5 定点控制电路原理框图 定点控制的硬件部分设计为一块PC机的插件，其通过I/O总线与微机相连，采用并行方式交换数据和信息.系统在开始采集时，首先由主机给出控制信号，打开视频同步信号的控制门，由场同步信号对行脉冲计数器(计数器一)和列脉冲计数器(计数器二)清零，并同时启动计数器一，开始计数.在行脉冲主数达到比较器一的预置值时，产生一级控制信号.一级控制信号对列脉冲计数器(计数器二)清零，随后启动计数器二.计数器二通过对10MHZ晶体振荡器产生脉冲的计数，可以将每一行视频信号划分为520个图像点.当计数器达到比较器二的预置值时，产生二级控制信号.二级控制信号一方面启动ADC，另一方面产生一个计算机中断服务，该服务将此时所采集的数据写入缓冲区.同时，二级控制信号将计数器二清零. 四、实验结果与误差分析 |
|
图6 模拟剂量场时所用的铅砖 图7和图8为将一号铅砖置于剂量场中，探测器对其模拟的场强分布进行数据测量，重建的三维图形.其中：图7为对测量的投影数据未经坪场修正，图10则为经过坪场修正后的处理结果.图9为将三号铅砖置于一号铅砖之上，对所测量的数据(经过坪场修正)，进行重建后该剂量场强度分布的三维图形.图10为将二号铅砖置于剂量场中，探测器在其下面进行数据测量，并对测量的数据经过坪场修正后，重建该剂量场强度分布的三维图形.图11和图12为将三号铅砖置于二号铅砖之上，探测器对构造的剂量场进行数据测量，所重建该剂量场强度分布的三维图形.其中图11未经坪场修正，图12则经过坪场修正. |
|
图7 剂量场重建图像一 |
图8 剂量场重建图像二 |
图9 剂量场重建图像三 |
|
图10 剂量场重建图像四 |
图11 剂量场重建图像五 |
图12 剂量场重建图像六 |
| 影响测量系统精度的主要因素有：前端探测系统的随机误差；不同滤波函数对重建图像质量的影响；闪烁光纤芯直径大小对重建精度的影响；数据采样速率所产生的影响；探测器旋转中心偏移产生的影响. 对于本文所设计的阵列式闪烁光纤探测器(有效探测区域100mm×100mm)，将其放置所构造的剂量场中进行实时成像测量，根据理论推导和实验结果的数据分析，可估算可能导致的各种误差［9］，以综合评估系统的性能.若输入的投影数据，其相对误差不超过±3%，则模拟实验和计算表明，重建误差可控制在3%左右.对于直径为1mm的光纤，获取的投影数据平均相对误差经折算约为0.5%，重建平均相对误差约为0.4%.投影方向数一般小于7个时，则完全不能重建.当方向数增加，则重建图像误差逐渐减小.采样频率的选取同样应满足Niquist定理，否则，会影响图像重建精度及空间分辨.对于一般的剂量场分布，若不考虑各种其他因素的影响，当投影方向数为60，采样间隔等于1mm时，图像重建精度是非常高的.通过模拟运算，其重建场平均相对误差非常小，约百分之零点几.中心偏移对重建图像质量的影响十分大.在制作阵列式闪烁光纤探测器时，一定要切实注意精确地确定其旋转中心位置，否则，会产生很大误差，甚至导致变形.若中心偏差控制在不超过0.1个象素点，则产生的重建误差可控制在2.0%以内. 综合各种因素，本系统的成像测量的总体平均相对误差可控制在5%以内，空间分辨率不低于1mm. 国家自然科学基金资助项目 |
工商网监
评论