核輻射劑量場實時成像測量系統(tǒng)的研究

本文利用閃爍體的吸收發(fā)光特點，并結合計算機斷層掃描技術，提出了對核輻射劑量場分布進行實時成像測量的新方法——陣列式吸收發(fā)光CT法，研制出閃爍光纖陣列構成的核探測器及其伺服控系統(tǒng).采用高靈敏度的電荷耦合器件(CCD)拾取探測器產生的微弱閃爍光信號，并用定點采集的方法對視頻信號進行數據的快速采集.在圖像重建方法上，提出了迭代濾波反投影重建方法和利用非完全投影進行數據修復，對獲取的投影數據進行變換和處理.通過選擇濾波函數及其適當的參數，獲得最佳濾波效果，以重建劑量場的二維場分布，從而建立了核輻照劑量場的實時、高精度的成像測量系統(tǒng).實驗及模擬實驗結果均很好地證實了，該測量方法的可行性和測量系統(tǒng)工作的可靠性.
　　關鍵詞:核輻射劑量場；閃爍體；斷層掃描；迭代濾波反投影；圖像重建

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射線起，人們就開始采用感光膠片進行劑量測定，以后發(fā)展了量熱劑量學、化學劑量學以及利用熱釋光現象的劑量測定技術.近年來，又研制出電離室探測器和微型半導體探測器.但從目前國內外所采用的劑量測量方法和儀器來看，各有許多不足之處，如都不能進行劑量場的強度分布及形態(tài)的實時測量；一些方法還需要另外的數據讀出設備，致使無法現場獲取結果；大多也僅能進行點測量，且材料的一致性差；除感光膠片法外，其他空間分辨都較差.因此，現在尚無能夠比較好地適合用于劑量場的強度分布及形態(tài)的實時測量系統(tǒng).對于劑量場的測量，包括束流的總強度及其橫截面上的強度分布及形態(tài)的實時成像測量，目前在國內外沿屬空白(未見報道).而對刀、x刀、醫(yī)學射線加速器等這些對射線有一定匯聚要求的輻照源的射線匯聚束斑空間形狀的強度分布的實時成像測量，具有特別重要和迫切的實用價值.

二、基本原理
　　對各種核輻劑量射場的探測，其原理跟它們與物質的相互作用是密切相關的.我們所提出的陣列式吸收發(fā)光CT法是一種將閃爍體發(fā)光特性與圖像重建技術相結合的核探測方法［1］.該方法依據射線與物質相互作用的機制，及閃爍體自身的發(fā)光特點進行能量變換，將劑量場上各點的強度轉換成與其成線性關系的閃爍光強度.借鑒計算機斷層掃描技術的基本思想，設計出陣列式探測器.對探測平面各象素點上形成的閃爍光，分別沿軸線方向線積分后接收，采用特定的測量方式，可獲得不同方向、不同位置的完備的投影數據.最后利用圖像重建技術，可有效地實現對劑量場進行實時成像測量.
　　圖像重建的主要方法一般有：直接反投影法、傅立葉變換重建法、卷積反投影重建法、代數迭代法等.目前所采用的各種重建方法，都還存在著一些不足的地方，而算法對重建圖像的質量與速度起著關鍵的作用.對不同的目標應用不同的算法與之相適應.在卷積反投影算法中，選擇不同的卷積函數，對重建圖像的質量影響是很大的，需根據不同的情況，作相應的調整［2］.本文對卷積反投影重建法進行延伸，使其不僅能在空間域進行卷積處理，而且能方便地選擇適當的濾波函數和參數在頻域進行頻譜修正，達到最佳的處理效果，從而使重建精度和空間分辨率都得到進一步的提高.綜合卷積反投影法和代數迭代法的長處，我們提出了一種迭代濾波反投影法，可更好地實現圖像重建［1，3］.迭代濾波反投影重建法是一種迭代優(yōu)化的過程：在每次迭代運算中，首先根據上次的重建結果，依次在每個投影方向上計算重建圖像的投影，再同實測的投影數據相比較，將差值再濾波反投影在圖像上，以修正重建結果，即完成一次迭代運算，并將該次的運算結果作為下一次迭代的初值.重復上述過程，直到投影誤差總和小于給定的閾值或設定的迭代數，從而結束重建過程.
　　另外，在實際應用中可根據需要考慮利用非完全投影重建法來進行適當的數據處理.假設在整個(s,θ)平面上，投影函數P(s,θ)是解析的，即使有部分數據丟失，可根據其解析特性，將所需要的數據有效地估算出.由于閃爍光纖直徑的限制，影響了空間分辨率的進一步提高.采樣的投影數據對于s變量是離散的，但也應注意到投影數據是投影方向上各點數值的積分，其隱含著該方向上各點數據的連續(xù)性.故采用非完全投影重建法對一些點的投影數據進行估計和補齊，補齊缺少的數據一般必須滿足三個條件［1，4］：
　　(1)在數據缺少區(qū)域，根據已測數據的變化規(guī)律，使欲補齊的數據與其保持連續(xù)、光滑，由被測場的解析特性決定.
　　(2)對于各方向的投影數據，保持其積分的相等.
　　(3)保持所有投影積分的相等性，這是雷當變換所要求的.

三、系統(tǒng)構成
　　整個系統(tǒng)分為四個部分：前端探測系統(tǒng)、機械旋轉掃描系統(tǒng)、定點數據采集系統(tǒng)、數據處理及圖像重建系統(tǒng)，見圖1.