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

    基于級聯(lián)隨機森林與活動輪廓的3D MR圖像分割

    2019-06-11 06:42:58馬超劉亞淑駱功寧王寬全
    自動化學(xué)報 2019年5期
    關(guān)鍵詞:體素級聯(lián)輪廓

    馬超 劉亞淑 駱功寧 王寬全

    精確分割磁共振(Magnetic resonance,MR)圖像中的組織結(jié)構(gòu)并進(jìn)行后續(xù)的量化分析,能夠有效揭示疾病的病理生理學(xué)特征及其對特定解剖結(jié)構(gòu)的影響,在臨床診斷、治療方案規(guī)劃以及組織生物物理模型構(gòu)建等方面發(fā)揮著重要作用[1?2].然而,在多模態(tài)復(fù)雜條件磁共振圖像體數(shù)據(jù)中實現(xiàn)對組織結(jié)構(gòu)可重復(fù)的精確分割仍然是一項極具挑戰(zhàn)的任務(wù).首先,由于受射頻場不均勻性、不同軟組織之間的差異性和部分容積效應(yīng)等影響,磁共振圖像質(zhì)量退化嚴(yán)重,表現(xiàn)為圖像灰度的不均勻性和不同組織間信號的混疊[3],增加了磁共振圖像精細(xì)結(jié)構(gòu)分割的難度[4].另外,圖像的低對比度、軟組織的弱邊緣和復(fù)雜形狀結(jié)構(gòu),也給分割帶來了困難[5].

    近十年來,大量磁共振圖像自動化分割方法被提出,其中,基于隨機森林(Random forest,RF)的方法和基于活動輪廓模型(Actives contour model,ACM)的方法是較為成功的兩類方法[6?7].

    隨機森林適合于處理大量的具有高維特征的多類別數(shù)據(jù),在許多醫(yī)學(xué)圖像分割任務(wù)中被證明具有較高的精確性和魯棒性[8?9],近年獲得了持續(xù)增加的關(guān)注[4,10?11].例如,Schneider等提出一種血管三維分割與中心線提取框架[12],該框架在不同尺度與方向上提取圖像局部特征,并分別利用傾斜隨機森林與多元霍夫投票實現(xiàn)血管的三維分割與中心線估計.Nicholas等利用多模態(tài)灰度、幾何以及不對稱特征驅(qū)動兩級隨機森林,實現(xiàn)了完整腦部以及腦腫瘤的分割[10].Chao等提出一種兩階段分割模型[13],該模型利用泛化的霍夫變換與三維活動形狀模型進(jìn)行腎皮質(zhì)定位,進(jìn)而利用改進(jìn)的隨機森林方法,實現(xiàn)了腹部三維CT圖像中腎臟多個組織結(jié)構(gòu)的全自動分割.近年由于受到深度學(xué)習(xí)在自然場景圖像分割中良好效果的啟發(fā)[14],一些研究者提出了隨機森林的級聯(lián)架構(gòu)[6,15],以實現(xiàn)有限的醫(yī)學(xué)圖像訓(xùn)練樣本上具有良好性能分割模型的學(xué)習(xí).盡管基于隨機森林的方法能夠獲得良好的分割效果,但由于訓(xùn)練數(shù)據(jù)的數(shù)量和質(zhì)量所引起的過擬合問題依然存在[11],并且典型的隨機森林的分割結(jié)果并不具有幾何約束[15].

    活動輪廓模型一般利用圖像的邊緣或區(qū)域信息驅(qū)動活動輪廓向目標(biāo)邊緣演化并收斂[16],在醫(yī)學(xué)圖像分割領(lǐng)域具有廣泛的應(yīng)用.例如,Avendi等[17]與Ngo等[18]將深度學(xué)習(xí)與活動輪廓模型整合,實現(xiàn)短軸心臟磁共振圖像中左心室的自動分割.Hoogi等[19]提出了一種自適應(yīng)局部窗方法,提升了活動輪廓模型在CT與磁共振圖像上對肝損傷的分割效果.Wang等[20]利用基于配準(zhǔn)的測地線活動輪廓模型,實現(xiàn)了腹主動脈瘤內(nèi)腔與外壁的分割.孫文燕等[21]將模糊聚類與活動輪廓模型結(jié)合,提高了局部分割方法對初始輪廓的魯棒性.然而,大多數(shù)典型的基于活動輪廓模型的方法都不是完全自動化的分割方法,其分割流程與分割結(jié)果都依賴特定的輪廓初始化,并且對于具有灰度非均勻和低組織對比度等復(fù)雜條件的圖像,較難獲得滿意的分割結(jié)果[22].

    為解決上述問題,本文在相關(guān)工作[5?6,15]基礎(chǔ)上,提出一種新的整合了級聯(lián)隨機森林(Concatenated random forests,cRFs)與活動輪廓模型的醫(yī)學(xué)磁共振圖像體數(shù)據(jù)三維分割方法.該方法從多模態(tài)磁共振圖像體數(shù)據(jù)中提取多尺度局部魯棒統(tǒng)計特征,并利用具有級聯(lián)架構(gòu)的隨機森林迭代的對訓(xùn)練樣本進(jìn)行學(xué)習(xí)得到組織分類器序列;在迭代過程中,每個層級產(chǎn)生的組織概率圖作為擴展的圖像源被輸入到下一層級的分類器訓(xùn)練中.利用訓(xùn)練好的級聯(lián)分類器可實現(xiàn)對測試樣本的初步分割,該分割結(jié)果進(jìn)一步作為初始輪廓和形狀先驗被整合進(jìn)一個尺度可調(diào)的活動輪廓模型,以完成最終的分割.與以往基于隨機森林的分割方法相比,本文方法能夠有效地在多尺度、多層級上整合多模態(tài)圖像信息,并且將隨機森林的體素分類任務(wù)轉(zhuǎn)化為輪廓曲線演化,對于組織結(jié)構(gòu)的精確、光滑的分割具有重要意義.此外,與以往基于活動輪廓模型的方法相比,本文方法可實現(xiàn)自動化的分割,并且對于復(fù)雜條件磁共振圖像具有更好的魯棒性.

    與本文相似的工作是文獻(xiàn)[6,15?16].本文方法與文獻(xiàn)[16]中方法均使用了局部高斯分布構(gòu)建了活動輪廓模型的能量泛函,主要區(qū)別在于本文方法是在多模式成像條件下構(gòu)建能量泛函并整合了形狀約束項,這對于在單一成像模式下呈現(xiàn)相似灰度的組織結(jié)構(gòu)的分割具有重要意義;同時本文方法是一種自動化的分割方法,而文獻(xiàn)[16]方法需要手工設(shè)定初始輪廓.本文方法與文獻(xiàn)[6]中方法以及文獻(xiàn)[15]中的LINKS方法均使用隨機森林實現(xiàn)了在深度學(xué)習(xí)中常見的級聯(lián)架構(gòu),以迭代對分割結(jié)果進(jìn)行精煉.而本文方法與LINKS方法使用了不同的特征提取方式;為避免過擬合問題,LINKS方法利用了隨機選擇樣本與特征的常規(guī)的隨機性策略,而本文進(jìn)一步提出了一種級聯(lián)隨機森林訓(xùn)練的Dropout策略,以提升模型在小規(guī)模且類別分布非均衡訓(xùn)練樣本(如腦腫瘤數(shù)據(jù),其98% 為健康組織體素)上的泛化能力;本文將體素分類與活動輪廓演化進(jìn)行了整合.相比于文獻(xiàn)[6]中單一成像模式下的兩相分割方法,本文方法為多成像模式下的多相分割方法,并且有著不同的輸入數(shù)據(jù)處理方式、訓(xùn)練策略以及先驗信息整合方式.

    本文剩余部分內(nèi)容安排如下:第1節(jié)詳細(xì)介紹了本文提出的基于級聯(lián)隨機森林與活動輪廓模型的磁共振圖像三維分割方法;第2節(jié)說明了算法實現(xiàn)過程和實驗結(jié)果;第3節(jié)對實驗結(jié)果進(jìn)行了討論;第4節(jié)對本文工作進(jìn)行了總結(jié).

    1 本文方法

    圖1為以腦腫瘤磁共振圖像的三維分割為例說明本文方法的流程,這里將磁共振體數(shù)據(jù)的分割問題轉(zhuǎn)化為組織體素分類和組織輪廓曲線演化的整合問題.具體的,當(dāng)有磁共振體數(shù)據(jù)輸入時,分割算法分三個階段來實施:

    1)從多模態(tài)磁共振體數(shù)據(jù)中提取隨機的魯棒統(tǒng)計特征以獲取圖像局部與環(huán)境信息,這些圖像特征作為輸入項用以驅(qū)動后續(xù)分割階段.

    2)利用級聯(lián)隨機森林對組織結(jié)構(gòu)進(jìn)行初步分割,即在級聯(lián)框架內(nèi)將隨機森林用作級聯(lián)分類器,以體素分類的方式產(chǎn)生組織結(jié)構(gòu)的組織概率圖序列.將磁共振體數(shù)據(jù)內(nèi)的每一個體素標(biāo)識為具有最大概率的組織類別,由此獲得組織結(jié)構(gòu)的初步分割.

    3)為精煉初步分割結(jié)果,體素分類被進(jìn)一步整合進(jìn)輪廓曲線演化方案中.通過驅(qū)動活動輪廓演化,并在理想的組織結(jié)構(gòu)輪廓處收斂,從而獲得最終的分割結(jié)果.

    圖1 分割框架流程圖Fig.1 Flowchart of the proposed segmentation framework

    分別對體素分類和輪廓曲線演化階段進(jìn)行離線訓(xùn)練后,該分割框架即可部署用于自動的組織結(jié)構(gòu)分割任務(wù).下面詳細(xì)說明分割框架的三個階段.

    1.1 隨機魯棒統(tǒng)計特征

    設(shè)I:?→R 為待分割圖像體數(shù)據(jù),其中,?∈R3為圖像域.設(shè)B(x)??為以體素x為中心的鄰域區(qū)域,Rcentric∈B(x)為中心體素x的尺度可調(diào)的局部區(qū)域,Rrandom∈B(x)為x鄰域內(nèi)隨機位置與尺度的局部區(qū)域.

    在對中心體素x的鄰域區(qū)域B(x)進(jìn)行正則化處理后[23],利用高斯核函數(shù)Kσ(u)控制局部區(qū)域尺度和區(qū)域內(nèi)權(quán)重,可計算關(guān)于x的局部魯棒統(tǒng)計特征量如下:

    其中,Ψ(x)為局部加權(quán)灰度均值;Q(x)為局部四分位差,即局部區(qū)域體素灰度值的上四分位數(shù)(Q3)與下四分位數(shù)(Q1)的差值;γ(x)為局部加權(quán)灰度標(biāo)準(zhǔn)差.

    由此,定義局部與環(huán)境魯棒統(tǒng)計特征向量分別為:

    則完整隨機魯棒統(tǒng)計特征向量為:其中,參數(shù)b∈{0,1}表明特征向量輸入體素分類階段還是輪廓演化階段,分別如圖2(a)和(b)所示.理論上,對于體數(shù)據(jù)中的某一體素,通過改變局部區(qū)域的位置和尺度,可以提取任意數(shù)量的隨機特征.在當(dāng)前工作中,局部區(qū)域最大尺度限定為5×5×5,搜索鄰域限定為31×31×31.

    圖2 隨機魯棒統(tǒng)計特征學(xué)習(xí)方案二維示意圖Fig.2 A two-dimensional illustration of the random robust statistics features learning scheme

    1.2 級聯(lián)架構(gòu)下的體素分類

    本文利用隨機森林在級聯(lián)架構(gòu)下確定體素x∈?所屬的類別c∈C,以實現(xiàn)對組織結(jié)構(gòu)的分割.圖3是所提出級聯(lián)體素分類框架的流程圖.以體數(shù)據(jù)的隨機魯棒統(tǒng)計特征f(x,I)為輸入,級聯(lián)分類器將在每一層級輸出組織概率圖FM,由某一層級輸出的組織概率圖提取的隨機魯棒統(tǒng)計特征f(x,FM)將作為增廣特征與原始輸入特征一同輸入到下一層級,由此實現(xiàn)圖像信息的迭代的融合和精煉.級聯(lián)架構(gòu)下的信息迭代融合,通過向后續(xù)迭代分類器引入具有空間約束和較好初始化的特征,提升了體素分類的精度.

    在級聯(lián)架構(gòu)中分類器的選擇可以很靈活,本文選擇隨機森林作為分類器,因其能夠高效地處理從體數(shù)據(jù)中提取的大量、高維度的圖像特征.隨機森林由眾多分類/回歸樹組成,作為一項有監(jiān)督學(xué)習(xí)技術(shù),通??煞钟?xùn)練和測試兩個階段對其進(jìn)行操作,下一節(jié)將詳細(xì)說明面向組織結(jié)構(gòu)分割的隨機森林的訓(xùn)練與標(biāo)識過程.

    圖3 級聯(lián)架構(gòu)下的體素分類流程圖Fig.3 Overview of the voxel-wise classi fication within the proposed concatenated scheme

    1.3 隨機森林訓(xùn)練與標(biāo)識

    在級聯(lián)隨機森林第一個層級的訓(xùn)練中,每棵決策樹t將只利用由原始圖像I獲取的隨機魯棒統(tǒng)計特征f(x,I)作為輸入,通過學(xué)習(xí)得到一個弱分類器pt(c|f(x,I)).學(xué)習(xí)過程為利用訓(xùn)練樣本的高維特征迭代的在每一個樹節(jié)點對訓(xùn)練體素進(jìn)行分類.為避免過擬合并提升模型的泛化能力,除在訓(xùn)練中實施隨機森林常規(guī)的隨機性策略外[15],本文還引入一種原本用于卷積神經(jīng)網(wǎng)絡(luò)的Dropout策略[24].具體為,當(dāng)對該層級隨機森林訓(xùn)練完成后,將某一訓(xùn)練圖像的全部體數(shù)據(jù)輸入該森林用以輸出下一個層級所需要的組織概率圖,此時,按照某一概率(例如:0.5)獨立地從森林中移除一些樹.由此,削弱了相鄰層級分類器之間的依賴,增強了每一層級的獨立分類能力.

    在級聯(lián)隨機森林后續(xù)層級的訓(xùn)練中,從原始圖像I以及上一層級組織概率圖FM中獲取的隨機魯棒統(tǒng)計特征f(x,I,FM)被用來訓(xùn)練分類器pt(c|f(x,I,FM)),其訓(xùn)練過程與第一個層級一致.

    在標(biāo)識階段,利用訓(xùn)練好的級聯(lián)隨機森林對目標(biāo)圖像I的每一體素x進(jìn)行標(biāo)識.具體為,將體素x的高維特征輸入級聯(lián)隨機森林,在級聯(lián)架構(gòu)的某一層級每一棵樹將產(chǎn)生體素x屬于類別c的一個估計pt(c|f(x,I,FM)).則體素x在該層級的類別估計可計算為層級內(nèi)全部nt棵樹所做估計的平均,即

    體素x最終標(biāo)識為在最終層級具有最大概率值的類別估計,即pcRFs(c|f(x,I,FM))=maxcp(c|f(x,I,FM)),由此實現(xiàn)對目標(biāo)圖像I的初步分割.

    1.4 體素分類與輪廓演化的整合

    級聯(lián)隨機森林是對每一個體素獨立的進(jìn)行分類,并未考慮體素之間的關(guān)聯(lián)性,因此其分割結(jié)果不是具有幾何約束的,如分割結(jié)果中存在的散點和空洞.而活動輪廓模型可以通過衡量局部區(qū)域體素間的相似性構(gòu)建能量泛函,并通過最小化能量泛函驅(qū)動活動輪廓演化,從而獲得閉合的目標(biāo)輪廓.但是,標(biāo)準(zhǔn)的活動輪廓模型對于存在灰度非均勻和信號混疊等復(fù)雜條件的磁共振圖像分割效果不夠理想,并且由于需要手工設(shè)定初始輪廓,導(dǎo)致其無法實現(xiàn)自動化的目標(biāo)分割.為此,本文將級聯(lián)隨機森林的體素分類與活動輪廓模型驅(qū)動的活動輪廓演化進(jìn)行整合.通過活動輪廓演化對體素分類結(jié)果進(jìn)行精煉;另一方面,級聯(lián)隨機森林的體素分類結(jié)果則為活動輪廓模型提供了初始輪廓和形狀先驗,使其能夠自動地分割具有復(fù)雜成像條件的磁共振圖像.具體整合過程如下文所述.

    設(shè)圖像體數(shù)據(jù)為?,輪廓曲線/曲面C將體數(shù)據(jù)分割成N個互不相交的子區(qū)域.由此,每一體素x可以通過整合x鄰域內(nèi)多模態(tài)概率分布函數(shù)(Probability distribution function,PDF)[5]與級聯(lián)隨機森林輸出的空間先驗信息進(jìn)行描述.則體素x屬于第i個組織類別可描述為:其中,M為圖像模態(tài)的數(shù)量,高斯核函數(shù)Kη通過尺度參數(shù)η控制體素x的鄰域的范圍和權(quán)重,概率密度函數(shù)p中的參數(shù)μj,i對第j個圖像模式下局部區(qū)域?i的圖像特征進(jìn)行估計,fj(y)為從第j個圖像模式獲取的體素y的局部魯棒統(tǒng)計特征向量.值得注意的是,文獻(xiàn)[6]中面向兩相分割的先驗信息整合方式此處并不適用,因為多相目標(biāo)的輪廓描述與距離度量將導(dǎo)致后續(xù)能量泛函的復(fù)雜的水平集表示形式.

    由此,對于給定的體素x∈?,其擬合能量可定義為:

    最終的能量泛函定義為:

    基于Heaviside函數(shù)H(·)可構(gòu)造多相水平集函數(shù)以表示圖像的多個分區(qū).設(shè)為水平集函數(shù)集合,水平集函數(shù)數(shù)量L由待分割目標(biāo)類別數(shù)量決定,利用水平集函數(shù)的組合M(Φ)即可實現(xiàn)對圖像不同分區(qū)的描述[25].則式(7)中的能量泛函可由水平集函數(shù)集合Φ描述為:

    2 實驗結(jié)果

    圖4 級聯(lián)隨機森林對多個目標(biāo)圖像在不同層級做出的組織概率圖的估計Fig.4 The tissue probability maps estimated from different levels of the concatenated random forests for several target subjects

    我們利用公開的磁共振圖像體數(shù)據(jù)集[4,27?28]對所提出模型進(jìn)行了驗證,其中左心房數(shù)據(jù)[28]為1.5T場強掃描儀獲取的bSSFP模式的磁共振圖像體數(shù)據(jù),分辨率為1.25mm×1.25mm×2.7mm;腦腫瘤數(shù)據(jù)[4]為1.5T與3T場強掃描儀獲取的 T1、T1c、T2和FLAIR 模式的磁共振圖像體數(shù)據(jù),分辨率為1mm×1mm×1mm;尾狀核數(shù)據(jù)[27]為1.5T場強掃描儀獲取的IRSPGR模式的磁共振圖像體數(shù)據(jù),分辨率為0.9375mm×0.9375mm×1.5mm.這里僅進(jìn)行了最低限度的數(shù)據(jù)預(yù)處理,即將每種成像模式下每個體數(shù)據(jù)內(nèi)的體素的灰度值規(guī)范化為0~255.實驗環(huán)境為Inte(R)Core(TM)3.4GHz CPU,48GB RAM,Visual studio community 2015.除個別圖像外,實驗中級聯(lián)隨機森林包含5個層級,每個層級訓(xùn)練60個分類決策樹,最大樹深度為50,葉節(jié)點最小允許樣本數(shù)為8.決策樹訓(xùn)練過程中,每個節(jié)點考慮10000個隨機采樣的特征及其20個隨機分布的閾值以確定該節(jié)點分裂函數(shù).在目標(biāo)輪廓精煉階段,活動輪廓模型按文獻(xiàn)[5,29]所述方法進(jìn)行設(shè)置以獲得最優(yōu)分割結(jié)果.我們將從定性和定量結(jié)果兩方面呈現(xiàn)本文方法的分割效果,并在第3節(jié)對隨機森林的級聯(lián)架構(gòu)、參數(shù)設(shè)置以及模型的不同組成部分對分割結(jié)果的影響等方面進(jìn)行討論,并將本文方法與相關(guān)方法在復(fù)雜條件圖像上的分割精度和魯棒性等方面進(jìn)行比較分析.

    2.1 定性結(jié)果

    圖4是具有級聯(lián)架構(gòu)的隨機森林在不同層級上對多個測試樣本的分割結(jié)果.圖4中第1~3列分別是尾狀核、左心房和腦腫瘤的磁共振圖像體數(shù)據(jù).圖4(a)和(b)分別為原始目標(biāo)圖像體數(shù)據(jù)的一個切片圖像及其分割金標(biāo)準(zhǔn).圖4(c)~(f)分別對應(yīng)級聯(lián)隨機森林在第1、2、3、5層級輸出的三維組織概率圖的一個切片圖像.為表述清晰,我們只用綠色突出強調(diào)了組織概率圖中概率值大于0.6的體素.

    為更好理解所提出的分割框架中不同階段對分割結(jié)果的影響,圖5(a)~(c)分別呈現(xiàn)了分割框架中沒有輪廓精煉的級聯(lián)隨機森林(第二階段)、沒有形狀約束的活動輪廓模型以及整合了級聯(lián)隨機森林和活動輪廓模型(最終階段)的對比分割結(jié)果,圖5(d)為分割金標(biāo)準(zhǔn).圖5中第1~4行分別對應(yīng)腦腫瘤完整結(jié)構(gòu)、尾狀核、左心房和腦腫瘤子結(jié)構(gòu)—水腫的三維分割結(jié)果.

    圖6進(jìn)一步清晰展示了所提出方法對復(fù)雜條件圖像進(jìn)行三維分割的精確性和魯棒性.圖6(a)~(c)分別是左心房、腦腫瘤子結(jié)構(gòu)—水腫和腦腫瘤子結(jié)構(gòu)—腫瘤核的分割結(jié)果在心臟磁共振圖像軸向、腦磁共振圖像冠向和軸向上的多個切片圖像.輪廓線的不同顏色—藍(lán)色、紅色和綠色分別對應(yīng)所提出算法的分割結(jié)果、分割金標(biāo)準(zhǔn)和算法分割結(jié)果與金標(biāo)準(zhǔn)相符合部分.

    2.2 定量結(jié)果

    圖7是在左心房與腦腫瘤體數(shù)據(jù)集上,利用DC(Dice coefficient)系數(shù)作為評價指標(biāo),對不同參數(shù)影響下級聯(lián)隨機森林分割效果的定量評價.圖中DC系數(shù)值為各自體數(shù)據(jù)集訓(xùn)練集上的十折交叉驗證結(jié)果的平均值.此實驗僅考慮影響的趨勢,因此當(dāng)調(diào)整某一個參數(shù)時,其他參數(shù)被設(shè)置為某一固定值而不是為獲得最優(yōu)分割效果進(jìn)行取值.圖7(a)~(d)分別為每層級決策樹數(shù)量、最大樹深度、每個葉節(jié)點最小樣本數(shù)以及級聯(lián)架構(gòu)層級數(shù)目對分割結(jié)果的影響.

    圖5 分割框架不同階段的分割結(jié)果對比Fig.5 Comparison of different components in the proposed segmentation framework

    表1是利用平均DC值在具有復(fù)雜成像結(jié)果的左心房數(shù)據(jù)集上,對不同分割方法的兩相分割結(jié)果進(jìn)行的對比定量評價.對比方法包括標(biāo)準(zhǔn)的基于活動輪廓(ACM)[26,30]和隨機森林(RF)[9,31]的方法,競爭輪廓方法(Competitive contours)[32],多尺度塊多模式圖譜方法(M3AS)[33],以及文獻(xiàn)[6]中提出的基于級聯(lián)隨機森林(Combined random forest,CRF)和整合隨機森林與活動輪廓(CRF-AC)的兩種分割方法.此外,Proposed1和Proposed2分別表示所提出分割框架中未包含輪廓精煉的第二階段方法和包含輪廓精煉的最終階段方法.表1中結(jié)果除M3AS方法使用了獨立獲取的數(shù)據(jù)集外,其余方法皆為各自文獻(xiàn)在相同的公開數(shù)據(jù)集上報告的結(jié)果.另外,在Surface-to-Surface這一指標(biāo)上Proposed2的值為0.64±0.39mm,是目前在該左心房數(shù)據(jù)上獲得的最優(yōu)分割結(jié)果.

    圖6 低質(zhì)量磁共振圖像體數(shù)據(jù)三維分割結(jié)果在不同視角下的多個切片圖像Fig.6 Multiple slices of the 3D segmentation results for low quality volumetric MR images in different views

    圖7 級聯(lián)隨機森林不同參數(shù)對分割結(jié)果的影響Fig.7 Impact of different parameters in the concatenated random forests

    為進(jìn)一步驗證所提出方法在多模式磁共振體數(shù)據(jù)上的多相分割效果,表2顯示了不同分割方法在BRATS15腦腫瘤體數(shù)據(jù)測試集上的在線評價結(jié)果,量化評價指標(biāo)為DC、PPV(Positive predictive value)和Sensitivity.對比方法為4種基于深度學(xué)習(xí)的方法,分別是FCNNs-CRFs[34]、INPUTCASCADECNN[35]、CNNS[36]和DeepMedic[3].表2中INPUTCASCADECNN,利用了BRATS13測試集,CNNs利用了BRATS15初始的包含53個數(shù)據(jù)的測試集,這兩個數(shù)據(jù)集都包含在當(dāng)前的BRATS15測試集中.FCNNs-CRFs、DeepMedic以及本文提出的Proposed2方法皆在包含110個數(shù)據(jù)的BRATS15測試集上進(jìn)行了評測.

    表1 不同分割方法在左心房數(shù)據(jù)集上的DC系數(shù)對比Table 1 Dice coefficients(DC)of different methods on the left atrial dataset

    表2 不同方法在BRATS15測試集上的對比Table 2 Comparison of different methods on the BRATS15 test set

    3 討論

    3.1 隨機森林級聯(lián)架構(gòu)對分割結(jié)果的影響

    如圖4所示,在級聯(lián)架構(gòu)下,估計出的組織概率圖隨著層級的增加其精度和清晰度也隨之增加.產(chǎn)生這一現(xiàn)象的原因為在初始層級只有圖像原始特征用于估計組織概率圖,導(dǎo)致許多虛假目標(biāo)邊緣,并且對真實目標(biāo)的體素標(biāo)注也不完全(圖4(c)).而在后續(xù)的層級中,級聯(lián)架構(gòu)迭代的將前一層級輸出的組織概率圖與圖像原始特征進(jìn)行整合,以此實現(xiàn)對目標(biāo)組織結(jié)構(gòu)的不斷精煉(圖4(d)~(f)).在每一層級,每一棵獨立的決策樹以組織概率圖的形式產(chǎn)生對目標(biāo)的一個弱分類結(jié)果,利用式(4)對層級內(nèi)全部決策樹的分類結(jié)果進(jìn)行融合,則一些決策樹產(chǎn)生的錯誤邊緣(即組織概率圖中具有錯誤高概率值估計的體素)將由其他決策樹產(chǎn)生的組織概率圖加以修正.隨著層級的增加,輸入特征不斷得到精煉,決策樹分類能力得到持續(xù)提升,可以看到,第5層級輸出的對目標(biāo)組織結(jié)構(gòu)邊緣的估計(圖4(f)),已經(jīng)較好地接近分割金標(biāo)準(zhǔn)(圖4(b)).

    3.2 級聯(lián)隨機森林與活動輪廓模型對分割結(jié)果的影響

    如圖5所示,對于完整腦腫瘤以及尾狀核兩種組織結(jié)構(gòu),未進(jìn)行輪廓精煉的級聯(lián)隨機森林(圖5(a))、未整合形狀先驗的活動輪廓模型(圖5(b))以及整合了級聯(lián)隨機森林與活動輪廓的分割方法(圖5(c)),參照于分割金標(biāo)準(zhǔn)(圖5(d)),都取得了較好的分割效果.盡管由于級聯(lián)隨機森林的體素分類方案是對每一體素獨立的進(jìn)行標(biāo)識,并且沒有進(jìn)行進(jìn)一步的活動輪廓的精煉,導(dǎo)致其對尾狀核的分割結(jié)果存在較多散點和空洞,但目標(biāo)輪廓仍然較好的符合分割金標(biāo)準(zhǔn),并且這些散點和空洞可以在分割結(jié)果精煉階段,通過活動輪廓的演化進(jìn)行填充.

    然而,對于左心房和腦腫瘤子結(jié)構(gòu)—水腫來說,盡管活動輪廓模型的初始種子點以及活動輪廓曲線演化次數(shù)已經(jīng)被精細(xì)地控制,但由于目標(biāo)組織結(jié)構(gòu)與周圍組織結(jié)構(gòu)的灰度存在混疊,導(dǎo)致其分割結(jié)果泄漏問題較為嚴(yán)重(圖5(b)).相比于活動輪廓模型,級聯(lián)隨機森林得益于級聯(lián)的學(xué)習(xí)與測試方案,分割結(jié)果的泄漏問題得到較好的控制,對這兩種組織結(jié)構(gòu)的分割效果較好,但分割結(jié)果中組織結(jié)構(gòu)內(nèi)部的空洞問題進(jìn)一步劣化(圖5(a)).

    最后,整合了級聯(lián)隨機森林與活動輪廓的分割方法對于全部的組織結(jié)構(gòu)均取得了更為精確和具有幾何約束的分割結(jié)果(圖5(c)).該分割框架中的形狀先驗約束項(式(5))有效控制了活動輪廓向相鄰組織結(jié)構(gòu)的泄漏,同時,活動輪廓的演化填充了級聯(lián)隨機森林分割結(jié)果的空洞,并對細(xì)節(jié)進(jìn)一步做了精煉.

    3.3 對復(fù)雜條件圖像的分割結(jié)果

    如圖6所示,盡管由于較低的空間分辨率和較強的灰度不均勻,分割這些低成像質(zhì)量的磁共振圖像十分困難,本文方法的分割結(jié)果仍然與分割金標(biāo)準(zhǔn)契合良好,驗證了本文方法對于復(fù)雜條件磁共振體數(shù)據(jù)的分割魯棒性.

    3.4 模型參數(shù)對分割結(jié)果的影響

    圖7呈現(xiàn)了不同參數(shù)對分割精度的影響.如圖7(a)所示,從使用1棵決策樹到使用2棵決策樹有著明顯的精度提升,DC值從0.67±0.07提升到0.74±0.06.隨著決策樹數(shù)量的增加精度持續(xù)提升,直至趨于穩(wěn)定.對這一現(xiàn)象的一個合理解釋是多棵決策樹的融合修正了部分決策樹的錯誤分類,而訓(xùn)練樣本的數(shù)量決定了隨機森林所能夠達(dá)到的最大精度,以及達(dá)到這一精度時所需的決策樹數(shù)量.如圖7(b)所示,分割精度從樹深度5~30逐漸提升.為使樹深度能夠適應(yīng)訓(xùn)練樣本規(guī)模以及樹中每個葉節(jié)點的最小樣本數(shù)量,本文中將最大樹深度設(shè)定為50.減小葉節(jié)點最小樣本數(shù)量(圖7(c))能夠提升分割精度,而該參數(shù)過小的取值容易引起模型過擬合.如圖7(d)所示,由于迭代的利用組織概率圖進(jìn)行體素分類,隨著級聯(lián)架構(gòu)層級數(shù)量的增加,分割精度持續(xù)提升并在一定層級數(shù)量后趨于穩(wěn)定,尤其是從第1層級到第2層級,此時組織概率圖被初次使用,作為增廣圖像輸入到第2層級,導(dǎo)致分割精度提升明顯.以上實驗結(jié)果進(jìn)一步證明了級聯(lián)分類架構(gòu)對于實現(xiàn)目標(biāo)分割的有效性.

    3.5 分割效果對比分析

    如表1所示,即使沒有進(jìn)行最終的輪廓精煉,所提出的方法(Proposed1)仍然取得了較高的分割精度,在左心房體數(shù)據(jù)訓(xùn)練集與測試集上的分割結(jié)果的DC值分別為0.78±0.06與0.74±0.07.Proposed1方法較低的DC值主要由所分割出目標(biāo)內(nèi)部的空洞引起,如圖4和圖5所示,即便如此,分割出的目標(biāo)外邊緣依然保持了與分割金標(biāo)準(zhǔn)較高的契合度.由此,Proposed1方法的分割結(jié)果能夠為后續(xù)的曲線演化提供較好的初始輪廓,Proposed2方法因此獲得了良好的分割精度,訓(xùn)練集與測試集上的DC值分別達(dá)到0.94±0.04與0.90±0.03.而ACM和RF等標(biāo)準(zhǔn)的分割方法未能利用圖像體數(shù)據(jù)的豐富環(huán)境信息引導(dǎo)分割過程,以獲得更好的分割精度.Proposed1為基于學(xué)習(xí)的分割方法,其在訓(xùn)練集和測試集上取得了相近的分割結(jié)果,進(jìn)一步證明了本文所提出的級聯(lián)隨機森林訓(xùn)練策略在抑制過擬合方面的優(yōu)勢.

    在對左心房的分割中,競爭輪廓方法取得了較高的分割精度,其DC值達(dá)到了0.92,考慮到該方法為半自動的分割方法,而本文方法為全自動的分割方法,因此本文方法仍然具有一定優(yōu)勢.M3AS方法利用了獨立獲取的數(shù)據(jù)集報告了分割結(jié)果,雖然在不同數(shù)據(jù)集上對不同方法進(jìn)行對比不完全可靠,但仍然可從數(shù)據(jù)對比中體現(xiàn)本文方法的良好精度.Proposed1與CRF方法獲得了相似的DC值,驗證了本文所提出的隨機森林級聯(lián)框架的有效性,進(jìn)一步來說,Proposed1方法有著更加精簡的級聯(lián)架構(gòu),并且不需要對訓(xùn)練數(shù)據(jù)集進(jìn)行手工擴展,因此有著更大的應(yīng)用潛力.在所有對比方法中,CRF-AC方法獲得了最高的0.93±0.05的DC值,然而該方法為單一成像模式下的兩相分割模型,而本文所提出的Proposed2能夠綜合利用多種成像模式進(jìn)行多相分割,實驗中利用的單一成像模式的左心房標(biāo)準(zhǔn)數(shù)據(jù)集并未充分體現(xiàn)本文方法的通用性.因此將在后續(xù)通過表2的實驗結(jié)果對本文方法的性能做進(jìn)一步分析.目前利用深度學(xué)習(xí)對左心房結(jié)構(gòu)進(jìn)行分割的研究還很少見,例如與本文方法進(jìn)行對比的競爭輪廓方法為一種基于形變模型的方法,而M3AS為基于多圖譜融合的方法.這主要是由于深度學(xué)習(xí)方法對數(shù)據(jù)量存在嚴(yán)重的依賴,而當(dāng)前公開的左心房數(shù)據(jù)集的數(shù)據(jù)量還比較有限,本文所用標(biāo)準(zhǔn)數(shù)據(jù)集以及M3AS方法所用的獨立數(shù)據(jù)集,其訓(xùn)練集均只包含10個磁共振體數(shù)據(jù).

    在表2所示的腦腫瘤分割結(jié)果中,由于INPUTCASCADECNN與CNNS使用了較小規(guī)模的測試集,直接的量化對比不十分可靠,但仍然可以看出Proposed2所具有的競爭力,在Complete和Core區(qū)域獲得了接近INPUTCASCADECNN和優(yōu)于CNNS方法的DC值.在Enhancing區(qū)域,INPUTCASCADECNN與CNNS方法的DC值均明顯高于另外三種方法,這主要是由于Enhancing區(qū)域在腦腫瘤中往往成彌散狀,即使手工分割仍然存在較為顯著的不一致性,而擴充后的完整的BRATS15測試集的分割金標(biāo)準(zhǔn)是融合了手工分割與算法分割的結(jié)果,因此分割金標(biāo)準(zhǔn)有可能會受到算法分割結(jié)果的系統(tǒng)性影響,這也是BRATS15數(shù)據(jù)集存在的一個局限.對于使用BRATS15測試集中全部110個數(shù)據(jù)進(jìn)行評測的Proposed2、FCNNs-CRFs和DeepMedic方法,Proposed2的優(yōu)勢較為明顯.Proposed2在三個評測區(qū)域均取得了最高的DC值,其PPV和Sensitivity值則與另外兩種方法較為接近,但是Proposed2方法在取得較高的PPV與Sensitivity值的同時,在同一評測區(qū)域該方法的PPV和Sensitivity取值更為均衡,驗證了本文方法在進(jìn)行多相分割時,在控制過分割和欠分割方面的良好的綜合性能.

    在對腦腫瘤的分割中,與本文方法進(jìn)行對比的是4種深度學(xué)習(xí)方法.雖然深度學(xué)習(xí)方法在自然場景圖像分割中體現(xiàn)了良好性能,但對于心臟和腦部磁共振圖像來說,不同患者間磁共振成像結(jié)果的多樣性遠(yuǎn)低于CIFAR和ImageNet等數(shù)據(jù)集中的自然場景圖像[35],因此深度學(xué)習(xí)模型龐大的模型規(guī)模所具有的優(yōu)勢此時并不明顯.此外標(biāo)準(zhǔn)數(shù)據(jù)集中較少的訓(xùn)練數(shù)據(jù)也限制了深度學(xué)習(xí)方法的性能.最后,本文方法有著與卷積神經(jīng)網(wǎng)絡(luò)相似的能夠?qū)Ψ指罱Y(jié)果進(jìn)行逐級精煉的層級架構(gòu),但參數(shù)規(guī)模卻遠(yuǎn)小于表2中對比方法的卷積神經(jīng)網(wǎng)絡(luò)的參數(shù)規(guī)模,因此本文方法能夠利用較少的訓(xùn)練數(shù)據(jù)獲得更高的分割精度.

    4 結(jié)論

    本文提出了一種整合了級聯(lián)隨機森林與活動輪廓模型的醫(yī)學(xué)磁共振圖像三維分割方法.與以往基于隨機森林的方法將分割問題定義為像素分類任務(wù)不同,本文方法通過輪廓曲線演化進(jìn)一步精煉了體素分類,由此獲得具有幾何約束的分割結(jié)果.相比于標(biāo)準(zhǔn)活動輪廓模型需要輸入初始輪廓且通常難以分割具有復(fù)雜條件的圖像的局限,本文方法可以自動地將演化曲線初始化為與目標(biāo)輪廓接近的曲線,并且通過有效的特征學(xué)習(xí)和形狀先驗整合方案提升了對復(fù)雜條件圖像的分割精度和魯棒性.實驗結(jié)果表明,本文方法能夠有效分割磁共振圖像體數(shù)據(jù)中具有復(fù)雜成像結(jié)構(gòu)的組織結(jié)構(gòu),左心房分割結(jié)果的DC值為0.90±0.03,Surfaceto-Surface值為0.64±0.39mm,腦腫瘤分割結(jié)果在Complete、Core和Enhancing三個評測區(qū)域的DC值分別達(dá)到0.87,0.74和0.64.測試時間在10s以內(nèi).多個數(shù)據(jù)集上的對比實驗結(jié)果顯示了本文方法相比于其他同類方法在分割精度和魯棒性上的優(yōu)勢.如前文所述,BRATS15數(shù)據(jù)集存在一定的局限性,因此在更廣泛數(shù)據(jù)集上對本文方法進(jìn)行驗證,并對隨機森林的不同整合方式進(jìn)行深入的探索將是我們下一步的工作內(nèi)容.

    猜你喜歡
    體素級聯(lián)輪廓
    基于超體素聚合的流式細(xì)胞術(shù)自動門控方法
    基于多級細(xì)分的彩色模型表面體素化算法
    OPENCV輪廓識別研究與實踐
    運用邊界狀態(tài)約束的表面體素加密細(xì)分算法
    基于實時輪廓誤差估算的數(shù)控系統(tǒng)輪廓控制
    基于體素格尺度不變特征變換的快速點云配準(zhǔn)方法
    級聯(lián)LDPC碼的STBC-OFDM系統(tǒng)
    電子制作(2016年15期)2017-01-15 13:39:09
    基于級聯(lián)MUSIC的面陣中的二維DOA估計算法
    在線學(xué)習(xí)機制下的Snake輪廓跟蹤
    計算機工程(2015年4期)2015-07-05 08:27:39
    LCL濾波器在6kV級聯(lián)STATCOM中的應(yīng)用
    電測與儀表(2014年1期)2014-04-04 12:00:34
    最新美女视频免费是黄的| 国产伦人伦偷精品视频| 亚洲美女黄片视频| 欧美黄色片欧美黄色片| 亚洲精品美女久久av网站| av有码第一页| 国产成人精品无人区| 正在播放国产对白刺激| 欧美乱码精品一区二区三区| 国产视频一区二区在线看| 少妇粗大呻吟视频| 波多野结衣av一区二区av| 国产一卡二卡三卡精品| 国产免费现黄频在线看| 久久性视频一级片| 在线观看免费午夜福利视频| av天堂久久9| 黑人操中国人逼视频| 国产精品国产高清国产av| 精品国产亚洲在线| 亚洲片人在线观看| 日韩欧美免费精品| 黄片大片在线免费观看| 99热国产这里只有精品6| 久久精品国产清高在天天线| 亚洲九九香蕉| 国产极品粉嫩免费观看在线| 高清在线国产一区| 日日摸夜夜添夜夜添小说| 亚洲午夜理论影院| 一区二区三区激情视频| 天天影视国产精品| 淫秽高清视频在线观看| 欧美激情极品国产一区二区三区| 最好的美女福利视频网| 亚洲成人精品中文字幕电影 | 日本五十路高清| 精品国内亚洲2022精品成人| 国产日韩一区二区三区精品不卡| 国产人伦9x9x在线观看| 亚洲国产中文字幕在线视频| 久久久国产成人精品二区 | 夫妻午夜视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产单亲对白刺激| 97超级碰碰碰精品色视频在线观看| 在线观看日韩欧美| 后天国语完整版免费观看| 香蕉丝袜av| 久久中文字幕一级| 精品久久久久久,| 成人国语在线视频| 深夜精品福利| 欧美亚洲日本最大视频资源| 亚洲三区欧美一区| av天堂久久9| 99热国产这里只有精品6| 桃红色精品国产亚洲av| 十八禁网站免费在线| 午夜福利影视在线免费观看| 国产成人精品无人区| 久久久国产精品麻豆| 亚洲色图av天堂| 亚洲人成电影免费在线| 成在线人永久免费视频| 国产成人精品在线电影| www.999成人在线观看| 黄片小视频在线播放| 成在线人永久免费视频| 免费少妇av软件| 欧美一级毛片孕妇| svipshipincom国产片| 欧美久久黑人一区二区| 国产亚洲精品一区二区www| 亚洲人成网站在线播放欧美日韩| 日韩成人在线观看一区二区三区| 午夜福利免费观看在线| 久9热在线精品视频| 人人澡人人妻人| 九色亚洲精品在线播放| 色在线成人网| 巨乳人妻的诱惑在线观看| 国产精品自产拍在线观看55亚洲| 久久久国产一区二区| 手机成人av网站| 一区福利在线观看| 一级片免费观看大全| 国产精品 欧美亚洲| 国产精品免费视频内射| 亚洲,欧美精品.| avwww免费| 欧美日韩中文字幕国产精品一区二区三区 | 欧美色视频一区免费| 国产av一区在线观看免费| 国产1区2区3区精品| 精品福利观看| 亚洲一区二区三区色噜噜 | 中亚洲国语对白在线视频| 精品欧美一区二区三区在线| 9191精品国产免费久久| a级毛片黄视频| 久久婷婷成人综合色麻豆| 久久欧美精品欧美久久欧美| 一区福利在线观看| 日韩三级视频一区二区三区| 亚洲精品粉嫩美女一区| 亚洲黑人精品在线| 人人妻人人澡人人看| 成人18禁高潮啪啪吃奶动态图| 午夜91福利影院| 精品少妇一区二区三区视频日本电影| 亚洲精品一区av在线观看| 欧美黄色片欧美黄色片| 男女下面插进去视频免费观看| 国产欧美日韩综合在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久精品亚洲精品国产色婷小说| 男人舔女人的私密视频| 亚洲熟妇中文字幕五十中出 | 午夜91福利影院| 亚洲成人精品中文字幕电影 | 日韩免费高清中文字幕av| 午夜福利欧美成人| 欧美色视频一区免费| 亚洲中文日韩欧美视频| 在线观看免费视频网站a站| 日韩高清综合在线| 男男h啪啪无遮挡| 亚洲精品一区av在线观看| 69av精品久久久久久| 国内久久婷婷六月综合欲色啪| 日韩中文字幕欧美一区二区| 啦啦啦 在线观看视频| e午夜精品久久久久久久| 欧美 亚洲 国产 日韩一| 一区福利在线观看| 99香蕉大伊视频| 99久久99久久久精品蜜桃| 巨乳人妻的诱惑在线观看| 久久精品成人免费网站| 国产成人精品久久二区二区免费| 久久天堂一区二区三区四区| 在线观看免费午夜福利视频| 亚洲va日本ⅴa欧美va伊人久久| 午夜a级毛片| 久久久久九九精品影院| 高清欧美精品videossex| 久久亚洲精品不卡| 国产熟女午夜一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲片人在线观看| 精品久久久久久电影网| 久久香蕉精品热| 日韩三级视频一区二区三区| 亚洲第一欧美日韩一区二区三区| 亚洲精品中文字幕在线视频| а√天堂www在线а√下载| 国产高清videossex| 欧美一区二区精品小视频在线| 天天添夜夜摸| 最近最新中文字幕大全电影3 | 窝窝影院91人妻| 国产男靠女视频免费网站| 夜夜爽天天搞| 在线国产一区二区在线| 波多野结衣av一区二区av| 久久久久国内视频| 亚洲自偷自拍图片 自拍| 国产精品国产av在线观看| 久9热在线精品视频| 欧美乱色亚洲激情| 亚洲av成人av| 国产97色在线日韩免费| 老汉色∧v一级毛片| 亚洲av日韩精品久久久久久密| 色综合欧美亚洲国产小说| 亚洲人成网站在线播放欧美日韩| 午夜福利免费观看在线| 夫妻午夜视频| 午夜影院日韩av| 免费看十八禁软件| 一本大道久久a久久精品| 丰满的人妻完整版| av片东京热男人的天堂| 老汉色∧v一级毛片| 黄频高清免费视频| 日韩大尺度精品在线看网址 | 性色av乱码一区二区三区2| 男女高潮啪啪啪动态图| 新久久久久国产一级毛片| 亚洲五月天丁香| 免费一级毛片在线播放高清视频 | 91av网站免费观看| 国产99久久九九免费精品| 村上凉子中文字幕在线| 亚洲色图 男人天堂 中文字幕| 国产熟女xx| 嫩草影视91久久| 午夜福利影视在线免费观看| 三级毛片av免费| 熟女少妇亚洲综合色aaa.| 人人妻人人添人人爽欧美一区卜| 侵犯人妻中文字幕一二三四区| 一a级毛片在线观看| 69av精品久久久久久| 国产亚洲欧美精品永久| 久久精品亚洲熟妇少妇任你| 久久国产乱子伦精品免费另类| 亚洲人成77777在线视频| 人人妻人人澡人人看| 可以在线观看毛片的网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲熟妇熟女久久| 在线观看日韩欧美| 热re99久久国产66热| 国产野战对白在线观看| 亚洲性夜色夜夜综合| 99riav亚洲国产免费| 露出奶头的视频| 亚洲一区高清亚洲精品| 老司机在亚洲福利影院| 欧美激情高清一区二区三区| 法律面前人人平等表现在哪些方面| 女人被躁到高潮嗷嗷叫费观| 亚洲五月婷婷丁香| 国产精品98久久久久久宅男小说| 极品教师在线免费播放| 亚洲男人的天堂狠狠| 又黄又爽又免费观看的视频| 香蕉久久夜色| 免费观看精品视频网站| www日本在线高清视频| 日韩免费av在线播放| 久久精品国产清高在天天线| 日韩人妻精品一区2区三区| 国产精品一区二区在线不卡| 超碰97精品在线观看| 中文字幕人妻丝袜一区二区| 大香蕉久久成人网| 伦理电影免费视频| 精品一区二区三区视频在线观看免费 | 91九色精品人成在线观看| 亚洲欧美日韩另类电影网站| 男人舔女人下体高潮全视频| 色综合站精品国产| 国产国语露脸激情在线看| 午夜91福利影院| 波多野结衣av一区二区av| 日日干狠狠操夜夜爽| 久久久国产欧美日韩av| 一二三四在线观看免费中文在| 91九色精品人成在线观看| 久久国产乱子伦精品免费另类| 久久精品国产亚洲av高清一级| 搡老乐熟女国产| 一进一出抽搐gif免费好疼 | 欧美中文日本在线观看视频| 亚洲一区中文字幕在线| 日韩成人在线观看一区二区三区| 久久天躁狠狠躁夜夜2o2o| 欧美日韩亚洲国产一区二区在线观看| 欧美精品啪啪一区二区三区| 日本欧美视频一区| 久久中文字幕人妻熟女| 精品午夜福利视频在线观看一区| 欧美性长视频在线观看| 日韩欧美一区视频在线观看| 嫁个100分男人电影在线观看| 人妻丰满熟妇av一区二区三区| 国产亚洲精品久久久久久毛片| 深夜精品福利| 男人的好看免费观看在线视频 | 91成年电影在线观看| 久久久精品国产亚洲av高清涩受| 国产1区2区3区精品| 国产在线观看jvid| 母亲3免费完整高清在线观看| 亚洲avbb在线观看| 久久精品91蜜桃| 国产一区二区在线av高清观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 日本五十路高清| 亚洲狠狠婷婷综合久久图片| 波多野结衣一区麻豆| 亚洲精品久久午夜乱码| 国产精品久久久av美女十八| 一级片免费观看大全| 91大片在线观看| 欧美一区二区精品小视频在线| 国产一区二区三区综合在线观看| 视频区图区小说| 欧美性长视频在线观看| 少妇被粗大的猛进出69影院| 麻豆国产av国片精品| 精品国内亚洲2022精品成人| 在线观看免费视频日本深夜| 久久中文看片网| 操出白浆在线播放| 99国产精品99久久久久| 亚洲片人在线观看| 国产极品粉嫩免费观看在线| 99久久99久久久精品蜜桃| av视频免费观看在线观看| 亚洲人成电影免费在线| 不卡av一区二区三区| 欧美午夜高清在线| 国产精品美女特级片免费视频播放器 | 国产黄a三级三级三级人| 97超级碰碰碰精品色视频在线观看| 黄色成人免费大全| 亚洲精品美女久久久久99蜜臀| 亚洲成国产人片在线观看| 妹子高潮喷水视频| 露出奶头的视频| 国产成人啪精品午夜网站| 色哟哟哟哟哟哟| 天堂中文最新版在线下载| 亚洲精品国产色婷婷电影| 在线观看免费视频日本深夜| 男女床上黄色一级片免费看| 免费搜索国产男女视频| 女人被狂操c到高潮| 精品卡一卡二卡四卡免费| 美女大奶头视频| 一级,二级,三级黄色视频| 久久热在线av| 好男人电影高清在线观看| 亚洲性夜色夜夜综合| 久久精品国产99精品国产亚洲性色 | 久久九九热精品免费| 咕卡用的链子| 亚洲第一青青草原| 国产成人av教育| 777久久人妻少妇嫩草av网站| 色综合站精品国产| 久久伊人香网站| 少妇粗大呻吟视频| 久久久久久久精品吃奶| 精品国内亚洲2022精品成人| 国产又爽黄色视频| 久久精品亚洲av国产电影网| netflix在线观看网站| 久久午夜亚洲精品久久| 窝窝影院91人妻| 18禁观看日本| 亚洲狠狠婷婷综合久久图片| 制服人妻中文乱码| 亚洲专区中文字幕在线| 国产精品 国内视频| 久久精品国产综合久久久| 国产蜜桃级精品一区二区三区| 很黄的视频免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成av片中文字幕在线观看| 免费在线观看黄色视频的| 香蕉久久夜色| 亚洲精品久久成人aⅴ小说| 一级作爱视频免费观看| 不卡av一区二区三区| 男女下面进入的视频免费午夜 | 成人18禁高潮啪啪吃奶动态图| 丰满饥渴人妻一区二区三| 国产午夜精品久久久久久| 久久这里只有精品19| 免费在线观看视频国产中文字幕亚洲| 搡老岳熟女国产| 桃红色精品国产亚洲av| 久久精品亚洲精品国产色婷小说| 另类亚洲欧美激情| 无人区码免费观看不卡| 欧美不卡视频在线免费观看 | 少妇裸体淫交视频免费看高清 | 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡| 久久精品91无色码中文字幕| 成人精品一区二区免费| 在线免费观看的www视频| 波多野结衣一区麻豆| 丝袜美足系列| 精品国产一区二区三区四区第35| 久久中文字幕一级| 9色porny在线观看| 国产精品影院久久| 久久九九热精品免费| 麻豆一二三区av精品| 国产真人三级小视频在线观看| 三级毛片av免费| 老司机深夜福利视频在线观看| 99热只有精品国产| www.www免费av| 精品欧美一区二区三区在线| 在线观看www视频免费| 久久久水蜜桃国产精品网| 桃色一区二区三区在线观看| 国产伦一二天堂av在线观看| 亚洲精品一区av在线观看| 久久人妻熟女aⅴ| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看视频国产中文字幕亚洲| 新久久久久国产一级毛片| 一级片免费观看大全| 淫妇啪啪啪对白视频| 不卡一级毛片| 叶爱在线成人免费视频播放| 亚洲人成电影免费在线| 一级黄色大片毛片| 亚洲精品粉嫩美女一区| 天堂影院成人在线观看| 老熟妇乱子伦视频在线观看| 午夜福利一区二区在线看| 母亲3免费完整高清在线观看| 黄片播放在线免费| avwww免费| 日本撒尿小便嘘嘘汇集6| 每晚都被弄得嗷嗷叫到高潮| 日韩人妻精品一区2区三区| 美国免费a级毛片| 久久精品影院6| 久久精品国产99精品国产亚洲性色 | 国产97色在线日韩免费| 午夜视频精品福利| 18美女黄网站色大片免费观看| 久久久久久久久中文| 日韩 欧美 亚洲 中文字幕| 夜夜夜夜夜久久久久| 可以在线观看毛片的网站| 欧美av亚洲av综合av国产av| 麻豆国产av国片精品| 韩国av一区二区三区四区| 神马国产精品三级电影在线观看 | 99久久精品国产亚洲精品| 激情在线观看视频在线高清| 啦啦啦 在线观看视频| 欧美大码av| 亚洲精品一区av在线观看| 一进一出抽搐动态| 亚洲精品在线美女| 亚洲国产欧美一区二区综合| 丝袜美足系列| 88av欧美| 黑人猛操日本美女一级片| 亚洲欧美日韩无卡精品| 久久久久九九精品影院| 岛国在线观看网站| 黄色视频不卡| 9色porny在线观看| 日韩 欧美 亚洲 中文字幕| 操出白浆在线播放| 一级黄色大片毛片| 香蕉丝袜av| 亚洲精品av麻豆狂野| 国产高清激情床上av| 久久午夜亚洲精品久久| 久久国产精品影院| 啦啦啦观看免费观看视频高清| 嫩草影视91久久| 给我免费播放毛片高清在线观看| 国产精品1区2区在线观看.| 最近视频中文字幕2019在线8| 男女那种视频在线观看| 最新中文字幕久久久久| 最新在线观看一区二区三区| 黄色日韩在线| 一个人免费在线观看电影| 亚洲一区二区三区不卡视频| 亚洲av成人av| 国产大屁股一区二区在线视频| 久久精品国产亚洲av香蕉五月| 麻豆av噜噜一区二区三区| 国产精品一区二区性色av| 天堂av国产一区二区熟女人妻| 黄色女人牲交| 久久草成人影院| 国产午夜精品久久久久久一区二区三区 | 好男人在线观看高清免费视频| 伦理电影大哥的女人| 国产亚洲欧美在线一区二区| 亚洲不卡免费看| 日本与韩国留学比较| 中亚洲国语对白在线视频| 久久久久久久久中文| 国产私拍福利视频在线观看| bbb黄色大片| 黄色日韩在线| 舔av片在线| 国产大屁股一区二区在线视频| 国内揄拍国产精品人妻在线| 在线免费观看不下载黄p国产 | 国产视频一区二区在线看| 丰满乱子伦码专区| 亚洲在线自拍视频| 日本 av在线| 欧美中文日本在线观看视频| 午夜福利高清视频| 日韩欧美精品免费久久 | 岛国在线免费视频观看| 人人妻人人看人人澡| 一二三四社区在线视频社区8| 免费看日本二区| 亚洲专区中文字幕在线| 国产精华一区二区三区| 白带黄色成豆腐渣| 久久久精品欧美日韩精品| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 久久国产精品人妻蜜桃| 久久亚洲精品不卡| 又爽又黄无遮挡网站| 国产精品影院久久| 亚洲人成伊人成综合网2020| 亚洲专区中文字幕在线| 亚洲最大成人手机在线| 婷婷亚洲欧美| 女人十人毛片免费观看3o分钟| 色5月婷婷丁香| 精品免费久久久久久久清纯| 少妇高潮的动态图| 免费人成视频x8x8入口观看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品影视一区二区三区av| 三级男女做爰猛烈吃奶摸视频| 性插视频无遮挡在线免费观看| 在线观看免费视频日本深夜| 亚洲av熟女| 9191精品国产免费久久| 天堂√8在线中文| 乱码一卡2卡4卡精品| 搡女人真爽免费视频火全软件 | 在线播放无遮挡| 波多野结衣高清作品| 蜜桃久久精品国产亚洲av| 国产三级在线视频| 能在线免费观看的黄片| av在线蜜桃| 最近在线观看免费完整版| 亚洲第一区二区三区不卡| 久久久色成人| 精品一区二区三区视频在线观看免费| 精品欧美国产一区二区三| 亚洲 国产 在线| 可以在线观看的亚洲视频| 久久久久久九九精品二区国产| 亚洲人成电影免费在线| 日韩人妻高清精品专区| 夜夜夜夜夜久久久久| 中文字幕人成人乱码亚洲影| 婷婷丁香在线五月| 久久婷婷人人爽人人干人人爱| 中文在线观看免费www的网站| 欧美黑人欧美精品刺激| 午夜精品一区二区三区免费看| 亚洲av成人av| 国产在线男女| 99国产综合亚洲精品| 麻豆一二三区av精品| 亚洲自偷自拍三级| 亚洲三级黄色毛片| 亚洲无线在线观看| 人人妻人人看人人澡| 亚洲av第一区精品v没综合| 欧美不卡视频在线免费观看| 最近中文字幕高清免费大全6 | 麻豆av噜噜一区二区三区| 亚洲成人精品中文字幕电影| 99久国产av精品| 国产不卡一卡二| 中文在线观看免费www的网站| 波多野结衣高清作品| 中文字幕av在线有码专区| 国产黄色小视频在线观看| 国产乱人伦免费视频| 午夜福利视频1000在线观看| 永久网站在线| 午夜精品在线福利| 麻豆久久精品国产亚洲av| 简卡轻食公司| 日本五十路高清| 精品99又大又爽又粗少妇毛片 | 精品一区二区免费观看| 久久精品国产亚洲av天美| 免费电影在线观看免费观看| 高潮久久久久久久久久久不卡| 国产一区二区在线av高清观看| 日本一本二区三区精品| 亚洲自偷自拍三级| 欧美一级a爱片免费观看看| 国产黄a三级三级三级人| 啦啦啦韩国在线观看视频| 国产成年人精品一区二区| 怎么达到女性高潮| 十八禁网站免费在线| 国产探花在线观看一区二区| 亚洲国产欧美人成| 欧美日本视频| 美女xxoo啪啪120秒动态图 | 亚洲性夜色夜夜综合| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 亚洲色图av天堂| 18禁裸乳无遮挡免费网站照片| 久久精品91蜜桃| 99精品久久久久人妻精品| 91午夜精品亚洲一区二区三区 | 国产中年淑女户外野战色| 亚洲成人精品中文字幕电影| 精品99又大又爽又粗少妇毛片 | 国产国拍精品亚洲av在线观看| 国产伦精品一区二区三区四那| 国产成人福利小说| 成人精品一区二区免费| 偷拍熟女少妇极品色| 老熟妇乱子伦视频在线观看| 真人一进一出gif抽搐免费| 嫩草影视91久久|