• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于深度序列特征學(xué)習(xí)的臨床感染性角膜炎圖像分類

    2021-12-30 06:52:44許葉圣孔鳴謝文加段潤平方鉦清林宇蕭朱強湯斯亮吳飛姚玉峰
    工程 2021年7期
    關(guān)鍵詞:眼科醫(yī)生子塊角膜炎

    許葉圣 ,孔鳴 ,謝文加,段潤平 ,方鉦清 ,林宇蕭 ,朱強 ,湯斯亮 ,吳飛 ,*,姚玉峰 ,*

    a Department of Ophthalmology, Sir Run Run Shaw Hospital, School of Medicine, Zhejiang University, Hangzhou 310016, China

    b College of Computer Science and Technology, Zhejiang University, Hangzhou 31002, China

    1. 引言

    通常來說,角膜疾病的診斷是由醫(yī)生通過觀察,結(jié)合個人的經(jīng)驗和知識得出的。近年來,使用深度卷積神經(jīng)網(wǎng)絡(luò)(CNN)的深度學(xué)習(xí)算法在醫(yī)療圖像處理領(lǐng)域大放異彩。疾病分類和診斷算法已在各醫(yī)學(xué)成像技術(shù)領(lǐng)域得到了應(yīng)用與驗證,包括計算機斷層掃描(CT)、磁共振成像(MRI)、眼底照相、光學(xué)相干斷層掃描(OCT)和病理圖像等[1]。這得益于醫(yī)學(xué)成像技術(shù)可以非常自然地導(dǎo)出豐富的圖像數(shù)據(jù),商業(yè)化的醫(yī)學(xué)成像技術(shù)對這些醫(yī)學(xué)圖像設(shè)定了統(tǒng)一的標(biāo)準(zhǔn),可以在短時間內(nèi)從單個或多個醫(yī)學(xué)中心收集這些圖像。

    然而大多數(shù)臨床疾病的診斷并不需要商業(yè)化的醫(yī)學(xué)成像技術(shù),因此許多醫(yī)療機構(gòu)在臨床診療中不進行圖像記錄保存,導(dǎo)致這類疾病圖像數(shù)據(jù)的收集主要靠各醫(yī)療中心的歷史積累,但針對這些臨床疾病開發(fā)機器學(xué)習(xí)診斷系統(tǒng)也同等重要。針對皮膚病變良惡性分類的研究[2]是在非傳統(tǒng)醫(yī)療圖像領(lǐng)域的一次成功嘗試。角膜疾病分類也可以效仿這種方式。角膜病是世界上最常見的致盲原因之一[3,4]。保守估計,世界上約有450萬人因罹患角膜病導(dǎo)致角膜混濁,從而蒙受中到重度的視力損傷[4]。感染性角膜炎是最常見的角膜病[5]。正常的角膜擁有獨有的透明性,感染性角膜炎最顯著的特點是病原體在角膜內(nèi)生長,導(dǎo)致角膜局灶性團塊混濁和粗糙,不可避免地使每種病原微生物在組織內(nèi)生長并呈現(xiàn)獨特的特點[6]。感染性角膜炎的診斷主要取決于眼科醫(yī)生對角膜感染性病變視覺特征的鑒別。臨床上,眼科醫(yī)生通常依靠裂隙燈顯微鏡來觀察角膜是否異常。除了用作觀察工具外,裂隙燈顯微鏡還可以用于拍攝照片并同時記錄每位患者的角膜狀態(tài),從而有助于開發(fā)出具有優(yōu)質(zhì)標(biāo)注的人工智能(AI)數(shù)據(jù)集,為感染性角膜炎的識別和分析打下基礎(chǔ)。

    從1998年開始,我們積累了一個大型且擁有正確標(biāo)注的裂隙燈顯微圖像數(shù)據(jù)集,該數(shù)據(jù)集來自10 609名角膜疾病患者,總共115 408張圖像?;谶@個數(shù)據(jù)集我們設(shè)計出一種基于深度學(xué)習(xí)的方法,以端到端的方式進行感染性角膜炎的診斷。為了直觀地模擬眼科醫(yī)生診斷感染性角膜炎的方式,我們提出了一種特征學(xué)習(xí)機制,通過學(xué)習(xí)序列特征來識別不同的視覺模式,其中蘊含著豐富的視覺信息。對于一張臨床圖像,從感染病灶區(qū)域的中心到邊緣的樣本子塊被重新排列,組成一個序列有序集(SOS),并輸入神經(jīng)網(wǎng)絡(luò)進行特征學(xué)習(xí)。我們提出的序列特征學(xué)習(xí)機制可以利用感染性病變區(qū)域子塊之間的空間關(guān)系,并且可以分離出數(shù)據(jù)樣本中蘊含的值得探索的差異性因素。此外,該數(shù)據(jù)集提供了一個潛在的策略以實現(xiàn)更可靠、有效、準(zhǔn)確的診斷。使用該數(shù)據(jù)集對我們的模型進行評估,獲得了比400名眼科醫(yī)生更高的診斷準(zhǔn)確率。

    2. 相關(guān)工作

    2.1. 醫(yī)療數(shù)據(jù)挖掘

    多年來,電子病歷(EMR)積累了大量的醫(yī)學(xué)數(shù)據(jù),這使研究人員從中發(fā)現(xiàn)了許多隱藏知識。 數(shù)據(jù)挖掘方法被廣泛應(yīng)用于醫(yī)學(xué)數(shù)據(jù),以發(fā)現(xiàn)隱藏的知識,并利用提取的知識來輔助各種有害疾病的預(yù)測、診斷和治療。

    疾病預(yù)測對于預(yù)防疾病的發(fā)生和減少疾病的損傷具有重要意義。Yang等[7]使用患者的健康記錄來預(yù)測糖尿病潛在并發(fā)癥,同時發(fā)現(xiàn)了并發(fā)癥與實驗室類型之間的隱含關(guān)系。He等[8]使用EMR數(shù)據(jù)集預(yù)測肺癌術(shù)后并發(fā)癥,并同時從數(shù)據(jù)集中提取關(guān)鍵變量。

    具有預(yù)測診斷標(biāo)簽和藥物信息的EMR可幫助自動助理預(yù)測疾病診斷并為醫(yī)生提供快速診斷參考。Nee等[9]使用大型EMR文本數(shù)據(jù)集對每種疾病的EMR上下文進行建模,并在EMR中執(zhí)行準(zhǔn)確的疾病診斷預(yù)測。Wright等[10]使用數(shù)據(jù)挖掘方法從醫(yī)學(xué)數(shù)據(jù)集中獲取有用的關(guān)系和規(guī)則集,以預(yù)測接下來要開哪些藥。

    2.2. 傳統(tǒng)模型在醫(yī)療圖像中的應(yīng)用

    傳統(tǒng)的醫(yī)學(xué)圖像使用手工構(gòu)造特征,通過淺模型進行分類與分割。Scott等[11]在2003年使用梯度取向、拐角和邊緣強度檢測雙能量X射線圖像中的椎骨。此外,區(qū)域劃分和合并也是基于區(qū)域的方法中的較為知名的技術(shù)。Manousakas等[12]應(yīng)用區(qū)域劃分和合并技術(shù),試圖克服在MRI中使用均勻性測量時遇到的困難。Zhao等[13]介紹了形態(tài)學(xué)的基本數(shù)學(xué)理論和操作,并提出了一種新穎的形態(tài)學(xué)邊緣檢測方法,以區(qū)分帶有椒鹽噪聲的CT圖像中的肺部邊緣。實驗結(jié)果表明,與2006年最佳邊緣檢測方法相比,該方法在醫(yī)學(xué)圖像降噪和邊緣檢測方面均更有效。除此之外,K-means聚類也曾被Kaus等[14]應(yīng)用在心臟MRI中進行自動左心室的分割。Cordes等[15]通過使用分層聚類來衡量MRI的連通性,該方法可以檢測到低頻波動的相似性,結(jié)果表明,類似于已知神經(jīng)元連接的層次聚類能夠檢測心臟連通模式。2006年,Pohl等[16]提出了一種將有符號距離圖嵌入線性對數(shù)優(yōu)勢空間的方法,可以解決建模問題。上述方法專注于區(qū)域、邊緣和聚類等簡單手段,所以它們在現(xiàn)實世界數(shù)據(jù)上往往性能有限[17]。

    2.3. 深度學(xué)習(xí)在醫(yī)療圖像中的應(yīng)用

    在計算機輔助診斷中,深度學(xué)習(xí)現(xiàn)已廣泛用于醫(yī)學(xué)圖像識別[18,19]。深度學(xué)習(xí)的基本結(jié)構(gòu)是CNN,它的結(jié)構(gòu)有三層,即卷積層、池化層和全連接層。為了開發(fā)基于CNN的強大AI算法,我們通常需要大量帶標(biāo)注的數(shù)據(jù)。

    醫(yī)學(xué)圖像的標(biāo)準(zhǔn)化收集并不像收集一般自然圖像那樣容易。但是,如今幾個公共醫(yī)學(xué)圖像數(shù)據(jù)庫和多中心數(shù)據(jù)收集可以幫助解決該問題。某些類型的醫(yī)學(xué)圖像數(shù)據(jù)可以被大量收集,如X射線圖像、CT、心電圖和病理圖像。通過使用這些大數(shù)據(jù),基于CNN的AI算法可以對CT圖像進行解剖結(jié)構(gòu)分割[20],對胸部X射線圖像正?;虍惓=Y(jié)果進行分類[21],對肺癌或乳腺癌進行篩查[22,23],檢測出顱腦CT掃描中的危險情況[24],使用基于生成對抗網(wǎng)絡(luò)(GAN)的模型對肝臟病變進行分類[25],進行心臟病篩查[26,27],以及在病理圖像中探測淋巴結(jié)轉(zhuǎn)移等[27,28]。

    在眼科領(lǐng)域,由于眼底照相和OCT圖像易于收集,因此基于CNN的AI算法主要應(yīng)用領(lǐng)域是探查視網(wǎng)膜疾病,如糖尿病性視網(wǎng)膜病變、年齡相關(guān)性黃斑變性和青光眼[29-31]。

    當(dāng)前,人工智能輔助的醫(yī)學(xué)診斷系統(tǒng)主要應(yīng)用于醫(yī)學(xué)成像領(lǐng)域。如果疾病的診斷需要依賴自然觀察,則其主要取決于醫(yī)生的個人經(jīng)驗。皮膚病損是一個例子,當(dāng)前的AI算法可以將數(shù)碼皮膚照片中的惡性黑色素瘤與良性病變區(qū)分開[2]。角膜病也是一個例子,眼科醫(yī)生可以使用裂隙燈顯微鏡來獲得正確的診斷,但迄今為止,還沒有有關(guān)利用AI來提高角膜病的診斷準(zhǔn)確性的研究。

    3. 方法

    3.1. 圖像數(shù)據(jù)集構(gòu)建

    經(jīng)機構(gòu)倫理委員會批準(zhǔn),此研究的圖像數(shù)據(jù)集包括了1998年5月至2018年在浙江大學(xué)醫(yī)學(xué)院附屬邵逸夫醫(yī)院眼科通過裂隙燈顯微鏡獲得的115 408張臨床數(shù)字圖像,這些圖像來自10 609例89種角膜病患者。臨床圖像由兩種類型的裂隙燈顯微鏡拍攝,即蔡司裂隙燈顯微鏡SL 130(德國卡爾·蔡司公司),集成了SL Cam for Imaging模塊,每個圖像的分辨率為1024×768像素;以及附帶數(shù)碼相機Unit DC-1的Topcon裂隙燈顯微鏡(日本拓普康公司),其圖像分辨率為1740×1536像素或2048×1536像素。

    在數(shù)據(jù)集中,拍攝的活躍期角膜感染的圖像,包括細(xì)菌性角膜炎(BK)、真菌性角膜炎(FK)和單純皰疹病毒基質(zhì)性角膜炎(HSK)的圖像,被選擇用于算法分類的訓(xùn)練或測試集。來自角膜感染患者的所有圖像均帶有明確的臨床診斷標(biāo)注,至少有以下證據(jù)中的兩種:①角膜感染的臨床表現(xiàn)如圖1(a)所示。②通過相關(guān)的診斷性單藥或聯(lián)合藥物治療,影響和終止了角膜感染的進展,從而最終治愈。③感染部位樣本的病原體鑒定:在細(xì)菌和真菌感染中,通過顯微鏡檢查涂片或生物培養(yǎng)確認(rèn)病原診斷;在病毒感染中,通過聚合酶鏈反應(yīng)(PCR)評估淚液或角膜刮片組織樣本進行病原診斷。除了角膜感染的類別外,患有具有相似視覺特征的其他角膜病患者的圖像被歸類為“其他”類別。此類別包括各種角膜變性、泡性角結(jié)膜炎、各種角膜腫瘤、角膜乳頭狀瘤、角膜退行性變,甚至還包括棘阿米巴角膜炎。每個類別的代表性圖像如圖1(a)所示。

    最終數(shù)據(jù)集包含來自867例患者的2284幅圖像。訓(xùn)練集包括來自747例患者的387張隨機選擇的BK圖像、519張FK圖像、488張HSK圖像和528張其他角膜病圖像。測試集由120例患者隨機選擇的86張BK圖像、97張FK圖像、51張HSK圖像和128張其他診斷圖像組成。為了評估眼科醫(yī)生的疾病分類表現(xiàn),選擇測試集中每位患者的首次診斷圖像以構(gòu)建評估眼科醫(yī)生診斷水平的數(shù)據(jù)集(即總共使用了120張圖像來評估眼科醫(yī)生的表現(xiàn))。

    圖1. 典型裂隙燈顯微圖像與我們提出的四類角膜疾病的SOS模型中嵌入特征的t-分布隨機鄰域嵌入(t-SNE)的表示可視化。(a)細(xì)菌性角膜炎(BK)、真菌性角膜炎(FK)、單純皰疹病毒基質(zhì)性角膜炎(HSK)和其他角膜炎(上述三類角膜疾病之外的角膜炎)的典型裂隙燈顯微圖像。不同類別的疾病或同一類別疾病的不同階段表現(xiàn)出不同的視覺特征。(b)SOS模型通過t-SNE嵌入每種疾病類別的二維空間獲得的深度序列特征。t-SNE用于可視化高維數(shù)據(jù),這些數(shù)據(jù)是經(jīng)過診斷驗證的攝影測試集的SOS模型中的特征表示(362幅圖像)。彩色的點云代表疾病的不同類別,顯示了算法如何將疾病分組為不同的簇。插圖顯示了與各個點相對應(yīng)的圖像。

    3.2. 基于深度序列特征的診斷模型

    如前所述,我們發(fā)明了一種使用圖像序列特征的感染性角膜炎分類模型訓(xùn)練方法。為了展示我們提出的這一方法的優(yōu)勢,我們將其與另兩類方法——基于完整圖像特征的學(xué)習(xí)方法和基于圖像子塊特征的方法——進行了對比。

    基于完整圖像特征的深度學(xué)習(xí)模型接收原始的無標(biāo)注疾病圖像,直接由CNN進行分類。為解決訓(xùn)練圖像有限的問題,我們采用了遷移學(xué)習(xí)的方法[32,33]。在實驗中,我們選擇了如下三種經(jīng)典的圖像分類模型結(jié)構(gòu):VGG-16(由英國牛津大學(xué)計算幾何小組提出)[34]、GoogleLeNet-v3(由谷歌公司修改的LeNet模型的第三個版本)[35]以及稠密卷積網(wǎng)絡(luò)(DenseNet)[36]。

    在基于圖像子塊特征深度學(xué)習(xí)模型當(dāng)中,眼前節(jié)圖像由手工進行初步標(biāo)注,分割出四塊:角膜感染病灶區(qū)、角膜感染灶旁區(qū)、結(jié)膜充血區(qū)和前房積膿區(qū)。在該類方法中我們同樣使用了三種遷移模型結(jié)構(gòu),即VGG-16、GoogleLeNet-v3和DenseNet。在對圖像的每個子塊完成分類后,通過多數(shù)投票的方法給出整張疾病圖像的分類結(jié)果。

    在基于序列的深度學(xué)習(xí)模型中,對于每張圖片,如果其中存在病變區(qū)域,模型的關(guān)注點就會被設(shè)置在病變區(qū)域上。病變區(qū)的最小外接圓被進一步劃分為由小到大的K個圓環(huán),如圖2所示。落在由內(nèi)到外的第i個圓環(huán)上的所有圖像子塊組成一個子塊集合Si,所有集合組成一個由內(nèi)而外的序列{S1,S2, ...,SK}。為解決標(biāo)注數(shù)據(jù)不足的問題,在訓(xùn)練階段,我們通過隨機失活機制來隨機剔除每個集合中的一些元素,這樣我們能夠獲得更多的序列來擴大訓(xùn)練集多樣性,使得模型的訓(xùn)練更加穩(wěn)定。

    圖2. 采樣以及如何將病變區(qū)域分成K個集合。圓圈代表每個集合的邊界,正方形代表采樣區(qū)域。要注意的是為避免圖片過度重疊,僅顯示了一半的采樣區(qū)域。

    我們利用一個深層的殘差CNN(如DenseNet),通過編碼器-解碼器框架[37-39]來提取每個集合中的每個子塊的表征。卷積編碼器能夠?qū)⒌趇個集合中的第j個子塊pij轉(zhuǎn)化為一個能夠表示其性質(zhì)的表征向量fij,得到圖像子塊的表征的集合{F1,F2, ...,FK}。對每個集合Fi,通過最大池化計算獲得整個子塊集合的表征f′i,表示每個集合的整體性質(zhì)。由于病變區(qū)由內(nèi)到外的子塊集合構(gòu)成了一個序列,我們可以使用長短時記憶(long short-term memory, LSTM)模型[37]——一種經(jīng)典的序列學(xué)習(xí)模型,來將表征序列{f′1,f′2, ...,f′K}轉(zhuǎn)化為最終用于分類的表征。最終的圖像特征可以通過一個全連接層解碼,通過softmax計算得到每個疾病類別的分類概率。圖1(b)展示了每種疾病的圖像表征在二維空間上的分布。圖3展示了整個系統(tǒng)的結(jié)構(gòu)。經(jīng)由獲得的概率分布與真實標(biāo)簽的損失函數(shù),通過反向傳播精調(diào)模型參數(shù)[40,41]。

    3.3. 眼科醫(yī)生基于圖像的診斷能力分析

    我們從全國各地招募眼科醫(yī)生,測試他們基于圖像進行診斷的能力,并將其與已開發(fā)的深度學(xué)習(xí)方法進行比較。呈現(xiàn)給眼科醫(yī)生的圖像和測試集中每位患者的經(jīng)過診斷驗證的圖像都是從初次就診時隨機選擇的(即總共120張圖像)。所招募的眼科醫(yī)生具有不同的學(xué)術(shù)頭銜(從住院醫(yī)生到高年資醫(yī)生,甚至到醫(yī)學(xué)院的臨床教授)、工作機構(gòu)(從大學(xué)醫(yī)學(xué)院的教學(xué)醫(yī)院到公立市級醫(yī)院再到社區(qū)診所)和專業(yè)經(jīng)驗(分為1~5年、6~10年、11~15年、16~20年以及20年以上)。我們總共招募了421名眼科醫(yī)生。

    圖3. 病變區(qū)域進行深度序列特征學(xué)習(xí)的過程。對于每張裂隙燈顯微圖像,病變區(qū)域外接圓被劃分為K個圓環(huán)(直觀起見,此處K = 3)。我們從每個圓環(huán)(從最內(nèi)層到最外層)上采樣,采樣的區(qū)域用于生成集合序列。隨后通過最大池化和長短期記憶(LSTM)模型進行序列特征學(xué)習(xí)。f:采樣區(qū)域的特征;f’:集合的序列特征。

    眼科醫(yī)生手動檢視圖像從而進行基于圖像的診斷,主要遵循兩步流程。第一步,眼科醫(yī)生僅依靠圖像進行診斷,將來自測試集中每位患者的首次診斷圖像中的四類角膜病圖像,即BK、FK、HSK和其他角膜病的圖像提供給眼科醫(yī)生,眼科醫(yī)生通過手動選擇為每個圖像做出診斷。第二步,為眼科醫(yī)生提供與每個圖像相關(guān)的其他標(biāo)準(zhǔn)化和結(jié)構(gòu)化的醫(yī)學(xué)信息,包括簡短的病史、起病時間、疼痛程度和復(fù)發(fā)情況(如果有)以及藥物使用史。第三步,要求眼科醫(yī)生通過手動檢視并考慮其他醫(yī)學(xué)信息,對每個圖像做出診斷決定。所有眼科醫(yī)生均獨立執(zhí)行此程序,沒有時間限制。

    3.4. 統(tǒng)計學(xué)分析

    鑒于不同的學(xué)術(shù)頭銜、工作機構(gòu)和專業(yè)經(jīng)驗導(dǎo)致的置信度不同,使用社會科學(xué)統(tǒng)計軟件包(SPSS 18.0版;美國Cary公司)對眼科醫(yī)生的診斷數(shù)據(jù)進行統(tǒng)計分析。眼科醫(yī)生診斷準(zhǔn)確度的平均水平以(均值±標(biāo)準(zhǔn)差)%表示。使用Kolmogorov-Smirnov檢驗驗證數(shù)據(jù)的正態(tài)性。根據(jù)數(shù)據(jù)正態(tài)性,使用單因素方差分析(ANOVA)分析了不同醫(yī)院級別和職稱組之間的診斷準(zhǔn)確性差異。最小顯著差異法用于參數(shù)變量的事后分析。使用皮爾遜相關(guān)系數(shù)測試了診斷準(zhǔn)確性和專業(yè)經(jīng)驗?zāi)曩Y之間的相關(guān)性。采用逐步法進行多元線性回歸分析,從學(xué)術(shù)頭銜、醫(yī)院層級和專業(yè)經(jīng)驗?zāi)曩Y等方面分析了統(tǒng)計學(xué)因素的影響。進行了配對t檢驗(針對正態(tài)分布的變量)和Wilcoxon符號秩檢驗(針對非正態(tài)分布的變量),以確定在有無附加醫(yī)學(xué)信息的情況下,醫(yī)生在診斷準(zhǔn)確性上是否存在顯著差異。所有檢驗的顯著性水平設(shè)置為0.05。

    4. 結(jié)果

    4.1. 不同深度模型的性能對比

    完整圖像疾病的深度模型目前在醫(yī)學(xué)圖像診斷任務(wù)中非常流行,此類方法直接將原始醫(yī)學(xué)圖像交由CNN處理。我們的工作中選擇了VGG-16、GoogleLeNet-v3和DenseNet三種經(jīng)典模型結(jié)構(gòu),并對比了其對BK、FK、HSK的診斷準(zhǔn)確率,如表1所示??紤]到直接將完整圖像傳入CNN可能會含有一些無關(guān)的信息,我們進一步基于這三種結(jié)構(gòu)設(shè)計了基于圖像子塊的模型[42,43]。在基于圖像子塊的模型中,圖像包含由人工標(biāo)注的角膜感染病灶區(qū)、角膜感染灶旁區(qū)、結(jié)膜充血區(qū)和前房積膿區(qū)等區(qū)域信息,我們并不傳入完整的原始圖像,而是使用包含這些區(qū)域的圖像子塊進行分類。三種基于子塊的模型對每個子塊的識別準(zhǔn)確率分別能夠達(dá)到49.62%、51.52%和60%。在每個圖像子塊分類后,通過多數(shù)投票的方法給出整張圖像的診斷結(jié)果?;趫D像子塊的模型對圖像的分類準(zhǔn)確率分別可以達(dá)到52.50%、55.52%和66.30%,如表1所示。

    最后,我們使用基于圖像序列特征的模型,該方法能夠保持醫(yī)學(xué)圖像中潛在的空間結(jié)構(gòu)信息。如前所述,子塊集合序列的表征是通過由內(nèi)而外的順序排列獲得的(即SOS),使用SOS提取特征能達(dá)到78.73%的分類準(zhǔn)確率。除了產(chǎn)生集合序列外,我們還生成了隨機排序的子塊(ROP)和順序排列的子塊(SOP)來獲得特征序列。ROP通過隨機選擇產(chǎn)生子塊序列,而SOP按照從內(nèi)到外的順序采樣序列(但是沒有將屬于同一圓環(huán)的子塊聚成一個集合)。實驗結(jié)果表面,ROP方法能夠達(dá)到74.23%的準(zhǔn)確率(BK 75.29%、FK 68.04%、HSK 82.35%)。SOP方法達(dá)到75.14%的準(zhǔn)確率。這些實驗結(jié)果表明我們基于序列的深度學(xué)習(xí)模型在基于圖像的角膜病分類診斷任務(wù)上能獲得最好的表現(xiàn)。

    表1 不同深度模型在測試集上的診斷能力對比

    4.2.與眼科醫(yī)生的診斷正確率對比

    我們使用數(shù)據(jù)集評估了本文中考慮的所有算法,以比較每種算法與眼科醫(yī)生的表現(xiàn)。表2列出了所有算法的準(zhǔn)確性以及該數(shù)據(jù)集(120幅圖像)上眼科醫(yī)生的平均表現(xiàn)。表3列出了眼科醫(yī)生在臨床圖像診斷中的表現(xiàn)。從全國各地共招募421名眼科醫(yī)生參加這項研究。在沒有參考任何額外醫(yī)療信息的眼科醫(yī)生的平均準(zhǔn)確度為(49.27 ± 11.5)%(范圍:20.00%~86.67%),遠(yuǎn)低于AI深度學(xué)習(xí)模型所獲得的準(zhǔn)確性。例如,SOS算法的診斷準(zhǔn)確度達(dá)到80.00%,其中BK、FK和HSK的準(zhǔn)確度分別為53.33%、83.33%和93.33%(表2)。圖4描繪了受試者工作特征(ROC)曲線、SOS模型的混淆矩陣以及眼科醫(yī)生的表現(xiàn)。ROC曲線是分類模型的一種可視化。曲線下的面積(AUC)是性能的衡量,最大值為1。如果眼科醫(yī)生的敏感度-特異性點位于分類模型的曲線下方,則該模型可達(dá)到優(yōu)于眼科醫(yī)生的性能。

    表2 在樣本量為120的測試集上眼科醫(yī)生與深度模型的對比

    表3 根據(jù)醫(yī)院級別、工作年限和眼科醫(yī)生的專業(yè)職稱的平均分類準(zhǔn)確率

    對于工作地點對眼科醫(yī)生診斷水平的影響,教學(xué)醫(yī)院的眼科醫(yī)生的表現(xiàn)比市級醫(yī)院和社區(qū)診所的要好(P< 0.001),而市級醫(yī)院和社區(qū)診所之間沒有統(tǒng)計差異(P= 0.226)。專業(yè)級別較高的眼科醫(yī)生在診斷臨床圖像時具有更高的準(zhǔn)確性,如高級職稱和中級職稱醫(yī)師的表現(xiàn)優(yōu)于住院醫(yī)師(分別為P< 0.001和P= 0.003),但高級職稱和中級職稱之間差異無統(tǒng)計學(xué)意義(P= 0.071)。在執(zhí)業(yè)年資和診斷準(zhǔn)確性之間未發(fā)現(xiàn)顯著相關(guān)性(P=0.084)。

    綜合考慮醫(yī)院等級和醫(yī)生等級的因素,教學(xué)醫(yī)院的高級職稱眼科醫(yī)生的表現(xiàn)[準(zhǔn)確性為(57.08 ± 12.02)%,范圍:33.33%~86.67%]好于社區(qū)診所的住院醫(yī)生[準(zhǔn)確性為(41.99 ± 10.51)%,范圍:22.50%~63.33%]。逐步多元回歸分析得出影響診斷準(zhǔn)確性的三個模型。模型1(R2= 0.062)僅具有醫(yī)院水平的因素(β= 0.254,P<0.001);模型2(R2= 0.100)具有醫(yī)院水平(β= 0.239,P< 0.001)和專業(yè)職稱(β= 0.200,P< 0.001)的因素;模型3(R2= 0.109)具有所有三個因素:醫(yī)院水平(β= 0.227,P< 0.001)、專業(yè)職稱(β= 0.326,P< 0.001)和工作年限(β= ?0.164,P= 0.024)。

    當(dāng)進一步為眼科醫(yī)生提供每張圖像的附加醫(yī)學(xué)信息時,包括簡短的病史、發(fā)病時間、疼痛程度和復(fù)發(fā)發(fā)作(如果有)以及用藥史,平均總診斷準(zhǔn)確性從49.27%提高到57.16%,差異具有統(tǒng)計學(xué)意義(Wilcoxon符號秩檢驗,P< 0.001)。詳細(xì)而言,附加醫(yī)學(xué)信息后,BK的診斷準(zhǔn)確性從46.55%提高到55.55%(P< 0.001),F(xiàn)K從45.56%提高到56.28%(P< 0.001),HSK從65.01%提高到73.25%(P< 0.001),404名醫(yī)生的平均總準(zhǔn)確性提高了8.28%,9名醫(yī)生的準(zhǔn)確性下降了2.13%,其他8名醫(yī)生的準(zhǔn)確性保持不變。

    圖4. SOS模型與眼科醫(yī)生表現(xiàn)的ROC曲線及混淆矩陣對比圖。(a)SOS模型的ROC曲線;(b)、(c)在用于診斷水平測試的數(shù)據(jù)集上眼科醫(yī)生與SOS模型的混淆矩陣;(d)~(f)對應(yīng)的BK、FK與HSK等各子類的ROC曲線。

    5. 討論

    一般而言,人類通過視覺、聽覺、觸覺、味覺和嗅覺對事物進行判斷,并以此對事物進行適當(dāng)?shù)姆诸怺44]。其中視覺感知發(fā)揮著最重要的作用[45],視覺知識可以描述空間形狀、大小和相關(guān)性以及顏色和紋理之間的關(guān)系[46]。醫(yī)生對疾病的診斷主要依靠觀察和推理。在所有人類疾病中,角膜病具有最直接、最顯著的視覺感知變化;因為健康角膜具有完全透明的獨特特征,這與病理狀況下角膜內(nèi)外圖像的變化形成鮮明對比。人類專家通過對圖像的理解和分析來進行角膜病的診斷決策,這可能是AI輔助人類最合適的任務(wù)。

    一般而言,深度學(xué)習(xí)是由大量帶標(biāo)注的數(shù)據(jù)驅(qū)動的[47,48]。然而,尚不清楚有多少臨床圖像訓(xùn)練數(shù)據(jù)足以被開發(fā)用于診斷臨床疾病的AI系統(tǒng)。我們中心收集和記錄具有臨床圖像的角膜病例已有20余年,但是,如果按照每種疾病類別對所有圖像進行標(biāo)注,那么在最常見的疾病類別中可以有數(shù)千張臨床圖像,而某些罕見疾病類別中只有少數(shù)臨床圖像。每個角膜病類別中帶標(biāo)注數(shù)據(jù)的不平衡性導(dǎo)致我們將注意力集中在最常見的疾?。ㄈ绺腥拘越悄げ。┥?,以初步開發(fā)本研究中的AI診斷系統(tǒng)。

    在本研究中,我們證明了通過裂隙燈顯微鏡拍攝的臨床圖像,通過CNN進行的深度學(xué)習(xí)可以用于角膜感染性疾病的臨床診斷。我們共評估了3組9種深度學(xué)習(xí)架構(gòu),以開發(fā)用于角膜感染性疾病的圖像診斷系統(tǒng)。從圖像級別和圖像子塊級別模型的結(jié)果可以看出,盡管只有4個類別,但這仍是一個困難的任務(wù),尤其是對于VGG-16和GoogLeNet-v3。這兩種結(jié)構(gòu)在圖像子塊分類中均表現(xiàn)不佳,圖像子塊之間的投票并未顯著改善其性能。相比之下,DenseNet的圖像子塊分類達(dá)到60%,投票后達(dá)到66.3%。研究表明,如果模型在圖像子塊分類中表現(xiàn)得足夠好,那么從感染性病變區(qū)域著眼于圖像子塊可以產(chǎn)生比著眼于圖片整體更高的性能。ROP方法可被視為除表決之外組合圖像子塊特征的另一種方法。結(jié)果表明,即使沒有空間信息,采用適當(dāng)?shù)慕M合方法也可以進一步改善圖像子塊級模型。我們發(fā)現(xiàn),總體而言SOS是用于角膜病中僅依賴圖像進行診斷的最有前景的方法。SOS比其他方法更好的一個可能原因是在這種深度學(xué)習(xí)模型中,直接實現(xiàn)了對臨床圖像空間結(jié)構(gòu)的適當(dāng)利用。SOP表現(xiàn)不佳,是因為它沒有考慮病變區(qū)域的環(huán)形結(jié)構(gòu)。據(jù)我們所知,本研究首次提出了一種角膜病分類深度學(xué)習(xí)模型,該模型在僅依靠圖像進行診斷時比人類眼科醫(yī)生具有更高的準(zhǔn)確性。在本研究中,一般眼科醫(yī)生在僅依賴圖像的角膜病診斷中的表現(xiàn)要比AI系統(tǒng)差。毫無疑問,錯誤的診斷可能造成長期使用不適當(dāng)?shù)乃幬?,?dǎo)致識別特征變得更為模糊,使醫(yī)生診斷更加困難[6]。我們研究中的多元回歸分析表明,就職稱、工作機構(gòu)和年資而言,這3個統(tǒng)計學(xué)因素對診斷性能會有影響,但3個模型中的確定系數(shù)較低。這表明上述因素可能無法如實全面地決定眼科醫(yī)生對角膜病的診斷準(zhǔn)確性,或者影響診斷性能的因素可能非常復(fù)雜,無法通過上述3個因素簡單地準(zhǔn)確總結(jié)。因此,如果AI可以幫助臨床醫(yī)生以更高的診斷準(zhǔn)確性顯著提高臨床能力,那么這將極大地使角膜病患者受益,同時節(jié)省了醫(yī)療資源并減輕社會負(fù)擔(dān)。目前全世界仍有450萬人正遭受由角膜病引起角膜混濁導(dǎo)致的中度至重度視力障礙,特別是在發(fā)展中國家。提高診斷準(zhǔn)確性可有兩種途徑,一是完善醫(yī)師培訓(xùn)體系,加強醫(yī)師專業(yè)教育和培訓(xùn);另一種即開發(fā)實用的AI系統(tǒng)以協(xié)助診斷。我們的研究表明,通過使用臨床圖像來開發(fā)AI系統(tǒng)以提高角膜病診斷準(zhǔn)確性是切實可行的。在測試眼科醫(yī)生的表現(xiàn)時,我們發(fā)現(xiàn)當(dāng)向醫(yī)生提供圖像和病史時,診斷準(zhǔn)確性比單純提供圖像有所提高(從49.27%增至57.16%,P< 0.001)。該結(jié)果表明,其他信息可以幫助進一步提高診斷表現(xiàn),AI診斷系統(tǒng)也可能如此。研究表明,將數(shù)據(jù)驅(qū)動的機器學(xué)習(xí)與人類知識相結(jié)合可以有效地開發(fā)出可解釋、強大且通用的AI [49];諸如病史之類的信息可能包含類似于人的常識,可以使模型利用有限的訓(xùn)練數(shù)據(jù)來解決許多不同的任務(wù)。為了改善我們的AI診斷系統(tǒng)以提高準(zhǔn)確性,在未來的工作中可能需要設(shè)計一種多模態(tài)學(xué)習(xí)模型(即視覺和非視覺信息的有效組合),或是一個更合適的序列學(xué)習(xí)模型。

    不可否認(rèn)的是,在現(xiàn)階段我們的AI診斷準(zhǔn)確性只能由我們所收集的有限圖像數(shù)據(jù)來證實[50],并通過與使用相同臨床圖像的眼科醫(yī)生進行對比。這種AI系統(tǒng)在輔助醫(yī)生臨床實踐中的實際應(yīng)用,需要在將來進行更大規(guī)模的深入臨床評估[51]。

    6. 結(jié)論

    傳染性角膜炎是最常見的眼科疾病,可能導(dǎo)致失明。眼科醫(yī)生通過觀察裂隙燈圖像診斷疾病,利用計算機輔助圖像分析算法方便診斷。在這項工作中,我們提出了一個序列水平的深部模型,端到端診斷傳染性角膜炎。具體地說,利用深度卷積網(wǎng)絡(luò)良好的特征提取性能,首先提取角膜區(qū)域的細(xì)節(jié)模式,然后將局部特征分成符合空間結(jié)構(gòu)的有序集合,學(xué)習(xí)角膜圖像的全局表示并進行診斷。我們收集了超過10 000名患者的110 000張圖像。在此基礎(chǔ)上,充分的實驗對比結(jié)果表明,該模型是一種更為可行的結(jié)構(gòu),比傳統(tǒng)的CNN具有更好的診斷性能。另外,通過與400多位專業(yè)眼科醫(yī)生的診斷結(jié)果進行對比發(fā)現(xiàn),我們的模式可以大大超過專業(yè)人士的平均水平,達(dá)到頂級眼科醫(yī)生的水平表現(xiàn)。據(jù)我們所知,這是第一個關(guān)于感染性角膜炎診斷的研究,我們的研究有力地證明了使用人工智能進行這些類型疾病的臨床輔助診斷的潛力。

    致謝

    本研究得到了浙江省衛(wèi)生健康委員會(WKJZJ-1905、2018ZD007)、浙江省重點研究開發(fā)項目(2018C03082)和國家自然科學(xué)基金(61625107)的支持。感謝張仲非教授的意見和建議。

    Compliance with ethics guidelines

    Yesheng Xu, Ming Kong, Wenjia Xie, Runping Duan,Zhengqing Fang, Yuxiao Lin, Qiang Zhu, Siliang Tang, Fei Wu, and Yu-Feng Yao declare that they have no conflict of interest or financial conflicts to disclose.

    猜你喜歡
    眼科醫(yī)生子塊角膜炎
    基于八叉樹的地震數(shù)據(jù)多級緩存方法
    基于八叉樹的地震數(shù)據(jù)分布式存儲方法研究
    基于特征值算法的圖像Copy-Move篡改的被動取證方案
    送別
    海燕(2020年3期)2020-11-22 08:56:59
    如何區(qū)分角膜炎和結(jié)膜炎
    基于波浪式矩陣置換的稀疏度均衡分塊壓縮感知算法
    耳眼科醫(yī)生
    益壽寶典(2018年36期)2018-01-27 19:59:04
    眼科醫(yī)生給小朋友的信
    蒙西醫(yī)結(jié)合治療病毒性角膜炎的臨床療效
    How Cats See The World
    香蕉av资源在线| 97超视频在线观看视频| 免费观看在线日韩| 99久国产av精品国产电影| 国内精品一区二区在线观看| 麻豆精品久久久久久蜜桃| 亚洲av第一区精品v没综合| 国产精品一及| 草草在线视频免费看| 国产熟女欧美一区二区| 精品福利观看| 国产午夜福利久久久久久| 久久人妻av系列| 可以在线观看的亚洲视频| 亚洲精品色激情综合| 香蕉av资源在线| 99久久精品热视频| 亚洲av美国av| 一a级毛片在线观看| 午夜a级毛片| 亚洲av美国av| 麻豆国产97在线/欧美| 丰满人妻一区二区三区视频av| 国产精品亚洲一级av第二区| 少妇高潮的动态图| 精品久久久久久久人妻蜜臀av| 亚洲久久久久久中文字幕| 成年av动漫网址| 色视频www国产| 一边摸一边抽搐一进一小说| 亚洲欧美日韩卡通动漫| 露出奶头的视频| 久久久久久伊人网av| 欧美日韩精品成人综合77777| 日韩欧美一区二区三区在线观看| 欧美bdsm另类| 黑人高潮一二区| 俄罗斯特黄特色一大片| 一个人免费在线观看电影| 婷婷六月久久综合丁香| 国产一区二区三区av在线 | 亚洲成人久久性| 免费大片18禁| 美女xxoo啪啪120秒动态图| 黑人高潮一二区| 我要搜黄色片| 高清毛片免费观看视频网站| 国产视频内射| 亚洲不卡免费看| 欧美潮喷喷水| 国产综合懂色| 亚洲久久久久久中文字幕| 欧美成人精品欧美一级黄| 日韩 亚洲 欧美在线| 欧美色欧美亚洲另类二区| 嫩草影院入口| 亚洲国产精品合色在线| 12—13女人毛片做爰片一| 99热网站在线观看| av在线天堂中文字幕| a级毛片a级免费在线| 香蕉av资源在线| 精品人妻偷拍中文字幕| 内射极品少妇av片p| 黄色欧美视频在线观看| 又爽又黄无遮挡网站| 一进一出抽搐动态| 一进一出抽搐动态| 欧美3d第一页| 亚洲三级黄色毛片| 色哟哟哟哟哟哟| 国产精品一区二区三区四区久久| 中国美女看黄片| 一级黄片播放器| 桃色一区二区三区在线观看| 国产中年淑女户外野战色| 一级毛片我不卡| 日本a在线网址| 噜噜噜噜噜久久久久久91| 露出奶头的视频| 国产真实伦视频高清在线观看| 亚洲真实伦在线观看| 国产一区二区在线观看日韩| 五月玫瑰六月丁香| 久久久色成人| 亚洲人成网站高清观看| 99久久久亚洲精品蜜臀av| 亚洲精品色激情综合| 91久久精品国产一区二区成人| 成人漫画全彩无遮挡| 中国国产av一级| 久久精品影院6| 高清日韩中文字幕在线| 噜噜噜噜噜久久久久久91| АⅤ资源中文在线天堂| 午夜亚洲福利在线播放| 欧美成人精品欧美一级黄| 少妇高潮的动态图| 99热网站在线观看| 给我免费播放毛片高清在线观看| 中文在线观看免费www的网站| 欧美日韩综合久久久久久| 五月伊人婷婷丁香| 毛片一级片免费看久久久久| aaaaa片日本免费| 在线播放国产精品三级| 久久精品国产亚洲av涩爱 | 午夜精品一区二区三区免费看| 午夜免费激情av| 国产亚洲精品久久久com| 久久99热6这里只有精品| av专区在线播放| 欧美激情国产日韩精品一区| 国产熟女欧美一区二区| 国产精品免费一区二区三区在线| 亚洲欧美日韩卡通动漫| 最新在线观看一区二区三区| 午夜视频国产福利| 久久久久国产精品人妻aⅴ院| 亚洲电影在线观看av| 欧美不卡视频在线免费观看| 特级一级黄色大片| 精品人妻视频免费看| 久久久精品欧美日韩精品| 亚洲国产色片| 国产综合懂色| 在线观看66精品国产| 久久久午夜欧美精品| 亚洲成人精品中文字幕电影| 久久人妻av系列| 亚洲av成人精品一区久久| 亚洲国产日韩欧美精品在线观看| 特级一级黄色大片| 免费av不卡在线播放| 99久久久亚洲精品蜜臀av| 一区福利在线观看| 综合色av麻豆| 18禁在线无遮挡免费观看视频 | 免费无遮挡裸体视频| 深夜a级毛片| 精品久久国产蜜桃| 亚洲婷婷狠狠爱综合网| 欧美成人精品欧美一级黄| 国产精品野战在线观看| 国产极品精品免费视频能看的| 国产高清激情床上av| 超碰av人人做人人爽久久| 日韩亚洲欧美综合| 日本黄色片子视频| 精品国内亚洲2022精品成人| 国产精品国产高清国产av| 午夜福利高清视频| 男女啪啪激烈高潮av片| 日韩精品青青久久久久久| 日韩欧美一区二区三区在线观看| 日韩欧美一区二区三区在线观看| 美女xxoo啪啪120秒动态图| 亚洲欧美日韩高清专用| 91精品国产九色| 亚洲av成人精品一区久久| 十八禁国产超污无遮挡网站| 亚洲不卡免费看| 赤兔流量卡办理| 插逼视频在线观看| 日产精品乱码卡一卡2卡三| 久久精品综合一区二区三区| 久久精品综合一区二区三区| 欧美成人免费av一区二区三区| h日本视频在线播放| 国产成年人精品一区二区| 国产成人aa在线观看| 搞女人的毛片| 日日撸夜夜添| 美女 人体艺术 gogo| 欧美最黄视频在线播放免费| 婷婷精品国产亚洲av| 国产三级在线视频| 国产一区二区三区在线臀色熟女| 乱人视频在线观看| 禁无遮挡网站| 亚洲天堂国产精品一区在线| 人妻丰满熟妇av一区二区三区| 插逼视频在线观看| 久久久久久伊人网av| 少妇裸体淫交视频免费看高清| 午夜免费激情av| 99热这里只有是精品在线观看| 国产精品电影一区二区三区| 欧美高清成人免费视频www| 国产高潮美女av| 国产乱人视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站高清观看| 日韩欧美精品v在线| 性色avwww在线观看| 亚洲人成网站在线观看播放| 久久久成人免费电影| 男插女下体视频免费在线播放| 少妇的逼水好多| 精品国产三级普通话版| 久99久视频精品免费| 成人综合一区亚洲| 亚洲人成网站在线播放欧美日韩| 国产高潮美女av| 欧美精品国产亚洲| 欧美一区二区国产精品久久精品| 别揉我奶头~嗯~啊~动态视频| 国产中年淑女户外野战色| 在现免费观看毛片| 偷拍熟女少妇极品色| 午夜激情福利司机影院| 97超级碰碰碰精品色视频在线观看| 白带黄色成豆腐渣| 成年女人永久免费观看视频| 国产久久久一区二区三区| 国产高清视频在线播放一区| 欧美zozozo另类| 联通29元200g的流量卡| avwww免费| 人妻少妇偷人精品九色| 精品福利观看| 亚洲最大成人av| 亚洲久久久久久中文字幕| 亚洲成人av在线免费| 欧美日韩在线观看h| av免费在线看不卡| 淫秽高清视频在线观看| 日本一二三区视频观看| 国产精品一二三区在线看| 黄色欧美视频在线观看| 亚洲人成网站在线播放欧美日韩| 国产精品一及| 久久久久国内视频| 精品久久久久久久人妻蜜臀av| 一区福利在线观看| 免费在线观看成人毛片| 午夜免费男女啪啪视频观看 | 午夜精品在线福利| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 真实男女啪啪啪动态图| 欧美日本亚洲视频在线播放| 国产av在哪里看| 丰满乱子伦码专区| 国产欧美日韩一区二区精品| 成人美女网站在线观看视频| 日本黄色视频三级网站网址| 日本 av在线| 免费搜索国产男女视频| 十八禁国产超污无遮挡网站| 亚洲国产精品久久男人天堂| 91午夜精品亚洲一区二区三区| 天天躁日日操中文字幕| 国产伦一二天堂av在线观看| 欧美精品国产亚洲| 国产精品永久免费网站| 人人妻人人澡欧美一区二区| 1024手机看黄色片| 欧美日韩综合久久久久久| 一级毛片久久久久久久久女| 国产在视频线在精品| 成人av一区二区三区在线看| 国产黄a三级三级三级人| 可以在线观看的亚洲视频| 99久久精品一区二区三区| 变态另类成人亚洲欧美熟女| 亚洲av一区综合| 五月玫瑰六月丁香| 亚洲美女视频黄频| 禁无遮挡网站| 日韩大尺度精品在线看网址| 欧美潮喷喷水| 国内少妇人妻偷人精品xxx网站| 免费在线观看影片大全网站| 久久久久久久久久成人| 成人特级黄色片久久久久久久| 国产三级在线视频| 99热只有精品国产| 亚洲欧美清纯卡通| 丝袜喷水一区| 日韩,欧美,国产一区二区三区 | 一夜夜www| 亚洲在线观看片| 3wmmmm亚洲av在线观看| 国产亚洲欧美98| 成人av一区二区三区在线看| 最近手机中文字幕大全| 99久久九九国产精品国产免费| 国产精品爽爽va在线观看网站| 观看免费一级毛片| 欧美又色又爽又黄视频| 精品一区二区三区视频在线| 波野结衣二区三区在线| 亚洲精品色激情综合| 韩国av在线不卡| 国产在视频线在精品| 18+在线观看网站| 国产精品福利在线免费观看| 男女之事视频高清在线观看| 蜜桃亚洲精品一区二区三区| 天天一区二区日本电影三级| 日韩精品有码人妻一区| 国产成人影院久久av| 国产成人a区在线观看| 91午夜精品亚洲一区二区三区| 观看美女的网站| 最近最新中文字幕大全电影3| 婷婷精品国产亚洲av在线| 国产片特级美女逼逼视频| av天堂中文字幕网| 国产精品久久久久久av不卡| 黄色一级大片看看| 男女之事视频高清在线观看| 国产成人一区二区在线| 中文字幕熟女人妻在线| 小蜜桃在线观看免费完整版高清| 丰满的人妻完整版| 最近在线观看免费完整版| 国产精品亚洲美女久久久| 又粗又爽又猛毛片免费看| 亚洲18禁久久av| 亚洲精品久久国产高清桃花| 色尼玛亚洲综合影院| 三级国产精品欧美在线观看| 国产精品女同一区二区软件| 久久精品夜色国产| 久久精品人妻少妇| 又爽又黄无遮挡网站| 婷婷亚洲欧美| 噜噜噜噜噜久久久久久91| 亚洲图色成人| 日韩av不卡免费在线播放| 18禁在线无遮挡免费观看视频 | 成人三级黄色视频| 精品人妻熟女av久视频| 成年版毛片免费区| 亚洲av一区综合| 免费人成在线观看视频色| 精品一区二区免费观看| aaaaa片日本免费| .国产精品久久| 久久久久久久午夜电影| 久久精品国产亚洲av天美| 久久欧美精品欧美久久欧美| 俺也久久电影网| 伊人久久精品亚洲午夜| 亚洲18禁久久av| 老司机午夜福利在线观看视频| 在线a可以看的网站| 丰满乱子伦码专区| 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看| 欧美激情国产日韩精品一区| 麻豆国产97在线/欧美| 日韩欧美国产在线观看| 一级a爱片免费观看的视频| 最近在线观看免费完整版| 成年版毛片免费区| 嫩草影院精品99| 亚洲天堂国产精品一区在线| 亚洲美女黄片视频| 免费观看的影片在线观看| 亚洲欧美日韩东京热| 2021天堂中文幕一二区在线观| 别揉我奶头~嗯~啊~动态视频| 寂寞人妻少妇视频99o| 禁无遮挡网站| 久久精品91蜜桃| 男人的好看免费观看在线视频| 国产高清不卡午夜福利| 一本久久中文字幕| 日本黄大片高清| 一本一本综合久久| 午夜福利高清视频| 99热全是精品| 日韩欧美精品v在线| 干丝袜人妻中文字幕| 精品久久久噜噜| 人人妻,人人澡人人爽秒播| 中文字幕精品亚洲无线码一区| 午夜福利高清视频| 韩国av在线不卡| av天堂在线播放| 一级黄片播放器| 欧美又色又爽又黄视频| 性欧美人与动物交配| 亚洲自拍偷在线| 99久国产av精品国产电影| 欧美在线一区亚洲| av福利片在线观看| 99久久精品国产国产毛片| 亚洲第一区二区三区不卡| 亚洲精品国产av成人精品 | 男人和女人高潮做爰伦理| 日本黄大片高清| 国产淫片久久久久久久久| 日韩欧美 国产精品| 久久久国产成人免费| 麻豆久久精品国产亚洲av| 简卡轻食公司| 免费看a级黄色片| 亚洲一区高清亚洲精品| 一a级毛片在线观看| 国产精品久久电影中文字幕| 国产精品一及| 国产精品嫩草影院av在线观看| 小说图片视频综合网站| 老司机影院成人| 国产爱豆传媒在线观看| 亚洲人与动物交配视频| 欧美+亚洲+日韩+国产| 免费在线观看成人毛片| 亚洲精品成人久久久久久| 午夜福利成人在线免费观看| 九九久久精品国产亚洲av麻豆| 尤物成人国产欧美一区二区三区| 亚洲国产欧洲综合997久久,| 亚洲人成网站在线观看播放| 男女之事视频高清在线观看| 婷婷精品国产亚洲av| 久久久久九九精品影院| ponron亚洲| 欧美xxxx性猛交bbbb| 亚洲精品在线观看二区| 久久草成人影院| 成人国产麻豆网| 精品久久久久久久久久久久久| 久久久久久久久久黄片| 国产av麻豆久久久久久久| 日韩强制内射视频| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩乱码在线| 欧美最新免费一区二区三区| 日韩亚洲欧美综合| 日本五十路高清| 中出人妻视频一区二区| 久久99热这里只有精品18| 中文在线观看免费www的网站| 波多野结衣高清作品| 国产欧美日韩精品一区二区| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| 成人无遮挡网站| 国产一级毛片七仙女欲春2| 国产精品一区www在线观看| 日日啪夜夜撸| 亚洲av电影不卡..在线观看| 国产精品免费一区二区三区在线| a级一级毛片免费在线观看| 美女xxoo啪啪120秒动态图| 黄色一级大片看看| 美女cb高潮喷水在线观看| 国产爱豆传媒在线观看| 国产精品伦人一区二区| 久久久久久伊人网av| 日韩大尺度精品在线看网址| 听说在线观看完整版免费高清| 熟妇人妻久久中文字幕3abv| videossex国产| 国产精品一区二区性色av| 精品无人区乱码1区二区| 亚洲自拍偷在线| 亚洲不卡免费看| 亚洲自偷自拍三级| 久久亚洲国产成人精品v| 69av精品久久久久久| 九九在线视频观看精品| 51国产日韩欧美| 美女免费视频网站| 最近最新中文字幕大全电影3| 国产成人福利小说| 日韩国内少妇激情av| 天天躁日日操中文字幕| 91av网一区二区| 欧美国产日韩亚洲一区| 国产成人福利小说| 日韩国内少妇激情av| 91久久精品国产一区二区三区| 欧美精品国产亚洲| 久久久久免费精品人妻一区二区| 美女内射精品一级片tv| 国产女主播在线喷水免费视频网站 | 中文亚洲av片在线观看爽| 99在线视频只有这里精品首页| 亚洲真实伦在线观看| 精品久久久久久成人av| 免费高清视频大片| 中国美白少妇内射xxxbb| 国产高清有码在线观看视频| 夜夜看夜夜爽夜夜摸| 人妻少妇偷人精品九色| 欧美日韩精品成人综合77777| 老师上课跳d突然被开到最大视频| 国产一区二区在线av高清观看| 精品少妇黑人巨大在线播放 | 男女边吃奶边做爰视频| 少妇熟女aⅴ在线视频| 性插视频无遮挡在线免费观看| 在线观看66精品国产| 久99久视频精品免费| 一级a爱片免费观看的视频| 老师上课跳d突然被开到最大视频| 亚洲无线在线观看| 日韩在线高清观看一区二区三区| av视频在线观看入口| 成人特级黄色片久久久久久久| 中文字幕精品亚洲无线码一区| 在线观看一区二区三区| 老司机午夜福利在线观看视频| 亚洲自拍偷在线| 成人二区视频| 欧美日韩在线观看h| 韩国av在线不卡| av专区在线播放| 一进一出好大好爽视频| 国内久久婷婷六月综合欲色啪| 国产成人91sexporn| 久久久久精品国产欧美久久久| 高清日韩中文字幕在线| 国产又黄又爽又无遮挡在线| 搡女人真爽免费视频火全软件 | 日韩在线高清观看一区二区三区| 淫秽高清视频在线观看| 联通29元200g的流量卡| 国产国拍精品亚洲av在线观看| 亚洲人成网站高清观看| 黑人高潮一二区| 欧美一级a爱片免费观看看| 亚洲熟妇中文字幕五十中出| 啦啦啦观看免费观看视频高清| 九色成人免费人妻av| 国产探花在线观看一区二区| 麻豆国产97在线/欧美| 日韩,欧美,国产一区二区三区 | 日韩成人av中文字幕在线观看 | 99热6这里只有精品| 久久精品国产亚洲av香蕉五月| 亚洲av中文字字幕乱码综合| av卡一久久| 91午夜精品亚洲一区二区三区| 精品久久久久久成人av| 久久精品综合一区二区三区| 日韩,欧美,国产一区二区三区 | 国产乱人视频| 一级av片app| 草草在线视频免费看| 十八禁国产超污无遮挡网站| 有码 亚洲区| 国产伦在线观看视频一区| 美女被艹到高潮喷水动态| 97在线视频观看| 国产激情偷乱视频一区二区| 长腿黑丝高跟| 99久久精品一区二区三区| 狂野欧美激情性xxxx在线观看| 日日摸夜夜添夜夜添av毛片| 18禁在线播放成人免费| 天天一区二区日本电影三级| 国产成人a∨麻豆精品| 在线天堂最新版资源| 亚洲欧美中文字幕日韩二区| 亚洲激情五月婷婷啪啪| a级毛色黄片| 精品久久久久久久久亚洲| 如何舔出高潮| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲av天美| 日韩欧美三级三区| 免费人成在线观看视频色| 天美传媒精品一区二区| av在线蜜桃| 人妻久久中文字幕网| 女生性感内裤真人,穿戴方法视频| 国产白丝娇喘喷水9色精品| 少妇人妻精品综合一区二区 | 精品乱码久久久久久99久播| 搞女人的毛片| 在线免费十八禁| 成人国产麻豆网| 18+在线观看网站| 一本久久中文字幕| 搡老岳熟女国产| 欧美潮喷喷水| 一区二区三区免费毛片| 午夜福利18| 联通29元200g的流量卡| 男女视频在线观看网站免费| av天堂在线播放| 国产成人a区在线观看| 久久6这里有精品| 97碰自拍视频| 丝袜喷水一区| 岛国在线免费视频观看| 又粗又爽又猛毛片免费看| 乱码一卡2卡4卡精品| 校园人妻丝袜中文字幕| 一进一出抽搐gif免费好疼| 国产精品久久电影中文字幕| 国产一区二区三区在线臀色熟女| 亚洲第一区二区三区不卡| 午夜福利成人在线免费观看| a级毛片a级免费在线| 欧美三级亚洲精品| 欧美另类亚洲清纯唯美| 国产精品三级大全| 成人亚洲欧美一区二区av| 女人被狂操c到高潮| 麻豆国产97在线/欧美| 小说图片视频综合网站| 成年女人毛片免费观看观看9| 男女边吃奶边做爰视频| 熟女电影av网| 午夜福利成人在线免费观看| 露出奶头的视频|