畢帆，涂文勇，石慧烽，馮琨，樊文慧，胡海生

上海交通大學(xué)醫(yī)學(xué)院附屬第九人民醫(yī)院 口腔頜面-頭頸腫瘤科，上海市，200011

0 引言

計(jì)算機(jī)斷層成像技術(shù)自出現(xiàn)以來，廣泛應(yīng)用于放射診斷和放射治療中?，F(xiàn)代精準(zhǔn)放射治療都是建立在CT斷層圖像上進(jìn)行靶區(qū)勾畫、計(jì)劃設(shè)計(jì)和劑量計(jì)算等工作，由于很多數(shù)患者為術(shù)后放療，部分患者有假牙植入或植入鈦板等金屬物，導(dǎo)致在CT模擬定位掃描中產(chǎn)生金屬偽影，周圍軟組織無法識(shí)別，造成靶區(qū)勾畫和劑量分布計(jì)算等方面存在誤差[1]，給腫瘤患者的精準(zhǔn)治療帶來困難。CT模擬定位中金屬偽影的出現(xiàn)主要是由于在掃描中X線穿過高密度金屬物時(shí)產(chǎn)生衰減，導(dǎo)致對(duì)應(yīng)金屬物在射線方向后投影數(shù)據(jù)缺失，周圍組織因部分容積效應(yīng)、光子不足信號(hào)喪失嚴(yán)重而產(chǎn)生金屬偽影，常表現(xiàn)為星條狀暗區(qū)或放射狀高密度影[2]。

隨著軟硬件技術(shù)的發(fā)展，多種去偽影算法相繼出現(xiàn)[3-5]，MAR（Metal Artifact Reduction）算法是GE公司推出的去金屬偽影算法，目的是揭示被金屬所掩蓋的解剖細(xì)節(jié)。該算法是基于二維投影函數(shù)算法，有助于減少體內(nèi)金屬造成的光子饑餓、光束硬化和條紋狀偽影。MAR算法[6-7]首先確定與金屬物相對(duì)應(yīng)的“毀損”的原始數(shù)據(jù)，由校正后數(shù)據(jù)代替金屬毀損數(shù)據(jù)得到無瑕疵的數(shù)據(jù)，具體是把圖像分類后正向投影，生成無瑕疵數(shù)據(jù)，與原始投影數(shù)據(jù)疊加結(jié)合，得到最終的投影數(shù)據(jù)，從而可顯示隱藏在偽影下的解剖細(xì)節(jié)。本研究通過原型模體外加鈦板金屬物作為CT掃描實(shí)驗(yàn)?zāi)ｓw，分別采用常規(guī)算法和MAR算法對(duì)模體進(jìn)行掃描，通過計(jì)算分析金屬影響層面的圖像質(zhì)量客觀評(píng)價(jià)指標(biāo)，探究金屬偽影對(duì)周圍區(qū)域像素 CT值影響的規(guī)律，探討MAR算法在CT模擬機(jī)金屬偽影去除中的能力。

1 材料與方法

1.1 材料

原型模體為有機(jī)玻璃制作的空心圓柱形模體，高141 mm，外直徑213 mm，模體內(nèi)為純凈水，該模體為GE公司所生產(chǎn)，型號(hào)為2206352。在模體表面用雙面膠粘貼鈦板，鈦板長為7.8 mm，寬為1.5 mm，中間為10個(gè)小孔。將鈦板粘貼于模體表面中間，如圖1所示。使用GE DISCOVERY RT大孔徑螺旋定位CT掃描系統(tǒng)，掃描參數(shù)：管電壓 120 KeV；管電流300 mA；FOV 25 cm；層厚2.5 mm；螺距1.375:1，分別采用常規(guī)算法和MAR算法對(duì)模體進(jìn)行掃描。

圖1 金屬鈦板和掃描模體Fig.1 Metal titanium and phantom

1.2 圖像質(zhì)量客觀評(píng)價(jià)指標(biāo)分析

根據(jù)常規(guī)掃描和MAR模式掃描的兩組圖像，鈦板偽影共影響34個(gè)層面，分別對(duì)上述層面進(jìn)行測量，圖像測量采用同一窗寬窗位（WW：300 Hu，WL：80 Hu）。分別在各影響層面畫出感興趣區(qū)（Region of Interest,ROI）、測量感興趣區(qū)的CT值和噪聲值（即SD值，代表感興趣區(qū)CT值的標(biāo)準(zhǔn)差），分別計(jì)算圖像的偽影指數(shù)（Artifact Index,AI）[8]和對(duì)比噪聲比（Contrast to Noise Ratio,CNR）[9]，計(jì)算方式如式（1）、（2）所示。每組圖像隨機(jī)選擇一個(gè)偽影最嚴(yán)重區(qū)域，記為ROI 1；另選擇兩個(gè)無偽影區(qū)域，記為ROI 2和ROI 3，ROI=76 mm2。

1.3 統(tǒng)計(jì)學(xué)分析

使用Microsoft Excel 2016建立數(shù)據(jù)庫，采用SPSS 21.0進(jìn)行統(tǒng)計(jì)學(xué)分析，將所得的數(shù)據(jù)進(jìn)行檢驗(yàn)，AI和CNR值等計(jì)量資料用（±s）表示，計(jì)量資料的統(tǒng)計(jì)用配對(duì)t檢驗(yàn)，P＜0.01為具有統(tǒng)計(jì)學(xué)差異。

2 結(jié)果

2.1 兩種算法下AI指數(shù)和CNR的對(duì)比

通過兩種算法分別對(duì)CT模體進(jìn)行掃描，每組掃描共獲得88層，其中受偽影影響有34層面，分別對(duì)兩組圖像受影響層面進(jìn)行測量，根據(jù)式（1）、式（2）分別算出偽影指數(shù)和對(duì)比噪聲比。根據(jù)偽影指數(shù)式（1）可得，AI指數(shù)越大，說明圖像受到金屬偽影的影響越大；反之越小。在兩種不同的掃描算法下，模體的AI指數(shù)和CNR值的統(tǒng)計(jì)結(jié)果如表1所示。在偽影指數(shù)指標(biāo)方面：常規(guī)算法的AI指數(shù)（20.65±5.04）要明顯高于MAR算法的AI指數(shù)（10.28±2.60），表明通過MAR算法，圖像的偽影影響明顯下降。對(duì)比噪聲比的影響受到CT值和SD值的影響。CNR指數(shù)是客觀評(píng)價(jià)圖像質(zhì)量標(biāo)準(zhǔn)中的一項(xiàng)重要指標(biāo)，可以有效反應(yīng)圖像中ROI的識(shí)別程度。常規(guī)算法的CNR值（5.61±1.36）低于MAR算法的CNR值（7.81±1.12），通過MAR算法可提高圖像的CNR值，表明各偽影影響層面，圖像質(zhì)量有一定程度提升。分別對(duì)常規(guī)算法和MAR算法下的兩個(gè)參數(shù)行配對(duì)t檢驗(yàn)，AI指數(shù)（t=12.14，P＜0.01）和CNR（t=-10.12，P＜0.01）統(tǒng)計(jì)學(xué)差異明顯，圖2表示了AI指數(shù)和CNR的兩種算法的箱式圖。

表1 兩種算法圖像質(zhì)量定量測量結(jié)果與比較（±s）Tab.1 Quantitative measurement results and comparison of image quality under two algorithms（±s）

表1 兩種算法圖像質(zhì)量定量測量結(jié)果與比較（±s）Tab.1 Quantitative measurement results and comparison of image quality under two algorithms（±s）

圖2 兩種算法AI指數(shù)和CNR指數(shù)對(duì)比箱式圖Fig.2 AI index and CNR index contrast box plot under two algorithms

2.2 不同層面下的去偽影前后的影響

通過上述掃描和計(jì)算，兩組算法下的掃描后受影響的CT圖像層面均為34層，為進(jìn)一步探究各層面影響大小，分別將兩組算法下的AI值制成折線圖，如圖3所示。從折線圖上來看，在常規(guī)算法的掃描下，各層面的AI指數(shù)均大于MAR算法下的AI指數(shù)，各相鄰層面的AI指數(shù)呈現(xiàn)階梯上下變化趨勢，這與金屬鈦板的形狀有一定關(guān)系，本實(shí)驗(yàn)研究使用的金屬鈦板為長條形，中間有相關(guān)空腔，故在斷層掃描中該層面AI指數(shù)較小。進(jìn)一步對(duì)比AI值下降程度，分別統(tǒng)計(jì)影響層面的AI值下降速率，達(dá)到21.72%～88.40%。AI值下降比例最大的為I 37.5層面，常規(guī)算法和MAR算法下的AI值為19.52和2.26，下降比率達(dá)到88.40%。

圖3 兩種算法不同層面?zhèn)斡坝绊懘笮ig.3 The influence of artifacts at different slices under two algorithms

3 討論

在CT模擬機(jī)對(duì)腫瘤定位時(shí)，金屬植入物會(huì)產(chǎn)生嚴(yán)重的偽影，影響鄰近組織結(jié)構(gòu)的觀察，導(dǎo)致在靶區(qū)勾畫、計(jì)劃設(shè)計(jì)、劑量計(jì)算等方面存在不確定性，影響精準(zhǔn)放療的實(shí)施。故在CT模擬定位中如何有效去除金屬偽影同時(shí)最大限度提高圖像質(zhì)量成為繼續(xù)解決的問題。金屬偽影產(chǎn)生的主要原因是人體組織與金屬的密度差異過大[10-11]。金屬偽影的產(chǎn)生與金屬自身的高衰減系數(shù)造成的光子饑餓效應(yīng)有關(guān)，這些偽影以條紋、放射狀等形式顯示出來，這取決于金屬的形狀和性質(zhì)。目前金屬偽影的解決方法主要集中在優(yōu)化掃描參數(shù)設(shè)置[12]、采用能譜CT掃描[13]以及圖像后處理方法上[14]。優(yōu)化掃描參數(shù)的方法是增加管電壓和管電流，應(yīng)用效果一般，且大大增加掃描劑量；且改變管電壓會(huì)影響CT模擬機(jī)的電子密度，從而造成TPS劑量計(jì)算的不準(zhǔn)確性。能譜CT嘗試從硬件角度解決該問題，但目前無法解決因硬化效應(yīng)造成的偽影。圖像后處理方法主要是從算法和技術(shù)上對(duì)偽影問題進(jìn)行處理，主要包含迭代重建法、投影插值法及混合方法等。如部分金屬偽影去除技術(shù)是利用正弦曲線修復(fù)方法將原始數(shù)據(jù)中受到金屬影響的部分替換，雖然該方法在一定程度上消除了金屬偽影，但由于原始數(shù)據(jù)的替換導(dǎo)致數(shù)據(jù)的丟失，造成圖像信噪比降低。

本研究通過原型模體外加鈦板金屬物作為CT掃描實(shí)驗(yàn)?zāi)ｓw，分別采用常規(guī)算法和MAR算法對(duì)模體進(jìn)行掃描，通過計(jì)算分析金屬影響層面的圖像質(zhì)量客觀評(píng)價(jià)指標(biāo)，探討MAR算法在CT模擬機(jī)金屬偽影去除中的能力，研究結(jié)果表明在圖像質(zhì)量方面，應(yīng)用MAR算法后，AI指數(shù)和CNR兩項(xiàng)指標(biāo)均有提升，兩算法的對(duì)比有統(tǒng)計(jì)學(xué)意義。在兩種算法下不同層面?zhèn)斡坝绊懖煌鲗用娴腁I指數(shù)影響大小達(dá)到21.72%～88.40%。MAR算法是GE公司新推出的去偽影算法[15]，該算法僅需一次采集即可得到校正后的圖像，在減輕偽影的同時(shí)減少了檢查時(shí)間，提升了患者檢查的舒適度；該技術(shù)可以用于人工關(guān)節(jié)置換、口腔植入物、內(nèi)固定金屬支架及放射性粒子等的偽影消除等。綜上所述，MAR算法可以明顯降低CT模擬機(jī)定位掃描時(shí)產(chǎn)生的射束硬化和光子饑餓金屬偽影，同時(shí)可整體提高CT圖像的圖像質(zhì)量。本研究結(jié)果可為CT模擬機(jī)金屬偽影去除、優(yōu)化掃描基礎(chǔ)提供參考，將進(jìn)一步指導(dǎo)醫(yī)學(xué)物理人員開展金屬偽影對(duì)靶區(qū)勾畫和劑量計(jì)算等方面的工作。

本研究也存在一定的局限性，使用單一的模體對(duì)圖像的各項(xiàng)參數(shù)進(jìn)行了評(píng)價(jià)，但在實(shí)際臨床應(yīng)用中，掃描部位、掃描參數(shù)、金屬材料和不同的重建方式均會(huì)產(chǎn)生不一樣的效果。本模體的制作是將金屬物放置于模體表面，無法完全仿真金屬物在體內(nèi)進(jìn)行CT掃描的影響。下一步研究將進(jìn)一步優(yōu)化掃描模體，使之更符合金屬在人體中的情況；同時(shí)可回顧性研究病人的圖像資料，以克服單一模體所得數(shù)據(jù)的片面性問題。