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

    基于自適應(yīng)符號函數(shù)的主動(dòng)輪廓模型?

    2019-10-26 18:05:44翁桂榮何志勇
    軟件學(xué)報(bào) 2019年12期
    關(guān)鍵詞:輪廓線輪廓邊緣

    翁桂榮,何志勇

    (蘇州大學(xué) 機(jī)電工程學(xué)院,江蘇 蘇州 215021)

    圖像分割技術(shù)在計(jì)算機(jī)視覺、場景識(shí)別、醫(yī)學(xué)圖像分析、移動(dòng)目標(biāo)跟蹤等領(lǐng)域具有非常重要的意義.基于曲線演化模型和偏微分方程等的分割方法在理論上已經(jīng)取得新的突破,基于水平集方法的主動(dòng)輪廓模型已成為圖像分割的熱點(diǎn)[1?4].1988年,Osher等人[5]在分析流體力學(xué)問題時(shí),創(chuàng)造性地提出水平集方法,主要思想是:把平面上的閉合曲線的演化過程上升為高一維的函數(shù)與水平面交集的演變過程,即零水平集的變化過程,通過求解零水平集函數(shù)來獲得運(yùn)動(dòng)的邊界,用于界面追蹤和形狀建模.但是傳統(tǒng)的水平集演化由于采用逆向差分的數(shù)值解法,計(jì)算效率不高.Chop[6]首先提出了窄帶方法,同時(shí)給出了重新初始化的概念,就是在數(shù)值演化過程中只更新零水平集附近帶狀區(qū)域里的水平集數(shù)值,計(jì)算量就大幅減少.但多次重新初始化,計(jì)算量也隨之增加,快速地生成符號距離函數(shù),對于提高水平集方法的速度至關(guān)重要.

    1989年,Mumford和Shah提出了一個(gè)經(jīng)典的基于能量最小化的分割模型,即Mumford-Shah(M-S)模型[7].該模型的能量泛函由保真項(xiàng)、平滑項(xiàng)和長度約束項(xiàng)組成.所謂保真項(xiàng)顯示原圖像與之分割結(jié)果的相似性;平滑項(xiàng)的作用是使得分割結(jié)果光滑;長度約束項(xiàng)用于約束模型的曲線長度.其中,保真項(xiàng)和平滑項(xiàng)利用區(qū)域信息特征,長度約束項(xiàng)基于邊緣梯度信息.

    邊緣型的方法適用于目標(biāo)邊緣清楚的圖像,對于沒有明顯梯度變化的邊緣則逼近效果不理想.最早的水平集模型是由Malladiet[8]和Caselles[9]等人提出的幾何活動(dòng)輪廓模型,該模型基于曲線的平均曲率運(yùn)動(dòng)和圖像的梯度信息,構(gòu)建了基于水平集方法的速度函數(shù).Siddiqis[10]提出在測地線活動(dòng)輪廓模型中增加面積項(xiàng)功能,目的是提高變形曲線跨越輪廓線中較小縫隙的能力.

    2001年,Chan-Vese提出了基于簡化M-S能量泛函的CV模型[11],該模型不依賴于圖像的梯度信息,對噪聲不敏感,較為適用于分割梯度無意義或邊緣模糊的圖像.Kimmel在模型[12]中集成了邊緣的方向性信息,并通過最小化邊緣和圖像梯度間夾角的方式對測地線活動(dòng)輪廓模型進(jìn)行了改良.

    保證演化方程處于穩(wěn)定狀態(tài)是幾何主動(dòng)輪廓模型初始化的標(biāo)準(zhǔn)方法,理想狀態(tài)下使得演化方程穩(wěn)定的狀態(tài)解是一個(gè)符號距離函數(shù).幾何主動(dòng)輪廓模型的問題主要有兩個(gè),一是通過重新初始化得到規(guī)則的水平集函數(shù)并不能保證零水平集函數(shù)可以運(yùn)動(dòng)至目標(biāo)邊界,二是在實(shí)際操作過程中如何選擇重新初始化的時(shí)間.

    Li針對上述問題,提出了解決重新初始化的方法[13],即距離規(guī)則的水平集(DRLSE)方法,通過正規(guī)項(xiàng)的作用,保證水平集函數(shù)是符號距離函數(shù).其構(gòu)造初始超平面也非常簡潔,提高了分割速度.為了進(jìn)一步提高模型的運(yùn)行速度,并克服單阱勢函數(shù)的缺點(diǎn),Li構(gòu)建了一個(gè)雙阱函數(shù),發(fā)展出快速距離規(guī)則水平集[14]以及區(qū)域擬合驅(qū)動(dòng)的主動(dòng)輪廓模型(RSF)[15,16],利用勢函數(shù)的特點(diǎn),主動(dòng)輪廓模型在演化過程中就不需要重新初始化.

    但是,DRLSE,RSF模型未能解決另一個(gè)核心問題,即對初始輪廓敏感.此外,DRLSE模型還存在下述不穩(wěn)定的問題,如抗噪聲能力差,收斂速度較慢,在演化過程中容易陷入虛假邊界,會(huì)從弱邊緣處泄露[17].

    Wang提出了一種基于邊緣的自適應(yīng)水平集演化方程[18],將面積項(xiàng)系數(shù)修改為可變符號函數(shù),達(dá)到改善對初始輪廓敏感的目的,克服了DRLSE單向運(yùn)動(dòng)的缺點(diǎn).該模型具有較好的操作性,分割速度大幅提高,特別是在邊界較清晰或背景噪聲低的情況下,能達(dá)到很好的分割效果.但是該模型有抗噪聲能力較差、容易陷入虛假邊界以及從弱邊緣處泄露等不穩(wěn)定問題[19,20].

    本文提出了一種致力于提高初始輪廓魯棒性的主動(dòng)輪廓模型,設(shè)計(jì)了一個(gè)自適應(yīng)符號函數(shù)結(jié)合自適應(yīng)邊緣指示函數(shù)的能量泛函,使得主動(dòng)輪廓演化能根據(jù)自適應(yīng)符號函數(shù)的方向從初始輪廓開始自動(dòng)進(jìn)行膨脹及收縮.本文模型提高了對初始輪廓的魯棒性,同時(shí)解決了主動(dòng)輪廓模型收斂速度慢以及在弱邊緣處分割錯(cuò)誤的問題.此外,為了演化的穩(wěn)定性,還改進(jìn)了距離規(guī)則項(xiàng).

    1 經(jīng)典的主動(dòng)輪廓模型

    1.1 DRLSE模型

    設(shè)置一個(gè)灰度圖像I:Ω→R,其圖像范圍為Ω?R2.

    距離規(guī)則水平集[14]模型在傳統(tǒng)的模型上添加了距離規(guī)則功能,其主要作用是解決水平集函數(shù)與符號距離函數(shù)之間的偏差,使得水平集函數(shù)在迭代過程中不再考慮重新初始化的問題.該模型的能量泛函定義為

    公式(1)中,p(?)為勢阱函數(shù),?為梯度算子,φ為水平集函數(shù),δε(?)是Dirac函數(shù),Hε(?)是Heaviside函數(shù),g(?)為邊緣指示函數(shù),μ,λ,α為常數(shù).

    其中,c0為正值常數(shù),R0代表圖像區(qū)域中的某塊區(qū)域.

    DRLSE模型為水平集方法提供了新的研究思路,其存在如下問題.

    1) 用于加快零水平集演化速率的面積項(xiàng)為單一值(正或負(fù)值),在能量泛函進(jìn)行最小化中,只能是從正(負(fù))趨向0的單向選擇,即沒有自我調(diào)整能力;

    2) 面積項(xiàng)及長度項(xiàng)高度依賴邊緣指示函數(shù),而邊緣指示函數(shù)由圖像的梯度構(gòu)建,當(dāng)梯度值很大時(shí),其函數(shù)值越接近0.經(jīng)過高斯濾波后的目標(biāo)邊界會(huì)變得模糊變寬,導(dǎo)致目標(biāo)邊界區(qū)域變寬.當(dāng)目標(biāo)之間的距離很近時(shí),目標(biāo)邊界可能會(huì)因高斯平滑而相互連接,導(dǎo)致分割失??;

    3) 在實(shí)驗(yàn)中常數(shù)α需要人為設(shè)置,其取值大小對分割的結(jié)果影響很大.α過小會(huì)使演化速度降低,從而陷入虛假邊界,導(dǎo)致演化停止;α過大,會(huì)使邊界處的分割出現(xiàn)泄漏,導(dǎo)致分割失??;

    4) DRLSE中的雙阱勢dp(s)是一個(gè)三段式函數(shù),當(dāng)分母s為0時(shí),無法計(jì)算;該雙阱勢函數(shù)在零點(diǎn)附近演化速度接近為1,為演化最大值,使得水平集迅速向前進(jìn)行演化,從而讓圖像里的弱邊界s急劇減小至0,導(dǎo)致演化曲線侵入被分割目標(biāo)內(nèi)部,分割可能失敗[17].另外,在|?φ|=1附近時(shí),其斜率小,演化速度慢.

    1.2 RSF模型

    Li等人提出了一種區(qū)域擬合驅(qū)動(dòng)的主動(dòng)輪廓(RSF)模型[16],也稱局部二值擬合(LBF)模型[15].RSF模型可以在較小鄰域內(nèi)精確找出細(xì)小目標(biāo),并加快大目標(biāo)在大鄰域處的演化速度.RSF模型的核心思想是利用圖像局部區(qū)域像素的灰度信息,用局部二值擬合能量項(xiàng)代替全局二值擬合能量項(xiàng),通過調(diào)節(jié)核函數(shù)的方差σ控制中心像素的鄰域.根據(jù)輪廓定義了一個(gè)數(shù)據(jù)擬合能量公式:

    其中,Kσ(?)為Gaussian函數(shù),y?x代表x和y像素點(diǎn)的距離,H(?)是Heaviside函數(shù),φ(?)是Lipschitz函數(shù),擬合值f1(?),f2(?)局部地逼近輪廓外、內(nèi)圖像灰度的平均值,I(?)是灰度值,λ1,λ2為常數(shù).

    公式(8)中,H(?)不是一個(gè)連續(xù)函數(shù),用H(?)代替[11].RSF模型的梯度下降流方程如下:

    其中,μ,ν為常數(shù),e1(?),e2(?)表達(dá)如下:

    其梯度下降流方程中:第1項(xiàng)表示擬合部分;第2項(xiàng)對零水平輪廓有縮短長度或者光滑的效果,對于維持輪廓的規(guī)則性是必不可少的;第3項(xiàng)是內(nèi)部能量項(xiàng),由單勢阱距離規(guī)則組成,主要作用是矯正水平集與符號距離間的誤差,避免曲線在演化過程中重復(fù)初始化.

    RSF模型對于灰度不均勻的圖像分割效果顯著,而對于亮度不均勻的圖像卻無法進(jìn)行分割.另外,如果初始輪廓設(shè)置不當(dāng),RSF模型往往會(huì)陷入局部極小.

    1.3 自適應(yīng)水平集演化模型

    Wang和He[18]提出了輪廓提取的自適應(yīng)水平集演化模型(ALSE)解決初始輪廓敏感性的方法,其水平集函數(shù)的偏微分方程即其梯度下降流方程如下[18]:

    分析文獻(xiàn)[18]我們可以發(fā)現(xiàn),ALSE模型引入由圖像信息調(diào)整的權(quán)系數(shù)v(I)代替常值權(quán)系數(shù)α.在DRLSE模型的基礎(chǔ)上增加了一個(gè)符號函數(shù),符號函數(shù)的方向由I(x,y)與輪廓線外、內(nèi)圖像的平均值的均值差值決定,克服了DRLSE模型中面積項(xiàng)為單一方向運(yùn)動(dòng)的問題.因此,偏微分方程能根據(jù)初始輪廓內(nèi)外的圖像灰度信息調(diào)整運(yùn)動(dòng)方向,改善模型對初始輪廓的魯棒性.

    該模型具有較好的操作性,分割速度大幅提高,特別是在邊界清晰或背景噪聲低的情況下,能發(fā)揮很好的分割效果.但是該模型僅僅修正了DRLSE模型中的方向,而DRLSE模型中的抗噪聲能力差、容易陷入虛假邊界以及從弱邊緣處泄露等不穩(wěn)定問題,沒有得到進(jìn)一步解決.

    2 自適應(yīng)符號函數(shù)的主動(dòng)輪廓模型

    針對DRLSE,RSF初始輪廓位置敏感性,本文設(shè)計(jì)了一個(gè)自適應(yīng)符號函數(shù)用于提高初始輪廓的魯棒性.為了進(jìn)一步克服ALSE及上述DRLSE的幾個(gè)缺點(diǎn),改進(jìn)了邊緣指示函數(shù),加快了演化速度.此外,還修正了距離正則項(xiàng),使得演化更穩(wěn)定.其梯度下降流方程為

    2.1 自適應(yīng)邊緣指示函數(shù)分析

    公式(18)中,圖像的標(biāo)準(zhǔn)差(std2)反映了圖像像素值與均值的離散程度[21].公式(19)中β的大小表達(dá)了目標(biāo)相距平均值的遠(yuǎn)近關(guān)系.

    公式(19)能根據(jù)目標(biāo)圖像信息自適應(yīng)調(diào)整gβ函數(shù)的曲線.公式(18)中的β如同公式(6)中α的作用:當(dāng)圖像目標(biāo)相距平均值較遠(yuǎn)時(shí),圖像的方差S較大,β,gβ對應(yīng)增大,面積項(xiàng)能量隨之提高,可以使演化速度變快,并解決了分割時(shí)演化曲線進(jìn)入虛假邊界的問題;圖像的方差S較小時(shí),邊緣處會(huì)模糊,β,gβ隨之減少,使得面積項(xiàng)能量降低,降低演化速度且不會(huì)出現(xiàn)邊界泄漏現(xiàn)象.圖1所示為兩種邊緣指示函數(shù)的對比圖(兩種邊緣指示函數(shù)的對比:藍(lán)色線為DRLSE方法;紅色為本文方法,共3條:β=0.75,β=1.0,β=1.25),自適應(yīng)邊緣指示函數(shù)對遠(yuǎn)離邊界(梯度值小)的區(qū)域給予較大的值,此區(qū)域由邊界閾值β自適應(yīng)根據(jù)圖像特征決定,有利于加快演化速度,同時(shí)避免陷入虛假邊界中;接近邊界(梯度值大)的區(qū)域給予較小的值,當(dāng)零水平集運(yùn)動(dòng)到目標(biāo)邊緣時(shí),演化速度迅速減少到0,避免出現(xiàn)邊界泄漏現(xiàn)象.該函數(shù)能夠使零水平集保持較好的光滑性,從而加快演化速度.β又叫邊界閾值,是調(diào)節(jié)gβ函數(shù)的重要參數(shù).

    Fig.1 Adaptive boundary indicator function圖1 自適應(yīng)邊緣指示函數(shù)

    為驗(yàn)證本文提出的自適應(yīng)邊緣指示函數(shù)的有效性,根據(jù)DRLSE模型,僅改變gβ函數(shù),得到如圖2(d)~圖2(f)中所示的分割結(jié)果,其中,圖2(a)~圖2(c)是DRLSE分割結(jié)果(圖2(a)~圖2(c)為DRLSE中的邊緣指示函數(shù)方法,圖2(d)~圖2(f)為自適應(yīng)邊緣指示函數(shù)方法,綠線為初始輪廓,青線為迭代過程輪廓線,紅線為分割結(jié)果).

    Fig.2 Segmentation results with two boundary indicator functions圖2 兩種邊緣指示函數(shù)方法對比分割圖

    每隔10次迭代計(jì)算一次分割面積,當(dāng)前后兩次面積相對誤差小于1%時(shí),自動(dòng)停止迭代,N為迭代次數(shù).DRLSE實(shí)驗(yàn)參數(shù):Δt=1,μ=0.2,λ=5,ε=1.0,α=1.5,c0=2,σ=1.5.實(shí)驗(yàn)結(jié)果表明:本文自適應(yīng)邊緣指示函數(shù)模型在保持較快演化速度的同時(shí),可以有效避免陷入虛假邊界中.

    2.2 雙阱勢演化速度dp(s)

    DRLSE的dp(s)必須考慮s=0的特殊情況,且在零點(diǎn)附近演化速度接近為1,達(dá)到演化最大值.如圖3所示,此情況下水平集迅速向前演化,有可能導(dǎo)致零水平集進(jìn)入目標(biāo)內(nèi)部,而本文的dp(s)降低了在零阱勢處的演化速度,提高了對弱邊界的魯棒性.同時(shí),從圖3看出:本文的雙阱勢在s=1的鄰域內(nèi)比DRLSE的雙阱勢更加陡峭,這就使得保持符號距離函數(shù)的作用更加靈敏.

    Fig.3 Comparison between theevolution speedof two double-well evolution speed functions圖3 兩種雙阱勢函數(shù)演化速度的對比

    對比實(shí)驗(yàn)條件相同,參數(shù)同上,用DRLSE的方法僅改變雙阱勢函數(shù),進(jìn)行分割比較,結(jié)果如圖4所示(第1排為原圖,第2排為DRLSE的雙阱勢函數(shù),第3排為本文的雙阱勢函數(shù);綠線為初始輪廓,青線為迭代過程輪廓線,紅線為分割結(jié)果).

    Fig.4 Segmentation results withtwo double-well evolution speed functions圖4 兩種雙阱勢函數(shù)演化對比效果圖

    第1列為了說明兩種雙阱勢函數(shù)的穿透邊界的情況,將α從1.5增加到2.5,加快演化速度,可以看到,DRLSE模型出現(xiàn)了邊界泄露的情況,而運(yùn)用了本文雙阱勢函數(shù)的模型則得到了理想的結(jié)果;第2列說明在正常演化速度(α=1.5)的情況下,避免了目標(biāo)遺漏;第3列說明本文方法的dp(s)可以加快演化速度.以上結(jié)果說明,本文的雙阱勢函數(shù)對于曲線的演化有較好的穩(wěn)定性.

    2.3 自適應(yīng)符號函數(shù)分析

    DRLSE模型及相關(guān)論文中涉及的模型,在初始輪廓線處的曲線迭代過程具有單向性[20,21].α>0時(shí)為收縮運(yùn)動(dòng),α<0時(shí)為膨脹運(yùn)動(dòng),直至梯度下降流方程為0時(shí)停止.因此,初始曲線必須完全設(shè)置在目標(biāo)內(nèi)部或目標(biāo)外部.為了克服DRLSE模型對初始輪廓較為敏感的缺點(diǎn),本文增加一個(gè)圖像自適應(yīng)符號函數(shù)ψ(Iσ,c1,c2),使模型能根據(jù)圖像目標(biāo)自適應(yīng)改變水平集的運(yùn)動(dòng)方向.(c1+c2)/2為輪廓線內(nèi)、外圖像均值的平均值,偏黑區(qū)域< (c1+c2)/2<偏白區(qū)域,Iσ?(c1+c2)/2表示偏黑的區(qū)域?yàn)樨?fù)值,而偏白的區(qū)域?yàn)檎?,這就產(chǎn)生了正負(fù)號,其絕對值大小反映了偏離均值的強(qiáng)度,但要對此強(qiáng)度值進(jìn)行進(jìn)一步處理以便加入模型中.如圖5所示,arctan((Iσ?(c1+c2)/2)/τ)函數(shù)平滑了各像素點(diǎn)偏離均值的強(qiáng)度,限定了最大幅度,使目標(biāo)邊界(零點(diǎn))兩側(cè)的符號相反,從而讓零水平集可以自適應(yīng)地選擇演化方向,當(dāng)輪廓線嵌入目標(biāo)內(nèi)時(shí),輪廓線開始進(jìn)行膨脹運(yùn)動(dòng);當(dāng)輪廓線在目標(biāo)外時(shí),輪廓線開始進(jìn)行收縮運(yùn)動(dòng).不僅如此,arctan((Iσ?(c1+c2)/2)/τ)在邊界(零點(diǎn))附近具有很強(qiáng)的靈敏度,可以較精確地表達(dá)邊界的能量變化.在遠(yuǎn)離邊界的區(qū)域,自適應(yīng)符號函數(shù)絕對值較大,可以加快演化速度;在接近邊界的區(qū)域,自適應(yīng)符號函數(shù)快速下降,演化速度迅速降低,避免出現(xiàn)邊界泄漏現(xiàn)象.其中,τ(tau)可以調(diào)整強(qiáng)度大小.

    Fig.5 Adaptive sign function圖5 自適應(yīng)符號函數(shù)

    3 實(shí)驗(yàn)分析

    自適應(yīng)符號函數(shù)的主動(dòng)輪廓模型算法如下.

    1) 參數(shù)常規(guī)設(shè)置如下:α=15,Δt=5,μ=0.008,ε=1.0,λ=5,σ=1.5,c0=2,τ=35.內(nèi)迭代次數(shù)=5,外迭代次數(shù)=80.另外,如果目標(biāo)太小,Gaussian濾波會(huì)讓邊界模糊,就可能丟失細(xì)小的目標(biāo),可以取σ=0.7;

    2) 實(shí)驗(yàn)圖像進(jìn)Gaussian濾波后,計(jì)算其邊緣指示函數(shù)gβ(I)、自適應(yīng)符號系數(shù)ψ(I)及β(I);

    3) 設(shè)置初始輪廓線參數(shù);

    4) 根據(jù)公式(16)進(jìn)行參數(shù)計(jì)算;

    5) 每迭代5次,計(jì)算水平集函數(shù)面積變化,如果變化小于1%,則停止迭代并輸出分割參數(shù);

    6) 檢測是否完成總迭代次數(shù)?是,輸出分割參數(shù).否則,返回到步驟4);

    第5)步僅僅提供一個(gè)客觀的評判標(biāo)準(zhǔn),并非必要的步驟.

    3.1 模型中各重要函數(shù)的實(shí)驗(yàn)分析

    公式(16)中,本文修改了面積項(xiàng)的能量公式,面積項(xiàng)的能量ψ(Iσ,c1,c2)?gβ的最小值由當(dāng)中的兩個(gè)函數(shù)共同決定.當(dāng)演化輪廓曲線C運(yùn)動(dòng)到邊界區(qū)域時(shí),自適應(yīng)符號函數(shù)ψ(Iσ,c1,c2)到達(dá)最小值.邊緣指示函數(shù)gβ與ψ(Iσ,c1,c2)同樣具有引導(dǎo)主動(dòng)輪廓曲線向同質(zhì)邊界運(yùn)動(dòng)的功能,在邊界處也達(dá)到最小值,提高了模型的魯棒性.同時(shí),由于ψ(Iσ,c1,c2)?gβ具有不同質(zhì)地區(qū)域符號相反,因此能改變運(yùn)動(dòng)方向,很好地提高了模型對初始輪廓的魯棒性.圖6用圖解的方式表達(dá)了重要函數(shù)的演化.圖6(a)為分析原圖像及初始輪廓線.圖6(b)為自適應(yīng)邊緣指示函數(shù)效果圖.圖6(c)為ψ(Iσ,c1,c2)在初始輪廓下的效果圖,白色區(qū)域?yàn)檎?,黑色區(qū)域?yàn)樨?fù)值,紅色線為正負(fù)區(qū)域的邊界.圖6(d)為ψ(Iσ,c1,c2)?gβ在初始輪廓下的效果圖,在目標(biāo)邊緣趨近于0,兩側(cè)符號相反.此系數(shù)在目標(biāo)邊緣附近比圖6(c)更平滑,紅色線為正負(fù)區(qū)域的邊界.圖6(e)為在初始輪廓下,第1次迭代后的面積項(xiàng)β(I)ψ(Iσ,c1,c2)gβδε(φ)的能量圖,紅色線為正負(fù)區(qū)域能量的邊界,綠線為初始輪廓線.初始輪廓線切割目標(biāo)區(qū)(負(fù)能量區(qū)域)的線段向外膨脹,以求能量趨于零方向,初始輪廓線在目標(biāo)區(qū)(負(fù)能量區(qū)域)外的線段為正能量,因此,此部分線段就做收縮運(yùn)動(dòng),向能量減小的方向運(yùn)動(dòng),以求能量趨于零方向.另外,由于δε(φ)函數(shù)的因素,僅僅考慮了初始輪廓線附近的能量,遠(yuǎn)離輪廓曲線的能量就大幅減弱.圖6(f)為第1次迭代分割結(jié)果.

    Fig.6 The segmentation results withparameters in the iteration圖6 迭代前后各參數(shù)的效果圖

    圖7為幾次迭代前后,輪廓曲線變化的效果圖(第1行為第2次迭代前后的效果圖;第2行為第5次迭代前后的效果圖;第3行為最后一次(25)迭代前后的效果圖).

    Fig.7 The segmentation process of the proposed model圖7 迭代效果圖

    圖7(a)為自適應(yīng)符號函數(shù)ψ(Iσ,c1,c2)迭代前后的正負(fù)區(qū)域的變化,綠色線為迭代前的,紅色線為迭代后;圖7(b)為第n次迭代后的面積項(xiàng)能量圖,紅色線為正負(fù)區(qū)域能量的邊界,綠色線為前一次迭代結(jié)果,同時(shí)作為下次迭代的初始輪廓線;圖7(c)為迭代分割結(jié)果,綠色線為初始輪廓,紅色線為分割結(jié)果.

    3.2 初始輪廓敏感性實(shí)驗(yàn)

    本文模型的首要目的就是要提高水平集模型對初始輪廓的魯棒性,如圖8所示,選取4副灰度圖像,每個(gè)灰度圖像設(shè)置3個(gè)不同形狀或位置的初始輪廓線,分別得到演化過程及分割結(jié)果(綠線為初始輪廓,青色線為迭代過程的輪廓線(子圖),紅線為分割結(jié)果,N為迭代次數(shù)).實(shí)驗(yàn)說明:本文模型可以選取不同的初始輪廓線,可以從不同的位置嵌入或包含目標(biāo)邊界,對于有空洞的目標(biāo),也能正確有效地分割.

    Fig.8 The evolution process and segmentation results with the proposed model under different initial contours圖8 使用本文模型,選取不同的初始輪廓、所得到的演化過程及分割結(jié)果

    3.3 評價(jià)實(shí)驗(yàn)

    為說明本文模型的速度及精確性,引入一個(gè)標(biāo)準(zhǔn)圖庫,分割評價(jià)數(shù)據(jù)庫[24].本文模型及ALSE模型采用固定的邊框初始輪廓.由于RSF模型的初始輪廓必須要精心選取才能正確分割,否則分割就會(huì)失敗,考慮到DRLSE模型的初始輪廓要包含目標(biāo)且初始輪廓線的位置會(huì)嚴(yán)重影響速度,故本節(jié)不參與比較.圖9所示是采用分割評價(jià)數(shù)據(jù)庫圖像選取固定的初始輪廓(95:M-70,100:N-80,M及N為圖像的大小)及分割結(jié)果.綠線為初始輪廓,紅線為分割結(jié)果.第1列為評價(jià)數(shù)據(jù)庫原圖;第2列為數(shù)據(jù)庫分割圖;第3列為RSF模型分割圖;第4列為ALSE模型分割圖;第5列為本文模型分割圖:第4列、第5列的初始輪廓為95:M-70,100:N-80(M及N為圖像的大小).

    Fig.9 The segmentation results with three different models圖9 3種模型評價(jià)數(shù)據(jù)庫的分割結(jié)果

    為了更加準(zhǔn)確地將本模型的分割結(jié)果與專家分割進(jìn)行定量[25]比較,本文采用基于面積的dice similarity coefficient(DSC)[23]來度量模型算法的性能,定義如下:

    其中,N1和N2分別代表得到的分割和真實(shí)的分割區(qū)域,S(Ni)是指分割的區(qū)域面積.DSC的值越接近1,其分割效果越好.表1顯示了分別用RSF模型、ALSE模型及本文模型對分割評價(jià)數(shù)據(jù)庫[25]中的圖像進(jìn)行分割所得結(jié)果的DSC計(jì)算以及分割時(shí)間記錄.表1數(shù)據(jù)表明:3種模型的準(zhǔn)確度均較高且基本相同,但在計(jì)算時(shí)間上有很大差別.本文模型速度上遠(yuǎn)遠(yuǎn)優(yōu)于RSF模型,也優(yōu)于ALSE模型.

    Table 1 The number of time spent and DSC in Fig.9 under the same initialcontour (s)表1 比較RSF、ALSE及本文模型的分割圖像相似系數(shù)值及分割時(shí)間(DSC/Time(s)) (秒)

    3.4 實(shí)驗(yàn)比較

    為了說明本文模型對亮度不均及弱邊界圖像的分割效果,采用固定初始輪廓線(均采用距邊框14個(gè)像素),選取了6副圖像對DRLSE,RSF,ALSE及本文模型進(jìn)行比較,如圖10所示(第1列:原圖;第2列:DRLSE模型;第3列:RSF模型;第4列:ALSE模型;第5列:本文模型.綠線為初始輪廓(均采用距離邊框14個(gè)像素為初始輪廓),青色線表示其迭代過程中輪廓線的演化(子圖),紅線為分割結(jié)果,N為迭代次數(shù)).

    Fig.10 Segmentation results with four models under the same initial contour圖10 4種模型在相同初始輪廓下的分割效果圖

    Fig.10 Segmentation results with four models under the same initial contour (Continued)圖10 4種模型在相同初始輪廓下的分割效果圖(續(xù))

    表2數(shù)據(jù)及圖10的實(shí)驗(yàn)結(jié)果表明:DRLSE模型分割速度慢,如果初始輪廓與目標(biāo)交叉,分割就失敗,但對弱邊界圖像的處理較好;RSF模型對初始輪廓及亮度不均的圖像敏感;ALSE模型分割速度比較快,提高了初始輪廓的魯棒性,對亮度不均及弱邊界的圖像比較敏感;本文模型分割速度較快,同時(shí)提高了對初始輪廓、亮度不均及弱邊界的圖像的魯棒性.

    Table 2 The number of time spent with four models (s)表2 DRLSE模型、RSF模型、ALSE模型及本文模型的分割速度 (秒)

    為了進(jìn)一步分析本文模型對亮度不均及弱邊界圖像的分割效果的魯棒性,如圖11所示,將本文模型與Canny,Sobel及Laplacian of a Gaussian這3種經(jīng)典的邊緣檢測模型進(jìn)行分割比較(第1列:原圖;第2列:Canny模型,參數(shù)T=0.3,sigma=1.5;第3列Sobel模型,采用系統(tǒng)默認(rèn)值;第4列:Laplacian of a Gaussian模型,參數(shù)T=0.3,sigma=2.8;第5列:本文模型.紅色線為分割結(jié)果).圖11表明:本文模型的分割結(jié)果具有較高的精確且光滑的曲線及較好的魯棒性,其分割效果優(yōu)于其他幾種經(jīng)典模型.

    Fig.11 Contrast experiments results of the proposed model and three classical edge-based models圖11 本文模型與3種經(jīng)典的邊緣檢測模型的分割效果圖

    Fig.11 Contrast experiments results of the proposed model and three classical edge-based models (Continued)圖11 本文模型與3種經(jīng)典的邊緣檢測模型的分割效果圖(續(xù))

    綜上所述可以發(fā)現(xiàn):本模型采用自適應(yīng)符號函數(shù)ψ(I),通過計(jì)算輪廓線內(nèi)、外的灰度平均值,從而大致計(jì)算出目標(biāo)區(qū)域,在迭代后不斷改善目標(biāo)區(qū)域.函數(shù)在目標(biāo)區(qū)域邊緣趨近于0,兩側(cè)符號相反,使得面積項(xiàng)能量有正負(fù)之值.因此,構(gòu)建的面積項(xiàng)自適應(yīng)地實(shí)現(xiàn)收縮及膨脹,控制零水平集曲線向ψ=0處演化,從而正確識(shí)別目標(biāo)的邊界.

    4 總結(jié)

    本文提出了一種邊緣型主動(dòng)輪廓模型.通過輪廓線內(nèi)、外的灰度平均值,構(gòu)建了一個(gè)自適應(yīng)符號函數(shù).該函數(shù)不僅在目標(biāo)兩側(cè)具有相反的符號,而且能夠反映出演化曲線距離目標(biāo)邊界的遠(yuǎn)近程度并自動(dòng)據(jù)此調(diào)節(jié)演化速度,可以使曲線在迭代過程中自適應(yīng)地調(diào)整運(yùn)動(dòng)的方向.因此,可以設(shè)置不同的初始輪廓線參數(shù).利用本文提出的自適應(yīng)邊緣指示函數(shù),進(jìn)一步使得演化曲線改善容易陷入虛假邊界及易在邊界泄漏的問題,同時(shí)加快了演化速度.改進(jìn)了距離規(guī)則項(xiàng),實(shí)現(xiàn)對弱邊界的準(zhǔn)確檢測,提高了演化的穩(wěn)定性.通過與分割評價(jià)數(shù)據(jù)庫進(jìn)行對比可以看出,本文模型有較高的分割精度.通過與3個(gè)著名的模型對比更進(jìn)一步說明:本文模型分割速度較快,提高了對初始輪廓的魯棒性,抗噪聲能力強(qiáng).

    猜你喜歡
    輪廓線輪廓邊緣
    OPENCV輪廓識(shí)別研究與實(shí)踐
    基于HTML5的凸輪廓線圖解法App教學(xué)軟件研究
    基于實(shí)時(shí)輪廓誤差估算的數(shù)控系統(tǒng)輪廓控制
    一張圖看懂邊緣計(jì)算
    節(jié)日帽
    在線學(xué)習(xí)機(jī)制下的Snake輪廓跟蹤
    多輪廓線的三維形體重構(gòu)技術(shù)研究與實(shí)現(xiàn)*
    基于鼻子下輪廓線的鼻尖定位法
    創(chuàng)造早秋新輪廓
    舒適廣告(2008年9期)2008-09-22 10:02:48
    在邊緣尋找自我
    雕塑(1999年2期)1999-06-28 05:01:42
    久久久久国产一级毛片高清牌| 婷婷六月久久综合丁香| 男女午夜视频在线观看| 在线永久观看黄色视频| av在线天堂中文字幕| 午夜免费激情av| 午夜福利成人在线免费观看| 国产精品久久视频播放| 免费看日本二区| 啦啦啦观看免费观看视频高清| 亚洲欧美精品综合一区二区三区| 久久婷婷人人爽人人干人人爱| 搡老妇女老女人老熟妇| 欧美中文日本在线观看视频| 国产午夜精品论理片| 亚洲中文字幕一区二区三区有码在线看 | 成人三级做爰电影| 精品高清国产在线一区| 一进一出抽搐gif免费好疼| 亚洲精品久久成人aⅴ小说| 十八禁网站免费在线| 在线观看www视频免费| 久久午夜综合久久蜜桃| 日日爽夜夜爽网站| 国产精品99久久99久久久不卡| 麻豆成人av在线观看| 在线观看www视频免费| 久久欧美精品欧美久久欧美| 在线观看一区二区三区| 99精品久久久久人妻精品| 久久草成人影院| 国产成人精品久久二区二区免费| www.自偷自拍.com| 亚洲av日韩精品久久久久久密| 国产不卡一卡二| 国产野战对白在线观看| 亚洲欧美一区二区三区黑人| 欧美成人免费av一区二区三区| 欧美绝顶高潮抽搐喷水| 最好的美女福利视频网| 搡老妇女老女人老熟妇| 老司机午夜福利在线观看视频| 国产熟女xx| 午夜亚洲福利在线播放| www.www免费av| 国产精品永久免费网站| 精品一区二区三区视频在线观看免费| 男女床上黄色一级片免费看| 欧美日韩一级在线毛片| 亚洲国产日韩欧美精品在线观看 | 午夜影院日韩av| 国产午夜精品久久久久久| 久久午夜亚洲精品久久| 黄色毛片三级朝国网站| 男女午夜视频在线观看| 性色av乱码一区二区三区2| 国产成人欧美在线观看| 国产精品99久久99久久久不卡| 亚洲片人在线观看| 高清毛片免费观看视频网站| 国产探花在线观看一区二区| 国内少妇人妻偷人精品xxx网站 | 久久这里只有精品19| 久99久视频精品免费| 精品熟女少妇八av免费久了| 久久伊人香网站| 一二三四在线观看免费中文在| 美女午夜性视频免费| av视频在线观看入口| 国产v大片淫在线免费观看| 国产免费男女视频| 男女做爰动态图高潮gif福利片| 在线国产一区二区在线| 99久久久亚洲精品蜜臀av| 欧美黄色片欧美黄色片| 国内精品久久久久久久电影| 一本大道久久a久久精品| 无遮挡黄片免费观看| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美激情综合另类| 国产精品一区二区三区四区免费观看 | 日韩欧美国产一区二区入口| 欧美黄色片欧美黄色片| 日韩精品青青久久久久久| 亚洲自拍偷在线| 天堂动漫精品| 日本一二三区视频观看| 午夜亚洲福利在线播放| 午夜福利在线在线| 天天添夜夜摸| 亚洲av电影不卡..在线观看| 日本精品一区二区三区蜜桃| 久久精品91无色码中文字幕| 精品国内亚洲2022精品成人| 国内少妇人妻偷人精品xxx网站 | 午夜福利免费观看在线| 女生性感内裤真人,穿戴方法视频| 久久久久久国产a免费观看| 成人手机av| 性欧美人与动物交配| 欧美久久黑人一区二区| 国产精品综合久久久久久久免费| 欧美日韩中文字幕国产精品一区二区三区| 欧美zozozo另类| 亚洲免费av在线视频| 日本熟妇午夜| 久久久久久久久久黄片| 国产单亲对白刺激| 日本黄色视频三级网站网址| 欧美一区二区国产精品久久精品 | 一级黄色大片毛片| 美女黄网站色视频| 久久久精品国产亚洲av高清涩受| 亚洲va日本ⅴa欧美va伊人久久| 婷婷六月久久综合丁香| 亚洲成av人片免费观看| 亚洲国产欧美网| 国内久久婷婷六月综合欲色啪| 日本黄色视频三级网站网址| 九色成人免费人妻av| 国产av又大| 久久 成人 亚洲| 欧美又色又爽又黄视频| 亚洲七黄色美女视频| svipshipincom国产片| 日本黄色视频三级网站网址| 色综合站精品国产| 美女扒开内裤让男人捅视频| 亚洲 欧美 日韩 在线 免费| 床上黄色一级片| 国内精品久久久久久久电影| 亚洲va日本ⅴa欧美va伊人久久| 草草在线视频免费看| 丝袜人妻中文字幕| av中文乱码字幕在线| 他把我摸到了高潮在线观看| 久久久国产成人免费| 中文字幕高清在线视频| 日韩 欧美 亚洲 中文字幕| www.www免费av| 啦啦啦韩国在线观看视频| 久久欧美精品欧美久久欧美| 亚洲aⅴ乱码一区二区在线播放 | 天天添夜夜摸| 91大片在线观看| 悠悠久久av| 嫩草影院精品99| 精品午夜福利视频在线观看一区| 动漫黄色视频在线观看| 在线观看美女被高潮喷水网站 | av天堂在线播放| 黄色a级毛片大全视频| 精品国产超薄肉色丝袜足j| 亚洲av电影在线进入| 精品无人区乱码1区二区| 国产精品乱码一区二三区的特点| 亚洲18禁久久av| 国产精品 国内视频| 最近最新中文字幕大全免费视频| 久久伊人香网站| 每晚都被弄得嗷嗷叫到高潮| 狠狠狠狠99中文字幕| 日本在线视频免费播放| 久久久久久人人人人人| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| av中文乱码字幕在线| 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 香蕉久久夜色| 美女黄网站色视频| 欧美一级毛片孕妇| 欧美日韩黄片免| 黄色视频,在线免费观看| 成人三级做爰电影| 欧美+亚洲+日韩+国产| 窝窝影院91人妻| 久久久久久九九精品二区国产 | 国产精品99久久99久久久不卡| 九色成人免费人妻av| 亚洲第一欧美日韩一区二区三区| 成人特级黄色片久久久久久久| 视频区欧美日本亚洲| 国产v大片淫在线免费观看| 国产精品99久久99久久久不卡| 国产一区在线观看成人免费| 精品国产乱码久久久久久男人| 这个男人来自地球电影免费观看| 国产97色在线日韩免费| 久久久久久久久免费视频了| 一个人免费在线观看的高清视频| 人成视频在线观看免费观看| 亚洲中文日韩欧美视频| 国产高清激情床上av| 国产欧美日韩一区二区精品| 国产乱人伦免费视频| 国产一区二区激情短视频| 人妻久久中文字幕网| 成人特级黄色片久久久久久久| 天天添夜夜摸| 91老司机精品| 亚洲专区字幕在线| 怎么达到女性高潮| 又黄又爽又免费观看的视频| 国产不卡一卡二| 精品免费久久久久久久清纯| 欧美丝袜亚洲另类 | 中文字幕人妻丝袜一区二区| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 亚洲精华国产精华精| 亚洲精品av麻豆狂野| 国产精品香港三级国产av潘金莲| 午夜老司机福利片| cao死你这个sao货| 国产精品 国内视频| 99精品在免费线老司机午夜| 一级a爱片免费观看的视频| a在线观看视频网站| 成人永久免费在线观看视频| 99在线人妻在线中文字幕| 国产69精品久久久久777片 | 哪里可以看免费的av片| 老司机福利观看| 一进一出抽搐gif免费好疼| 12—13女人毛片做爰片一| АⅤ资源中文在线天堂| 亚洲欧洲精品一区二区精品久久久| 国内少妇人妻偷人精品xxx网站 | 91大片在线观看| 日韩精品青青久久久久久| 欧美在线一区亚洲| 国产av一区在线观看免费| 亚洲精品中文字幕在线视频| 亚洲av第一区精品v没综合| 午夜精品一区二区三区免费看| 18禁黄网站禁片免费观看直播| 在线观看免费午夜福利视频| av免费在线观看网站| 国产三级黄色录像| 久久久久精品国产欧美久久久| 国产精品久久久久久亚洲av鲁大| 最近在线观看免费完整版| 九色成人免费人妻av| 亚洲成av人片免费观看| 精品少妇一区二区三区视频日本电影| 在线观看美女被高潮喷水网站 | 国产99白浆流出| 在线观看美女被高潮喷水网站 | 日韩免费av在线播放| 欧美黑人巨大hd| 女人爽到高潮嗷嗷叫在线视频| 深夜精品福利| 国产精品亚洲av一区麻豆| 三级男女做爰猛烈吃奶摸视频| 亚洲av美国av| 舔av片在线| 最近最新中文字幕大全电影3| 亚洲乱码一区二区免费版| 无限看片的www在线观看| 两个人免费观看高清视频| 久久久久久亚洲精品国产蜜桃av| 欧美一区二区精品小视频在线| 黄频高清免费视频| 亚洲av中文字字幕乱码综合| 亚洲18禁久久av| 国产精品 欧美亚洲| 欧美日本亚洲视频在线播放| 一进一出抽搐动态| 日本黄大片高清| 色播亚洲综合网| 香蕉久久夜色| 黄色视频,在线免费观看| 一级毛片女人18水好多| 黄片小视频在线播放| 色av中文字幕| 国产一区在线观看成人免费| 免费看日本二区| 日韩精品青青久久久久久| 两性夫妻黄色片| 亚洲精品美女久久久久99蜜臀| 在线观看免费午夜福利视频| 亚洲成人久久性| 麻豆av在线久日| 国产成+人综合+亚洲专区| 亚洲精品美女久久av网站| 久久伊人香网站| 国产成人系列免费观看| 国产成人欧美在线观看| 久久国产乱子伦精品免费另类| 狂野欧美白嫩少妇大欣赏| 999久久久精品免费观看国产| 亚洲国产欧美一区二区综合| 国产伦在线观看视频一区| 国产一区二区激情短视频| 国产精品一区二区免费欧美| 男女做爰动态图高潮gif福利片| 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| 男人的好看免费观看在线视频 | 19禁男女啪啪无遮挡网站| 久久精品国产99精品国产亚洲性色| 成人av一区二区三区在线看| 国产激情久久老熟女| 亚洲 国产 在线| 国产亚洲av嫩草精品影院| 国产精品免费视频内射| 一个人观看的视频www高清免费观看 | 久久国产精品人妻蜜桃| 美女免费视频网站| 俄罗斯特黄特色一大片| 国产成人精品久久二区二区免费| 欧美日韩中文字幕国产精品一区二区三区| 高潮久久久久久久久久久不卡| 国产蜜桃级精品一区二区三区| 成人亚洲精品av一区二区| 久久国产乱子伦精品免费另类| 午夜激情福利司机影院| 男人舔女人的私密视频| 热99re8久久精品国产| 欧美日韩黄片免| 嫁个100分男人电影在线观看| 久久这里只有精品中国| x7x7x7水蜜桃| 日韩精品青青久久久久久| 在线观看日韩欧美| 少妇粗大呻吟视频| 女警被强在线播放| 嫩草影视91久久| 精品欧美国产一区二区三| 亚洲成人免费电影在线观看| 禁无遮挡网站| 亚洲国产高清在线一区二区三| 在线观看午夜福利视频| 国产亚洲精品久久久久久毛片| 亚洲国产精品成人综合色| 精品国产美女av久久久久小说| 色av中文字幕| 男女做爰动态图高潮gif福利片| 色av中文字幕| 婷婷六月久久综合丁香| 级片在线观看| 久久久久久久久久黄片| 高潮久久久久久久久久久不卡| 久久久久久久久久黄片| 日韩大尺度精品在线看网址| 哪里可以看免费的av片| 精品久久久久久久人妻蜜臀av| 国产欧美日韩精品亚洲av| 50天的宝宝边吃奶边哭怎么回事| 国产高清视频在线观看网站| 日日夜夜操网爽| 亚洲国产欧洲综合997久久,| 日韩大尺度精品在线看网址| 国产三级黄色录像| 亚洲午夜理论影院| 成在线人永久免费视频| 搞女人的毛片| 久久久久久久久中文| 精华霜和精华液先用哪个| 国产97色在线日韩免费| 亚洲成人中文字幕在线播放| 国产精品久久电影中文字幕| 国产av一区二区精品久久| 一卡2卡三卡四卡精品乱码亚洲| 国产一区二区在线观看日韩 | 男女之事视频高清在线观看| 成人高潮视频无遮挡免费网站| 亚洲七黄色美女视频| 在线十欧美十亚洲十日本专区| 午夜精品一区二区三区免费看| 国产精品久久久久久精品电影| svipshipincom国产片| 久久99热这里只有精品18| 欧美中文综合在线视频| 久久99热这里只有精品18| 国产av一区在线观看免费| 他把我摸到了高潮在线观看| 黄色a级毛片大全视频| 欧美激情久久久久久爽电影| 69av精品久久久久久| 无限看片的www在线观看| 亚洲精品美女久久久久99蜜臀| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 欧美成人性av电影在线观看| 欧美最黄视频在线播放免费| 久久精品人妻少妇| av超薄肉色丝袜交足视频| 真人做人爱边吃奶动态| 少妇人妻一区二区三区视频| 俺也久久电影网| 韩国av一区二区三区四区| 国内揄拍国产精品人妻在线| 很黄的视频免费| 一边摸一边做爽爽视频免费| 黄色视频,在线免费观看| 亚洲人成网站高清观看| 亚洲国产精品999在线| 久久久国产成人免费| 国产精品99久久99久久久不卡| 最近视频中文字幕2019在线8| 91老司机精品| 国产又色又爽无遮挡免费看| 亚洲人与动物交配视频| 久久香蕉精品热| 日日干狠狠操夜夜爽| 欧美高清成人免费视频www| 久久久久久久久久黄片| 欧美另类亚洲清纯唯美| 国产三级中文精品| 久久国产乱子伦精品免费另类| 18禁黄网站禁片午夜丰满| 又粗又爽又猛毛片免费看| 国产视频内射| av国产免费在线观看| 精品国产亚洲在线| 欧美色视频一区免费| 国内少妇人妻偷人精品xxx网站 | 久久久精品大字幕| 国产精品爽爽va在线观看网站| 波多野结衣巨乳人妻| 俄罗斯特黄特色一大片| 老熟妇乱子伦视频在线观看| 禁无遮挡网站| 嫩草影视91久久| 国产成人欧美在线观看| 男人舔女人下体高潮全视频| 日韩有码中文字幕| 精品人妻1区二区| 亚洲一区高清亚洲精品| 久久这里只有精品中国| 午夜精品一区二区三区免费看| 久久久久久国产a免费观看| 搡老熟女国产l中国老女人| 欧美成人性av电影在线观看| 亚洲中文字幕日韩| 在线观看免费午夜福利视频| 亚洲成人精品中文字幕电影| 精品午夜福利视频在线观看一区| 亚洲成人国产一区在线观看| 精品久久久久久久人妻蜜臀av| 日韩欧美免费精品| 午夜福利免费观看在线| 黑人操中国人逼视频| 日本a在线网址| 国产精品免费视频内射| 别揉我奶头~嗯~啊~动态视频| 午夜福利在线在线| 亚洲国产中文字幕在线视频| 精品乱码久久久久久99久播| 男男h啪啪无遮挡| 久久人妻福利社区极品人妻图片| 黄色成人免费大全| 亚洲欧美日韩高清专用| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 国产欧美日韩一区二区三| 中文字幕久久专区| 91av网站免费观看| 亚洲专区国产一区二区| 99在线视频只有这里精品首页| 亚洲真实伦在线观看| 婷婷丁香在线五月| 男女床上黄色一级片免费看| 久久中文字幕人妻熟女| 日韩大尺度精品在线看网址| 国产欧美日韩精品亚洲av| 欧美日本视频| 久久午夜亚洲精品久久| 我的老师免费观看完整版| 亚洲成人久久爱视频| 国产高清视频在线播放一区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲人成网站在线播放欧美日韩| 免费人成视频x8x8入口观看| 中文资源天堂在线| 亚洲av日韩精品久久久久久密| 国产精品 欧美亚洲| 午夜福利高清视频| 真人一进一出gif抽搐免费| 在线观看免费视频日本深夜| 日韩大码丰满熟妇| 久久久久久免费高清国产稀缺| 听说在线观看完整版免费高清| 老熟妇乱子伦视频在线观看| 99在线视频只有这里精品首页| 国产男靠女视频免费网站| 日韩欧美在线二视频| www.999成人在线观看| 色综合婷婷激情| 亚洲aⅴ乱码一区二区在线播放 | 午夜激情福利司机影院| 草草在线视频免费看| 国产激情欧美一区二区| 搡老妇女老女人老熟妇| 999久久久精品免费观看国产| 国产视频内射| 18禁国产床啪视频网站| 欧美日本视频| 99riav亚洲国产免费| 亚洲 国产 在线| 亚洲avbb在线观看| 亚洲精品粉嫩美女一区| 日本熟妇午夜| 亚洲欧美日韩高清专用| 18禁黄网站禁片免费观看直播| 国产精品一区二区免费欧美| 我要搜黄色片| 国产午夜福利久久久久久| 国产区一区二久久| 我要搜黄色片| 动漫黄色视频在线观看| 国产69精品久久久久777片 | 午夜精品在线福利| 亚洲精品一区av在线观看| 亚洲人成伊人成综合网2020| 亚洲最大成人中文| 一本一本综合久久| 精品久久久久久,| 午夜久久久久精精品| 视频区欧美日本亚洲| 午夜激情av网站| 亚洲人成电影免费在线| 在线观看一区二区三区| 一级毛片精品| 最近最新中文字幕大全免费视频| 99久久久亚洲精品蜜臀av| 国产爱豆传媒在线观看 | 午夜a级毛片| 国产av麻豆久久久久久久| 亚洲aⅴ乱码一区二区在线播放 | 人妻久久中文字幕网| 亚洲av成人av| 岛国在线观看网站| 久99久视频精品免费| 久久亚洲精品不卡| 国产精品野战在线观看| 亚洲精华国产精华精| 波多野结衣高清无吗| 伦理电影免费视频| 久久中文字幕人妻熟女| 在线观看66精品国产| 精品熟女少妇八av免费久了| 亚洲国产看品久久| 国产精品一区二区免费欧美| 最近视频中文字幕2019在线8| 亚洲av五月六月丁香网| 久9热在线精品视频| 国产亚洲精品第一综合不卡| 法律面前人人平等表现在哪些方面| 一二三四社区在线视频社区8| 亚洲av片天天在线观看| 亚洲五月天丁香| 午夜视频精品福利| 亚洲精品国产一区二区精华液| 一本久久中文字幕| 中文字幕精品亚洲无线码一区| 午夜福利欧美成人| 精品久久久久久久久久免费视频| 俄罗斯特黄特色一大片| 黑人巨大精品欧美一区二区mp4| 别揉我奶头~嗯~啊~动态视频| 啦啦啦免费观看视频1| 久久热在线av| 亚洲熟妇中文字幕五十中出| 黑人操中国人逼视频| www.www免费av| 97碰自拍视频| 天堂影院成人在线观看| 国产精品日韩av在线免费观看| 成人国产综合亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦免费观看视频1| 国产成人av激情在线播放| 免费看美女性在线毛片视频| 亚洲中文字幕日韩| e午夜精品久久久久久久| 婷婷六月久久综合丁香| 亚洲av中文字字幕乱码综合| 一个人观看的视频www高清免费观看 | www.精华液| 国产一区二区在线av高清观看| 老司机深夜福利视频在线观看| 毛片女人毛片| 2021天堂中文幕一二区在线观| 成人18禁在线播放| 露出奶头的视频| 免费无遮挡裸体视频| 国产熟女午夜一区二区三区| 亚洲欧美日韩无卡精品| 国产又黄又爽又无遮挡在线| 50天的宝宝边吃奶边哭怎么回事| 90打野战视频偷拍视频| 国产精品免费视频内射| 99热6这里只有精品| 国产高清激情床上av| 一区二区三区国产精品乱码| 亚洲国产精品成人综合色| 男人舔奶头视频| 欧美中文日本在线观看视频| 日本黄色视频三级网站网址| 精品久久蜜臀av无| 美女黄网站色视频| 一本综合久久免费| 一本久久中文字幕| 巨乳人妻的诱惑在线观看| 亚洲欧美日韩高清在线视频| av视频在线观看入口| 老熟妇仑乱视频hdxx| 成人三级做爰电影| 女生性感内裤真人,穿戴方法视频| 在线观看66精品国产|