核輻射劑量場(chǎng)實(shí)時(shí)成像測(cè)量系統(tǒng)的研究
掃描二維碼
隨時(shí)隨地手機(jī)看文章
本文利用閃爍體的吸收發(fā)光特點(diǎn)，并結(jié)合計(jì)算機(jī)斷層掃描技術(shù)，提出了對(duì)核輻射劑量場(chǎng)分布進(jìn)行實(shí)時(shí)成像測(cè)量的新方法——陣列式吸收發(fā)光CT法，研制出閃爍光纖陣列構(gòu)成的核探測(cè)器及其伺服控系統(tǒng).采用高靈敏度的電荷耦合器件(CCD)拾取探測(cè)器產(chǎn)生的微弱閃爍光信號(hào)，并用定點(diǎn)采集的方法對(duì)視頻信號(hào)進(jìn)行數(shù)據(jù)的快速采集.在圖像重建方法上，提出了迭代濾波反投影重建方法和利用非完全投影進(jìn)行數(shù)據(jù)修復(fù)，對(duì)獲取的投影數(shù)據(jù)進(jìn)行變換和處理.通過(guò)選擇濾波函數(shù)及其適當(dāng)?shù)膮?shù)，獲得最佳濾波效果，以重建劑量場(chǎng)的二維場(chǎng)分布，從而建立了核輻照劑量場(chǎng)的實(shí)時(shí)、高精度的成像測(cè)量系統(tǒng).實(shí)驗(yàn)及模擬實(shí)驗(yàn)結(jié)果均很好地證實(shí)了，該測(cè)量方法的可行性和測(cè)量系統(tǒng)工作的可靠性.
關(guān)鍵詞:核輻射劑量場(chǎ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
一、引 言
輻射劑量學(xué)是試圖探討射線能量的傳遞及生物組織對(duì)其能量的吸收，并用實(shí)驗(yàn)的方法測(cè)定輻射量值.從早期使用X射線起，人們就開(kāi)始采用感光膠片進(jìn)行劑量測(cè)定，以后發(fā)展了量熱劑量學(xué)、化學(xué)劑量學(xué)以及利用熱釋光現(xiàn)象的劑量測(cè)定技術(shù).近年來(lái)，又研制出電離室探測(cè)器和微型半導(dǎo)體探測(cè)器.但從目前國(guó)內(nèi)外所采用的劑量測(cè)量方法和儀器來(lái)看，各有許多不足之處，如都不能進(jìn)行劑量場(chǎng)的強(qiáng)度分布及形態(tài)的實(shí)時(shí)測(cè)量；一些方法還需要另外的數(shù)據(jù)讀出設(shè)備，致使無(wú)法現(xiàn)場(chǎng)獲取結(jié)果；大多也僅能進(jìn)行點(diǎn)測(cè)量，且材料的一致性差；除感光膠片法外，其他空間分辨都較差.因此，現(xiàn)在尚無(wú)能夠比較好地適合用于劑量場(chǎng)的強(qiáng)度分布及形態(tài)的實(shí)時(shí)測(cè)量系統(tǒng).對(duì)于劑量場(chǎng)的測(cè)量，包括束流的總強(qiáng)度及其橫截面上的強(qiáng)度分布及形態(tài)的實(shí)時(shí)成像測(cè)量，目前在國(guó)內(nèi)外沿屬空白(未見(jiàn)報(bào)道).而對(duì)刀、x刀、醫(yī)學(xué)射線加速器等這些對(duì)射線有一定匯聚要求的輻照源的射線匯聚束斑空間形狀的強(qiáng)度分布的實(shí)時(shí)成像測(cè)量，具有特別重要和迫切的實(shí)用價(jià)值.
二、基本原理
對(duì)各種核輻劑量射場(chǎng)的探測(cè)，其原理跟它們與物質(zhì)的相互作用是密切相關(guān)的.我們所提出的陣列式吸收發(fā)光CT法是一種將閃爍體發(fā)光特性與圖像重建技術(shù)相結(jié)合的核探測(cè)方法［1］.該方法依據(jù)射線與物質(zhì)相互作用的機(jī)制，及閃爍體自身的發(fā)光特點(diǎn)進(jìn)行能量變換，將劑量場(chǎng)上各點(diǎn)的強(qiáng)度轉(zhuǎn)換成與其成線性關(guān)系的閃爍光強(qiáng)度.借鑒計(jì)算機(jī)斷層掃描技術(shù)的基本思想，設(shè)計(jì)出陣列式探測(cè)器.對(duì)探測(cè)平面各象素點(diǎn)上形成的閃爍光，分別沿軸線方向線積分后接收，采用特定的測(cè)量方式，可獲得不同方向、不同位置的完備的投影數(shù)據(jù).最后利用圖像重建技術(shù)，可有效地實(shí)現(xiàn)對(duì)劑量場(chǎng)進(jìn)行實(shí)時(shí)成像測(cè)量.
圖像重建的主要方法一般有：直接反投影法、傅立葉變換重建法、卷積反投影重建法、代數(shù)迭代法等.目前所采用的各種重建方法，都還存在著一些不足的地方，而算法對(duì)重建圖像的質(zhì)量與速度起著關(guān)鍵的作用.對(duì)不同的目標(biāo)應(yīng)用不同的算法與之相適應(yīng).在卷積反投影算法中，選擇不同的卷積函數(shù)，對(duì)重建圖像的質(zhì)量影響是很大的，需根據(jù)不同的情況，作相應(yīng)的調(diào)整［2］.本文對(duì)卷積反投影重建法進(jìn)行延伸，使其不僅能在空間域進(jìn)行卷積處理，而且能方便地選擇適當(dāng)?shù)臑V波函數(shù)和參數(shù)在頻域進(jìn)行頻譜修正，達(dá)到最佳的處理效果，從而使重建精度和空間分辨率都得到進(jìn)一步的提高.綜合卷積反投影法和代數(shù)迭代法的長(zhǎng)處，我們提出了一種迭代濾波反投影法，可更好地實(shí)現(xiàn)圖像重建［1，3］.迭代濾波反投影重建法是一種迭代優(yōu)化的過(guò)程：在每次迭代運(yùn)算中，首先根據(jù)上次的重建結(jié)果，依次在每個(gè)投影方向上計(jì)算重建圖像的投影，再同實(shí)測(cè)的投影數(shù)據(jù)相比較，將差值再濾波反投影在圖像上，以修正重建結(jié)果，即完成一次迭代運(yùn)算，并將該次的運(yùn)算結(jié)果作為下一次迭代的初值.重復(fù)上述過(guò)程，直到投影誤差總和小于給定的閾值或設(shè)定的迭代數(shù)，從而結(jié)束重建過(guò)程.
另外，在實(shí)際應(yīng)用中可根據(jù)需要考慮利用非完全投影重建法來(lái)進(jìn)行適當(dāng)?shù)臄?shù)據(jù)處理.假設(shè)在整個(gè)(s,θ)平面上，投影函數(shù)P(s,θ)是解析的，即使有部分?jǐn)?shù)據(jù)丟失，可根據(jù)其解析特性，將所需要的數(shù)據(jù)有效地估算出.由于閃爍光纖直徑的限制，影響了空間分辨率的進(jìn)一步提高.采樣的投影數(shù)據(jù)對(duì)于s變量是離散的，但也應(yīng)注意到投影數(shù)據(jù)是投影方向上各點(diǎn)數(shù)值的積分，其隱含著該方向上各點(diǎn)數(shù)據(jù)的連續(xù)性.故采用非完全投影重建法對(duì)一些點(diǎn)的投影數(shù)據(jù)進(jìn)行估計(jì)和補(bǔ)齊，補(bǔ)齊缺少的數(shù)據(jù)一般必須滿足三個(gè)條件［1，4］：
(1)在數(shù)據(jù)缺少區(qū)域，根據(jù)已測(cè)數(shù)據(jù)的變化規(guī)律，使欲補(bǔ)齊的數(shù)據(jù)與其保持連續(xù)、光滑，由被測(cè)場(chǎng)的解析特性決定.
(2)對(duì)于各方向的投影數(shù)據(jù)，保持其積分的相等.
(3)保持所有投影積分的相等性，這是雷當(dāng)變換所要求的.
三、系統(tǒng)構(gòu)成
整個(gè)系統(tǒng)分為四個(gè)部分：前端探測(cè)系統(tǒng)、機(jī)械旋轉(zhuǎn)掃描系統(tǒng)、定點(diǎn)數(shù)據(jù)采集系統(tǒng)、數(shù)據(jù)處理及圖像重建系統(tǒng)，見(jiàn)圖1.
|
圖1 系統(tǒng)總體框圖 1.前端探測(cè)系統(tǒng) |
|
圖2 閃爍光纖陣列構(gòu)成的核探測(cè)器及其與光接收器連接[!--empirenews.page--] 設(shè)計(jì)的前端探測(cè)系統(tǒng)是由多根閃爍光纖構(gòu)成的，其產(chǎn)生的閃爍光通過(guò)光導(dǎo)光纖耦合至CCD進(jìn)行數(shù)據(jù)讀取.因此，由于各根光纖性能的不一致，端面處理及反光特性的差異，傳輸效率和光耦合效率的不同，以及可能受到的損傷而引起性能的改變等等，必將會(huì)致使相同的輸入，有不相同的輸出響應(yīng).同時(shí)，還有光學(xué)成像系統(tǒng)的光損失及CCD光敏元的不均勻性等.為此，構(gòu)造的陣列式前端探測(cè)器在實(shí)際應(yīng)用中一般還須進(jìn)行坪場(chǎng)修正.所謂坪場(chǎng)修正，就是對(duì)敏根光纖在CCD上獲取的輸出響應(yīng)數(shù)據(jù)乘以一修正因子，使它們各自的綜合性能保持一致，即具有相同的場(chǎng)強(qiáng)與電信號(hào)的轉(zhuǎn)換特性. |
|
圖3 步進(jìn)電機(jī)控制電路的原理框圖 3.定點(diǎn)數(shù)據(jù)采集系統(tǒng) |
|
圖4 定點(diǎn)采集系統(tǒng)結(jié)構(gòu)框圖 定點(diǎn)控制是根據(jù)監(jiān)視器屏幕二維空間上某點(diǎn)的位置，確定與其相對(duì)應(yīng)的一維視頻信號(hào)中該點(diǎn)的時(shí)刻.定點(diǎn)控制電路的原理框圖如圖5所示. |
|
圖5 定點(diǎn)控制電路原理框圖[!--empirenews.page--] 定點(diǎn)控制的硬件部分設(shè)計(jì)為一塊PC機(jī)的插件，其通過(guò)I/O總線與微機(jī)相連，采用并行方式交換數(shù)據(jù)和信息.系統(tǒng)在開(kāi)始采集時(shí)，首先由主機(jī)給出控制信號(hào)，打開(kāi)視頻同步信號(hào)的控制門(mén)，由場(chǎng)同步信號(hào)對(duì)行脈沖計(jì)數(shù)器(計(jì)數(shù)器一)和列脈沖計(jì)數(shù)器(計(jì)數(shù)器二)清零，并同時(shí)啟動(dòng)計(jì)數(shù)器一，開(kāi)始計(jì)數(shù).在行脈沖主數(shù)達(dá)到比較器一的預(yù)置值時(shí)，產(chǎn)生一級(jí)控制信號(hào).一級(jí)控制信號(hào)對(duì)列脈沖計(jì)數(shù)器(計(jì)數(shù)器二)清零，隨后啟動(dòng)計(jì)數(shù)器二.計(jì)數(shù)器二通過(guò)對(duì)10MHZ晶體振蕩器產(chǎn)生脈沖的計(jì)數(shù)，可以將每一行視頻信號(hào)劃分為520個(gè)圖像點(diǎn).當(dāng)計(jì)數(shù)器達(dá)到比較器二的預(yù)置值時(shí)，產(chǎn)生二級(jí)控制信號(hào).二級(jí)控制信號(hào)一方面啟動(dòng)ADC，另一方面產(chǎn)生一個(gè)計(jì)算機(jī)中斷服務(wù)，該服務(wù)將此時(shí)所采集的數(shù)據(jù)寫(xiě)入緩沖區(qū).同時(shí)，二級(jí)控制信號(hào)將計(jì)數(shù)器二清零. 四、實(shí)驗(yàn)結(jié)果與誤差分析 |
|
圖6 模擬劑量場(chǎng)時(shí)所用的鉛磚 圖7和圖8為將一號(hào)鉛磚置于劑量場(chǎng)中，探測(cè)器對(duì)其模擬的場(chǎng)強(qiáng)分布進(jìn)行數(shù)據(jù)測(cè)量，重建的三維圖形.其中：圖7為對(duì)測(cè)量的投影數(shù)據(jù)未經(jīng)坪場(chǎng)修正，圖10則為經(jīng)過(guò)坪場(chǎng)修正后的處理結(jié)果.圖9為將三號(hào)鉛磚置于一號(hào)鉛磚之上，對(duì)所測(cè)量的數(shù)據(jù)(經(jīng)過(guò)坪場(chǎng)修正)，進(jìn)行重建后該劑量場(chǎng)強(qiáng)度分布的三維圖形.圖10為將二號(hào)鉛磚置于劑量場(chǎng)中，探測(cè)器在其下面進(jìn)行數(shù)據(jù)測(cè)量，并對(duì)測(cè)量的數(shù)據(jù)經(jīng)過(guò)坪場(chǎng)修正后，重建該劑量場(chǎng)強(qiáng)度分布的三維圖形.圖11和圖12為將三號(hào)鉛磚置于二號(hào)鉛磚之上，探測(cè)器對(duì)構(gòu)造的劑量場(chǎng)進(jìn)行數(shù)據(jù)測(cè)量，所重建該劑量場(chǎng)強(qiáng)度分布的三維圖形.其中圖11未經(jīng)坪場(chǎng)修正，圖12則經(jīng)過(guò)坪場(chǎng)修正. |
|
圖7 劑量場(chǎng)重建圖像一 |
圖8 劑量場(chǎng)重建圖像二 |
圖9 劑量場(chǎng)重建圖像三 |
|
圖10 劑量場(chǎng)重建圖像四 |
圖11 劑量場(chǎng)重建圖像五 |
圖12 劑量場(chǎng)重建圖像六 |
| 影響測(cè)量系統(tǒng)精度的主要因素有：前端探測(cè)系統(tǒng)的隨機(jī)誤差；不同濾波函數(shù)對(duì)重建圖像質(zhì)量的影響；閃爍光纖芯直徑大小對(duì)重建精度的影響；數(shù)據(jù)采樣速率所產(chǎn)生的影響；探測(cè)器旋轉(zhuǎn)中心偏移產(chǎn)生的影響. 對(duì)于本文所設(shè)計(jì)的陣列式閃爍光纖探測(cè)器(有效探測(cè)區(qū)域100mm×100mm)，將其放置所構(gòu)造的劑量場(chǎng)中進(jìn)行實(shí)時(shí)成像測(cè)量，根據(jù)理論推導(dǎo)和實(shí)驗(yàn)結(jié)果的數(shù)據(jù)分析，可估算可能導(dǎo)致的各種誤差［9］，以綜合評(píng)估系統(tǒng)的性能.若輸入的投影數(shù)據(jù)，其相對(duì)誤差不超過(guò)±3%，則模擬實(shí)驗(yàn)和計(jì)算表明，重建誤差可控制在3%左右.對(duì)于直徑為1mm的光纖，獲取的投影數(shù)據(jù)平均相對(duì)誤差經(jīng)折算約為0.5%，重建平均相對(duì)誤差約為0.4%.投影方向數(shù)一般小于7個(gè)時(shí)，則完全不能重建.當(dāng)方向數(shù)增加，則重建圖像誤差逐漸減小.采樣頻率的選取同樣應(yīng)滿足Niquist定理，否則，會(huì)影響圖像重建精度及空間分辨.對(duì)于一般的劑量場(chǎng)分布，若不考慮各種其他因素的影響，當(dāng)投影方向數(shù)為60，采樣間隔等于1mm時(shí)，圖像重建精度是非常高的.通過(guò)模擬運(yùn)算，其重建場(chǎng)平均相對(duì)誤差非常小，約百分之零點(diǎn)幾.中心偏移對(duì)重建圖像質(zhì)量的影響十分大.在制作陣列式閃爍光纖探測(cè)器時(shí)，一定要切實(shí)注意精確地確定其旋轉(zhuǎn)中心位置，否則，會(huì)產(chǎn)生很大誤差，甚至導(dǎo)致變形.若中心偏差控制在不超過(guò)0.1個(gè)象素點(diǎn)，則產(chǎn)生的重建誤差可控制在2.0%以內(nèi). 綜合各種因素，本系統(tǒng)的成像測(cè)量的總體平均相對(duì)誤差可控制在5%以內(nèi)，空間分辨率不低于1mm. |
下載文檔
