基于解剖標(biāo)志點(diǎn)約束的腦部PET/CT圖像圖譜配準(zhǔn)與腦區(qū)分割

陳 朝 峰，邱 天 爽，馮 洪 波，張 延 軍，王 洪 凱*

(1.大連理工大學(xué) 電子信息與電氣工程學(xué)部，遼寧 大連 116024；2.鹽城工學(xué)院 信息工程學(xué)院，江蘇 鹽城 224051；3.大連醫(yī)科大學(xué)附屬第一醫(yī)院 核醫(yī)學(xué)科，遼寧 大連 116011 )

0 引 言

正電子發(fā)射斷層掃描(positron emission tomography，PET)是一種非侵入性核功能成像技術(shù)，它利用放射性核素在分子水平上實(shí)現(xiàn)可視化和測量人腦的代謝情況．因此，在精神病學(xué)、神經(jīng)外科和神經(jīng)病學(xué)等醫(yī)學(xué)領(lǐng)域，PET在人腦疾病的診斷、分期、治療規(guī)劃和隨訪等方面具有明顯的優(yōu)勢[1-6]．在臨床前期和臨床應(yīng)用中，可使用腦部PET成像來確定組織的功能特征．準(zhǔn)確的定量測量使醫(yī)生能夠評估治療期間和治療后病變部位的變化，從而更好地評估腫瘤灌注、通透性、血容量和治療反應(yīng)．盡管PET的功能性遠(yuǎn)高于常規(guī)結(jié)構(gòu)成像，但仍需要計(jì)算機(jī)斷層掃描(computed tomography，CT)或磁共振(magnetic resonance，MR)成像的解剖學(xué)信息來準(zhǔn)確定位放射性示蹤劑．因此，PET越來越多地和CT或者M(jìn)R整合在一起組成多模態(tài)醫(yī)學(xué)影像系統(tǒng)，以結(jié)合高分辨率的解剖結(jié)構(gòu)信息與放射性示蹤劑代謝高靈敏度的功能成像信息．當(dāng)PET與CT或MR結(jié)合時(shí)，利用功能信息(PET)和結(jié)構(gòu)信息(CT或MR)可取得比單獨(dú)使用任何一種方法都更高的靈敏度和特異性．現(xiàn)在，臨床實(shí)踐中的新規(guī)范是獲取一體化的PET/CT或PET/MRI多模態(tài)圖像來代替單一的PET，綜合利用功能和結(jié)構(gòu)信息，可以準(zhǔn)確標(biāo)記和識別代謝異常位置，為患者提供更準(zhǔn)確的診斷．

對于腦部PET圖像的腦區(qū)量化，經(jīng)由受過神經(jīng)解剖訓(xùn)練的專家根據(jù)個(gè)體的MR圖像手動繪制每個(gè)區(qū)域，仍是定義解剖感興趣體積(volume of interest，VOI)的金標(biāo)準(zhǔn)．然而，手工繪制VOI非常耗時(shí)并且過度依賴于操作者，容易出現(xiàn)評分員偏差甚至可能出現(xiàn)低的再現(xiàn)性．在PET/CT圖像中雖然包含結(jié)構(gòu)成像CT，但是CT相對MR成像的密度分辨率較低，在軟組織特別是在組織間密度變化微小的區(qū)域，CT成像質(zhì)量變差，CT成像組織密度對比范圍小的劣勢就凸顯出來．因此，近年來腦部PET圖像的計(jì)算機(jī)輔助分割從特定腦區(qū)分割到全腦分割方法[7-9]得到了廣泛研究．一些研究專注于基于馬爾可夫場(Markov random field，MRF)[9-10]，利用腦部PET圖像的空間特征進(jìn)行分割．然而，基于MRF的PET圖像分割方法只能粗略地將腦部PET圖像劃分為灰質(zhì)(gray matter，GM)、白質(zhì)(white matter，WM)和腦脊液(cerebrospinal fluid，CSF)．當(dāng)前更加主流的做法是基于圖譜的分割方法[11-18]來劃分腦部的每個(gè)VOI．典型的軟件工具箱有SPM[19]、PMOD[11]、FreeSurfer[20]、MIMneuro[11]等．這些工具箱可以直接為個(gè)體的腦部PET圖像提供詳細(xì)的全腦分割．從實(shí)用角度來看，通過這些工具包可以方便快捷地獲得相對準(zhǔn)確的全腦分區(qū)結(jié)果．然而，這些分割方法潛在局限性是它們的基本原理主要基于全局灰度匹配優(yōu)化策略．對于局部腦區(qū)的解剖結(jié)構(gòu)特別是在深層腦結(jié)構(gòu)部分，這類方法的VOI勾畫精度有限．此外，這類方法忽略了腦解剖結(jié)構(gòu)中不同腦區(qū)間的空間位置信息，而這些空間位置信息對提高分割的準(zhǔn)確率是特別關(guān)鍵的．

鑒于上述腦部PET圖像分割的困難和獨(dú)特挑戰(zhàn)，人們?nèi)匀恍枰芯扛泳_的腦部PET圖像分割方法．本文的思想最初是由用于醫(yī)學(xué)圖像分析和手術(shù)規(guī)劃的灰度和解剖標(biāo)志點(diǎn)聯(lián)合配準(zhǔn)的研究所啟發(fā)[21-22]．根據(jù)文獻(xiàn)檢索情況，目前對于腦部PET/CT圖像分割還沒有將基于解剖標(biāo)志點(diǎn)和基于灰度相結(jié)合的相關(guān)研究．本文將這一思想推廣到腦部PET/CT圖像分割領(lǐng)域，提出一種結(jié)合解剖形態(tài)和解剖標(biāo)志點(diǎn)的雙模態(tài)PET/CT圖像腦區(qū)分割方法．為了充分利用人體解剖形態(tài)的先驗(yàn)知識，將中國人3D結(jié)構(gòu)腦部圖譜參照PET/CT圖像的灰度值填充成相應(yīng)的三維圖像，作為配準(zhǔn)到個(gè)體的參考模板．然后，使用人工指定的解剖標(biāo)志點(diǎn)作為引導(dǎo)，實(shí)現(xiàn)兼顧空間位置信息和灰度信息的配準(zhǔn)工作．

1 實(shí)驗(yàn)方法

圖1展示了本文基于解剖標(biāo)志點(diǎn)和圖像灰度腦部配準(zhǔn)方法的工作流程．本文采用公開的高質(zhì)量頭部解剖圖譜[23-24]作為參考圖譜，并將參考圖譜按照PET/CT圖像的灰度值填充之后配準(zhǔn)到每個(gè)受試者頭部PET/CT圖像．

圖1 頭部PET/CT圖像配準(zhǔn)方法的流程Fig.1 The workflow of the method for the head PET/CT images registration

1.1 頭部圖譜的預(yù)處理

為了利用腦部精確解剖結(jié)構(gòu)的先驗(yàn)知識，增強(qiáng)配準(zhǔn)的魯棒性，本文將已有的人體頭部形狀模型配準(zhǔn)到測試個(gè)體圖像，以獲得每個(gè)個(gè)體腦部的配準(zhǔn)結(jié)果．具體地，本文使用通用頭部數(shù)字模型BodyParts3D[23]，并用中國人腦部數(shù)字圖譜(Chinese2020)[24]替換BodyParts3D的大腦．Chinese2020 是結(jié)構(gòu)精細(xì)的中國人標(biāo)準(zhǔn)腦模板，可以更好地表征中國人群的腦部形態(tài)學(xué)特征．這里的頭部圖譜是用三維三角形網(wǎng)格繪制的．因此，根據(jù)個(gè)體PET/CT圖像中對應(yīng)位置的灰度值，對三角形網(wǎng)格繪制的頭部圖譜進(jìn)行填充，各腦區(qū)填充的像素值取自實(shí)際PET/CT圖像中的像素值，從而實(shí)現(xiàn)由曲面網(wǎng)格數(shù)據(jù)轉(zhuǎn)換為三維圖像數(shù)據(jù)．

1.2 基于解剖標(biāo)志點(diǎn)的配準(zhǔn)

本文解剖標(biāo)志點(diǎn)的概念是指具有解剖學(xué)意義的特征點(diǎn)并指示特定的解剖位置．大腦中常用的解剖標(biāo)志點(diǎn)有前連合點(diǎn)(anterior commissure，AC)、后連合點(diǎn)(posterior commissure，PC)和各個(gè)腦區(qū)中心等．通常，通過解剖標(biāo)志點(diǎn)來確保在同一物種內(nèi)不同個(gè)體之間的對應(yīng)關(guān)系．基于解剖標(biāo)志點(diǎn)的配準(zhǔn)方法是解決參考模板和個(gè)體樣本之間的對應(yīng)和轉(zhuǎn)換問題．采用薄板樣條(thin-plate-spline，TPS)方法[21-22，25-26]實(shí)現(xiàn)解剖標(biāo)志點(diǎn)間配準(zhǔn)，通過TPS插值方式求解參考圖譜結(jié)構(gòu)映射到個(gè)體圖像的非線性空間變換Fl．通過變換Fl，整個(gè)模板網(wǎng)格模型被映射到單個(gè)圖像空間．基于解剖標(biāo)志點(diǎn)的配準(zhǔn)方法利用映射模板拓?fù)浣Y(jié)構(gòu)到個(gè)體圖像來表示個(gè)體的解剖結(jié)構(gòu)，從而實(shí)現(xiàn)對個(gè)體的分割．

1.3 基于圖像灰度的配準(zhǔn)

基于圖像灰度的非線性配準(zhǔn)方法[27]采用互信息(mutual information，MI)作為相似性度量，自適應(yīng)隨機(jī)梯度下降(adaptive stochastic gradient descent，ASGD)作為優(yōu)化器，通過B樣條空間變換把參考圖譜配準(zhǔn)到個(gè)體圖像上．經(jīng)過對不同的相似性測度進(jìn)行測試，發(fā)現(xiàn)MI對于本文的PET/CT測試數(shù)據(jù)是最魯棒的．基于圖像灰度的配準(zhǔn)方法可以看作估計(jì)源圖像到目標(biāo)圖像一致性的變換，使參考圖譜(源圖像)和訓(xùn)練圖像(目標(biāo)圖像)之間的相似性測度最小化．參考圖譜配準(zhǔn)到個(gè)體之后，也就是將參考圖譜的解剖結(jié)構(gòu)分別映射到個(gè)體圖像的相似解剖位置，得到非線性變形場Fi．

1.4 灰度信息與解剖標(biāo)志點(diǎn)相結(jié)合的圖譜配準(zhǔn)

上面提到的兩種配準(zhǔn)方法，其中基于解剖標(biāo)志點(diǎn)的配準(zhǔn)方法實(shí)現(xiàn)了關(guān)鍵解剖特征位置的精確對齊，但是忽略了整體結(jié)構(gòu)的匹配；與之相反，基于圖像灰度的配準(zhǔn)方法不考慮局部特征，而是對配準(zhǔn)過程進(jìn)行全局優(yōu)化．因此，本文提出將這兩種配準(zhǔn)方法相互融合，從而實(shí)現(xiàn)兩種方法的優(yōu)勢互補(bǔ)，兼顧圖像的空間解剖位置信息和灰度信息以獲得更加準(zhǔn)確的配準(zhǔn)結(jié)果．

基于組合的配準(zhǔn)方法是對上述兩種空間變換函數(shù)線性加權(quán)組合，即：

Fc(x)=W(x)·Fl(x)+(1-W(x))·Fi(x)

(1)

其中x是三維空間中任意點(diǎn)的空間位置．Fl(x)和Fi(x)分別是基于解剖標(biāo)志點(diǎn)和基于圖像灰度的配準(zhǔn)方法得到的x處的變形向量．Fc(x)是線性加權(quán)之后的變形向量．變形向量定義為位移矢量，即F(x)∶=(ΔxΔyΔz)．

理論上，當(dāng)x接近解剖標(biāo)志點(diǎn)時(shí)，加權(quán)函數(shù)W(x)應(yīng)接近1；當(dāng)x遠(yuǎn)離解剖標(biāo)志點(diǎn)時(shí)，加權(quán)函數(shù)W(x)應(yīng)接近0．這樣Fc(x)在解剖標(biāo)志點(diǎn)附近時(shí)相當(dāng)于Fl(x)，在遠(yuǎn)離解剖標(biāo)志點(diǎn)時(shí)相當(dāng)于Fi(x)．因此，將d(x)定義為從x到其最近的解剖標(biāo)志點(diǎn)的距離，即：

(2)

其中Li是第i個(gè)解剖標(biāo)志點(diǎn)的坐標(biāo)．因此，將函數(shù)W(x)轉(zhuǎn)換為關(guān)于d(x)的函數(shù)，即W(d(x))，簡記為W(d)．W(d)的取值范圍是[0，1]，并且滿足單調(diào)遞減．另外，W(d)也應(yīng)該滿足凹面的，這樣當(dāng)d接近無窮大時(shí)，它的下降速度會變慢．

綜上，可以得到函數(shù)W(d)[21]為

(3)

其中參數(shù)δ控制W(d)下降速度，經(jīng)過測試發(fā)現(xiàn)實(shí)驗(yàn)結(jié)果對這個(gè)參數(shù)不敏感．δ在4～8 mm時(shí)，結(jié)果相似．對于式(3)取經(jīng)驗(yàn)參數(shù)δ=6 mm．利用式(3)定義的加權(quán)函數(shù)，計(jì)算組合變形場Fc(x)并將其應(yīng)用于參考圖譜，得到每個(gè)測試對象對應(yīng)的變形模板圖譜．

2 實(shí)驗(yàn)與結(jié)果

2.1 測試數(shù)據(jù)

作為回顧性研究，從醫(yī)院數(shù)據(jù)庫中獲取了健康體檢中20例無癥狀測試者的頭部PET/CT圖像數(shù)據(jù)．本文所用數(shù)據(jù)采集自西門子Biograph64型PET/CT掃描儀(52個(gè)硅酸镥(LSO)晶體環(huán)的PET和64排螺旋CT)．放射化學(xué)純度皆大于95%．患者禁食4～6 h以上，靜脈注射顯像劑[18F]-FDG 7.4×106Bq/kg，安靜避光環(huán)境平臥60 min后，以治療體位行CT和PET頭部斷層顯像，圖像同機(jī)融合，利用4.0 mm高斯卷積核濾波后的反投影重建圖像．PET采集圖像數(shù)據(jù)切片厚度為1.5 mm，單個(gè)像素大小1.0 mm，采用336×336矩陣．CT圖像的采集參數(shù)為管電壓100～140 kV，電流28～298 mA，像素尺寸為0.59～1.37 mm，層間距為1.25～3.00 mm，采用512×512矩陣．考慮到年齡是影響腦部示蹤劑代謝的重要因素，因此本文的測試數(shù)據(jù)涵蓋了不同的性別和年齡段，以驗(yàn)證方法對于不同性別和年齡測試數(shù)據(jù)具有較好的普適性．表1顯示了無癥狀受試者的年齡和性別分布．

表1 PET/CT測試圖像不同年齡和性別的分布Tab.1 Distribution of collected PET/CT images of different ages and genders

2.2 實(shí) 驗(yàn)

由于現(xiàn)有基于圖譜的腦部PET分割方法大多采用全局灰度匹配的優(yōu)化策略，對于局部腦區(qū)的解剖結(jié)構(gòu)特別是在深層腦結(jié)構(gòu)部分的VOI勾畫精度有限．因此本文選擇了6個(gè)深層腦區(qū)作為腦部PET圖像的VOI，包括左右丘腦(left and right thalamus，THA_L and THA_R)、左右尾狀核(left and right caudate nucleus，CAU_L and CAU_R)、左右側(cè)殼核(left and right putamen，PUT_L and PUT_R)．相應(yīng)地，選擇每個(gè)VOI的中心作為解剖標(biāo)志點(diǎn)，由臨床醫(yī)學(xué)專家在PET圖像中手工標(biāo)注．基于全局灰度的匹配優(yōu)化策略在深層腦區(qū)配準(zhǔn)效果不理想，為了改善配準(zhǔn)效果，在這些區(qū)域引入解剖標(biāo)志點(diǎn)引導(dǎo)和約束配準(zhǔn)過程，實(shí)現(xiàn)深層腦區(qū)的精確配準(zhǔn)．而對于腦皮層區(qū)域，一方面，基于灰度配準(zhǔn)良好的全局優(yōu)化性能，配準(zhǔn)效果已經(jīng)相當(dāng)精確；另一方面，腦皮層的環(huán)形結(jié)構(gòu)使得不容易選擇子腦區(qū)的中心點(diǎn)．因此，本文解剖標(biāo)志點(diǎn)的指定都集中于核心腦區(qū)．

本文方法是通過把預(yù)處理的頭部圖譜配準(zhǔn)到個(gè)體頭部的PET/CT圖像，進(jìn)而得到分割好的個(gè)體大腦．通過與傳統(tǒng)的基于圖像灰度的配準(zhǔn)(intensity-based registration，IR)方法對比，對本文方法性能進(jìn)行評價(jià)．基于圖像灰度的配準(zhǔn)方法將預(yù)處理的頭部圖譜配準(zhǔn)到每個(gè)測試圖像中．本文方法是基于同樣的圖像配準(zhǔn)策略，增加關(guān)鍵的解剖標(biāo)志點(diǎn)約束配準(zhǔn)過程．

2.3 定性結(jié)果

圖2分別展示了個(gè)體融合的頭部PET/CT圖像通過基于傳統(tǒng)灰度的配準(zhǔn)方法與本文基于圖像灰度和解剖標(biāo)志點(diǎn)聯(lián)合的配準(zhǔn)方法，個(gè)體全腦分割結(jié)果在矢狀面、冠狀面和軸向面3個(gè)切面上效果圖．

從圖2可以明顯看出，本文所提出的配準(zhǔn)方法可以準(zhǔn)確地刻畫大腦的各個(gè)腦區(qū)，獲得了良好的配準(zhǔn)結(jié)果．基于圖像灰部形狀的全局對齊，在頭部的高密度區(qū)域(如顱骨)可以實(shí)現(xiàn)很好的配準(zhǔn)，但是在頭部的低密度區(qū)域(如腦室)，會出現(xiàn)嚴(yán)重的配準(zhǔn)誤差．相比之下，在基于圖像灰度配準(zhǔn)結(jié)果精度不高的腦區(qū)，解剖標(biāo)志點(diǎn)的介入可以約束配準(zhǔn)過程，提高配準(zhǔn)精度．在配準(zhǔn)過程中引入解剖標(biāo)志點(diǎn)，也就是利用深度腦區(qū)的空間解剖位置信息優(yōu)化配準(zhǔn)，從而提高了腦區(qū)VOI的配準(zhǔn)精度，如丘腦、尾狀核、殼核等．圖2為本文方法配準(zhǔn)結(jié)果與基于傳統(tǒng)灰度的配準(zhǔn)結(jié)果對比，第一行在矢狀面、冠狀面和軸向面用黑色圈起來的區(qū)域，解剖標(biāo)志點(diǎn)的引入將VOI的配準(zhǔn)限制在正確的解剖位置，即相應(yīng)腦區(qū)的真實(shí)解剖位置，本文方法明顯改善了腦區(qū)VOI的配準(zhǔn)結(jié)果．

圖2 3個(gè)切面下基于灰度和基于聯(lián)合的配準(zhǔn)方法個(gè)體全腦分割結(jié)果Fig.2 Whole brain parcellation results based on intensity and joint registration methods in three sections

2.4 定量結(jié)果

為了定量評估腦區(qū)分割結(jié)果的準(zhǔn)確性，對基于圖像灰度的分割結(jié)果與基于圖像灰度和解剖標(biāo)志點(diǎn)聯(lián)合的分割結(jié)果中的6個(gè)深層腦區(qū)CAU_R、CAU_L、PUT_R、PUT_L、THA_R和THA_L作為腦部PET圖像的VOI與醫(yī)學(xué)專家提供的金標(biāo)準(zhǔn)進(jìn)行比較．使用相似系數(shù)S(dice similarity coefficient，DSC)、平均表面距離D(average surface distance，ASD)和器官體積恢復(fù)系數(shù)R(recovery coefficient，RC)3個(gè)評價(jià)指標(biāo)來進(jìn)行定量分析，具體定義如下：

(4)

(5)

(6)

式中：S測量分割腦區(qū)Rr與專家分割金標(biāo)準(zhǔn)對應(yīng)區(qū)域Rt的重疊率，|·|表示感興趣區(qū)域的體素?cái)?shù)量，∩表示對兩個(gè)感興趣區(qū)域取交集．D測量分割腦區(qū)與專家分割金標(biāo)準(zhǔn)對應(yīng)區(qū)域之間的平均表面距離，nr和nt分別代表分割腦區(qū)和專家分割金標(biāo)準(zhǔn)對應(yīng)區(qū)域中的曲面網(wǎng)格頂點(diǎn)數(shù)，di是從分割腦區(qū)的第i個(gè)頂點(diǎn)到專家分割金標(biāo)準(zhǔn)對應(yīng)區(qū)域表面的距離，dj是專家分割金標(biāo)準(zhǔn)的第j個(gè)頂點(diǎn)到分割腦區(qū)對應(yīng)區(qū)域表面的距離．R測量分割腦區(qū)與專家分割金標(biāo)準(zhǔn)對應(yīng)區(qū)域體積比．

以上3個(gè)評價(jià)指標(biāo)可以全面定量地反映腦區(qū)分割結(jié)果的準(zhǔn)確性．DSC代表分割結(jié)果與金標(biāo)準(zhǔn)的重疊率，從整體上反映了分割的VOI形狀、大小和位置的匹配程度，其值越接近1，配準(zhǔn)精度越高．ASD表示分割結(jié)果曲面的各個(gè)頂點(diǎn)到金標(biāo)準(zhǔn)的最小距離的均值，體現(xiàn)了分割的VOI表面距離誤差，其值越小，相應(yīng)的分割精度越高．RC表示分割的VOI與金標(biāo)準(zhǔn)的體積比，反映了分割的VOI體積的精度，其值越接近1，對應(yīng)的分割精度越高．

圖3比較了基于圖像灰度的分割方法與基于圖像灰度和解剖標(biāo)志點(diǎn)聯(lián)合的分割方法在DSC、ASD和RC 3項(xiàng)評價(jià)指標(biāo)下6個(gè)大腦VOI的分割結(jié)果．對于選定的6個(gè)深層腦區(qū)CAU_R、CAU_L、PUT_R、PUT_L、THA_R和THA_L，本文方法顯示了很好的配準(zhǔn)性能，并且這3個(gè)評價(jià)指標(biāo)的分布都比基于圖像灰度的分割方法更集中．無論是從DSC、ASD還是RC上進(jìn)行對比，本文方法都要明顯優(yōu)于基于圖像灰度的分割方法，本文方法得到的DSC更大，ASD更小，RC更接近1．

(a)相似系數(shù)

DSC平均值為79.91%，尤其是右丘腦的DSC值為87.39%．令人鼓舞的是，本文方法得到的平均ASD為0.95 mm，小于像素分辨率(1.37 mm)．本文方法的器官體積配準(zhǔn)恢復(fù)精度同樣優(yōu)于對比方法，并且RC更接近1．將解剖標(biāo)志點(diǎn)和圖像灰度融合到同一個(gè)配準(zhǔn)框架中，進(jìn)一步提高了配準(zhǔn)精度，實(shí)現(xiàn)了DSC最大、ASD最小、RC最接近1．

3 討 論

本文提出了一種新的基于PET/CT圖像的全腦腦區(qū)分割方法．該方法結(jié)合了基于解剖標(biāo)志點(diǎn)和基于圖像灰度配準(zhǔn)方法各自的優(yōu)點(diǎn)，采用BodyParts3D人體三維模型和中國人3D結(jié)構(gòu)腦部圖譜Chinese2020作為參考模板，提高了腦區(qū)分割的精度．通過DSC、ASD和RC 3項(xiàng)評價(jià)指標(biāo)對腦部PET/CT圖像分割效果評估，表明該方法在形狀匹配方面比基于圖像灰度的方法具有更高的精度．

本文方法的主要優(yōu)點(diǎn)是充分利用解剖標(biāo)志點(diǎn)引導(dǎo)局部腦區(qū)的精確配準(zhǔn)，提出了結(jié)合解剖標(biāo)志點(diǎn)和圖像灰度兩種信息的配準(zhǔn)策略，提高了局部腦區(qū)特別是深層腦區(qū)的配準(zhǔn)精度，快速準(zhǔn)確地對腦部PET/CT圖像實(shí)現(xiàn)全腦腦區(qū)分割．相比之下，在傳統(tǒng)腦部圖譜方法中，基于單模態(tài)的全局灰度信息的配準(zhǔn)方法則不能確保局部的、像素值較低的腦區(qū)的配準(zhǔn)精度．對于這一局限，本文引入這些區(qū)域的解剖標(biāo)志點(diǎn)作為配準(zhǔn)過程的約束條件，改善這些區(qū)域的配準(zhǔn)精度．解剖標(biāo)志點(diǎn)一般選擇需要改善區(qū)域的中心或是器官上具有解剖意義的標(biāo)志點(diǎn)．基于解剖標(biāo)志點(diǎn)的配準(zhǔn)方法實(shí)現(xiàn)了關(guān)鍵解剖特征位置的精確對齊，但是忽略了整體結(jié)構(gòu)的匹配；與之相反，基于圖像灰度的配準(zhǔn)方法不考慮局部特征，而是對配準(zhǔn)過程進(jìn)行全局優(yōu)化．因此，本文提出將這兩種配準(zhǔn)方法相互融合，從而實(shí)現(xiàn)兩種方法的優(yōu)勢互補(bǔ)，兼顧圖像的空間解剖位置信息和灰度信息以獲得更加準(zhǔn)確的配準(zhǔn)結(jié)果．在配準(zhǔn)過程中首先對要處理的像素進(jìn)行判斷，當(dāng)待配準(zhǔn)像素靠近解剖標(biāo)志點(diǎn)時(shí)，形變場以基于解剖標(biāo)志點(diǎn)的配準(zhǔn)方法為主，而當(dāng)待配準(zhǔn)像素遠(yuǎn)離解剖標(biāo)志點(diǎn)時(shí)，形變場以基于圖像灰度的配準(zhǔn)方法為主，從而保證了配準(zhǔn)的效果，提高了配準(zhǔn)精度．本文所提的聯(lián)合配準(zhǔn)方法，結(jié)合兩種配準(zhǔn)方法的各自優(yōu)勢，充分利用了圖像的灰度信息和空間解剖位置信息，既考慮全局優(yōu)化，又通過解剖標(biāo)志點(diǎn)實(shí)現(xiàn)局部精確對齊．從1.4節(jié)灰度信息和解剖標(biāo)志點(diǎn)相結(jié)合的圖譜配準(zhǔn)可以知道，影響配準(zhǔn)結(jié)果的主要因素是形變場的計(jì)算．本文綜合考慮實(shí)驗(yàn)的可行性和準(zhǔn)確性，選取了e指數(shù)作為加權(quán)函數(shù)．為了對基于網(wǎng)格的腦部圖譜模板與PET/CT體素圖像匹配，本文提出將已有的腦部圖譜網(wǎng)格模板填充成體素圖像，然后基于圖像灰度進(jìn)行配準(zhǔn)．值得一提的是，本文方法不受人種種族限制．如果要分割其他種族或群體的個(gè)體腦區(qū)，選擇相應(yīng)種族或群體的頭部模板代替本文中國人腦部圖譜進(jìn)行配準(zhǔn)即可．

4 結(jié) 語

本文提出了一種聯(lián)合解剖標(biāo)志點(diǎn)和圖像灰度的配準(zhǔn)方法，能夠?qū)崿F(xiàn)對腦部PET/CT圖像全腦腦區(qū)的精確分割．為了充分利用腦部圖譜的解剖形態(tài)先驗(yàn)知識和腦區(qū)之間的解剖空間位置信息，將已有腦部圖譜網(wǎng)格填充為與配準(zhǔn)個(gè)體PET/CT灰度值相應(yīng)的體素圖像，通過在深度腦區(qū)選定解剖標(biāo)志點(diǎn)來引導(dǎo)配準(zhǔn)過程．利用解剖標(biāo)志點(diǎn)的解剖位置信息來保證關(guān)鍵解剖區(qū)域的局部精確性．與傳統(tǒng)的配準(zhǔn)方法相比，該方法不僅具有更好的全局形狀匹配精度，而且兼顧了局部區(qū)域的配準(zhǔn)．本文提出方法的全腦腦區(qū)分割結(jié)果得到了核醫(yī)學(xué)臨床醫(yī)生的充分肯定．今后的工作重點(diǎn)是實(shí)現(xiàn)關(guān)鍵解剖標(biāo)志點(diǎn)的自動檢測，更好地實(shí)現(xiàn)腦部PET/CT圖像的全腦腦區(qū)分割．