基于卷積神經(jīng)網(wǎng)絡(luò)的瞼板腺形態(tài)人工智能分析系統(tǒng)的構(gòu)建

朱敏穎 林曉蕾 張祖輝 戴琦

瞼板腺功能障礙（meibomian gland dysfunction，MGD）是一種慢性、彌漫性瞼板腺異常，可引起淚膜異常、眼部刺激癥狀、炎癥反應(yīng)及眼表損害[1]。目前MGD的病因及發(fā)病機制尚不明確，但有研究發(fā)現(xiàn)瞼板腺形態(tài)和功能存在密切關(guān)聯(lián)。Pult等[2]發(fā)現(xiàn)瞼板腺形態(tài)缺失與淚膜脂質(zhì)層厚度和穩(wěn)定性密切相關(guān)。Yin等[3]研究認為瞼板腺形態(tài)上不均勻缺失與MGD密切相關(guān)。Ban等[4]認為上眼瞼瞼板腺缺失與瞼板腺分泌分數(shù)及淚膜脂質(zhì)層相關(guān)。目前Arita等[5]的分級法是使用較多的評價瞼板腺形態(tài)缺失的指標，但其人工標注工作量大、效率低、誤差大?；谏疃葘W(xué)習(xí)的人工智能（artificial intelligence，AI）系統(tǒng)在分類及圖像識別方面的效率和精度都遠高于圖像處理算法和早期的機器學(xué)習(xí)算法，且成本大幅降低，是目前醫(yī)學(xué)診斷當中最火熱和最具有前景的發(fā)展方向。深度學(xué)習(xí)經(jīng)過多年的發(fā)展，人們提出了許多深度神經(jīng)網(wǎng)絡(luò)，其中，卷積神經(jīng)網(wǎng)絡(luò)（convolutional neural network，CNN）已被成熟應(yīng)用于人臉識別、文字識別、圖像檢索等多個領(lǐng)域[6-7]。對于物體識別、圖片分類等，CNN的準確率可以接近甚至優(yōu)于人工分析方式[8-9]。Wang等[10]和Zhou等[11]使用CNN對瞼板腺萎縮進行分級，結(jié)果證明AI技術(shù)對瞼板腺萎縮的自動評價是有效的。但目前研究主要關(guān)注整個瞼板區(qū)域和萎縮區(qū)域的識別，并沒有對每根瞼板腺腺體進行檢測和評價。筆者團隊前期研究提出瞼板腺彎曲度及瞼板腺密度兩項定量評價指標[12]，可以成為AI系統(tǒng)的進階評價指標。本研究建立的AI系統(tǒng)是基于CNN的U-Net網(wǎng)絡(luò)構(gòu)建的，該系統(tǒng)能夠自動處理瞼板腺圖片，識別出瞼板腺的單根腺體，并使用前期提出的瞼板腺量化評價指標來準確評估瞼板腺的各項形態(tài)參數(shù)。

1 對象和方法

1.1 對象 選取2020年12月至2021年3月在溫州醫(yī)科大學(xué)附屬眼視光醫(yī)院杭州院區(qū)就診的60例受試者作為研究對象。其中阻塞性MGD患者32例（MGD組），健康志愿者28例（正常組）。MGD組男16例，女16例；年齡 22～50（33.28±9.28）歲；正常組男 11 例，女 17 例；年齡8～51（25.25±11.19）歲。兩組受試者性別比較差異無統(tǒng)計學(xué)意義（P＞0.05），但年齡比較差異有統(tǒng)計學(xué)意義（P＜0.01）。排除標準：（1）患有眼部疾病或已知影響眼前節(jié)解剖結(jié)構(gòu)情況的患者，如有眼部炎癥、眼部手術(shù)史、佩戴隱形眼鏡和（或）外傷史；（2）有影響瞼板腺功能的全身用藥史的患者；（3）患有任何其他已知會影響淚膜的眼部疾病或全身性疾病者。所有患者均由兩位經(jīng)驗豐富的眼科醫(yī)生根據(jù)眼部癥狀、瞼緣異常和瞼板腺評分診斷為阻塞性MGD。如果滿足以下3個標準中的任何一個，則患者被診斷為阻塞性MGD：（1）眼部癥狀評分≥3分；（2）瞼緣評分≥2分；和（或）（3）瞼板腺評分≥3分[13]。只測量受試者的右眼。本研究經(jīng)溫州醫(yī)科大學(xué)附屬眼視光醫(yī)院醫(yī)學(xué)倫理委員會審批通過，所有受試者均簽署知情同意書。

1.2 AI模型的訓(xùn)練和測試

1.2.1 圖像預(yù)處理 使用眼表分析儀（德國 OCULUS Optikgerate GmbH歐科路光學(xué)儀器有限責(zé)任公司，型號：Keratography 5M）拍攝60例受試者右眼上下眼瞼的瞼板腺圖像，并將其分為40幅訓(xùn)練圖像和20幅測試圖像。訓(xùn)練圖像中的瞼板腺由兩位高級職稱醫(yī)師手動標注。在訓(xùn)練前，研究小組對40幅訓(xùn)練圖像進行了如下的預(yù)處理。

（1）圖像轉(zhuǎn)換為灰度，如下所示：

Igray=Ir×0.299+Ig×0.587+Ib×0.114

其中，Igray是灰度值，Ir是RGB圖像的R通道的值，Ig是RGB圖像的G通道的值，Ib是RGB圖像的B通道的值。

（2）標準化。使用以下公式對圖像進行標準化：

其中，Xstd是像素標準值，Xgray是像素灰度值，μ是像素平均值，σ是像素標準偏差。

（3）歸一化。使用以下公式對圖像進行歸一化：

其中，Xnorm是像素歸一化值。

（4）直方圖用限制對比度自適應(yīng)直方圖均衡調(diào)整。瞼板腺圖像不均勻，有些圖像對比度低，有些曝光過度。因此，采用對比度受限的自適應(yīng)直方圖均衡化算法來提高全局對比度，防止在相對均勻的圖像區(qū)域過度放大噪聲。

（5）伽馬校正。每個像素的亮度由以下公式調(diào)整：

其中，X[范圍從（0，1）]是圖像中的原始像素值，Y是調(diào)整后的像素值。

原始的圖像（a）和預(yù)處理完成的圖像（b）見圖1。

圖1 原始圖像（a）和預(yù)處理圖像（b）

由于數(shù)據(jù)集只有40幅圖像，并且太小而無法訓(xùn)練一個CNN，因此筆者通過從每個訓(xùn)練圖像中隨機選擇子圖像塊來進行數(shù)據(jù)增強，然后，使用80%的子圖像用于培訓(xùn)，剩下的20%用于驗證。圖2顯示了從40萬個子圖像中獲得的40幅圖像。

圖2 子圖像示例

1.2.2 系統(tǒng)設(shè)計與訓(xùn)練 本研究所使用的CNN結(jié)構(gòu)是基于U-Net網(wǎng)絡(luò)構(gòu)建的。卷積層的濾波器數(shù)量為32、64、128、64和32。交叉熵用于損失函數(shù)，Adam用于學(xué)習(xí)速率優(yōu)化。每個卷積層以整流線性單元（ReLU）作為激活函數(shù)，在兩個連續(xù)的卷積層之間使用0.2的丟失率來防止模型過度擬合。CNN模型的網(wǎng)絡(luò)結(jié)構(gòu)見圖3。

圖3 卷積神經(jīng)網(wǎng)絡(luò)（CNN）模型的網(wǎng)絡(luò)結(jié)構(gòu)

在GTX 1070 8 G GPU上進行120代計算，并使用4個小批量補丁，對迷你U-Net的訓(xùn)練持續(xù)大約15 h后，驗證集獲得了最高0.895的交并比（intersection-ofunion，IoU），表明該模型已經(jīng)過訓(xùn)練完成，可以從圖像中提瞼板腺。

預(yù)測圖像的大小通常＞64×64像素，因此研究人員需要將預(yù)測圖像分成多個64×64像素的子畫面。最后，將所有子圖的結(jié)果合并，得到最終的預(yù)測結(jié)果。為了減少預(yù)測誤差，研究人員使用5個像素作為水平和垂直步長來預(yù)測每個子畫面。這意味著每個子畫面將產(chǎn)生25個預(yù)測結(jié)果，研究人員將這25個預(yù)測結(jié)果的平均值作為子畫面的預(yù)測結(jié)果。

研究人員通過傳統(tǒng)的圖像處理方法對圖像進行預(yù)測，得到預(yù)測結(jié)果。提取步驟如下：（1）對圖像進行中值濾波；（2）將圖像二值化；（3）在圖像中找到瞼板腺的輪廓。如果兩個等高線之間的垂直距離很小，則認為這些等高線屬于同一個瞼板腺并且是連接在一起的。提取的瞼板腺在圖4中以灰白色顯示。

圖4 提取的瞼板腺

1.3 評價指標

1.3.1 瞼板腺功能指標

1.3.1.1 眼表疾病指數(shù)（ocular surface disease index，OSDI）[14]所有受試者均接受OSDI問卷和眼部癥狀問卷調(diào)查，并回答是否存在14種MGD相關(guān)的眼部癥狀：眼疲勞、分泌物、異物感、干燥、不適感、黏稠感、疼痛、溢淚、瘙癢、發(fā)紅、沉重感、眩光、過度眨眼和瞼板腺囊腫或瞼腺炎。根據(jù)出現(xiàn)的這些癥狀的數(shù)量，對癥狀從0～14進行評分，分數(shù)越高，說明癥狀越明顯。

1.3.1.2 淚河高度（tear meniscus height，TMH） 使用眼表分析儀測量淚河高度并進行瞼板腺紅外拍照，TMH測量方法采用眨眼5 s后測量下眼瞼中央TMH，正常值在0.20 mm左右。

1.3.1.3 淚膜破裂時間（tear film break-up time，TBUT）和角膜熒光素染色（corneal fluorescein staining，CFS）評分[13]滴入熒光素后進行CFS評分和TBUT測定。TBUT測定3次，取其平均值，正常值為＞10 s。CFS評分標準如下：在5個區(qū)域（中央、顳、鼻、上、下）評估角膜熒光素染色（0個點=0，1～5個點=1，6～15個點=2，16～30個點=3，＞30個點=4）；如果發(fā)現(xiàn)融合染色，染料擴散到周圍的上皮或基質(zhì)，則根據(jù)融合區(qū)域的數(shù)量或角膜絲狀物的存在進行評分（1個融合區(qū)域，加1分；2個或2個以上融合區(qū)域，加2分），染色評分范圍為0～20分，分數(shù)越高，說明角膜上皮損傷越重。

1.3.1.4 瞼緣異常評分 根據(jù)黏膜交界處前后移位、血管充盈、瞼板腺孔堵塞、眼瞼邊緣不規(guī)則4項參數(shù)對眼瞼邊緣異常進行評分，每項評分0～4分，分數(shù)越高說明瞼緣異常程度越大。

1.3.1.5 瞼板腺表達能力評分[15]評分標準如下：0分，清亮的瞼脂容易排出；1分，輕微壓力下排出混濁的瞼脂；2分，超過中等壓力下才能排出混濁的瞼脂；3分，即使很大的壓力下也無瞼脂排出。每個眼瞼瞼板腺表達能力評分范圍為0～45分，分數(shù)越低說明瞼板腺分泌功能越接近正常。

1.3.1.6 瞼板腺評分[16]評分標準如下：0分，無萎縮；1分，萎縮面積＜1/3的瞼板總面積；2分，萎縮面積介于1/3～2/3的瞼板總面積；3分，萎縮面積＞2/3的瞼板總面積。將上下瞼的瞼板腺評分相加，得到每只眼睛的瞼板腺評分范圍為0～6，分數(shù)越高，說明瞼板腺萎縮越嚴重。

1.3.2 瞼板腺形態(tài)指標 提取瞼板腺腺體后，計算腺體的以下指標。

（1）瞼板腺密度。計算每根瞼板腺的面積（Smg）（以像素為單位），然后用每個瞼板上所有腺體面積總和與瞼板的面積（St）作對比，公式如下:

（2）平均眼瞼瞼板腺寬度。所有瞼板腺寬度總和除以腺體數(shù)量，包括平均上眼瞼瞼板腺寬度、平均下眼瞼瞼板腺寬度、平均全眼瞼瞼板腺寬度。

（3）平均眼瞼瞼板腺長度。所有瞼板腺長度總和除以腺體數(shù)量，包括平均上眼瞼瞼板腺長度、平均下眼瞼瞼板腺長度、平均全眼瞼瞼板腺長度。

（4）平均眼瞼瞼板腺彎曲度。所有瞼板腺彎曲度總和除以腺體數(shù)量，包括平均上眼瞼瞼板腺彎曲度、平均下眼瞼瞼板腺彎曲度、平均全眼瞼瞼板腺彎曲度。

2 結(jié)果

2.1 兩組受試者瞼板腺功能指標比較 與正常組比較，MGD組患者TBUT減少，瞼板腺表達能力評分降低，OSDI、瞼緣異常評分和瞼板腺評分均升高，差異均有統(tǒng)計學(xué)意義（均 P＜0.01），見表 1。

表1 兩組受試者瞼板腺功能指標比較

2.2 兩組受試者瞼板腺形態(tài)指標比較 與正常組比較，MGD患者平均上眼瞼瞼板腺長度、平均下眼瞼瞼板腺長度和平均全眼瞼瞼板腺長度均明顯縮短，上眼瞼瞼板腺密度、下眼瞼瞼板腺密度和全眼瞼瞼板腺密度均明顯減小，差異均有統(tǒng)計學(xué)意義（均P＜0.01）。兩組受試者平均眼瞼瞼板腺寬度和彎曲度比較差異均無統(tǒng)計學(xué)意義（均 P＞0.05），見表 2。

表2 兩組受試者瞼板腺形態(tài)指標比較

2.3 60例受試者瞼板腺形態(tài)指標與功能指標的相關(guān)性分析 平均上眼瞼瞼板腺長度與TBUT、瞼板腺表達能力評分均呈正相關(guān)（均P＜0.05），與OSDI、瞼板腺評分均呈負相關(guān)（均P＜0.05）；平均上眼瞼瞼板腺寬度與瞼緣異常評分、瞼板腺評分均呈正相關(guān)（均P＜0.05）；平均上眼瞼瞼板腺彎曲度與瞼板腺表達能力評分呈負相關(guān)（P＜0.05）；瞼板腺密度與TBUT、瞼板腺表達能力評分均呈正相關(guān)（均P＜0.05），與CFS評分、瞼板腺評分均呈負相關(guān)（均P＜0.05）。平均下眼瞼瞼板腺長度與TBUT、瞼板腺表達能力評分均呈正相關(guān)（均P＜0.05），與瞼緣異常評分、瞼板腺評分均呈負相關(guān)（均P＜0.05）；平均下眼瞼瞼板腺寬度與TBUT呈正相關(guān)（P＜0.05），與瞼板腺評分呈負相關(guān)（P＜0.05）；瞼板腺密度與TBUT、瞼板腺表達能力評分均呈正相關(guān)（P＜0.05），與CFS評分、瞼板腺評分均呈負相關(guān)（均P＜0.05）。平均全眼瞼瞼板腺長度與TBUT、瞼板腺表達能力評分均呈正相關(guān)（均P＜0.05），與 OSDI、瞼板腺評分均呈負相關(guān)（均 P＜0.05）；瞼板腺密度與TBUT、瞼板腺表達能力評分均呈正相關(guān)（均 P＜0.05），與 OSDI、CFS評分、瞼緣異常評分、瞼板腺評分均呈負相關(guān)（均P＜0.05），見表3。

表3 瞼板腺形態(tài)指標與功能指標的相關(guān)性分析（r值）

3 討論

本研究介紹了一個基于CNN的AI系統(tǒng)，該系統(tǒng)使用了計算機視覺領(lǐng)域的最新技術(shù)[17-18]，可以自動分析瞼板腺的形態(tài)，包括瞼板腺長度、寬度、密度和彎曲度。計算機會在不斷的訓(xùn)練中自動提取圖片的高維度特征并建立這些特征與結(jié)果間的內(nèi)在聯(lián)系。這就在一定程度上排除了人為主觀因素可能對結(jié)果產(chǎn)生的干擾。自動獲取瞼板腺腺體并計算其形態(tài)特征是一個巨大的挑戰(zhàn)，以往的研究對輸入圖像有嚴格的要求，如果眼瞼區(qū)域出現(xiàn)皺褶、反光等干擾因素時，會導(dǎo)致腺體的形狀從規(guī)則的條狀變形為片狀或分散狀，檢測結(jié)果并不理想[19-20]。為了解決這一問題，本研究的算法對圖像識別進行了優(yōu)化，并將這些算法與傳統(tǒng)的數(shù)字圖像處理算法（如對比度調(diào)整、低通濾波、腐蝕、膨脹等）相結(jié)合，提高了相同背景下瞼板腺的識別率。此外，由于使用了先進的圖像分割模型和U-Net系統(tǒng)，即使在本研究中受試者人數(shù)較少，檢測瞼板腺形態(tài)的準確性也非常高（IoU達到0.895，重復(fù)率達到100%）。目前筆者團隊僅使用GTX 1070 8 G GPU已經(jīng)可以在100 ms內(nèi)處理一張圖片，隨著算力的提高，使得數(shù)萬幅圖像在毫秒內(nèi)同時完成分析成為可能。計算機輔助的瞼板腺分析系統(tǒng)可以大大提高分析的準確性和效率，降低分析的成本，克服人工標記主觀上的誤差。瞼板腺形態(tài)的詳細分析是確定MGD嚴重程度的重要指標[21-22]。本研究比較了正常組和MGD組瞼板腺長度、寬度、彎曲度和密度的差異。結(jié)果顯示，MGD組患者平均瞼板腺長度較正常組明顯縮短，瞼板腺密度也較正常組明顯減小，但兩組平均瞼板腺寬度和彎曲度無明顯差異。MGD患者的瞼板腺分泌物排出阻力增加，瞼板腺腺體的壓力過大，眼表情況的進一步惡化，導(dǎo)致瞼板腺病變的加重，瞼板腺長度縮短，并且出現(xiàn)腺體的萎縮，導(dǎo)致瞼板腺密度開始降低。如果病情得不到控制，瞼板腺功能障礙進一步加重，瞼板腺密度會進一步降低，直至接近0%。與之前的研究[23-26]相似，本研究中AI系統(tǒng)顯示瞼板腺的腺體長度和密度與癥狀、淚膜穩(wěn)定性、瞼緣異常評分和瞼板腺表達能力均相關(guān)，說明瞼板腺的腺體長度和密度可以作為計算機輔助的MGD診斷的有效指標。與以往的研究[5，26-28]不同，本研究采用高性能、高精度的AI系統(tǒng)，可以提取每個腺體并分析其長度、寬度、彎曲度等，每根腺體提取積累后在計算腺體面積和整個瞼板面積的比例，精確度較傳統(tǒng)的腺體面積計算方法大大提高。結(jié)果表明，AI系統(tǒng)能夠準確地得到瞼板腺的缺失面積，減少了由于分級過渡區(qū)附近的誤差而導(dǎo)致的分級系統(tǒng)的不一致性和可變性。如果想要在人工分析中采用筆者團隊的這種瞼板腺分析方法是不現(xiàn)實的，因為這種分析方法需要大量的人工工作量，不可能作為常規(guī)的評估指標，但是AI系統(tǒng)的特點恰恰適合這種需要大量分析工作的任務(wù)。

本研究也有局限性。研究中使用的算法是一種通用算法。如果將該算法進一步發(fā)展成為專門針對瞼板腺圖像的一種特殊算法，則瞼板腺形態(tài)識別有望獲得更高的準確率。此外，本研究的樣本量也較小，未來的研究中，筆者團隊將招募更多的受試者來訓(xùn)練和測試AI系統(tǒng)。

綜上所述，基于CNN的AI系統(tǒng)只需幾十個訓(xùn)練實例就可以自動檢測和識別瞼板腺。它是一個準確、高效的瞼板腺形態(tài)學(xué)評價系統(tǒng)，在臨床上，該系統(tǒng)可以大大提高MGD診斷的靈敏度和準確度，減少觀察者主觀方面導(dǎo)致的誤差，降低醫(yī)務(wù)人員的工作量。在未來的研究中，將通過AI系統(tǒng)得到瞼板腺形態(tài)各項參數(shù)，并進行一種更加細致的量化評分，從而提出一種新的客觀的MGD形態(tài)學(xué)診斷標準，這具有非常大的市場前景和臨床意義。