基于計(jì)算機(jī)視覺(jué)的糖尿病視網(wǎng)膜病變自動(dòng)篩查系統(tǒng)

朱江兵 柯 鑫 劉 暢 楊金奎 何建國(guó)*

(1.北京大恒圖像視覺(jué)有限公司，北京100085;2.首都醫(yī)科大學(xué)附屬北京同仁醫(yī)院內(nèi)分泌科，北京100730)

糖尿病性視網(wǎng)膜病變(diabetic retinopathy，DR，以下簡(jiǎn)稱(chēng)糖網(wǎng))，是美國(guó)和歐洲工作人口中最主要的致盲眼病。據(jù)世界衛(wèi)生組織(World Health Organization，WHO)預(yù)測(cè)，到2030年，全球DR患者人數(shù)將增加到3.66億，糖網(wǎng)病防治將成為一個(gè)更為嚴(yán)重的世界性問(wèn)題［1］。研究［2-4］表明，對(duì) DR 患者進(jìn)行早期診斷和治療能有效防止視覺(jué)的損失以及失明，而防治的關(guān)鍵則是通過(guò)眼底照相檢查，定期隨訪(fǎng)發(fā)現(xiàn)病情的進(jìn)展，及時(shí)進(jìn)行激光干預(yù)治療。但是目前世界上超過(guò)50%的糖尿病患者沒(méi)有接受任何形式的眼部檢查，基于眼底影像的糖尿病視網(wǎng)膜病變的篩查工作基本還是依靠眼科醫(yī)生的肉眼觀察進(jìn)行。但當(dāng)面臨大規(guī)模篩查時(shí)，需要醫(yī)生分析和處理的數(shù)據(jù)量非常大，人工判讀方法既費(fèi)時(shí)又費(fèi)力，而且主觀性強(qiáng)，數(shù)據(jù)分析復(fù)雜并且難以量化，很難做到定量隨訪(fǎng)。所以上述問(wèn)題已成為實(shí)施大規(guī)模糖尿病視網(wǎng)膜病變篩查的主要障礙，臨床上迫切需要一種客觀、準(zhǔn)確的方法，來(lái)輔助眼科醫(yī)生對(duì)糖尿病患者眼底照相的結(jié)果進(jìn)行快速分析［5-6］。

近年來(lái)，隨著計(jì)算機(jī)輔助診斷技術(shù)的發(fā)展，基于計(jì)算機(jī)視覺(jué)的相關(guān)技術(shù)已經(jīng)在肝臟疾病、呼吸系統(tǒng)疾病的影像診斷中得到開(kāi)發(fā)和應(yīng)用［7-9］。基于計(jì)算機(jī)視覺(jué)的糖網(wǎng)自動(dòng)篩查系統(tǒng)正是基于大規(guī)模人口篩查的需求，通過(guò)計(jì)算機(jī)視覺(jué)的相關(guān)技術(shù)來(lái)輔助醫(yī)生判別和精確提取影像的病變指標(biāo)，智能分析患者的疾病信息，從而實(shí)現(xiàn)糖網(wǎng)病的大規(guī)模自動(dòng)篩查。由于計(jì)算機(jī)可以全面利用影像信息進(jìn)行精確的定量計(jì)算，不僅可以大幅提高醫(yī)生的診斷效率，讓醫(yī)生從繁重的閱片工作中解脫出來(lái)，而且能去除人的主觀性，避免因個(gè)人知識(shí)和經(jīng)驗(yàn)差異引起診斷誤差，對(duì)于建立客觀準(zhǔn)確的臨床診斷標(biāo)準(zhǔn)很有價(jià)值。本研究結(jié)合國(guó)際Messidor數(shù)據(jù)庫(kù)［6］，基于計(jì)算機(jī)視覺(jué)相關(guān)技術(shù)設(shè)計(jì)糖網(wǎng)病自動(dòng)檢測(cè)算法提取包括出血、滲出、微血管瘤等各類(lèi)糖網(wǎng)病變特征，再根據(jù)DR的臨床診斷標(biāo)準(zhǔn)，對(duì)眼底影像進(jìn)行自動(dòng)分級(jí)診斷，實(shí)現(xiàn)了全新的糖網(wǎng)病自動(dòng)篩查系統(tǒng)，并與經(jīng)過(guò)專(zhuān)家驗(yàn)證的篩查結(jié)果進(jìn)行了對(duì)比，研究表明，基于計(jì)算機(jī)視覺(jué)算法開(kāi)發(fā)的糖網(wǎng)自動(dòng)篩查系統(tǒng)能準(zhǔn)確、高效的完成眼科影像的糖網(wǎng)篩查工作，節(jié)省了醫(yī)療成本，便于未來(lái)建立客觀準(zhǔn)確的大規(guī)模健康數(shù)據(jù)管理系統(tǒng)，展現(xiàn)出很好的臨床應(yīng)用前景。

1 材料與方法

1.1 數(shù)據(jù)

國(guó)際Messidor數(shù)據(jù)庫(kù)［6］建立的目的即是便于對(duì)DR作計(jì)算機(jī)輔助診斷的研究。Messidor數(shù)據(jù)庫(kù)是當(dāng)前已公開(kāi)的包含眼底圖像數(shù)目最多的數(shù)據(jù)庫(kù)，包括2001幅彩色眼底圖片，其中包括800張散瞳的圖片和400張未散瞳的圖片。這些圖片都是在常規(guī)的臨床眼底檢查中使用眼底照相機(jī)拍攝的。眼底圖像數(shù)據(jù)庫(kù)的所有圖片分為3組，每組400張，每組再分成4個(gè)圖片集，共12個(gè)圖片集，每個(gè)數(shù)據(jù)集為100張圖，眼底圖片的分辨率分別為:1 400×960、2 240×1488、2 304×1 536，所有圖像保存為未壓縮的TIFF格式。對(duì)于每張圖像，Messidor數(shù)據(jù)庫(kù)提供了2個(gè)診斷結(jié)果，包括視網(wǎng)膜病變的等級(jí)和黃斑水腫的風(fēng)險(xiǎn)等級(jí)。這些診斷結(jié)果由醫(yī)學(xué)專(zhuān)家根據(jù)表1中所示的分級(jí)方案得到，也被作為對(duì)本文篩查系統(tǒng)性能分析的參考標(biāo)準(zhǔn)。根據(jù)這個(gè)參考標(biāo)準(zhǔn)，共有546張圖像被判定為正常，654張圖像被判定為DR，具體來(lái)說(shuō)，糖網(wǎng)1級(jí)共有153張，糖網(wǎng)2級(jí)共有247張，糖網(wǎng)3級(jí)有254張。黃斑水腫不在本文討論范圍。此外，有974張圖像沒(méi)有黃斑水腫的風(fēng)險(xiǎn)等級(jí);另分別有75張圖像風(fēng)險(xiǎn)等級(jí)為1和151張圖像風(fēng)險(xiǎn)等級(jí)為2。患者信息被抹去以確保患者的隱私。

表1 公開(kāi)的糖網(wǎng)病以及黃斑水腫的分級(jí)標(biāo)準(zhǔn)Tab.1 Grading schemes proposed for retino pathy grade and risk of macular edema

1.2 方法

本文提出的糖網(wǎng)自動(dòng)篩查系統(tǒng)通過(guò)分析患者的眼底檢查圖像中與DR相關(guān)的病變來(lái)提供診斷參考。如果檢查圖像中包含DR的病變或圖像質(zhì)量較低，那么被認(rèn)為很可能患有糖網(wǎng)病。為了做到這一點(diǎn)，本系統(tǒng)包含以下幾個(gè)模塊:

1.2.1 圖像預(yù)處理

在檢測(cè)眼底圖像中的解剖結(jié)構(gòu)和病變之前，眼底圖像的視場(chǎng)(region of interest，ROI)需要通過(guò)自適應(yīng)閾值分割和模板匹配的方法預(yù)先找到，圖像的黑色背景得以去除。

1.2.2 圖像質(zhì)量校正

為了從千差萬(wàn)別的眼底圖像中獲得穩(wěn)定的信息，對(duì)圖像進(jìn)行顏色、亮度以及曝光分布的歸一化處理十分必要。本文先根據(jù)圖像的曝光分布對(duì)圖像的亮度分布進(jìn)行反補(bǔ)償，然后選定一副質(zhì)量較好的圖像作為參考圖像，把被處理圖像的色調(diào)、飽和度和亮度信息都?xì)w一化到參考圖像的參考值上，從而保證了被處理圖像的穩(wěn)定性。

1.2.3 分割血管

血管系統(tǒng)是視網(wǎng)膜圖像中最重要的解剖結(jié)構(gòu)之一。血管的準(zhǔn)確分割對(duì)于從紅色病變中區(qū)分毛細(xì)血管以及其他解剖結(jié)構(gòu)，比如視盤(pán)的檢測(cè)都非常關(guān)鍵［10］。通過(guò)亮度閾值分割以及基于顏色特征的監(jiān)督分類(lèi)方法能有效分割出血管邊界。

1.2.4 分割視盤(pán)

視盤(pán)是眼底另一個(gè)很重要的解剖結(jié)構(gòu)。視盤(pán)的正確定位對(duì)于后期減少亮度病變的錯(cuò)判非常有用。視盤(pán)的檢測(cè)通過(guò)其亮度、形狀以及血管的方向等特征進(jìn)行定位。

1.2.5 紅色病變檢測(cè)

紅色病變，包括微血管瘤、出血和增生血管，是DR的重要特征［11］，因此，這些病變的正確檢測(cè)對(duì)DR篩查系統(tǒng)而言非常重要。小的紅色病灶的定位基于數(shù)學(xué)形態(tài)學(xué)來(lái)獲得候選區(qū)域，大的紅色病灶的定位則通過(guò)基于顏色特征的監(jiān)督分類(lèi)來(lái)完成［12］。使用有監(jiān)督的分類(lèi)器和一組描述的候選區(qū)域形狀、結(jié)構(gòu)、顏色和對(duì)比度的特征，可計(jì)算出檢測(cè)出的候選區(qū)域是一個(gè)真正紅色病變的概率。

1.2.6 亮度病變檢測(cè)

亮度病變，如滲出、棉絮斑、玻璃膜疣，經(jīng)常出現(xiàn)在DR人群篩查中。只有前兩種病變與DR相關(guān)。與紅色病變的檢出相類(lèi)似，像素級(jí)的有監(jiān)督分類(lèi)用于獲得候選的亮度病變區(qū)域［13］。每個(gè)候選區(qū)域是真正亮度病變的概率通過(guò)使用一組候選區(qū)域特征的監(jiān)督分類(lèi)平均值獲得，這些特征包括形狀、對(duì)比度、色彩和與最近紅色病變的距離［14］。

以上這些不同模塊的輸出必須被組合在一起以獲得關(guān)于患者檢查的一個(gè)最終結(jié)果。為了實(shí)現(xiàn)這個(gè)目標(biāo)，就需要計(jì)算上述模塊的各個(gè)輸出的特征，用監(jiān)督分類(lèi)器來(lái)判定病變的可能性［15］，從而為進(jìn)一步的診斷標(biāo)準(zhǔn)提供依據(jù)，實(shí)現(xiàn)對(duì)DR的自動(dòng)篩查。

參考醫(yī)學(xué)統(tǒng)計(jì)學(xué)的方法，本文選用靈敏度(sensitivity)及特異度(specificity)這兩個(gè)常用指標(biāo)來(lái)說(shuō)明篩查系統(tǒng)的檢測(cè)性能。其中靈敏度又稱(chēng)真陽(yáng)性率，即患者被診斷為陽(yáng)性的概率，計(jì)算公式是 :真陽(yáng)性/(真陽(yáng)性+假陰性)×100% ，此值越大，說(shuō)明診斷試驗(yàn)越靈敏，漏診率越低。特異度又稱(chēng)真陰性率，即實(shí)際上未患病的人被診斷為陰性的概率，計(jì)算公式是:真陰性/(真陰性+假陽(yáng)性)×100% ，此值越大，說(shuō)明診斷試驗(yàn)越精確，誤診的比例越低。

1.3 軟件及工作平臺(tái)

本文圖形圖像算法部分的代碼用Matlab實(shí)現(xiàn)，在Matlab2008a上編寫(xiě)，軟件系統(tǒng)部分的代碼用C++實(shí)現(xiàn)，在 VS2010上編寫(xiě)。程序在 CPU2.4GHZ、內(nèi)存4GB的Windows7的平臺(tái)下運(yùn)行。糖網(wǎng)自動(dòng)篩查系統(tǒng)流程示例如圖1。

圖1 糖網(wǎng)自動(dòng)篩查系統(tǒng)的流程示例Fig.1 Examples of the outputs of automated DR screening system

2 結(jié)果

為了驗(yàn)證本文提出的糖網(wǎng)自動(dòng)篩查系統(tǒng)的性能，一位具有4至20年DR篩選經(jīng)驗(yàn)的視網(wǎng)膜專(zhuān)家分別手動(dòng)分析了Messidor數(shù)據(jù)庫(kù)中的圖像，專(zhuān)家根據(jù)表1中提供的方案評(píng)定視網(wǎng)膜病變等級(jí)和黃斑水腫的風(fēng)險(xiǎn)等級(jí)，然后按照等級(jí)歸類(lèi)存儲(chǔ)。本文把人工檢測(cè)，篩查系統(tǒng)自動(dòng)檢測(cè)以及Messidor數(shù)據(jù)給定的結(jié)果做了對(duì)比，詳見(jiàn)表2。

表2 專(zhuān)家檢測(cè)、篩查系統(tǒng)自動(dòng)檢測(cè)以及Messidor數(shù)據(jù)參考標(biāo)準(zhǔn)的對(duì)比Tab.2 Outputs comparison of experts，CAD system and messidor data reference standard

本文進(jìn)一步使用ROCKIT軟件來(lái)分析人工判讀和機(jī)器判讀的檢測(cè)結(jié)果。筆者把原始統(tǒng)計(jì)結(jié)果作為輸入數(shù)據(jù)，ROCKIT軟件采用最大似然估計(jì)法擬合出一條二次接收者操作特征(receiver operator characteristics，ROC)曲線(xiàn)，ROC曲線(xiàn)下的面積Az反映人工判讀或者機(jī)器判讀的檢測(cè)性能。其中眼科專(zhuān)家和機(jī)器判定圖片正常異常與否這一實(shí)驗(yàn)中，采用了Messidor數(shù)據(jù)的全部1 200張圖，其對(duì)比結(jié)果如圖2所示。

分析表2中的數(shù)據(jù)可知，總共546張正常的眼底圖像中，自動(dòng)篩查系統(tǒng)患者正常的結(jié)果圖為512張，其中錯(cuò)把34張正常圖誤認(rèn)為異常，對(duì)比而言，人工錯(cuò)判的張數(shù)稍多，為44張，在總共654張非正常的眼底圖像中，機(jī)器系統(tǒng)的正確檢出結(jié)果為618張，漏檢36張異常圖片為正常，人工漏檢為32張。按照靈敏度和特異度的計(jì)算公式［11］，可知機(jī)器系統(tǒng)靈敏度即正確判斷患者的比率為93.8%，系統(tǒng)特異度即正確判斷非患者的比率為94.5%。相比而言，人工判讀的靈敏度為91.6%，特異度為94.9%。

通過(guò)ROC曲線(xiàn)的比較可知，計(jì)算機(jī)自動(dòng)判讀圖像正常異常的精度和人工判讀的精度是相當(dāng)?shù)摹＿@說(shuō)明本文開(kāi)發(fā)的糖網(wǎng)自動(dòng)篩查系統(tǒng)能準(zhǔn)確的完成眼科影像的糖網(wǎng)篩查工作。

圖2 正常異常圖片判斷的ROC曲線(xiàn)Fig.2 Normal and abnormal ROC curves

本文還進(jìn)一步測(cè)算了人工判定糖網(wǎng)和計(jì)算機(jī)自動(dòng)判定糖網(wǎng)的運(yùn)行時(shí)間，結(jié)果如表3所示。

表3 人工檢測(cè)與自動(dòng)運(yùn)行時(shí)長(zhǎng)對(duì)比Tab.3 Comparison of runs time of DR screening

由表3可知，糖網(wǎng)自動(dòng)篩查系統(tǒng)檢測(cè)一幅圖片的時(shí)間不到10s，而醫(yī)生如果要做到相對(duì)準(zhǔn)確的判斷，需要花費(fèi)更長(zhǎng)的時(shí)間來(lái)分析，尤其對(duì)于微血管瘤的判斷，而且，人工查找圖片并歸類(lèi)存儲(chǔ)，需要花費(fèi)相當(dāng)一部分時(shí)間，而對(duì)于計(jì)算機(jī)而言，完全可以自動(dòng)并且快速地完成這些流程，大大提高了閱片的效率。相對(duì)于人工閱片，糖網(wǎng)自動(dòng)篩查系統(tǒng)在時(shí)間上有明顯優(yōu)勢(shì)。計(jì)算機(jī)另一個(gè)巨大的優(yōu)勢(shì)是可以24 h不間斷工作，這是閱片醫(yī)生無(wú)法做到的。面對(duì)當(dāng)前大規(guī)模篩查的需求，糖網(wǎng)自動(dòng)篩查系統(tǒng)展現(xiàn)了非常好的臨床應(yīng)用前景。

3 討論

大規(guī)模糖網(wǎng)篩查是糖網(wǎng)病防治的關(guān)鍵環(huán)節(jié)，而傳統(tǒng)依靠專(zhuān)業(yè)醫(yī)生人工判讀的方式已經(jīng)很難滿(mǎn)足當(dāng)前大規(guī)模篩查的需求。本研究旨在利用計(jì)算機(jī)視覺(jué)相關(guān)技術(shù)，開(kāi)發(fā)出能夠用于輔助DR篩查的計(jì)算機(jī)自動(dòng)篩查系統(tǒng)。利用數(shù)學(xué)形態(tài)學(xué)、區(qū)域分割、支撐向量機(jī)(support vector machine，SVM)分類(lèi)等計(jì)算機(jī)視覺(jué)技術(shù)設(shè)計(jì)出檢測(cè)DR包括出血、滲出、微血管瘤等各類(lèi)病變特征的算法，再根據(jù)DR的臨床診斷標(biāo)準(zhǔn)，完成對(duì)眼底影像進(jìn)行自動(dòng)分級(jí)診斷，實(shí)現(xiàn)自動(dòng)篩查。實(shí)驗(yàn)表明，糖網(wǎng)自動(dòng)篩查系統(tǒng)判定靈敏度為93.8%，特異度為94.5%，與人工判定精度相當(dāng)，但檢測(cè)時(shí)間卻大大減少，僅為人工閱片時(shí)間的1/30。本文研究表明，基于計(jì)算機(jī)視覺(jué)算法開(kāi)發(fā)的糖網(wǎng)自動(dòng)篩查系統(tǒng)能準(zhǔn)確、高效的完成眼科影像的糖網(wǎng)篩查工作，大幅減少閱片醫(yī)生的工作量和人為的主觀性，節(jié)省了醫(yī)療成本，便于未來(lái)建立客觀準(zhǔn)確的大規(guī)模健康數(shù)據(jù)管理系統(tǒng)，具有很好的臨床應(yīng)用前景和社會(huì)效益。

［1］ Diabetes:World Health Organization report［EB/OL］.(2015-01-01)［2015-10-21］.http://www.who.int/mediacentre/facts heets/fs312/en.

［2］ Kinyoun J，Barton F，F(xiàn)isher M，et al.Detection of diabetic macular edema:ophthalmoscopy versus photography-Early Treatment Diabetic Retino-pathy Study Report Number 5.The ETDRS Research Group［J］.Ophthalmology，1989，96(6):746-750.

［3］ ［NO authors listed］.Early photocoagulation for diabetic retinopathy.ETDRS report number 9.Early Treatment Diabetic Retinopathy Study Research Group［J］.Ophthalmology，1991，98(5 suppl):766-785.

［4］ Bresnick G H，Mukamel D B，Dickinson J C，et al.A screening approach to the surveillance of patients with diabetes for the presence of vision-threatening retinopathy［J］.Ophthalmology，2000，107(1):19-24.

［5］ Abra`moff M D，Niemeijer M，Suttorp-Schulten M S，et al.Evaluation of a system for automatic detection of diabetic retinopathy from color fundus photographs in a large population of patients with diabetes［J］.Diabetes Care，2008，31(2):193-198.

［6］ Methods to evaluate segmentation and indexing techniques in the field of retinal ophthalmology［EB/OL］.(2010-09-06)［2015-10-21］.http://messidor.crihan.fr/.

［7］ 蔡斌，彭謹(jǐn)，王剛，等.基于模式識(shí)別建立的動(dòng)脈鈣化影像計(jì)算機(jī)輔助診斷系統(tǒng)［J］.四川醫(yī)學(xué)，2012，33(6):913-914，915.

［8］ 袁克虹，向蘭茜.用于計(jì)算機(jī)輔助診斷的肺實(shí)質(zhì)自動(dòng)分割方法［J］.清華大學(xué)學(xué)報(bào)(自然科學(xué)版)，2011，51(1):90-95.

［9］ 孫旭輝，田啟川，李臨生，等.面向CT影像的肺結(jié)節(jié)計(jì)算機(jī)輔助診斷算法［J］.計(jì)算機(jī)測(cè)量與控制，2011，19(2):295-302.

［10］Niemeijer M，Staal J J，van Ginneken B，et al.Comparative study of retinal vessel segmentation methods on a new publicly available database，in medical imaging［J］.Proc SPIE，2004，5370:648-656.

［11］Niemeijer M，Abra`moff M D，van Ginneken B.Information fusion for diabetic retinopathy CAD in digital color fundus photographs［J］.IEEE Tran Med Imaging，2009，28(5):775-785.

［12］ Philip S，F(xiàn)leming A D，Goatman K A.The efficacy of automated“disease/no disease”grading for diabetic retinopathy in a syste-matic screening programme［J］.Br J Ophthalmol，2007，91(11):1512-1517.

［13］ Abramoff M D，Suttorp-Schulten M S.Web-Based Screening for Diabetic Retinopa-thy in a Primary Care Population:The Eye-Check Project［J］.Telemed J E Health，2005，11(6):668-674.

［14］Tomi K，Valentina K，Iiris S，et al.DiaretDB1 V2.1-Diabetic retinopathy database and evaluation protocol［EB/OL］.(2009-05-06).［2015-10-21］.http://www2.it.lut.fi/project/imageret/diaretdb1_V2_1.

［15］ Zhang X，Thibault G，Decencière E，et al.Exudate detection in color retinal images for mass screening of diabetic retinopathy［J］.Med Image Anal，2014，18(7):1026-1043.