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

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

0 引 言

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

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

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

1 實驗方法

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

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

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

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

1.2 基于解剖標志點的配準

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

1.3 基于圖像灰度的配準

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

1.4 灰度信息與解剖標志點相結(jié)合的圖譜配準

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

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

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

(1)

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

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

(2)

其中Li是第i個解剖標志點的坐標．因此，將函數(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接近無窮大時，它的下降速度會變慢．

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

(3)

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

2 實驗與結(jié)果

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

作為回顧性研究，從醫(yī)院數(shù)據(jù)庫中獲取了健康體檢中20例無癥狀測試者的頭部PET/CT圖像數(shù)據(jù)．本文所用數(shù)據(jù)采集自西門子Biograph64型PET/CT掃描儀(52個硅酸镥(LSO)晶體環(huán)的PET和64排螺旋CT)．放射化學純度皆大于95%．患者禁食4～6 h以上，靜脈注射顯像劑[18F]-FDG 7.4×106Bq/kg，安靜避光環(huán)境平臥60 min后，以治療體位行CT和PET頭部斷層顯像，圖像同機融合，利用4.0 mm高斯卷積核濾波后的反投影重建圖像．PET采集圖像數(shù)據(jù)切片厚度為1.5 mm，單個像素大小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ù)涵蓋了不同的性別和年齡段，以驗證方法對于不同性別和年齡測試數(shù)據(jù)具有較好的普適性．表1顯示了無癥狀受試者的年齡和性別分布．

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

2.2 實 驗

由于現(xiàn)有基于圖譜的腦部PET分割方法大多采用全局灰度匹配的優(yōu)化策略，對于局部腦區(qū)的解剖結(jié)構(gòu)特別是在深層腦結(jié)構(gòu)部分的VOI勾畫精度有限．因此本文選擇了6個深層腦區(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)地，選擇每個VOI的中心作為解剖標志點，由臨床醫(yī)學專家在PET圖像中手工標注．基于全局灰度的匹配優(yōu)化策略在深層腦區(qū)配準效果不理想，為了改善配準效果，在這些區(qū)域引入解剖標志點引導(dǎo)和約束配準過程，實現(xiàn)深層腦區(qū)的精確配準．而對于腦皮層區(qū)域，一方面，基于灰度配準良好的全局優(yōu)化性能，配準效果已經(jīng)相當精確；另一方面，腦皮層的環(huán)形結(jié)構(gòu)使得不容易選擇子腦區(qū)的中心點．因此，本文解剖標志點的指定都集中于核心腦區(qū)．

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

2.3 定性結(jié)果

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

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

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

2.4 定量結(jié)果

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

(4)

(5)

(6)

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

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

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

(a)相似系數(shù)

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

3 討 論

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

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

4 結(jié) 語

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