康普頓相機的GEANT4模擬與反投影圖像重建

丁長驥 毛本將 袁永剛 姜志剛 楊朝文

1（四川大學(xué) 物理科學(xué)與技術(shù)學(xué)院 輻射物理及技術(shù)教育部重點實驗室 成都 610065）2（中國工程物理研究院 核物理與化學(xué)研究所 綿陽 621900）

康普頓相機的GEANT4模擬與反投影圖像重建

丁長驥1,2毛本將2袁永剛2姜志剛2楊朝文1

1（四川大學(xué) 物理科學(xué)與技術(shù)學(xué)院 輻射物理及技術(shù)教育部重點實驗室 成都 610065）
2（中國工程物理研究院 核物理與化學(xué)研究所 綿陽 621900）

模擬構(gòu)建了雙層康普頓相機，利用康普頓散射原理，獲取γ光子在兩層探測器發(fā)生散射前后沉積的能量和位置信息，通過一定算法進行圖像重建，獲得放射源位置信息。該相機具有靈敏度高、體積小、應(yīng)用范圍廣等優(yōu)點。通過蒙特卡羅方法模擬了具有雙層條狀結(jié)構(gòu)位置靈敏探測器（散射探測器、吸收探測器）組成的康普頓相機，其兩層探測器分別由Si和Ge材料構(gòu)成。利用反投影圖像重建算法實現(xiàn)了放射源圖像重建，當(dāng)單點源與探測器距離為40 mm時，成像的效率為0.38%，位置分辨率達到8.0 mm，角分辨率達到3°。對于不同位置和不同個數(shù)點源的情況下，檢驗了反投影重建算法在康普頓相機放射源空間重建中的效果。

康普頓相機，蒙特卡羅模擬，反投影圖像重建算法

20世紀(jì)70年代，Todd[1]首次提出了康普頓相機的概念，但由于當(dāng)時探測器材料和電子學(xué)數(shù)據(jù)采集技術(shù)的限制，該技術(shù)未得到深度發(fā)展。而同時期能對放射性物質(zhì)進行二維定位的單孔和編碼孔成像技術(shù)卻得到了較好發(fā)展。但相比之下，其成像效率和探測面積都不如康普頓相機，且體積更大，重量也更重。在放射性核素應(yīng)用越來越廣泛的今天，出于各種應(yīng)用和安全的考慮，對于放射源的檢測更為重要，需要精確的探知放射源所在具體地點，從而可有針對性地處理。

γ相機在天體物理、核設(shè)施退役、特殊核材料結(jié)構(gòu)分析、放射性核素示蹤等各種重要領(lǐng)域都有重要應(yīng)用。所以作為指標(biāo)優(yōu)良的康普頓相機應(yīng)用前景也是十分廣泛。在核材料檢測方面，康普頓相機有十分重要的應(yīng)用。

2007年，利弗莫爾實驗室成功地開發(fā)出一款康普頓相機SPEIR (SPEctroscopic Imager for γ-Rays)，該相機采用兩塊陣列（直徑70 mm）式Ge探測材料，每塊Ge材料的尺寸為2 cm×2 cm×2cm；探測成像的能量為150 keV?4 MeV，由于SPEIR系統(tǒng)采用Ge探測器，采用制冷的方式，系統(tǒng)成本較高[2]。之后也有許多由不同探測器、探測材料組成的康普頓相機問世。目前國內(nèi)對其研究較少，未制作實際應(yīng)用的康普頓相機，也并沒有系統(tǒng)地對其進行模擬研究。本文通過對其蒙特卡羅模擬獲得了理論數(shù)據(jù)，重建算法的構(gòu)建為實現(xiàn)算法優(yōu)化提供了支撐。

1 康普頓相機

康普頓散射中散射光子能量由式(1)決定：

式中，Ec為散射光子能量；E0為初始光子能量；θ為散射角；mec2為電子的靜止能量。

若其初始入射的光子為單能的話，當(dāng)散射角唯一時，則散射光子的能量也唯一確定。反之若已知初始能量E0與散射光子能量Ec，則散射角θ也是唯一確定的，康普頓相機正是利用此原理對放射源進行定位。在具體重建算法中會考慮散射材料對光子的多普勒效應(yīng)。

基本的康普頓相機由兩層探測器組成，如圖1所示。第一層探測器稱為散射探測器(Scattering detector)，均勻分割成多個單元，構(gòu)成位置和能量皆靈敏的探測器陣列。來自樣品的γ射線在第一個探測器中發(fā)生康普頓散射，在其中沉積能量為E1，散射位置為x1。第二層探測器稱為吸收探測器(Absorbing detector)，康普頓散射產(chǎn)生的散射光子被吸收探測器完全吸收，沉積能量E2，吸收位置為x2。γ在兩個探測器中的相互作用位置x1、x2和能量沉積E1、E2均可直接由兩個位置靈敏探測器測量，γ射線的能量E0可由E1加E2得到。兩層探測器多采用Si、Ge半導(dǎo)體探測器，以此獲得高的能量分辨率與精確的位置信息。

圖1 康普頓相機原理示意圖Fig.1 Schematic representation of Compton camera.

由于康普頓散射公式?jīng)]有給出任何方位角φ的信息，所以γ射線的發(fā)射位置一定處于一個圓椎表面上的某一點，這個圓椎的頂點為x1、中軸線為x1x2、半角為θ。經(jīng)過多次符合條件的事件，可以得到很多圓錐，理論上放射源所在位置就是這些圓錐所交匯的部分。圖像重建時在成像空間內(nèi)，每個像素或體素通過一定的算法得到其權(quán)重，經(jīng)過閾值篩選，確定放射源位置。

根據(jù)康普頓相機的原理，探測部分有散射探測器與吸收探測器兩部分，且同時對位置和能量靈敏。散射探測器記錄γ光子在散射探測器中發(fā)生康普頓效應(yīng)所沉積的能量與發(fā)生散射的位置；吸收探測器記錄發(fā)生康普頓效應(yīng)之后產(chǎn)生的散射光子在吸收探測器中能量沉積與發(fā)生能量沉積的位置，在探測器中沉積的能量即是γ光子轉(zhuǎn)移給電子的能量。所以理想的探測器為陣列式，可方便清晰地記錄和得到每次反應(yīng)所在位置與在此位置上所產(chǎn)生的能量沉積。

當(dāng)陣列式康普頓相機所包含探測單元較多時，所需電子學(xué)線路也較多，因此對符合電路與其他電子學(xué)系統(tǒng)等要求很高。因此基本的康普頓相機散射、吸收探測器由雙層正交放置的位置靈敏條狀半導(dǎo)體探測單元構(gòu)成，如圖2所示。

一方面雙層條狀探測單元可以對γ光子反應(yīng)處的橫縱坐標(biāo)同時輸出，獲得作用位置與其所沉積的能量。另一方面，將普通的陣列式輸出變?yōu)橛蓷l狀探測器分別輸出，簡化了電子學(xué)線路。

圖2 典型的康普頓相機結(jié)構(gòu)[3]Fig.2 Typical structure of Compton camera[3].

2 康普頓相機GEANT4模擬

在康普頓相機的研究過程中，使用GEANT4的標(biāo)準(zhǔn)模型作為工具進行計算機模擬，不僅可以模擬構(gòu)建出相機模型，獲得放射源圖像重建所需數(shù)據(jù)，構(gòu)建和驗證重建算法，而且能夠獲得探測器設(shè)計理論依據(jù)和數(shù)據(jù)來源，同時驗證相機性能，對探測器的設(shè)計提出改進。

康普頓相機中對能量和位置同時靈敏的探測單元尺寸應(yīng)當(dāng)較小，如果探測器尺寸過大，會導(dǎo)致輸出電路路數(shù)過多，且需要對大量事件進行一個符合甄別，在電子學(xué)實現(xiàn)上有難度；同時如果兩個探測器間距過大，會導(dǎo)致從散射探測器出來的散射光子進入吸收探測器的幾率變小，使康普頓相機效率降低。經(jīng)過多次初步模擬，確定康普頓相機結(jié)構(gòu)如圖3所示，散射探測器、吸收探測器均由雙層條狀探測單元構(gòu)成。模擬的γ光子能量為137Cs源0.662MeV。

圖3 康普頓相機模擬效果圖Fig.3 Compton camera simulation representation.

左側(cè)部分為散射探測器，由Si材料組成。獲得的沉積能量分辨率高，散射處位置分辨率好，同時對于0.662 MeV能量的γ光子散射效果較好。經(jīng)過多次初步模擬確定每個條狀單元尺寸為48 mm× 10mm×3 mm，第一層16個探測單元平行于x軸方

向放置；第二層16個探測單元平行于y軸方向放置，使得兩層探測單元正交。散射探測器整體尺寸為48mm×20 mm×48 mm。右側(cè)的吸收探測器由雙層條狀Ge材料組成，對γ光子沉積能量與作用位置探測靈敏，其探測單元結(jié)構(gòu)、放置與尺寸都與散射探測器相同。散射探測器與吸收探測器相距40 mm。

一個值得注意的問題是從放射源發(fā)出的大量γ光子，每個γ光子從產(chǎn)生到消失（能量在探測器中沉積完全或者逃逸出模擬空間）稱為一個事件。根據(jù)康普頓相機原理，符合條件的事件應(yīng)當(dāng)滿足：一個γ光子從放射源發(fā)出后，在散射探測器第一層發(fā)生一次康普頓散射，在第一層中產(chǎn)生次級電子，次級電子在第二層對應(yīng)單元產(chǎn)生信號，通過兩層分別產(chǎn)生信號的探測單元得到γ光子散射位置。之后的散射γ光子在吸收探測器中沉積完所有能量，意味著在吸收探測器中發(fā)生的過程可以是一次光電效應(yīng)，也可以是經(jīng)過康普頓效應(yīng)后再發(fā)生光電效應(yīng)等?？傊?，需要把散射γ光子所有的能量都沉積消耗在吸收探測器中。但在蒙特卡羅模擬中，大部分γ光子所經(jīng)歷的過程產(chǎn)生的事件，并不是我們所需要的。例如一個γ光子直接在散射探測器中發(fā)生光電效應(yīng)；或者在散射探測器中發(fā)生一次康普頓效應(yīng)后，產(chǎn)生的散射γ光子沒有進入吸收探測器，或者進入吸收探測器但并沒有將能量完全地沉積消耗。為此，模擬得到的大量事件需要經(jīng)過一定的篩選，得到符合要求的事件進行重建。

模擬過程中的甄別根據(jù)能量來進行：第一層探測器只發(fā)生一次康普頓散射并且γ光子在兩層探測器將能量都沉積完全。在實際探測器中，時間甄別與能量甄別都是需要的。時間甄別判斷不同能量沉積是否屬于同一個γ光子，然后再進行能量甄別，判斷是否滿足能量沉積完全這個條件。如果存在能譜展寬時，γ光子的能量甄別就變?yōu)橐粋€范圍，當(dāng)一個γ光子所有的沉積能量在此范圍內(nèi)，那么就判斷其符合能量甄別條件。

3 反投影圖像重建算法

根據(jù)康普頓相機對放射源探測的原理，需要采用適合康普頓相機的重建算法對獲得的數(shù)據(jù)進行計算，從而得到放射源所在位置空間信息。在反投影圖像重建算法中（后文簡稱反投影法），針對康普頓相機，我們可以把放射源所在的空間區(qū)域即成像空間分隔為立方體網(wǎng)格狀空間，每個立方體空間都作為一個體素。很明顯，在成像空間范圍一定的時候，分隔成的體素數(shù)量決定著體素的大小，體素的大小影響著對放射源重建之后的位置分辨率。

康普頓相機中，放射源所產(chǎn)生的γ射線距離散射探測器的軌跡為一條直線，但根據(jù)其探測原理，實際所確定出來的只有以直線軌跡為中心軸的圓錐面而已，即放射源在這樣的一個圓錐面上，而無法確定其精確位置。在康普頓相機的反投影法中，利用大量事件所產(chǎn)生的圓錐面的交集的位置，就是放射源所在位置概率最大的地方。

反投影示意圖見圖4。反投影法簡單地通過基于體素的權(quán)重進行投影從而重建圖像，各個體素在每次事件產(chǎn)生的圓錐經(jīng)過計算都可以獲得一定權(quán)重。權(quán)重作為對應(yīng)于特定的體素和圓錐的事件的一個近似估算的概率。其中每一個事件都可以獲得這次事件中每個體素的一個權(quán)重，然后統(tǒng)計所有事件得到每個體素所估計的強度，來確定分布在成像空間的γ射線發(fā)射核素的位置。放射源在每個體素的概率即是與這個體素所有事件權(quán)重的總和成正比。

圖4 反投影示意圖Fig.4 Backprojection diagram.

估計的權(quán)重計算通過整合概率密度函數(shù)在事件體素圓錐表面計算得到。

式中，ωij是在投影事件中，事件j在體素i所對應(yīng)的估計權(quán)重；r代表成像空間中某一點的位置坐標(biāo)。權(quán)重的計算可以由pdf (probability density function)函數(shù)計算各個體素得到。

式中，θoj為jth圓錐的孔徑大??；θ?θoj是從對應(yīng)視野的點到j(luò)th圓錐面的最小角距離；δθj是圓錐孔徑的角誤差。由于實際散射角與計算得到的圓錐角誤差不會完全相同，所以權(quán)重的具體計算公式在不同范圍內(nèi)有所不同[4]。

4 反投影法圖像重建

4.1 單點源條件

散射探測器靠近放射源的一面為XY平面，探測器表面中心點垂直向外為Z軸。在GEANT4中構(gòu)造一個點源，發(fā)射角度一定，坐標(biāo)(0,0,40)，位于Z軸上，與散射探測器表面距離為40 mm。散射探測器外側(cè)定義為一立方體成像空間，成像網(wǎng)格空間體積為60 mm×60 mm×60 mm，一共分割為121×121× 61個體素，每個體素在XY平面間隔0.5 mm。因此，成像空間在Z軸的范圍是10?70，每層切面相距1mm。

計算得到大量模擬數(shù)據(jù)，經(jīng)過篩選保留符合康普頓相機成像條件的事件。利用反投影法計算得到平行于XY平面不同距離處，即不同Z軸坐標(biāo)處切面的體素的權(quán)重，獲得放射源重建圖像，其中考慮了Si材料對散射光子的多普勒效應(yīng)。距離散射探測器由遠及近切面的成像效果如圖5所示。

圖5 (a)?(i)依次表示沿Z軸坐標(biāo)由大變小，即成像切面與散射探測器距離由遠及近時各切面的成像效果。根據(jù)康普頓相機成像的原理，平行于XY平面的切面與每次事件所形成的圓錐面相交形成一個橢圓，大量事件產(chǎn)生大量橢圓，而橢圓相交處即是放射源所在位置概率最大處。圖5中，隨著成像空間的切面Z軸坐標(biāo)由大變小，切面向探測器靠近，在圖5(g)中Z=17處達到最高分辨率，之后隨著距離越靠近得到的分辨率也在下降。圖5(g)的XZ面權(quán)重峰形如圖6所示。

圖6 單源圖5(g)的XZ面峰形Fig.6 XZ-palne peak shape of single source in Fig.5(g).

圖7 單源圖5(g)的XZ面峰形權(quán)重Fig.7 XZ-palne peak weight of single source in Fig.5(g).

將一維位置譜的半高寬(Full Width at Half Maximum, FWHM)定義為康普頓成像裝置的空間分辨率，圖7中可以得到在此條件下位置分辨率FWHM=8.0 mm，點源產(chǎn)生的γ光子為200000，篩選符合的個數(shù)為760，定義其成像效率為0.38%，角分辨率為3°，優(yōu)于通常采用的機械準(zhǔn)直器的成像系統(tǒng)空間位置的分辨率和成像效率。隨著放射源離探測器距離的增大，空間位置分辨率和效率都會明顯降低[5]。

4.2 雙點源條件

當(dāng)模擬放射源為雙點源時，探測器距離、探測器、成像的條件和算法都不發(fā)生改變。兩點源位于平行于XY平面的同一平面上，坐標(biāo)為(10,10,40)、(?10, ?10,40)。兩個點源發(fā)射γ光子的概率，發(fā)射角度完全相同。與單源相同，圖8(a)?(i)依次表示沿Z軸坐標(biāo)由大變小時各切面的成像效果，如圖8所示。

圖8 雙點源不同距離處切面的成像Fig.8 Imaging of facet with different distance by double point.

圖8 中顯示了雙點源情況下由遠及近不同切面處的重建圖像，在圖8(g)中可以得到清晰的雙點源重建圖像。圖9代表雙點源其中一個點源的峰形，可以得到點源空間分辨率為FWHM=10.0 mm， 200000個γ光子中符合要求的個數(shù)為713，成像效率為0.36%。與單源相比，其空間分辨率、角分辨率與成像效率都有所降低，原因是雙源所產(chǎn)生的γ光子相互干擾，降低了成像的質(zhì)量。

圖9 雙源圖8(g)中左側(cè)點源XZ面峰形權(quán)重Fig.9 XZ-palne peak weight of left single one of double source in Fig.8(g).

隨放射源數(shù)目增加，雖然分辨率下降明顯，但當(dāng)在反投影重建算法中，確定了放射源所在空間平面與成像空間平面的關(guān)系時，即通過重建算法可以確定放射源與探測器平面之間距離時，是可以對面源進行成像的。

5 結(jié)語

通過蒙特卡羅模擬與反投影重建算法可知，康普頓相機的γ光子的利用效率與位置分辨率都優(yōu)于單孔、編碼孔成像γ相機，空間分辨率可以達到FWHM=8.0 mm，角分辨率達到3°，γ光子成像效率達到0.38%，但反投影重建算法存在數(shù)據(jù)量大、計算時間長等不足。

通過蒙特卡羅模擬獲得了設(shè)計、制作、改進康普頓相機的理論數(shù)據(jù)。實際探測器中材料厚度更薄，所以之后工作中需要進一步改進康普頓相機的模擬，研究其性能參數(shù)，同時通過采用極大似然法優(yōu)化重建算法，提高康普頓相機的位置分辨率，進一步為制作生產(chǎn)提供依據(jù)。

1 Tdoo R W, Nightingale J M, Everett D B. A proposed γ camera[J]. Nature, 1974, 251: 124?132

2 Mihailescu L, Vetter K M, Burks M T, et al. SPEIR: a Ge Compton camera[R]. Lawrence Livermore National Laboratory, UCRL-TR-202620, March 1, 2004

3 David Peter Thomas Scraggs. Digital signal processing techniques for semiconductor Compton cameras[D]. University of Liverpool, November 2007

4 Mattafirri S. On Compton imaging[D]. UC Berkeley, 2010

5 楊永峰, Gono Y, Motomura S, 等. 一個用于多重示蹤技術(shù)的康普頓成像裝置[J]. 同位素, 2001, 14(3?4): 155?160 YANG Yongfeng, Gono Y, Motomura S, et al. A Compton camera for multitracer technique[J]. Journal of Isotopes, 2001, 14(3?4): 155?160

GEANT4 simulation and imaging reconstruction by backprojection algorithm for Compton camera

DING Changji1,2MAO Benjiang2YUAN Yonggang2JIANG Zhigang2YANG Chaowen1

1(College of Physical Science and Technology, Key Laboratory of Radiation Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, China) 2(Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics, Mianyang 621900, China)

Background: Due to its imaging principle advantages and the development of detectors, nuclear electronics, Compton camera has once again put forward. Purpose: This study aims to build a Compton camera via Monte Carlo simulation and test its backprojection algorithm. Methods: By using Monte Carlo simulation software GEANT4, a Compton camera with stripes structure is designed and data for reconstruction image of source is read out. Backprojection imaging reconstruction algorithm is studied for getting the reconstruction of source. Results: Using backprojection algorithm and simulated Compton camera, when the distance between source and the scattering detector is 40 mm, the spatial resolution is FWHM=8.0 mm, angular resolution is 3° and γ photon imaging efficiency is 0.38%, which are better than most of the pinhole and coded aperture γ camera imaging. Conclusion: Because of the existence of large amount of data, Compton camera needs long calculation time for measuring data. The next step is improving simulated Compton camera to study its performance parameters, while the maximum likelihood algorithm could improve the spatial resolution of Compton camera.

Compton camera, Monte Carlo simulation, Backprojection imaging reconstruction algorithm

TL814，R445.5

10.11889/j.0253-3219.2015.hjs.38.110402

丁長驥，男，1991年出生，2013年畢業(yè)于四川大學(xué)，現(xiàn)為碩士研究生，主要從事康普頓相機模擬與圖像重建

袁永剛，E-mail: yyg_2000@yeah.net

2015-06-24，

2015-09-11

CLC TL814 , R445.5