李振鵬，陳國(guó)明，孫艷寧，魏小花，李云蓉

（廣東第二師范學(xué)院計(jì)算機(jī)學(xué)院，廣州 510303）

0 引言

青光眼是一種常見(jiàn)的眼科疾病，患者通常表現(xiàn)為不同程度的視力下降和視野缺失，嚴(yán)重者甚至?xí)适б暳?。?jù)統(tǒng)計(jì)，2020 年全世界患青光眼的中老年人數(shù)約為7600 萬(wàn)，在將來(lái)還會(huì)持續(xù)增長(zhǎng)。青光眼的早期癥狀不明顯，有的患者表現(xiàn)為無(wú)癥狀，甚至發(fā)展到中后期才被確診。因此早發(fā)現(xiàn)、早診斷、早治療是治療青光眼的關(guān)鍵。

眼底照相是青光眼早期檢查的主要項(xiàng)目，通過(guò)觀察視盤(pán)與視杯的比例，即杯盤(pán)比是否有異常增大作為診斷的重要指標(biāo)之一。這種診斷方式主觀性較強(qiáng)，通常依靠眼科醫(yī)生的經(jīng)驗(yàn)判斷，如果杯盤(pán)比增大不明顯，則有可能造成誤診或漏診的情況。隨著深度學(xué)習(xí)的發(fā)展，計(jì)算機(jī)輔助診斷（CAD）逐漸應(yīng)用于臨床，有助于醫(yī)生對(duì)患者病情做出更加正確的判斷，提高工作質(zhì)量與效率。

為了降低青光眼早期誤診漏診的可能性，主要思路是提取青光眼病變特征，將正常與青光眼的眼底圖正確分類(lèi)。卷積神經(jīng)網(wǎng)絡(luò)（CNN）是CAD 系統(tǒng)中被廣泛應(yīng)用的深度學(xué)習(xí)方法，在進(jìn)行肺結(jié)節(jié)與前列腺癌的診斷中有著較高的準(zhǔn)確率。Harangi提出了一個(gè)深度卷積神經(jīng)網(wǎng)絡(luò)，結(jié)合另外四種卷積神經(jīng)網(wǎng)絡(luò)作為分類(lèi)層的輸出，實(shí)現(xiàn)了將皮膚鏡圖像分類(lèi)為黑素瘤、痣和脂溢性角化病三類(lèi)。Ozawa等基于GoogleNet構(gòu)建了一個(gè)基于CNN 的CAD 系統(tǒng)，并利用該系統(tǒng)對(duì)結(jié)腸鏡圖像進(jìn)行訓(xùn)練，在分析潰瘍性結(jié)腸炎患者的疾病嚴(yán)重程度上表現(xiàn)穩(wěn)健，有望于投入臨床幫助醫(yī)生進(jìn)行疾病的診斷。Chan 等對(duì)三種CNN網(wǎng)絡(luò)進(jìn)行遷移學(xué)習(xí)以分類(lèi)甲狀腺腫瘤，并進(jìn)行了訓(xùn)練和測(cè)試。結(jié)果表明，CNN 網(wǎng)絡(luò)的診斷效果比醫(yī)生要好，訓(xùn)練的CNN 網(wǎng)絡(luò)能提高大多數(shù)甲狀腺癌的診斷準(zhǔn)確性。Chougrad等提出了一種基于CNN 的乳腺癌分類(lèi)器，通過(guò)遷移學(xué)習(xí)克服了數(shù)據(jù)量偏少所帶來(lái)的過(guò)擬合問(wèn)題，達(dá)到了98.23%的準(zhǔn)確率，并能夠判斷腫瘤是惡性還是良性。Patel 等利用兩種胃腸息肉的視頻序列訓(xùn)練CNN 模型，在判斷息肉類(lèi)型是腺瘤性的還是增生性的結(jié)果與專(zhuān)科醫(yī)生的報(bào)告吻合較好。

醫(yī)學(xué)影像數(shù)據(jù)集是多種多樣的，相對(duì)于其他數(shù)據(jù)集而言，尚缺乏特定病癥的大型數(shù)據(jù)集，因此可用于深度學(xué)習(xí)的數(shù)據(jù)量較少。且圖像中存在著一定程度的噪聲干擾，不確定性較高，如何抵抗噪聲干擾并從較少的圖像數(shù)據(jù)中提取有效的隱藏信息成為提高青光眼識(shí)別準(zhǔn)確率的關(guān)鍵。

Maheshwari 等提出了一種基于位平面分解的青光眼診斷方法，對(duì)圖像進(jìn)行RGB 三通道分離并進(jìn)行位平面分解，再利用局部二值模式對(duì)每個(gè)位平面提取紋理特征，再將這些特征分別輸入到三個(gè)不同的支持向量機(jī)進(jìn)行分類(lèi)，實(shí)驗(yàn)表明在區(qū)分正常和青光眼病例方面是有效的。為了評(píng)估PET/CT 圖像區(qū)分非小細(xì)胞肺癌（NSCLC）的能力，Han 等比較了10 種機(jī)器學(xué)習(xí)模型和VGG16 深度學(xué)習(xí)方法，發(fā)現(xiàn)VGG16 算法優(yōu)于這些傳統(tǒng)機(jī)器學(xué)習(xí)方法。由于VGG16 網(wǎng)絡(luò)有較高性能，Montaha 等提出了一種基于VGG16 的網(wǎng)絡(luò)消融研究方法。以預(yù)處理后的乳腺X 光圖像作為數(shù)據(jù)集，評(píng)估該網(wǎng)絡(luò)與另外五種網(wǎng)絡(luò)的健壯性，研究表明，該網(wǎng)絡(luò)模型表現(xiàn)最佳。

本文主要提出了結(jié)合位平面分解和提取特征子圖兩種方法提取眼底圖特征，通過(guò)CNN 分類(lèi)器實(shí)現(xiàn)青光眼與正常眼的分類(lèi)。為了驗(yàn)證提取圖像特征的健壯性和魯棒性，采用FGSM對(duì)原始眼底圖像進(jìn)行對(duì)抗攻擊，再對(duì)攻擊后的圖像提取特征子圖和位平面分解，將位平面圖像分別輸入CNN 分類(lèi)器中以抵御對(duì)抗攻擊，并比較攻擊后的原始圖像與對(duì)應(yīng)特征子圖位平面圖像之間的分類(lèi)效果。

1 卷積神經(jīng)網(wǎng)絡(luò)分類(lèi)器

卷積神經(jīng)網(wǎng)絡(luò)（CNN）是一種包含卷積計(jì)算且具有深度結(jié)構(gòu)的前饋神經(jīng)網(wǎng)絡(luò)，是深度學(xué)習(xí)的代表算法之一。圖1是卷積神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖，該網(wǎng)絡(luò)簡(jiǎn)化了對(duì)圖像的復(fù)雜前期預(yù)處理，可以直接將圖片作為網(wǎng)絡(luò)的輸入，自動(dòng)對(duì)圖像提取特征，并且對(duì)圖片的變形具有高度不變性。實(shí)驗(yàn)使用的卷積神經(jīng)網(wǎng)絡(luò)分類(lèi)器共有三層卷積層，卷積核大小均為3x3，激活函數(shù)使用relu 函數(shù)，池化層采用最大池化策略，大小為2x2，全連接層激活函數(shù)使用softmax函數(shù)。

圖1 CNN網(wǎng)絡(luò)結(jié)構(gòu)圖

實(shí)驗(yàn)采用公共醫(yī)學(xué)圖像眼底影像作為數(shù)據(jù)集，總體分為兩類(lèi)，正常眼眼底與青光眼病變眼底。其中，正常樣本數(shù)為483，青光眼樣本數(shù)為167，把圖像分為訓(xùn)練集和測(cè)試集，比例為1 ∶1。將數(shù)據(jù)集輸入分類(lèi)器中，經(jīng)過(guò)80次迭代得到識(shí)別準(zhǔn)確率為77.6%，損失率為47.5%，如圖2、圖3 所示。實(shí)驗(yàn)結(jié)果顯示，原始圖像的分類(lèi)準(zhǔn)確率較低，準(zhǔn)確率曲線(xiàn)呈上升趨勢(shì)但未達(dá)到收斂狀態(tài)，并不適宜投入臨床試驗(yàn)，分類(lèi)性能還有一定的提升空間。CNN 的分類(lèi)性能受到迭代次數(shù)、網(wǎng)絡(luò)層數(shù)以及數(shù)據(jù)集的影響。本文以同一CNN 分類(lèi)器作對(duì)比實(shí)驗(yàn)，為了提高圖像識(shí)別率，采用位平面分解和生成特征子圖兩種手段對(duì)原始圖像進(jìn)行特征提取。

圖2 原始眼底影像分類(lèi)準(zhǔn)確率

圖3 原始眼底影像分類(lèi)損失率

2 快速梯度符號(hào)方法

Goodfellow 等提出了快速梯度符號(hào)（Fast Gradient Sign Method，F(xiàn)GSM）圖像攻擊算法。該文指出，一張能被常規(guī)分類(lèi)器分類(lèi)為熊貓的圖片，通過(guò)添加由網(wǎng)絡(luò)梯度生成的攻擊噪聲后，看起來(lái)仍然是熊貓，但分類(lèi)器卻將其分類(lèi)為長(zhǎng)臂猿。為了模擬醫(yī)學(xué)影像數(shù)據(jù)集遭遇噪聲攻擊的情況，分別利用高斯噪聲和FGSM對(duì)眼底圖像進(jìn)行攻擊，攻擊效果如圖4所示。

圖4 對(duì)原始圖像進(jìn)行攻擊

由圖4看出，經(jīng)過(guò)高斯攻擊后的圖像相對(duì)于原始圖像有著明顯差異，清晰度有所下降，部分圖像特征丟失。而經(jīng)過(guò)FGSM攻擊后的眼底圖像與原始圖像相似，仍能夠從圖像中觀測(cè)出原圖的整體特征。將兩種攻擊后的圖像輸入分類(lèi)器中，分類(lèi)準(zhǔn)確率和損失率如圖5、圖6、圖7和圖8所示。

圖5 高斯攻擊圖像分類(lèi)準(zhǔn)確率

圖6 高斯攻擊圖像分類(lèi)損失率

圖7 FGSM攻擊圖像分類(lèi)準(zhǔn)確率

圖8 FGSM攻擊圖像分類(lèi)損失率

經(jīng)過(guò)高斯噪聲攻擊的圖像分類(lèi)準(zhǔn)確率為66.7%，與原始圖像相比下降了11.0%；經(jīng)過(guò)FGSM 攻擊后的圖像分類(lèi)準(zhǔn)確率為72.9%，下降了4.7%。結(jié)果表明，高斯噪聲攻擊會(huì)降低圖像質(zhì)量，干擾噪聲會(huì)影響圖像分類(lèi)準(zhǔn)確率。FGSM相較于高斯噪聲攻擊而言，其隱蔽性更強(qiáng)，能在肉眼無(wú)法分辨的情況下對(duì)圖像進(jìn)行攻擊，更容易被惡意利用，在該種情況下，分類(lèi)器的分類(lèi)性能與人工觀測(cè)判斷相比有著明顯劣勢(shì)。如何在受到噪聲干擾的圖像數(shù)據(jù)集中提取出有效的圖像特征，成為抵御攻擊、提高圖像識(shí)別準(zhǔn)確率的關(guān)鍵。

3 特征子圖

VGG16網(wǎng)絡(luò)是2014 年牛津大學(xué)和Google DeepMind 公司聯(lián)合研發(fā)的深度網(wǎng)絡(luò)模型，該網(wǎng)絡(luò)一共有16 個(gè)訓(xùn)練參數(shù)，由13 層卷積層和3 層全連接層組成。VGG 網(wǎng)絡(luò)能將輸入的圖像提取出25 張?zhí)卣髯訄D，以及25 張?zhí)卣髯訄D融合后的一張整體特征圖，每張圖像共生成26 張圖像，形成了一個(gè)新的數(shù)據(jù)集。VGG16 網(wǎng)絡(luò)及其生成的特征子圖分別如圖9、圖10所示。

圖9 VGG16網(wǎng)絡(luò)

從圖10 可以看出，26 種圖像反映的圖像特征貢獻(xiàn)不一樣，某些子圖能夠很好反映出原始圖像的結(jié)構(gòu)特征，但某些子圖只能看到眼底輪廓，并且包含了大部分的黑暗區(qū)域，與原圖差異較大。我們通常認(rèn)為這些圖像是不適宜進(jìn)行分類(lèi)的，應(yīng)該抽取那些特征明顯或者更能反映杯盤(pán)比的圖像特征作為輸入圖像。為了驗(yàn)證假設(shè)，分別將所有特征圖像輸入分類(lèi)器，比較它們的識(shí)別率，準(zhǔn)確率與損失率如表1所示。

圖10 VGG網(wǎng)絡(luò)提取的特征子圖

表1 青光眼特征子圖的分類(lèi)準(zhǔn)確率和損失率

結(jié)果顯示，提取圖像特征后，圖像分類(lèi)準(zhǔn)確率提升顯著，平均識(shí)別率為79.23%，最高可達(dá)89.94%。這表明運(yùn)用VGG16 網(wǎng)絡(luò)對(duì)圖像提取特征子圖是有效的，能夠總體提高圖像的準(zhǔn)確率，但與此同時(shí)，特征子圖還有待改進(jìn)的空間，那些我們認(rèn)為應(yīng)該從數(shù)據(jù)集中剔除的特征子圖中有可能隱藏了一些肉眼不可見(jiàn)的特征，也可能是因?yàn)榉诸?lèi)器分類(lèi)失敗所導(dǎo)致的，與網(wǎng)絡(luò)梯度攻擊后的圖像分類(lèi)結(jié)果類(lèi)似。為了消除這種不確定性，我們將有可能影響分類(lèi)性能的特征子圖摒棄。

4 位平面分解

26 個(gè)特征圖像的最高識(shí)別率為89.9%，為了進(jìn)一步提高識(shí)別率，突破特征子圖所達(dá)到的準(zhǔn)確率瓶頸，我們對(duì)篩選出來(lái)的特征子圖進(jìn)行位平面分解，進(jìn)一步蒸餾圖像特征。用圖像作為底面，表示像素亮度大小的8位二進(jìn)制數(shù)作為高度，可形成一個(gè)立體直方圖，各像素位置相同的位形成的一個(gè)平面成為圖像的位平面，這個(gè)過(guò)程稱(chēng)為位平面分解。先將提取出的特征子圖進(jìn)行位平面分解，然后采用格雷碼進(jìn)行位平面置換，通常可以減小灰度變化所帶來(lái)的影響。圖11 顯示了青光眼圖像特征子圖的位平面分解結(jié)果。

圖11 位平面分解后進(jìn)行格雷碼置換

圖11（0）為被分解前的原圖，從圖11（1）到圖11（8）為最低位平面到最高位平面。通常情況下，位平面越高，圖像整體特征越明顯，更接近原圖，說(shuō)明高位平面圖保留了較多的原圖信息，在原圖組成上更為關(guān)鍵。低位平面圖像能反映局部特征，但整體與原圖有著一定差異。從圖11 看出，由于圖11（7）、圖11（8）位平面強(qiáng)調(diào)原圖信息，反而丟失了圖像的特征，導(dǎo)致圖像中出現(xiàn)了較多的無(wú)效區(qū)域。圖11（1）、圖11（2）位平面作為低位平面，更容易受到噪聲干擾，有可能承載了更多的噪聲信息，因此通過(guò)實(shí)驗(yàn)分析哪些位平面層更適合分類(lèi)顯得尤其關(guān)鍵。將篩選出的特征圖位平面圖像輸入CNN 分類(lèi)器中，得出最高的分類(lèi)準(zhǔn)確率達(dá)到97.41%，準(zhǔn)確率與損失率分別如圖12、圖13所示。

圖12 特征子圖最高分類(lèi)準(zhǔn)確率

圖13 對(duì)應(yīng)特征子圖分類(lèi)損失率

經(jīng)過(guò)80 次迭代，準(zhǔn)確率曲線(xiàn)在97.0%~97.5%之間呈震蕩形態(tài)，震蕩區(qū)間為0.5%，而損失率曲線(xiàn)在47%～41%之間呈震蕩形態(tài)，震蕩區(qū)間為6%。由于準(zhǔn)確率曲線(xiàn)的震蕩區(qū)間較小，在宏觀上看，準(zhǔn)確率曲線(xiàn)收斂于97.3%，平均損失率為44%左右，分類(lèi)性能良好，準(zhǔn)確率比原始圖像提高了19.7%，比高斯噪聲攻擊后的圖像提高了30.6%，比FGSM 攻擊后的圖像提高了24.4%。綜上所述，遭受噪聲攻擊的圖像雖然會(huì)造成分類(lèi)性能的下降，但通過(guò)對(duì)被攻擊的圖像生成特征子圖，再結(jié)合格雷碼位平面置換的方法，仍能夠從圖像中提取出關(guān)鍵特征，提高分類(lèi)器分類(lèi)性能，達(dá)到較高的準(zhǔn)確率與較低的損失率。

5 結(jié)語(yǔ)

本文針對(duì)青光眼臨床診斷依賴(lài)人工經(jīng)驗(yàn)判斷容易引發(fā)錯(cuò)誤的問(wèn)題，提出了一種基于CNN的青光眼分類(lèi)器。結(jié)合特征子圖和位平面分解提取青光眼圖像特征，通過(guò)實(shí)驗(yàn)結(jié)果說(shuō)明該方法能夠抵御噪聲干擾并提升分類(lèi)準(zhǔn)確率，有助于青光眼的早期診斷與篩查。該分類(lèi)器簡(jiǎn)單，會(huì)對(duì)分類(lèi)準(zhǔn)確率有所制約，將來(lái)的研究工作主要包括優(yōu)化分類(lèi)器網(wǎng)絡(luò)結(jié)構(gòu)，用不同網(wǎng)絡(luò)提高分類(lèi)器的分類(lèi)效能，降低損失率。