劉志華，卡迪力亞·庫(kù)爾班，李豐軍，嚴(yán)傳波

新疆醫(yī)科大學(xué) a. 公共衛(wèi)生學(xué)院；b. 醫(yī)學(xué)工程技術(shù)學(xué)院，新疆 烏魯木齊 830011

引言

肝包蟲病又稱作肝棘球蚴病，是一種因棘球絳蟲的幼蟲寄生于人體肝臟而導(dǎo)致的寄生蟲病[1-2]。肝包蟲病是牧區(qū)常見的人畜共患寄生蟲疾病，在我國(guó)多見于新疆、西藏、青海、甘肅、寧夏、內(nèi)蒙古、四川等西北農(nóng)牧區(qū), 嚴(yán)重危害人體健康和畜牧業(yè)的發(fā)展, 是嚴(yán)重的公共衛(wèi)生問(wèn)題[3]。肝包蟲病患者在患病初期無(wú)特異性的癥狀及體征，隨著包囊的生長(zhǎng)，患者出現(xiàn)腹部不適、疼痛、惡心、嘔吐等臨床體征表現(xiàn)，當(dāng)囊腫破裂時(shí)，將會(huì)引發(fā)機(jī)體感染甚至發(fā)生過(guò)敏反應(yīng)等并發(fā)癥，對(duì)人體健康造成極大危害。如今，治療肝包蟲病的首選方案是手術(shù)，但無(wú)論何種手術(shù)方式其術(shù)后并發(fā)癥的發(fā)生率都達(dá)到了60%[4-5]，對(duì)人體產(chǎn)生極大危害，然而通過(guò)對(duì)醫(yī)學(xué)圖像進(jìn)行病灶檢測(cè)可以協(xié)助醫(yī)生診斷疾病，有效減少疾病的死亡率。在病灶檢測(cè)方面，國(guó)外研究學(xué)者做了一些研究，提出一個(gè)深度目標(biāo)檢測(cè)網(wǎng)絡(luò)定位視盤和視杯來(lái)進(jìn)行青光眼的篩查和診斷[6]。有研究提出一種改進(jìn)的自動(dòng)腦腫瘤分割和識(shí)別方法，利用人工神經(jīng)網(wǎng)絡(luò)和灰度特征對(duì)腦瘤進(jìn)行分析診斷，區(qū)分良惡性病例[7]。Yahalomi等[8]訓(xùn)練了一種更快的用于目標(biāo)檢測(cè)的機(jī)器視覺(jué)神經(jīng)網(wǎng)絡(luò)，來(lái)識(shí)別和定位前后X射線圖像中的橈骨遠(yuǎn)端骨折，對(duì)骨折的識(shí)別正確率為96%。除此之外，國(guó)內(nèi)許多學(xué)者在肝包蟲病分類[9-10]、分割[11-13]方面做了一些研究，但存在超參數(shù)調(diào)整困難、模型預(yù)測(cè)存在很多假陽(yáng)性、誤分較多、需要人工交互、相鄰器官之間邊界輪廓分割不精確等問(wèn)題。以上僅僅是對(duì)肝包蟲病病灶分類、分割方面的研究，并沒(méi)有對(duì)肝包蟲病在肝臟發(fā)病的位置進(jìn)行統(tǒng)計(jì)學(xué)意義的研究，然而統(tǒng)計(jì)常發(fā)病的位置對(duì)區(qū)分鑒別肝硬化、肝囊腫等病變具有指導(dǎo)意義。因此，本文提出使用深度學(xué)習(xí)的方法對(duì)肝包蟲CT圖像進(jìn)行病灶檢測(cè)，然后統(tǒng)計(jì)并分析病灶區(qū)域在肝臟左右葉發(fā)生的位置，確定了常發(fā)病的位置，有助于輔助醫(yī)生快速診斷疾病，增強(qiáng)診斷準(zhǔn)確性和減少漏診率，為后續(xù)對(duì)病灶區(qū)域進(jìn)行紋理分析提供幫助。

1 方法

1.1 數(shù)據(jù)處理

在放射科醫(yī)師的指導(dǎo)下，本文選取單囊型和多囊型肝包蟲病CT圖像1000張，實(shí)驗(yàn)數(shù)據(jù)均來(lái)自新疆醫(yī)科大學(xué)各附屬醫(yī)院放射科。在獲取肝包蟲病CT圖像后，需要對(duì)圖像進(jìn)行預(yù)處理操作，使圖像盡可能地反映診斷時(shí)所需的圖像信息。本實(shí)驗(yàn)對(duì)圖像進(jìn)行以下預(yù)處理操作：① 對(duì)圖像進(jìn)行灰度歸一化，減少后續(xù)計(jì)算量；② 使用改進(jìn)的中值濾波算法對(duì)圖像進(jìn)行去噪；③ 常用的數(shù)據(jù)增強(qiáng)方法有添加噪聲、旋轉(zhuǎn)、翻轉(zhuǎn)、平移等[14]。圖像翻轉(zhuǎn)是在幾乎不改變圖像數(shù)據(jù)分布的情況下，數(shù)據(jù)翻倍。本文使用垂直翻轉(zhuǎn)、水平翻轉(zhuǎn)的方法[15]對(duì)圖像進(jìn)行數(shù)據(jù)增強(qiáng)，增強(qiáng)后的圖像共3000張。以單囊型肝包蟲病CT圖像為例，圖像預(yù)處理后的效果及圖像增強(qiáng)后的效果如圖1所示。

圖1 圖像預(yù)處理效果

1.2 目標(biāo)檢測(cè)模型

卷積神經(jīng)網(wǎng)絡(luò)基本結(jié)構(gòu)主要包括卷積層、池化層、全連接層及輸出層。卷積神經(jīng)網(wǎng)絡(luò)根據(jù)圖像數(shù)據(jù)逐層抽象提取原始醫(yī)學(xué)圖像中隱藏的數(shù)據(jù)信息，隨著特征提取網(wǎng)絡(luò)層數(shù)越深，提取出來(lái)的特征越豐富，提取特征的效果就越好。Faster RCNN常用的特征提取網(wǎng)絡(luò)主要有ZFNet、VGGNet、ResNet等。ResNet網(wǎng)絡(luò)在2015年首次提出，在ImageNet比賽分類任務(wù)上獲得第一名，之后許多方法基于ResNet網(wǎng)絡(luò)在圖像檢測(cè)、分割等領(lǐng)域得到廣泛應(yīng)用。ResNet網(wǎng)絡(luò)常被用于骨干網(wǎng)絡(luò)對(duì)輸入圖像進(jìn)行特征提取，為解決深度網(wǎng)絡(luò)訓(xùn)練過(guò)程中的退化問(wèn)題，其將殘差思想引入網(wǎng)絡(luò)結(jié)構(gòu)中，在深層網(wǎng)絡(luò)與淺層網(wǎng)絡(luò)間引入跳躍連接增加網(wǎng)絡(luò)層[16]。ResNet101采用bottleneck結(jié)構(gòu)搭建更深的網(wǎng)絡(luò)層，增強(qiáng)學(xué)習(xí)能力，利用1×1卷積減少通道數(shù)以降低計(jì)算量，采用恒等快捷連接實(shí)現(xiàn)跳躍連接，加快收斂速度。因此，本文將ResNet101網(wǎng)絡(luò)代替原始的VGG16網(wǎng)絡(luò)作為肝包蟲病CT圖像檢測(cè)的特征提取網(wǎng)絡(luò)，對(duì)肝包蟲病CT圖像進(jìn)行病灶檢測(cè)。bottleneck結(jié)構(gòu)圖如圖2所示。

圖2 殘差塊結(jié)構(gòu)圖

1.3 建立肝臟分區(qū)標(biāo)準(zhǔn)

在影像科醫(yī)師的指導(dǎo)下，找到肝包蟲病CT圖像的八個(gè)解剖標(biāo)志：下腔靜脈、肝右靜脈、肝中靜脈、肝左靜脈等，以肝左靜脈和下腔靜脈連線作為縱軸，同時(shí)建立坐標(biāo)系，以肝臟八段劃分法標(biāo)準(zhǔn)建立肝臟分區(qū)，并找到每一區(qū)域的質(zhì)心點(diǎn)。圖中紅色縱軸以肝左靜脈為界限分為Ⅰ、Ⅱ、Ⅲ段和Ⅳ段，紫色射線以肝中靜脈為界限分為Ⅳ段和Ⅴ/Ⅷ段，綠色射線以門靜脈為界限分為Ⅴ/Ⅷ段和Ⅵ/Ⅶ段，坐標(biāo)系中紅點(diǎn)為所在區(qū)域質(zhì)心點(diǎn)。肝臟分區(qū)圖如圖3所示，建立坐標(biāo)系分區(qū)圖如圖4所示，肝臟分區(qū)標(biāo)準(zhǔn)信息如表1所示。

圖3 肝臟分區(qū)圖

圖4 坐標(biāo)系圖

表1 肝臟分區(qū)標(biāo)準(zhǔn)信息

1.4 SCP方法

關(guān)于目標(biāo)檢測(cè)模型的定位問(wèn)題，不少的研究者進(jìn)行過(guò)相關(guān)研究，如一些學(xué)者首先使用Canny算子識(shí)別目標(biāo)區(qū)域的輪廓，再使用Hough梯度圓變換算法、最小二乘法、形心偏差法等方法定位目標(biāo)，該方法取得了不錯(cuò)的效果，但定位結(jié)果存在一定的偏差[17-19]。蔣蕓等[20]利用YOLO算法對(duì)眼底圖像視盤定位，該方法檢測(cè)效果較好，但對(duì)目標(biāo)的標(biāo)注花費(fèi)大量時(shí)間。李傳朋[21]通過(guò)計(jì)算目標(biāo)檢測(cè)框中心的方法實(shí)現(xiàn)定位，然而該方法對(duì)位置的統(tǒng)計(jì)存在一定的誤差。由于肝包蟲病病灶大小不一，拍攝CT圖像時(shí)每個(gè)人的體位不同，僅僅通過(guò)計(jì)算目標(biāo)檢測(cè)框的中心不足以說(shuō)明病灶在肝臟每一區(qū)域的位置。針對(duì)上述問(wèn)題本文提出一種Sobel算子和連通區(qū)域與極坐標(biāo)融合的方法（Sobel Operator and Connected Component and Polar Coordinate Fusion Method，SCP）對(duì)病灶在肝臟區(qū)域的位置進(jìn)行統(tǒng)計(jì)。本文基于CT圖像器官輪廓的質(zhì)心點(diǎn)和病灶的質(zhì)心點(diǎn)構(gòu)建極坐標(biāo)系，以水平方向?yàn)閱挝皇枷蛄?，利用空間余弦定理計(jì)算目標(biāo)向量和單位始向量的夾角，單位始向量的夾角即為極角，目標(biāo)向量的模即極徑。

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

2.1 目標(biāo)檢測(cè)結(jié)果

本實(shí)驗(yàn)?zāi)繕?biāo)檢測(cè)模型的檢測(cè)結(jié)果使用平均精度、均值平均精度等評(píng)價(jià)指標(biāo)進(jìn)行模型性能評(píng)價(jià)，模型性能評(píng)價(jià)結(jié)果如表2所示。圖5所示為模型對(duì)兩類肝包蟲病CT圖像的檢測(cè)結(jié)果。圖5a為單囊型肝包蟲病CT圖像檢測(cè)結(jié)果，圖5b多囊型肝包蟲病CT圖像檢測(cè)結(jié)果。

表2 模型對(duì)肝包蟲病CT圖像的性能評(píng)價(jià)

圖5 目標(biāo)檢測(cè)模型檢測(cè)結(jié)果

由表2可知，對(duì)于兩種類型的肝包蟲病CT圖像病灶檢測(cè)，基于ResNet101的特征提取網(wǎng)絡(luò)比基于VGG16的特征提取網(wǎng)絡(luò)在AP、mAP指標(biāo)上都取得了較好的效果，說(shuō)明更深層的網(wǎng)絡(luò)在肝包蟲病檢測(cè)任務(wù)上具有更優(yōu)秀的特征表達(dá)能力?；赗esNet101特征提取網(wǎng)絡(luò)的兩種類型肝包蟲病CT圖像的平均識(shí)別率均有所提高，然而對(duì)于多囊型肝包蟲病CT圖像的平均識(shí)別率提高了3.5%，由0.870提高至0.905，說(shuō)明基于ResNet101的特征提取網(wǎng)絡(luò)對(duì)多囊型肝包蟲病CT圖像的檢測(cè)識(shí)別效果較好。由表2可知，基于ResNet101特征提取網(wǎng)絡(luò)的漏檢率、錯(cuò)檢率均低于基于VGG16的特征提取網(wǎng)絡(luò)，說(shuō)明基于ResNet101的特征提取網(wǎng)絡(luò)病灶檢測(cè)結(jié)果較好，對(duì)病灶檢測(cè)時(shí)存在的錯(cuò)檢及漏檢率較小，使用深度學(xué)習(xí)方法對(duì)肝包蟲病的病灶檢測(cè)具有一定的可行性，在測(cè)試集上取得了良好的效果。

2.2 SCP方法結(jié)果

使用計(jì)算機(jī)估計(jì)某病常發(fā)位置，必須以器官CT輪廓為參考建立坐標(biāo)系，常見的坐標(biāo)系有直角坐標(biāo)系、極坐標(biāo)系，又由于從醫(yī)學(xué)角度對(duì)器官分區(qū)時(shí)某分區(qū)可能跨直角坐標(biāo)系的兩個(gè)象限，不易估計(jì)病灶位置。因此，本實(shí)驗(yàn)采用極坐標(biāo)系來(lái)描述病灶在器官中的位置。為了構(gòu)建極參考系，必須確定極點(diǎn)和極徑，本實(shí)驗(yàn)采用器官CT輪廓質(zhì)心點(diǎn)作為極點(diǎn)，病灶的質(zhì)心點(diǎn)作為終點(diǎn)。本實(shí)驗(yàn)基于這兩個(gè)質(zhì)心點(diǎn)構(gòu)建目標(biāo)向量，以水平方向?yàn)閱挝皇枷蛄?，利用空間余弦定理計(jì)算目標(biāo)向量和單位始向量的夾角，即為極角，目標(biāo)向量模即極徑。

本實(shí)驗(yàn)使用Sobel算子提取CT圖像器官輪廓，并使用連通域方法得到CT圖像器官輪廓的質(zhì)心點(diǎn)，圖6所示為Sobel算子提取CT圖像器官輪廓邊緣的結(jié)果，圖7所示為使用連通域方法得到的CT圖像器官輪廓質(zhì)心點(diǎn)位置，圖8所示為病灶質(zhì)心點(diǎn)與CT圖像器官輪廓質(zhì)心點(diǎn)在CT圖像中的位置，表3所示為部分圖像器官輪廓質(zhì)心點(diǎn)坐標(biāo)。

圖6 Sobel算子邊緣提取結(jié)果

圖7 CT圖像器官輪廓質(zhì)心點(diǎn)

圖8 質(zhì)心點(diǎn)位置圖

表3 多例圖片名稱對(duì)應(yīng)的圖像質(zhì)心點(diǎn)

基于兩個(gè)質(zhì)心點(diǎn)建立極坐標(biāo)系，如圖9所示為建立的極徑極角圖，根據(jù)肝臟分區(qū)標(biāo)準(zhǔn)信息，可以得到各區(qū)域圖像數(shù)量以及極徑極角范圍。本研究使用秩和檢驗(yàn)方法對(duì)不同區(qū)域的極徑和極角進(jìn)行統(tǒng)計(jì)學(xué)意義的研究，表4所示為肝包蟲病CT圖像在各區(qū)域的極坐標(biāo)信息統(tǒng)計(jì)結(jié)果。

圖9 極徑極角圖

表4 CT圖像在各區(qū)域的極坐標(biāo)信息統(tǒng)計(jì)結(jié)果

本研究對(duì)測(cè)試集中的數(shù)據(jù)使用秩和檢驗(yàn)方法統(tǒng)計(jì)得到了表4所示的單囊型和多囊型肝包蟲病的極徑和極角信息在肝臟分區(qū)中的統(tǒng)計(jì)分析結(jié)果，從表4可知，H=402.25，P＜0.001，四個(gè)區(qū)域的極角之間具有顯著性差異，四個(gè)區(qū)域的極角大小不全相同；H=39.69，P＜0.001，四個(gè)區(qū)域的極徑之間具有顯著性差異，四個(gè)區(qū)域的極徑大小不全相同。

2.3 統(tǒng)計(jì)肝包蟲病常發(fā)位置

基于肝包蟲病CT圖像在各區(qū)域的極徑極角范圍，可以得到單囊型和多囊型肝包蟲病CT圖像在各區(qū)域的數(shù)量，從而可以得到肝包蟲病在各區(qū)域的百分比。根據(jù)在各區(qū)域占比大小可以得出肝包蟲病在肝臟區(qū)域常發(fā)位置。表5所示為肝包蟲病CT圖像各區(qū)域數(shù)量一般統(tǒng)計(jì)結(jié)果，表6所示為肝包蟲病CT圖像各區(qū)域百分比統(tǒng)計(jì)結(jié)果。

表5 肝包蟲病CT圖像各區(qū)域數(shù)量一般統(tǒng)計(jì)結(jié)果[n（%）]

表6 肝包蟲病CT圖像各區(qū)域占比統(tǒng)計(jì)結(jié)果

本研究使用卡方檢驗(yàn)方法對(duì)得到的極坐標(biāo)進(jìn)行統(tǒng)計(jì)，得出表5所示的單囊型和多囊型肝包蟲病CT圖像在肝臟各區(qū)域的數(shù)量，由表5可知，單囊型肝包蟲病在Ⅱ區(qū)、Ⅴ區(qū)/Ⅷ區(qū)的數(shù)量多于多囊型肝包蟲病，而多囊型肝包蟲病在Ⅳ區(qū)、Ⅵ/Ⅶ區(qū)的數(shù)量多于單囊型肝包蟲病?；诒?統(tǒng)計(jì)結(jié)果，本文使用卡方檢驗(yàn)方法得到表6所示的單囊型和多囊型肝包蟲病CT圖像在肝臟各個(gè)區(qū)域的占比，由表6可知，χ2=20.896，P＜0.001，差異有統(tǒng)計(jì)學(xué)意義，單囊型和多囊型肝包蟲病在肝臟各個(gè)區(qū)域的占比不全相同。經(jīng)過(guò)兩兩比較后，單囊型肝包蟲病在Ⅳ區(qū)和Ⅴ/Ⅷ區(qū)、Ⅴ/Ⅷ區(qū)和Ⅵ/Ⅶ區(qū)的占比有差異，在Ⅱ區(qū)和Ⅳ區(qū)、Ⅴ/Ⅷ區(qū)、Ⅵ/Ⅶ區(qū)的占比無(wú)差異；多囊型肝包蟲病在Ⅳ區(qū)和Ⅴ/Ⅷ區(qū)、Ⅴ/Ⅷ區(qū)和Ⅵ/Ⅶ區(qū)的占比有差異，在Ⅱ區(qū)和Ⅳ區(qū)、Ⅴ/Ⅷ區(qū)、Ⅵ/Ⅶ區(qū)的占比無(wú)差異，且單囊型和多囊型肝包蟲病在肝臟的Ⅳ區(qū)占比較高，Ⅵ/Ⅶ區(qū)次之。又由于單囊型肝包蟲病囊壁發(fā)生破裂后，會(huì)逐步增長(zhǎng)為多囊型肝包蟲病，因此在單囊型肝包蟲病囊壁無(wú)破裂的情況下，單囊型和多囊型肝包蟲病常發(fā)病位置分別在肝臟的Ⅳ區(qū)、Ⅵ/Ⅶ區(qū)，Ⅴ區(qū)/Ⅷ區(qū)次之，在單囊型肝包蟲病破裂的情況下，肝包蟲病的常發(fā)位置主要在Ⅳ區(qū)、Ⅵ/Ⅶ區(qū)。

3 討論

本研究使用Faster RCNN模型和SCP方法統(tǒng)計(jì)肝包蟲病在肝臟區(qū)域常發(fā)病位置，基于ResNet101的目標(biāo)檢測(cè)模型對(duì)肝包蟲病病灶的檢測(cè)準(zhǔn)確率達(dá)到了89.6%，說(shuō)明該模型能夠有效地檢測(cè)出病灶從而可以得到病灶的質(zhì)心點(diǎn)坐標(biāo)；然后使用SCP方法統(tǒng)計(jì)出肝包蟲病常發(fā)病位置主要發(fā)生在肝臟的Ⅳ區(qū)、Ⅵ/Ⅶ區(qū)。一些學(xué)者在目標(biāo)定位方面使用邊緣檢測(cè)、目標(biāo)檢測(cè)等方法來(lái)確定目標(biāo)位置，然而這些研究方法存在一些問(wèn)題，如耗費(fèi)時(shí)間、計(jì)算量大、定位結(jié)果有誤差等，因此，以往的研究方法對(duì)目標(biāo)位置的定位具有一定的局限性。本研究將目標(biāo)檢測(cè)方法和極坐標(biāo)系相結(jié)合來(lái)確定病灶的發(fā)病位置，相比于使用直角坐標(biāo)系中的坐標(biāo)點(diǎn)確定常發(fā)位置更具有可靠性。除此之外，該方法使用極徑和極角確定病灶常發(fā)病位置結(jié)果更加精確，并且具有統(tǒng)計(jì)學(xué)意義。因此，本文方法相比于前人的研究具有一定的創(chuàng)新性。

4 結(jié)論

本文將Faster RCNN網(wǎng)絡(luò)應(yīng)用于肝包蟲病CT圖像的目標(biāo)檢測(cè)，并使用SCP方法統(tǒng)計(jì)肝包蟲病在肝臟的發(fā)病位置，該方法可以輔助醫(yī)生診斷疾病，減少漏檢、錯(cuò)檢的發(fā)生，同時(shí)能夠?yàn)獒t(yī)生診斷提供導(dǎo)向，更高效地診斷疾病，做到早發(fā)現(xiàn)、早診斷、早治療。本研究使用基于ResNet101網(wǎng)絡(luò)的目標(biāo)檢測(cè)模型能夠較好地檢測(cè)肝包蟲病病灶，模型檢測(cè)準(zhǔn)確率達(dá)到了89.6%，并使用 SCP 方法統(tǒng)計(jì)出肝包蟲病常發(fā)位置在肝臟的Ⅳ區(qū)、Ⅵ區(qū)/Ⅶ區(qū)。然而該方法對(duì)肝包蟲病CT圖像器官輪廓的質(zhì)心點(diǎn)位置確定存在一定的偏差，本文后續(xù)將繼續(xù)研究對(duì)器官輪廓質(zhì)心點(diǎn)的精確位置定位，通過(guò)大量數(shù)據(jù)對(duì)CT圖像器官輪廓的質(zhì)心點(diǎn)進(jìn)行聚類，以求得精確的質(zhì)心點(diǎn)位置。除此之外，該方法具有一定的普適性，可以將肝包蟲病常發(fā)病位置估計(jì)方法應(yīng)用于其他疾病的研究，通過(guò)確定疾病常發(fā)病的位置，對(duì)早期發(fā)現(xiàn)疾病具有一定的臨床價(jià)值，即能減輕患者病痛，又能節(jié)省醫(yī)療資源。