童云坤，陶 玲，李韙韜，黃海鵬，李怡燃

（南京航空航天大學(xué)自動化學(xué)院，南京210016）

0 引言

伽馬相機(jī)是用來給釋放伽馬光子的放射性元素進(jìn)行成像的設(shè)備[1]。在藥物開發(fā)與核醫(yī)學(xué)成像中，伽馬相機(jī)可以用來分析生物體中放射性標(biāo)記物的分布[2]，它不僅能反映生物系統(tǒng)的結(jié)構(gòu)特點(diǎn)[3]，還能夠顯示相應(yīng)的代謝功能信息，是一種常用的放射性核素醫(yī)學(xué)成像設(shè)備，已經(jīng)在臨床醫(yī)學(xué)中被廣泛應(yīng)用[4-5]。隨著動物模型在臨床藥物研究中的廣泛應(yīng)用，伽馬射線探測技術(shù)也從人體探測拓展至小動物的活體探測。由于小動物體內(nèi)組織的尺寸要遠(yuǎn)遠(yuǎn)小于人體內(nèi)組織的尺寸，原本適用于人體的伽馬成像系統(tǒng)的分辨力遠(yuǎn)遠(yuǎn)達(dá)不到小動物成像的需求，成像效果極不理想。而空間分辨力和檢測靈敏度呈反比，要求較高的空間分辨力就必須以犧牲相應(yīng)的靈敏度為代價。美國Cleveland Clinic 的小動物SPECT[6]系統(tǒng)采用25～36針孔數(shù)的鎢針孔加臨床雙探頭，空間分辨力可達(dá)0.7 mm。在國內(nèi)，清華大學(xué)的小動物SPECT 系統(tǒng)采用7 針孔平板鎢制針孔準(zhǔn)直器加臨床SPECT 探頭，可達(dá)到0.5 mm 的空間分辨力[7]。而SPECT 系統(tǒng)的成本構(gòu)造是傳統(tǒng)伽馬相機(jī)的數(shù)倍，同時，超高的空間高分辨力也意味著伽馬成像系統(tǒng)探測效率的降低和靈敏度的下降。因此如何平衡伽馬射線的檢測靈敏度和空間分辨力，實現(xiàn)在小視場下的高空間分辨力成像是小動物伽馬成像系統(tǒng)研制的一個最主要的挑戰(zhàn)。本文設(shè)計了一種伽馬射線檢測系統(tǒng)，通過配置定制的NaI（Tl）晶體以及相應(yīng)匹配的光電倍增管（photomultiplier tube，PMT），同時對軟件算法的圖像成幀技術(shù)方面進(jìn)行改善，從而平衡檢測靈敏度和空間分辨力性能指標(biāo)，最終對自然環(huán)境和特殊環(huán)境下伽馬射線進(jìn)行檢測，驗證對伽馬射線的測試性能，為后續(xù)研制成熟的小動物活體核素成像伽馬相機(jī)提供技術(shù)支持和實踐基礎(chǔ)。

1 伽馬射線檢測系統(tǒng)原理

目前伽馬射線檢測主要基于蓋格計數(shù)、閃爍體計數(shù)和半導(dǎo)體計數(shù)技術(shù)，國內(nèi)外在利用伽馬射線做圖像檢測時，大多采用閃爍體計數(shù)的方式。本文設(shè)計的方案也應(yīng)用了閃爍體計數(shù)的方式。整個檢測系統(tǒng)的原理如圖1 所示。系統(tǒng)主要由一塊NaI（Tl）晶體、一塊對應(yīng)的PMT 陣列、位置總和電路以及上位機(jī)組成，傳統(tǒng)小型伽馬相機(jī)的組成大多如此[8]。NaI（Tl）閃爍晶體接收外界的伽馬射線發(fā)出熒光，熒光進(jìn)入PMT，PMT 輸出電脈沖的幅值和它接收到的閃爍光的強(qiáng)度成正比，輸出的電壓經(jīng)過位置總和電路可以得到相應(yīng)的位置（X+、X-，Y+、Y-）及能量（E）信息，這些信息傳入到上位機(jī)被歸納與整理。

圖1 伽馬射線檢測系統(tǒng)原理圖

2 系統(tǒng)設(shè)計與實現(xiàn)

2.1 硬件設(shè)計方案

本文所設(shè)計的用于小動物成像系統(tǒng)的伽馬射線檢測系統(tǒng)在傳統(tǒng)伽馬相機(jī)硬件技術(shù)的基礎(chǔ)上做了一些考慮和調(diào)整[9]。就選擇NaI（Tl）晶體而言[10]，它的最大發(fā)射波長為415 nm，且在此波段沒有明顯自吸收，其發(fā)光效率也是同類晶體中最高的，因此可以和PMT 有很好的耦合。就選擇晶體的尺寸、形狀而言，如果晶體太薄，伽馬光子會穿透晶體到達(dá)PMT，雖然伽馬射線的波長太短，不會響應(yīng)PMT，但是透過晶體的這一部分伽馬光子無法在晶體內(nèi)部退激產(chǎn)生熒光，導(dǎo)致轉(zhuǎn)化的效率低下，宏觀上來看，這就造成了采集信息的時間變得更長；如果晶體太厚，雖然說絕大多數(shù)的伽馬光子會在晶體內(nèi)部退激產(chǎn)生相應(yīng)的熒光，但是可能導(dǎo)致伽馬光子在晶體內(nèi)多次作用之后才被PMT 探測，這個時候的閃爍點(diǎn)定位模糊，造成的空間分辨力下降，同時，一般NaI（Tl）晶體的能量分辨率為7%～10%，且對不同能量的分辨也不同，但是隨著晶體厚度的增加，晶體的長徑比會增加，它的能量分辨率就會下降，而最終形成的圖像是某個能量大小區(qū)間對應(yīng)的圖像，所以能量分辨率的下降會直接導(dǎo)致圖像質(zhì)量的下降。

考慮到研制的是小動物活體成像的伽馬射線檢測系統(tǒng)，因此不必選擇轉(zhuǎn)化效率很高的晶體，稍薄一些的晶體雖然轉(zhuǎn)化效率不會太高，但是已經(jīng)滿足在小鼠麻醉期間獲取相應(yīng)信息的要求，而且本文將要使用的锝同位素99mTc 釋放射線的能量大約為141 keV，屬于較低的能量[11]，它釋放的伽馬射線與NaI（Tl）晶體一般作用在表面的2～5 mm，所需的晶體厚度自然也不需要太厚。尤其關(guān)鍵的是，選取薄厚適宜的晶體可以使伽馬射線進(jìn)入晶體只經(jīng)過一次作用就產(chǎn)生熒光，由此產(chǎn)生的閃爍點(diǎn)的定位更為準(zhǔn)確，空間分辨力更高。如果使用較厚的晶體，檢測靈敏度并不會明顯改善，但是空間分辨力會急劇下降。測量人體的晶體厚度大約為8 mm，測量小動物的晶體比測量人體的晶體要更薄一些。綜上所述，本文最終定制了厚度為6 mm 的扁平狀的NaI（Tl）晶體進(jìn)行檢測，如圖2所示。

圖2 定制的NaI（Tl）閃爍晶體

本文采用的PMT 是日本濱松公司生產(chǎn)的位置靈敏PMT（R8520-00-C12），其單個有效探測面積為35 mm×35 mm，為多陽極位置靈敏PMT，其有效探測單元矩陣為10×10，增益約為0.9×109。它的噪聲小，放大倍數(shù)也較為穩(wěn)定，并且在415 nm 附近有很好的光譜響應(yīng)，R8520-00-C12 的光譜響應(yīng)如圖3 所示。探頭的固有空間分辨力取決于晶體陣列像素單元的尺寸和PMT 可探測的最小位置距離。經(jīng)過測試發(fā)現(xiàn)位置靈敏PMT 可探測的最小位置距離要小于閃爍晶體陣列像素單元，所以選取此款PMT 并不會對探頭的固有分辨力造成影響。同時本文設(shè)計在閃爍晶體和PMT 之間加入光導(dǎo)，以減少熒光在運(yùn)輸過程之間的損失。晶體、光導(dǎo)、PMT 之間涂有光學(xué)硅脂，用來排除空氣，防止系統(tǒng)的空間分辨力和檢測靈敏度的下降。對應(yīng)每個入射伽馬光子與晶體作用產(chǎn)生的熒光，PMT 分別輸出位置信號和能量信號，針對本實驗所用的多陽極位置靈敏PMT，本系統(tǒng)采用阻抗電荷分配法，即采用電阻網(wǎng)絡(luò)進(jìn)行陽極輸出信號的讀取。經(jīng)過相應(yīng)的放大電路和位置轉(zhuǎn)化電路轉(zhuǎn)化為X+、X-、Y+、Y-的4 路位置信息以及相應(yīng)的能量信息，將信息存入數(shù)據(jù)存儲器，然后通過傳輸控制協(xié)議（transmission control protocol，TCP）與上位機(jī)進(jìn)行通信。上位機(jī)對接收到的數(shù)據(jù)進(jìn)行解析、存儲、分析繪圖、實時顯示等操作。

圖3 R8520-00-C12 的光譜響應(yīng)

硬件的安裝組合方式如圖4 所示，黑色機(jī)箱部分為700 V 高壓電源和12 V 低壓電源，左側(cè)黑色方塊為NaI（Tl）晶體和PMT 陣列管，也就是探頭部分。由于PMT 受環(huán)境光的影響較大，所以需要將之密封，使其接收到的光信號全部來自NaI（Tl）吸收伽馬射線產(chǎn)生的熒光，這樣產(chǎn)生的圖像才會更為精準(zhǔn)。因為系統(tǒng)只做檢測使用，故沒有在探頭前安裝準(zhǔn)直器。若在做小動物活體成像實驗時，為實現(xiàn)系統(tǒng)對藥物的實時監(jiān)測，在保證空間分辨力的同時，也需要考慮靈敏度。盡管針孔準(zhǔn)直器擁有更高的空間分辨力，但是針孔準(zhǔn)直器視野內(nèi)不同區(qū)域的均勻性不同，無法同時兼顧到靈敏度問題，故一般直接采用高能通用平行孔準(zhǔn)直器，它的空間分辨力隨著距離的增加而降低，但是靈敏度隨著距離的增加卻沒有明顯變化，所以將實驗小動物安置得足夠近，就可以在保證空間分辨力的同時兼顧到檢測靈敏度，這樣產(chǎn)生的圖像與真實圖像的大小是一比一的。同時，為反映小動物體內(nèi)的放射性核素信息，此類平行孔準(zhǔn)直器配套的閃爍晶體探測器與常規(guī)探測器有所不同，一般會設(shè)計得比單個實驗小動物本身更大，來獲取整體動物核素信息。黑色電源箱上側(cè)是位置總和電路，采集到的數(shù)據(jù)存儲在相應(yīng)的數(shù)據(jù)存儲器當(dāng)中，等待上位機(jī)的傳輸調(diào)用，右側(cè)網(wǎng)線端口與上位機(jī)（PC）相連，組成一個完整的系統(tǒng)，兩者之間以TCP進(jìn)行交互。

圖4 硬件實物圖

2.2 軟件設(shè)計方案

上位機(jī)控制軟件是本系統(tǒng)的重要組成部分，主要包括與硬件的連接，數(shù)據(jù)的處理、顯示與存儲。本文基于LabWindows 設(shè)計了相應(yīng)的系統(tǒng)控制及顯示界面（如圖5 所示）。界面簡潔明了，主要包括一個“Canvas”控件，用來顯示一定能量范圍內(nèi)的512×512 散點(diǎn)圖；一個“Graph”控件，用來顯示伽馬射線的能譜圖，即相同能量的伽馬射線的數(shù)目，其中縱坐標(biāo)是伽馬射線的數(shù)目，橫坐標(biāo)是射線能量大小，根據(jù)自定義算法分為0～1 024 格，經(jīng)測試，絕大多數(shù)實驗所需能量都在其中。圖中可見的一條藍(lán)色和黃色的豎線是用來選取相應(yīng)的能寬，能寬中的伽馬射線是目標(biāo)射線，最終所成的圖像即是能寬中射線的位置坐標(biāo)，例如目標(biāo)元素99mT c 處于520 格位置，可選取510～530 作為能寬，而選取不同大小的能寬會對圖像的對比度、分辨力產(chǎn)生一定的影響。軟件還設(shè)計有相應(yīng)的與下位機(jī)進(jìn)行TCP 通信的連接，開始采集、停止采集等功能以及對圖像的采集時間、數(shù)目進(jìn)行標(biāo)注的功能，也包括圖像的凍結(jié)和存儲等功能。系統(tǒng)要實現(xiàn)的目標(biāo)都在軟件界面上得到了相應(yīng)的體現(xiàn)。

圖5 軟件界面

上位機(jī)和下位機(jī)的通信是通過TCP 完成的，這里采用的是基于TCP 的socket 編程，通過約定相應(yīng)的Winsock 版本，綁定相應(yīng)的地址來實現(xiàn)通信，協(xié)議指定相應(yīng)的START、GET、STOP 命令來控制數(shù)據(jù)的讀取與暫停。通過相應(yīng)的超級終端，可以觀測指定命令是否發(fā)送，來判斷上位機(jī)和下位機(jī)的通信連接是否成功。通過下位機(jī)傳輸上來的是5 組16 bit short型的數(shù)據(jù)，每一數(shù)據(jù)對應(yīng)相應(yīng)的位置和能量信息。下位機(jī)每次傳輸200 B 的數(shù)據(jù)，通過上位機(jī)指令一直循環(huán)發(fā)送。傳輸?shù)臄?shù)據(jù)中除去高6 位的校驗信息，剩下的低10 位即為傳輸?shù)挠行?shù)據(jù)，即X+、X-、Y+、Y-、E，這5 組數(shù)據(jù)值均在0～1 024，此處得到的E 即是能量的大小，而另外4 路位置信息通過重心法原理處理，得到的單位化（X，Y）坐標(biāo)在（-1，1）之間。將（X，Y）進(jìn)行拉伸、平移等操作可以得到相應(yīng)的512×512 測試圖像。

整個軟件打開之后需要點(diǎn)擊“連接”按鈕確定與下位機(jī)的TCP 通信正常。當(dāng)設(shè)備連接正確之后，點(diǎn)擊“start”即可開始實時采集相應(yīng)數(shù)據(jù)，點(diǎn)擊“stop”即為停止采集。在采集過程中，軟件還設(shè)計了采集時長和計數(shù)功能，當(dāng)計數(shù)累計到一定程度，即整個圖像畫面比較豐富或者能譜圖已經(jīng)產(chǎn)生了較為明顯的能峰時，可以選擇保存相應(yīng)的實驗數(shù)據(jù)。本系統(tǒng)對軟件的內(nèi)存和性能都進(jìn)行了相應(yīng)的優(yōu)化，且在實驗過程中運(yùn)行流暢。

3 性能測試實驗

在本系統(tǒng)搭建完成之后，為進(jìn)行性能測試，采集了2 組數(shù)據(jù)進(jìn)行對比。第一組數(shù)據(jù)是在普通的自然環(huán)境下，直接采集外界的伽馬信號。在沒有放射性源的環(huán)境下所接收到的伽馬信號應(yīng)該是微弱的、均勻的，且其在各個能量之間分布大致是相同的，沒有特殊的能峰存在，如圖6（a）所示。而在有99mTc 放射性元素的照射下，如圖6（b）所示，能譜圖上出現(xiàn)了明顯的能峰，能量大致集中在某一位置。綜合圖6（a）和（b）來看，本系統(tǒng)能夠?qū)崿F(xiàn)檢測相關(guān)環(huán)境中的伽馬射線。在檢測的過程中，形成明顯能峰的時間在30 s之內(nèi)，說明選取的晶體尺寸薄厚適宜，后續(xù)能夠在遠(yuǎn)小于小動物麻醉期的時間內(nèi)獲得想要的信息。同時，對采集到的能量進(jìn)行相應(yīng)的成像，由于自然環(huán)境中的伽馬射線能量大小分布廣、信號弱，故不取某一特定能量范圍進(jìn)行成像，而選取全范圍的成像，主要觀測圖像的分布是否均勻，如圖6（c）所示，所成的散點(diǎn)圖中每一個區(qū)域分布的伽馬射線能量大致是相同的，這與預(yù)期的一致。且沒有經(jīng)過任何矯正的散點(diǎn)陣列圖的線性度和均勻性遠(yuǎn)比常見的伽馬射線散點(diǎn)圖效果更好，從側(cè)面也驗證了系統(tǒng)的空間分辨力和檢測靈敏度得到了一定的平衡。從圖6 來看，本系統(tǒng)可以檢測到應(yīng)用環(huán)境中的伽馬射線，而且測量的準(zhǔn)確性也達(dá)到了需求，滿足了相關(guān)的檢測使用要求。

圖6 伽馬射線檢測系統(tǒng)采集的圖像

4 結(jié)語

本文基于LabWindows 這一開發(fā)平臺，利用定制的NaI（Tl）晶體和配套的PMT 等重要部件研制了一套可以實時檢測環(huán)境中伽馬射線的完整系統(tǒng)。系統(tǒng)結(jié)構(gòu)簡潔、操作方便，軟件控制程序有穩(wěn)定的運(yùn)行環(huán)境。通過驗證自然環(huán)境和特殊環(huán)境下的伽馬射線能量分布的不同，表明檢測系統(tǒng)可以對環(huán)境中的伽馬射線進(jìn)行有效檢測，為研制用于小動物活體核素成像的伽馬相機(jī)打下了基礎(chǔ)。由于設(shè)計目標(biāo)是對小動物成像，對系統(tǒng)的空間分辨力要求較高。通過探測器硬件設(shè)計，能夠在理論上達(dá)到固有分辨力需求，但是在實際采集過程中，由于PMT 與閃爍晶體以及電子學(xué)電路的耦合等方面達(dá)不到理想情況，所產(chǎn)生的誤差會對系統(tǒng)的性能造成一定影響。后續(xù)將通過改進(jìn)晶體和PMT 數(shù)目以及增添準(zhǔn)直器，同時對采集到的圖像進(jìn)行二次線性、均勻性及能量校準(zhǔn)，更進(jìn)一步提高系統(tǒng)分辨力，最終形成完整的小動物活體核素檢測伽馬相機(jī)，來檢測載體表面放射性核素99mTc 標(biāo)記放出的伽馬射線，形成相應(yīng)的動態(tài)二維圖像，直觀形象地來描述活體動物對相關(guān)藥物的處置行為，以及藥物在小動物體內(nèi)各個組織器官，特別是藥效部位、吸收部位、代謝消除部位以及重要的功能臟器部位的相應(yīng)動態(tài)分布過程，這將對評價藥物的分布特征、吸收特征具有積極的指導(dǎo)意義[12]，也將輔助指導(dǎo)和決策相關(guān)的臨床試驗設(shè)計。