周錫明 王禹鑫 張 弛 任 燁 蔡 彪

醫(yī)用X射線計(jì)算機(jī)斷層掃描(X-ray computed tomography,CT)系統(tǒng)從20世紀(jì)70年代發(fā)明至今,經(jīng)歷了多個(gè)快速發(fā)展的階段,從單排步進(jìn)式頭部專用CT進(jìn)化為多排螺旋全身CT,掃描時(shí)間從數(shù)分鐘進(jìn)化到數(shù)秒鐘[1].CT的硬件包括球管、探測(cè)器、動(dòng)力系統(tǒng)、控制系統(tǒng)、數(shù)據(jù)傳輸和重建等均有了巨大的技術(shù)進(jìn)步,其算法和軟件也在不停的發(fā)展之中.重建算法經(jīng)歷了代數(shù)重建技術(shù)(algebraic reconstruction technique,ART)、濾波反投影(filter back projection,FBP)算法和迭代重建(iterative reconstruction,IR)算法,到目前的人工智能(artificial intelligence,AI)重建算法[2-4].掃描方式上也有多種創(chuàng)新發(fā)明,從最初的單排步進(jìn)軸掃方式,發(fā)展到現(xiàn)在的大錐束螺旋、心臟門控掃描及雙能掃描等[5].在數(shù)據(jù)處理的算法上,也針對(duì)不同的圖像偽影的成因,發(fā)展了各種偽影消除算法,包括骨偽影、散射偽影、金屬偽影以及低信號(hào)偽影校正算法等等[6-7].然而,CT圖像質(zhì)量和系統(tǒng)性能的評(píng)定標(biāo)準(zhǔn)卻相對(duì)更新緩慢.

目前,國家標(biāo)準(zhǔn)"X射線計(jì)算機(jī)體層攝影設(shè)備通用技術(shù)條件"(YY/T 0310-2015)[8]中的圖像質(zhì)量測(cè)定方法只是針對(duì)CT掃描系統(tǒng)的基本圖像參數(shù)進(jìn)行了定義,并未針對(duì)CT掃描中的高級(jí)功能進(jìn)行圖像質(zhì)量的評(píng)估.在國家標(biāo)準(zhǔn)"X射線計(jì)算機(jī)體層攝影設(shè)備通用技術(shù)條件"中,對(duì)圖像質(zhì)量檢測(cè)部分中推薦的模體,包括業(yè)界廣泛使用的Catphan模體均為圓柱體,無法用于評(píng)估CT掃描系統(tǒng)中自動(dòng)曝光控制(automatic exposure control,AEC)的有效性.對(duì)于低劑量掃描,目前的研究多集中在對(duì)圖像質(zhì)量評(píng)價(jià)的方法上,而使用的模體大都屬于小型模體[9].在CT的輻射劑量低到一定程度后,圖像中會(huì)出現(xiàn)光子不足偽影,因此小型模體并不能模擬CT系統(tǒng)在低劑量下掃描患者的真實(shí)情況.為此,本研究針對(duì)AEC、低劑量和雙能掃描的3個(gè)高級(jí)掃描方式設(shè)計(jì)一種橢圓模體,為CT用戶和監(jiān)管部門對(duì)CT系統(tǒng)的高級(jí)掃描方式的性能評(píng)估提供參考.

1 CT系統(tǒng)高級(jí)掃描技術(shù)

1.1 CT掃描系統(tǒng)AEC

AEC是針對(duì)患者的體型,在掃描過程中動(dòng)態(tài)改變曝光電流的一種技術(shù).如果被掃描部位是類似于橢圓形狀,如肩膀或者胯部,在最終的圖像噪聲中,來自于長軸方向的掃描數(shù)據(jù)中的噪聲將會(huì)成為主要來源.被掃描物體在病床移動(dòng)方向上(Z軸方向)具有較大的尺寸變化,當(dāng)頭頸部或胸腹部連續(xù)掃描時(shí),如果球管的曝光電流保持恒定,則在Z軸方向上圖像的噪聲將出現(xiàn)較大的不均勻性.AEC是根據(jù)被掃描物體形狀(一般通過定位像來估計(jì))來計(jì)算出在各個(gè)方向上的衰減長度,然后根據(jù)長度來決定對(duì)應(yīng)方向上的電流強(qiáng)度.在衰減長度大的方向增加電流強(qiáng)度;相反再衰減長度小的方向降低電流強(qiáng)度.

1.2 低劑量掃描

CT等醫(yī)療放射成像設(shè)備目前已經(jīng)成為公眾受到的輻射的主要來源,約有50%的輻射劑量來自于醫(yī)療成像設(shè)備,并且有逐年上升的趨勢(shì)[10].醫(yī)療成像設(shè)備對(duì)于臨床診斷具有其他設(shè)備無法替代的意義,因此在醫(yī)院診療中,尤其在兒科中,影像設(shè)備中盡可能的低劑量原則(as low as reasonably achievable,ALARA)的概念也越來越受到重視[11].隨著新的圖像重建技術(shù)的出現(xiàn),包括AI在醫(yī)療診斷上的進(jìn)展,在肺部篩查等特殊的應(yīng)用領(lǐng)域中超低劑量的掃描方式開始出現(xiàn)[12-14].這些掃描方式不僅僅是單純地把掃描協(xié)議中的劑量降低,而是與高級(jí)的預(yù)處理、圖像重建以及圖像識(shí)別的算法相結(jié)合,在給患者盡可能低的輻射劑量的前提下,針對(duì)特定的應(yīng)用提供有用的診斷信息.

1.3 雙能掃描

雙能CT掃描的概念在CT被發(fā)明后不久被提出[15].雙能CT利用不同物質(zhì)對(duì)不同能量的X射線的吸收系數(shù)有差別原理,對(duì)同一掃描物體進(jìn)行不同X射線能量下的掃描,從而能達(dá)到區(qū)分被掃描物體中不同物質(zhì)目的.目前,各主要品牌的雙能掃描實(shí)現(xiàn)方式有很多種,包括西門子的雙源雙探測(cè)器方式、GE的快速切換方式和飛利浦的雙層探測(cè)器方式等[16-18].各品牌的雙能掃描方式各有優(yōu)缺點(diǎn),雙能CT掃描的最終圖像一般都會(huì)重建為單能圖,物質(zhì)分解圖,或者原子序數(shù)圖等表現(xiàn)形式.這些雙能圖像的準(zhǔn)確性目前并無統(tǒng)一的國家標(biāo)準(zhǔn)來進(jìn)行評(píng)估.雖然有些模體的供應(yīng)商提供一些已知化學(xué)成分的模體,如Catphan700中的CTP682就提供了水、有機(jī)玻璃、骨、特氟龍等材料,但是這些物質(zhì)的化學(xué)成分和一般雙能重建使用的基物質(zhì)對(duì)(水和碘)相差較大,導(dǎo)致理論值的計(jì)算與產(chǎn)品中采用的掃描技術(shù)相關(guān).針對(duì)雙能掃描,提出一種簡單易行的辦法對(duì)其物質(zhì)分解的準(zhǔn)確性進(jìn)行評(píng)估.

2 橢圓模體設(shè)計(jì)

2.1 橢圓模體結(jié)構(gòu)

根據(jù)自動(dòng)曝光控制和低劑量掃描方式的特點(diǎn),設(shè)計(jì)橢圓形的模體.模體的長軸為35 cm,短軸為25 cm,材料為有機(jī)玻璃聚甲基丙烯酸甲酯(polymethyl methacrylate,PMMA).橢圓模體的中心為20 cm空心圓柱,便于放置廠商自帶的質(zhì)量保證(quality assurance,QA)模體或Catphan模體進(jìn)行圖像質(zhì)量指標(biāo)的測(cè)量.在長軸的兩端,各設(shè)置2 cm直徑的圓孔,放置各種自定義材料做成的圓棒,如鋁棒、鐵棒及特氟龍棒等,以方便放入裝有各種溶液的試管.圖像質(zhì)量指標(biāo)的測(cè)量由Catphan500模體完成.橢圓模體設(shè)計(jì)見圖1.

圖1 橢圓模體設(shè)計(jì)示意圖

2.2 橢圓模體性能

模體設(shè)計(jì)成橢圓形狀,其目的是為了使AEC在掃描中起作用,且橢圓模體更符合臨床上患者的體型特征,能更好地模擬低劑量掃描下對(duì)圖像質(zhì)量的實(shí)際影響.

對(duì)于雙能掃描的測(cè)試,通過使用配制的1%、2%、4%和10%的4種不同質(zhì)量濃度的碘化鉀溶液進(jìn)來行測(cè)量,而不是臨床上的CT增強(qiáng)劑溶液,其原因是:①CT增強(qiáng)劑雖然主要成分的分子式是已知的,但是溶液中其他物質(zhì)的成分往往生產(chǎn)廠商并不提供;②商業(yè)化CT增強(qiáng)劑產(chǎn)品中,有效成分的濃度一般都有誤差,有的多達(dá)10%.使用CT增強(qiáng)劑溶液研究雙能CT掃描的準(zhǔn)確性時(shí),無法準(zhǔn)確計(jì)算理論值.而碘化鉀并非管制化學(xué)品,易溶于水,純度一般都在99%以上,且碘化鉀溶液的密度和溫度的關(guān)系也容易計(jì)算[19].使用碘化鉀來配置含碘溶液,其誤差來源很少.若使用一次配制的碘化鉀溶液>100 ml,很容易將碘化鉀溶液的精度控制在<0.001 g/ml.計(jì)算4種溶液在溫度t(℃)下的質(zhì)量密度ρ(g/cm3)為公式1:

式中m為每單位重量水中的物質(zhì)摩爾濃度,單位為mol/kg.

在以水和碘這兩種物質(zhì)為基物質(zhì)的4種濃度溶液,通過物質(zhì)分解得到的理論密度見表1.

表1 不同質(zhì)量濃度下碘化鉀水溶液25 ℃下的質(zhì)量密度和水基碘基物質(zhì)分解密度(mg/cm3)

2.3 橢圓模體測(cè)試方法

選擇3個(gè)不同品牌的CT產(chǎn)品(品牌A、B、C)進(jìn)行圖像質(zhì)量測(cè)試,分析在AEC、低劑量和雙能掃描3種高級(jí)掃描模式下橢圓模體對(duì)圖像質(zhì)量的影響.3個(gè)不同品牌CT分別覆蓋高中低端的CT產(chǎn)品,品牌A為高端256排的CT,品牌B為中端32排CT,品牌C為低端16排CT.其中,品牌A和B的CT產(chǎn)品中有雙能功能,品牌C的16排CT不具備雙能.

測(cè)試的圖像質(zhì)量指標(biāo)選用國家標(biāo)準(zhǔn)"X射線計(jì)算機(jī)體層攝影設(shè)備通用技術(shù)條件"(YY/T 0310-2015)中規(guī)定的CT值的準(zhǔn)確性、均勻性以及層厚,測(cè)試方法均遵守標(biāo)準(zhǔn)中的規(guī)定:①水的CT值均勻性≤±4 HU;②空氣的準(zhǔn)確性為(-1000±10)HU,水的準(zhǔn)確性為(0±4)HU;③切片層厚>2 mm時(shí)為1.0 mm,1~2 mm時(shí)為±50%,<1 mm時(shí)為±0.5 mm.對(duì)于雙能掃描,使用碘化鉀溶液來測(cè)試品牌A和B的雙能掃描中經(jīng)過物質(zhì)分解以后的水基和碘基圖的密度值,并與理論值進(jìn)行比較.

3 橢圓模體測(cè)試結(jié)果

3.1 AEC測(cè)試

自動(dòng)曝光測(cè)試采用常規(guī)的胸部掃描協(xié)議,管電壓120 kV,管電流130 mA,轉(zhuǎn)速0.75 s/360°,螺距1.0.每臺(tái)CT測(cè)試4次并取平均值.在測(cè)量過程中空氣的測(cè)量選取遠(yuǎn)離模體的感興趣區(qū)域(region of interest,ROI),水模的中心和3、6、9、12點(diǎn)位置的ROI的選取遵循國家標(biāo)準(zhǔn)(YY/T 0310-2015)中的規(guī)定,即在距離模體邊緣大約1 cm處,選取一個(gè)直徑大約是2 cm的圓形ROI進(jìn)行測(cè)量,中心ROI與周圍4個(gè)ROI的CT值之差的最大值為均勻性結(jié)果(見圖2、圖3和圖4).

實(shí)驗(yàn)結(jié)果顯示,在AEC打開并使用橢圓模體的情況下,各品牌的CT值準(zhǔn)確性、均勻性和切片層厚均在國家標(biāo)準(zhǔn)[8]的規(guī)定范圍內(nèi).需要指出的是雖然品牌C的水CT值在3、6、9點(diǎn)位置超出了4 HU的規(guī)定范圍,但是水的CT值準(zhǔn)確性是在圖像中心測(cè)量的,因此品牌C的結(jié)果同樣符合標(biāo)準(zhǔn).

圖2 各ROI選取的位置示意圖

圖3 不同品牌的CT系統(tǒng)在AEC打開時(shí)的CT值準(zhǔn)確性和均勻性測(cè)試結(jié)果

圖4 不同品牌的CT系統(tǒng)在AEC打開時(shí)的層厚測(cè)試結(jié)果

3.2 低劑量測(cè)試

低劑量圖像質(zhì)量測(cè)試中使用了肺部篩查協(xié)議中常用的掃描協(xié)議:管電壓120 kV,管電流30 mA,轉(zhuǎn)速0.8 s/360°,螺距1.35[13];同樣每臺(tái)CT測(cè)試4次并取平均數(shù)據(jù).

在低劑量打開并使用橢圓模體的情況下,各品牌產(chǎn)品的CT值準(zhǔn)確性、均勻性和切片層厚均在國家標(biāo)準(zhǔn)[8]的規(guī)定范圍內(nèi).因此,可以表明橢圓模體在低劑量掃描模式下,對(duì)圖像質(zhì)量研究的影響不大,見圖5和圖6.

圖5 不同品牌的CT產(chǎn)品在低劑量掃描時(shí)的CT值準(zhǔn)確性和均勻性測(cè)試結(jié)果

圖6 不同品牌的CT產(chǎn)品在低劑量掃描時(shí)的層厚測(cè)試結(jié)果

3.3 雙能測(cè)試結(jié)果

雙能掃描協(xié)議選取了品牌A和B提供的頭部協(xié)議,同樣每臺(tái)CT測(cè)試4次并取平均.在使用橢圓模體時(shí),兩個(gè)品牌的雙能掃描仍然能滿足圖像準(zhǔn)確性、均勻性和層厚的要求.

在兩個(gè)品牌的CT產(chǎn)品上分別測(cè)量的碘化鉀溶液的物質(zhì)分解密度值.兩個(gè)CT產(chǎn)品的水/碘密度和理論值都存在一定程度的偏差,品牌A的結(jié)果在低濃度時(shí)碘密度過高,而在高濃度時(shí)水密度過高,表明其物質(zhì)分解和碘化鉀濃度相關(guān);品牌C的結(jié)果與理論值更吻合,但是碘密度總體偏高,存在一個(gè)系統(tǒng)偏差.測(cè)試結(jié)果證明,不同品牌的CT產(chǎn)品之間的物質(zhì)分解結(jié)果存在差異,應(yīng)該與硬件掃描方式和算法有關(guān)系.品牌A和品牌C的物質(zhì)分解結(jié)果和理論值的誤差均在±20%以內(nèi)(見圖7、圖8和圖9).

圖7 不同品牌的CT產(chǎn)品在雙能掃描時(shí)的CT值準(zhǔn)確性和均勻性測(cè)試結(jié)果

圖8 不同品牌的CT產(chǎn)品在雙能掃描時(shí)的層厚測(cè)試結(jié)果

圖9 品牌A和品牌C對(duì)不同質(zhì)量密度的碘化鉀物質(zhì)分解結(jié)果

4 討論

本研究測(cè)試結(jié)果顯示,在標(biāo)準(zhǔn)測(cè)試模體的基礎(chǔ)上添加橢圓形模體,CT值的準(zhǔn)確性、均勻性及層厚的3項(xiàng)指標(biāo),在AEC、低劑量或者雙能掃描下3個(gè)不同品牌的CT產(chǎn)品仍能滿足國家標(biāo)準(zhǔn).使用低劑量時(shí),3個(gè)CT產(chǎn)品的層厚指標(biāo)都有不同程度的下降,表明在低劑量掃描時(shí)各產(chǎn)品的數(shù)據(jù)處理算法中均涉及探測(cè)器Z方向(沿病床的方向)的低通濾波,可能是為了降低低劑量掃描時(shí)的噪聲和偽影.

以雙能掃描物質(zhì)分解準(zhǔn)確度測(cè)量為例,本研究設(shè)計(jì)的橢圓體?？梢蕴峁?duì)高級(jí)掃描功能的性能指標(biāo)在不同CT產(chǎn)品之間的橫向比較提供一個(gè)標(biāo)尺.利用橢圓體?？梢苑謩e測(cè)量在關(guān)閉和打開AEC時(shí),達(dá)到相同的圖像質(zhì)量(相同的信噪比)時(shí)需要的輻射劑量,從而在不同的品牌和不同的產(chǎn)品間提供固定的參考依據(jù).

5 結(jié)論

本研究設(shè)計(jì)的一種新橢圓模體作為目前CT系統(tǒng)圖像質(zhì)量測(cè)試模體的補(bǔ)充,以方便研究CT系統(tǒng)在AEC、低劑量或者雙能掃描中的圖像質(zhì)量.在不同掃描技術(shù)下的實(shí)驗(yàn)驗(yàn)證了該橢圓模體并不會(huì)造成CT圖像質(zhì)量指標(biāo)超出國家硬性規(guī)定的標(biāo)準(zhǔn)范圍.簡單易行的雙能CT物質(zhì)分解準(zhǔn)確性的測(cè)量辦法通過實(shí)驗(yàn)進(jìn)行了驗(yàn)證.

對(duì)于各產(chǎn)品中廠商自行規(guī)定的噪聲,高低對(duì)比度分辨率其他圖像質(zhì)量指標(biāo)有待進(jìn)一步研究,這些指標(biāo)會(huì)隨著掃描方式的變化而改變.本研究設(shè)計(jì)的橢圓模體,可以通過與Catphan等模體的配合應(yīng)用于對(duì)各種高級(jí)掃描功能中的各種圖像質(zhì)量參數(shù)測(cè)量.