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

    視覺顯著區(qū)域和主動輪廓結(jié)合的圖像分割算法

    2022-05-17 06:02:08何亞茹葛洪偉
    計算機與生活 2022年5期
    關(guān)鍵詞:測度輪廓灰度

    何亞茹,葛洪偉+

    1.江南大學(xué) 江蘇省模式識別與計算機智能工程實驗室,江蘇 無錫214122

    2.江南大學(xué) 人工智能與計算機學(xué)院,江蘇 無錫214122

    圖像分割是實時圖像處理的基礎(chǔ)研究內(nèi)容,作為目標(biāo)提取、識別和跟蹤的預(yù)處理操作,分割結(jié)果的好壞直接影響后續(xù)步驟的判斷。圖像采集過程中,易受設(shè)備及外部環(huán)境的影響,獲取的圖像往往存在灰度分布不均且邊界模糊的問題,圖像分割面臨巨大的挑戰(zhàn)。Kass等人提出主動輪廓模型(active contour model,ACM)并首次應(yīng)用于圖像分割中,利用能量最小化原理,為目標(biāo)分割提供平滑的輪廓。由于此類模型在處理具有亞像素邊緣的圖像中表現(xiàn)突出,憑借允許嵌入各種先驗知識實現(xiàn)魯棒分割的優(yōu)勢,在計算機視覺和圖像分割中得到了廣泛的應(yīng)用。

    目前主要存在基于邊緣和區(qū)域的主動輪廓模型。前者主要依賴圖像梯度信息定義邊緣停止函數(shù),在梯度大的位置曲線停止演化,捕捉目標(biāo)邊界。測地主動輪廓(geodesic active contour,GAC)模型為具有代表性的基于邊界主動輪廓模型,該模型在高質(zhì)量圖像中取得精確的分割結(jié)果,但對噪聲敏感且不適合弱邊緣的提取。后者主要利用區(qū)域統(tǒng)計信息構(gòu)造驅(qū)動力,引導(dǎo)曲線演化至目標(biāo)邊界。根據(jù)區(qū)域統(tǒng)計信息類型的不同分為基于全局、局部的主動輪廓模型?;谌值闹鲃虞喞P椭芯哂写硇缘臑镃V(Chan-Vese)模型,該模型基于圖像由兩個同質(zhì)區(qū)域組成的假設(shè),能夠有效分割同質(zhì)圖像。Zhang等人在CV 模型的基礎(chǔ)上提出符號壓力函數(shù)(signed pressure function,SPF),該函數(shù)能夠有效阻止輪廓曲線在弱邊緣和模糊邊緣處的運動,檢測目標(biāo)邊界。然而,上述兩種模型,難以準(zhǔn)確捕捉灰度不均圖像的目標(biāo)邊界。為了解決這一問題,Li 等人引入高斯核函數(shù)來統(tǒng)計圖像局部信息,提出RSF(region-scalable fitting)模型,從而克服全局模型無法分割灰度呈非均勻變化圖像的缺點,但該模型對初始化的設(shè)置極其敏感,易出現(xiàn)錯誤的分割結(jié)果。鑒于對基于全局和局部兩種模型優(yōu)缺點的分析,許多研究學(xué)者嘗試將兩種模型結(jié)合,提出一系列混合模型。Wang 等人引入全局能量項,吸引演化曲線向目標(biāo)附近靠攏。同時,利用局部能量項進一步驅(qū)動演化曲線到達目標(biāo)邊界,該模型在灰度分布較為簡單的圖像中表現(xiàn)較好,但由于模型中的全局項無法在灰度呈分布不均勻圖像中有效識別目標(biāo)邊界,因此即便調(diào)節(jié)全局項權(quán)重系數(shù)也難以得到滿意的結(jié)果。趙杰等人通過設(shè)計雙核函數(shù)來提取局部信息,并使用改進的全局項加以輔助,提出的混合區(qū)域主動輪廓模型從模糊C 均值聚類算法(fuzzy C-means algorithm,F(xiàn)CM)中獲取目標(biāo)先驗信息,自動設(shè)置初始化輪廓。Ding 等人提出優(yōu)化LoG(Laplacian of Gaussian)能量項,作為一種全局檢測方法,與RSF 模型相結(jié)合,降低模型陷入局部極值的可能。然而,此模型仍需人工設(shè)置初始化輪廓,當(dāng)輪廓曲線遠離目標(biāo)邊界時,演化曲線易受背景干擾,收斂速度緩慢,且局部項的設(shè)計僅考慮局部空間信息,無法全面衡量局部鄰域點對中心點的影響程度,導(dǎo)致模型在待分割目標(biāo)弱邊緣處易發(fā)生泄漏,降低模型的分割精度。

    復(fù)雜場景中,僅使用主動輪廓模型往往難以得到理想的分割結(jié)果。陶永鵬等人通過提取CT 圖像的SIFT(scale invariant feature transform)特征,再利用加權(quán)隨機森林回歸算法獲取椎骨中心點位置,將其作為CV 分割模型的初始輪廓位置。該方法在自動定位方面具有一定的優(yōu)勢,但無法獲取待分割目標(biāo)形狀先驗信息。Li 等人提出了一種將模糊C 均值聚類與水平集方法相結(jié)合的算法。其中,水平集函數(shù)的初始化直接從模糊C 均值聚類中獲取,并且從模糊聚類的結(jié)果中估計水平集演化的參數(shù),減少了人工操作。Yang 等人使用類似的思想,將模糊聚類的結(jié)果作為初始約束信息,通過動態(tài)約束項將水平集方法與模糊聚類相結(jié)合,提出新的能量泛函。但由于聚類算法在復(fù)雜自然圖像中表現(xiàn)欠佳,導(dǎo)致主動輪廓模型獲取的初始值具有不準(zhǔn)確性,從而影響模型的分割結(jié)果。視覺顯著區(qū)域檢測算法是當(dāng)前計算機視覺領(lǐng)域的一大研究熱點。Fu 等人提出基于簇的視覺顯著性檢測算法,通過定義各簇對應(yīng)的對比度測度和空間測度,來衡量各簇的顯著值,并以非線性方式將兩個測度融合,得到最終顯著圖。該模型簡單、高效,獲取的目標(biāo)邊界清晰,能夠在灰度呈不均勻分布的醫(yī)學(xué)圖像或自然圖像中較好地提取視覺顯著區(qū)域。Yuan 等人結(jié)合回歸校正和正則化隨機游走排序算法,輸出更平滑、準(zhǔn)確的顯著圖。但此類模型僅能定位目標(biāo)顯著區(qū)域,難以將目標(biāo)完全分割出來。

    針對主動輪廓模型用于圖像分割時,演化曲線受背景干擾,演化速度緩慢,且易在弱邊緣處發(fā)生泄漏,難以保留目標(biāo)完整性的問題,本文根據(jù)每個簇的緊湊性定義緊湊對比度測度,與空間測度相結(jié)合,利用改進顯著性檢測算法對原始圖像進行預(yù)處理操作,提取目標(biāo)候選區(qū)域,自動設(shè)置初始化輪廓,并將獲取的目標(biāo)先驗信息用于優(yōu)化LoG 能量項的符號選取,自適應(yīng)決定曲線演化方向;其次,引入局部窗口內(nèi)鄰域灰度信息的變化,設(shè)計局部灰度測度,與局部核函數(shù)相結(jié)合,改進局部能量項,更為全面地衡量局部窗口內(nèi)鄰域像素點對中心點的影響權(quán)重,提高模型在目標(biāo)弱邊緣處的敏感程度。此外,本文模型增加長度項和距離正則項,用于平滑輪廓曲線,避免重新初始化。該模型無需人工設(shè)置初始化輪廓,并且提升了在目標(biāo)弱邊緣處的捕獲能力,取得了更為理想的分割效果。

    1 主動輪廓模型簡介

    1.1 RSF 模型

    Li 等人通過定義局部擬合項用于分割灰度分布不均勻圖像,提出局部尺度可控的RSF 模型。對于圖像域內(nèi)任意像素點,定義基于水平集函數(shù)的局部擬合項,由式(1)給出:

    式中,f()為像素點處的局部外(內(nèi))擬合均值;、為權(quán)重系數(shù),均為正數(shù);δ為高斯核函數(shù),控制局部窗口內(nèi)鄰域像素點處灰度值()對中心像素點處的影響程度;(·)為亥維賽函數(shù)。

    RSF 模型能夠有效分割灰度呈非均勻變化的圖像,但由于該模型的局部化特性,導(dǎo)致其對初始化輪廓的設(shè)置敏感,曲線演化易陷入局部極值,無法準(zhǔn)確定位目標(biāo)邊界,極大程度上限制了該模型的應(yīng)用。

    1.2 ACML 模型

    傳統(tǒng)LoG 算子對灰度變化十分敏感,在含噪圖像中,零交叉點易被判為錯誤邊緣點,構(gòu)建的能量項用于分割時,易出現(xiàn)“孤島”區(qū)域。為了解決此問題,Ding 等人提出基于圖像的優(yōu)化LoG 能量項,用于平滑同質(zhì)區(qū)域,保持邊緣信息,與僅由高斯核函數(shù)定義的局部能量項相結(jié)合,構(gòu)建混合模型,本文稱其為ACML(active contour model based on optimized LoG and RSF model)模型,降低演化曲線陷入局部極值的可能性。

    首先,使用變量來表示優(yōu)化后的LoG 算子,定義能量泛函,如式(2)所示:

    其次,假設(shè)待分割目標(biāo)灰度均值大于背景灰度均值,并使用Ω表示演化曲線包圍的內(nèi)部區(qū)域,并用優(yōu)化后的變量代替文獻[20]中傳統(tǒng)LoG 算子,構(gòu)建能量泛函:

    最后,使用θ表示全局項系數(shù),結(jié)合式(1)、式(3),定義外部能量約束項:

    該模型一定程度上避免曲線陷入局部極值。然而,該模型仍需人工設(shè)置初始化輪廓曲線,當(dāng)初始化輪廓設(shè)置遠離真實目標(biāo)邊界時,曲線容易受背景干擾,演化速度緩慢,且由于局部能量項僅考慮空間信息,使ACML 模型易在弱邊緣處發(fā)生泄漏,影響模型分割精度和效率。

    2 基于視覺顯著區(qū)域的初始輪廓提取算法

    本文結(jié)合視覺顯著性檢測方法獲取目標(biāo)形狀先驗信息,作為形狀約束引入主動輪廓模型中,提出一種基于視覺顯著區(qū)域的初始輪廓提取算法(initial contour based on visual salient regions,ISR),自動設(shè)置初始化輪廓曲線,消除人工干預(yù),初始輪廓曲線于目標(biāo)邊界附近開始演化,避免背景干擾,提高模型分割的精度和效率。

    2.1 基于簇的視覺顯著區(qū)域提取

    本文采用基于簇的視覺顯著性檢測算法獲取目標(biāo)先驗信息,并根據(jù)各簇的緊湊性對該算法進行改進。首先,使用聚類算法將圖像中像素點聚類為個簇,通過統(tǒng)計簇內(nèi)各像素點到圖像中心的空間距離,并比較簇C與其他簇(非簇C)之間的差異性,同時考慮每個簇的空間分布,分別定義空間測度ω(C)和緊湊對比度測度ω(C)。其次,以非線性方式融合兩個測度,獲取各簇對應(yīng)的顯著值。而后根據(jù)每個簇的顯著值得到簇中各像素點對應(yīng)的顯著值。

    分別定義空間測度ω(C) 和對比度測度ω(C)來表示各簇到圖像中心的空間距離以及各簇之間的差異性,如式(5)、式(6)所示:

    為了更好地抑制背景干擾,獲取圖像顯著區(qū)域,文獻[21]指出顯著目標(biāo)通常具有緊湊的空間分布,而背景區(qū)域較分散。本文將此概念延伸到簇,利用簇內(nèi)的各像素點到該簇中心位置的空間距離定義緊湊對比度測度,有效抑制背景干擾。首先計算圖像域中基于簇的對比度測度顯著值均值,如下所示:

    其中,表示聚類簇數(shù);使用來表示簇內(nèi)各像素點到該簇的中心位置的空間距離,并將其歸一化為[0,1]。定義如下:

    式中,<,本文所有實驗=0,=1。

    采取非線性方式融合兩種測度,獲得更為精確的視覺顯著區(qū)域,如式(10)所示:

    式中,(C)表示簇C顯著測度值,每個像素點屬于簇C的顯著性似然滿足高斯分布,定義如下:

    式中,v表示像素的特征向量;μ表示簇C的聚類中心;σ表示標(biāo)準(zhǔn)差。最后獲取每個像素點對應(yīng)的顯著值,如式(12)所示:

    圖1(a)中給出一幅X-Ray 血管原始圖像,圖1(b)、(c)給出對比度測度與緊湊對比度測度圖,再分別與空間測度以非線性方式相結(jié)合,獲取顯著圖與緊湊顯著圖,圖1(e)、(f)分別給出其對應(yīng)的二值圖與緊湊二值圖,再進一步利用形態(tài)學(xué)操作濾除孤立區(qū)域,獲取更為準(zhǔn)確的目標(biāo)形狀先驗信息。不難看出,使用改進后的對比度測度獲取的目標(biāo)先驗信息更加準(zhǔn)確,依此設(shè)置的初始化輪廓曲線能更好定位目標(biāo)位置,如圖1(d)中綠色曲線所示。

    圖1 對比度測度改進前后效果Fig.1 Effect of contrast measure before and after improvement

    2.2 初始輪廓提取算法

    本文結(jié)合改進視覺顯著性檢測算法,獲取目標(biāo)形狀先驗信息,提出基于視覺顯著區(qū)域的初始化輪廓提取算法ISR,步驟如下:

    1.使用聚類算法將圖像聚類為個簇;

    2.for←1 todo

    根據(jù)式(5)計算簇的空間測度ω(C);

    根據(jù)式(6)~(9)計算簇C的緊湊對比度測度ω(C);

    end

    3.根據(jù)式(10)融合兩個測度,得到各簇對應(yīng)的顯著值;

    4.for←1 todo

    根據(jù)式(12)獲得圖像每個像素點對應(yīng)的顯著值;

    end

    5.根據(jù)步驟4 獲取待分割圖像的視覺顯著區(qū)域,再利用閾值將獲取的顯著圖分割為二值化圖像,當(dāng)部分目標(biāo)錯分為背景或部分背景錯分為目標(biāo)時,構(gòu)成圖像的兩部分(目標(biāo)和背景)灰度差別越小,反之越大。使用最大化類間方差法來計算最優(yōu)閾值。將顯著值大于閾值的,標(biāo)簽記為1;將顯著值小于此閾值的,標(biāo)簽記為0,得到二值圖像,對二值化圖像進行形態(tài)學(xué)操作,濾除噪聲,將處理后得到的最終二值化圖像分界線作為初始演化輪廓曲線。

    視覺顯著區(qū)域的有效提取,能夠提供準(zhǔn)確的目標(biāo)先驗信息,為后續(xù)操作步驟提供可靠依據(jù)。本文在提取視覺顯著區(qū)域特征時,不局限于某一特定的聚類算法,通過實驗測試并分析means++、FCM、均值漂移(Mean shift)算法三種較為經(jīng)典的聚類算法在提取視覺顯著區(qū)域時的優(yōu)劣性。

    圖2(a)~(c)分別給出嵌入means++、FCM、Mean shift聚類算法的視覺顯著區(qū)域檢測算法所獲取的初始化輪廓(綠色曲線)。前兩種聚類算法均需人工給出聚類簇數(shù),means++算法高效率,易實現(xiàn),且在改進對比度測度的作用下,降低了對聚類簇數(shù)的敏感程度。由圖2(a)可以看出,在醫(yī)學(xué)和自然圖像中均表現(xiàn)較好。FCM 聚類算法以隸屬度來衡量每個數(shù)據(jù)點屬于各簇的可能性,在背景分布相對簡單的醫(yī)學(xué)圖像中獲取較為精確的視覺顯著區(qū)域,而在背景相對復(fù)雜且與目標(biāo)強度值相近的自然圖像中,難以將像素點準(zhǔn)確歸類,如該算法在小船圖像中的分割結(jié)果所示。Mean shift聚類算法雖無需人工給出聚類簇數(shù),但引入了帶寬參數(shù),當(dāng)帶寬大小設(shè)置不合理,收斂速度慢,無法準(zhǔn)確地提取視覺顯著區(qū)域。綜合上述分析,本文實驗均選取means++算法來獲取視覺顯著區(qū)域,其中簇數(shù)為5。

    圖2 使用不同聚類算法得到的初始化輪廓Fig.2 Initial contour obtained by different clustering algorithms

    3 基于混合區(qū)域的主動輪廓分割模型

    本章首先引入優(yōu)化變量來代替?zhèn)鹘y(tǒng)LoG 算子,構(gòu)建基于圖像全局信息的優(yōu)化LoG 能量項,改善傳統(tǒng)LoG 算子對灰度變化敏感的問題,并結(jié)合由顯著性檢測獲取的目標(biāo)先驗信息,設(shè)計自適應(yīng)符號函數(shù),對優(yōu)化LoG 能量項進行加權(quán),自行追蹤目標(biāo)邊緣方向;其次,設(shè)計局部灰度測度,與高斯核函數(shù)相結(jié)合,改進局部能量項,更加充分地衡量局部鄰域像素點對中心點的影響權(quán)重,提高模型對目標(biāo)弱邊緣的敏感程度。此外,引入長度項與距離正則項,用于平滑輪廓,避免水平集重新初始化,提出新的基于區(qū)域的混合主動輪廓模型(region based hybrid active contour model,RBHM)。

    3.1 自適應(yīng)符號函數(shù)加權(quán)的優(yōu)化LoG 能量項

    假設(shè)圖像由兩部分區(qū)域組成,且待分割目標(biāo)灰度均值大于背景灰度均值??紤]到傳統(tǒng)LoG 算子對灰度變化敏感問題,同文獻[12],本文使用優(yōu)化LoG算子來代替Δ,用于平滑同質(zhì)區(qū)域,保持目標(biāo)邊緣,使用表示圖像域,定義如下能量泛函:

    其中,為權(quán)重系數(shù),(·)為亥維賽函數(shù)。此外,本文利用目標(biāo)先驗信息判斷待分割目標(biāo)相對于圖像的亮暗程度,設(shè)計自適應(yīng)符號函數(shù),對優(yōu)化LoG 能量項進行加權(quán),自行判定目標(biāo)邊緣追蹤方向,增強模型自適應(yīng)能力。

    首先,利用改進的視覺顯著性檢測算法獲取目標(biāo)候選區(qū)域,使用大津法對顯著圖進行二值化操作,并利用圓盤半徑為2 的模板對獲得的二值化顯著圖依次進行腐蝕、膨脹操作。類似于文獻[24-25],選取腐蝕后的二值掩碼圖中最大連通域內(nèi)的像素點判為目標(biāo)種子點,標(biāo)簽記為1,其余為非目標(biāo)種子點,標(biāo)簽記為0,得到目標(biāo)標(biāo)簽矩陣。其次,對于彩色圖像,同文獻[26],選取原始圖像(R,G,B)通道中具有最大對比度的位圖Img(:,:,cmax),與標(biāo)簽矩陣以非線性方式結(jié)合,獲取具有最大對比度位圖中目標(biāo)種子強度均值與非目標(biāo)種子的強度均值,設(shè)計自適應(yīng)符號函數(shù)sgn(,)。具體操作步驟如下:

    (1)使用(,)表示使用視覺顯著性檢測算法獲取的顯著圖中像素點,標(biāo)簽矩陣表示二值化顯著圖經(jīng)過腐蝕操作后具有最大連通域的二值掩碼圖,定義如下:

    (2)同文獻[26],選取原始分割圖像R、G、B 通道中具有最大對比度的位圖Img(:,:,cmax),并使用(,)表示位圖Img(:,:,cmax)中任意像素點的灰度值,分別計算目標(biāo)種子點與非目標(biāo)種子點對應(yīng)的灰度均值,如下所示:

    (3)比較目標(biāo)種子點與非目標(biāo)種子點的均值大小,設(shè)計自適應(yīng)符號函數(shù),定義如下:

    結(jié)合式(13)、式(17),本文提出自適應(yīng)符號函數(shù)的優(yōu)化LoG 能量項,定義如下:

    圖3(a)中分別給出含有亮、暗目標(biāo)的兩幅合成圖像,圖3(b)分別為其對應(yīng)的顯著圖,圖3(c)中黃色封閉曲線所包圍的像素點為圖像經(jīng)預(yù)處理操作分別獲取的目標(biāo)種子,依此設(shè)計自適應(yīng)符號函數(shù),對LoG能量項進行加權(quán),自行決定曲線演化方向。

    圖3 目標(biāo)種子點Fig.3 Seeds of target

    不同于基于圖像全局信息的CV 模型,自適應(yīng)符號加權(quán)的優(yōu)化LoG 能量項作為一種全局檢測方法,能夠在灰度分布不均勻圖像中有效識別與背景強度值相差較大的目標(biāo)邊界,自適應(yīng)決定輪廓曲線的演化方向,避免模型陷入局部極值。然而,僅使用該模型容易出現(xiàn)過分割現(xiàn)象。

    3.2 改進局部能量項

    本文通過設(shè)計局部灰度測度δ,與局部核函數(shù)δ相結(jié)合,充分考慮局部窗口內(nèi)鄰域像素點對中心點的影響權(quán)重,構(gòu)建新的局部能量項。

    其次,設(shè)計局部灰度測度統(tǒng)計局部窗口內(nèi)中心像素點處的灰度信息變化,當(dāng)局部窗口位于同質(zhì)或弱邊緣區(qū)域,灰度變化較小,應(yīng)減弱空域中局部鄰域點對中心點的影響;當(dāng)局部窗口鄰域點為邊緣點或噪聲點時,局部窗口灰度變化較大,應(yīng)增強空域中鄰域點對中心點的作用。以此提高模型在弱邊緣處的敏感程度,防止錯誤邊緣點出現(xiàn)。

    對于圖像中的任一局部窗口N,使用ω表示中心像素點所在局部窗口的灰度分布均勻程度,定義如式(20):

    本文采取負指數(shù)函數(shù)來設(shè)計局部窗口內(nèi)中心像素點處的局部灰度測度,定義如下:

    不難看出,當(dāng)像素點對應(yīng)的ω值較小時,δ取得較大值;反之,亦然。結(jié)合局部灰度變化信息與空間信息改進局部能量項,重新定義式(1):

    式中,、均為正常數(shù);、表示局部擬合強度均值;()為亥維賽函數(shù),使用正則化函數(shù)H()代替計算,其定義如式(23)所示。

    改進的局部能量項在局部灰度統(tǒng)計信息與空間信息共同作用下,能夠有效提取目標(biāo)弱邊緣,極大程度上保持了目標(biāo)輪廓的連貫性與完整性,但對初始化輪廓設(shè)置仍然較敏感,且曲線演化易陷入局部極值,無法準(zhǔn)確定位目標(biāo)邊界。

    3.3 RBHM 模型的提出

    鑒于對上述兩種能量項的優(yōu)缺點分析,本文以線性方式將兩者結(jié)合起來,提出的RBHM 模型在分割過程中能夠有效保持目標(biāo)的完整性,具有較強的自適應(yīng)能力,能量泛函定義如下:

    采用同文獻[7]中對RSF模型的實現(xiàn)方式,局部均值擬合函數(shù)()、()分別通過以下卷積運算獲?。?/p>

    固定()、(),能量泛函(25)關(guān)于水平集函數(shù)變分,由梯度下降流得到演化曲線,如下所示:

    式中,=1,2,?表示梯度算子;div(·)為散度算子。使用ISR 算法獲取的初始化輪廓能夠有效識別視覺顯著區(qū)域,并將提取的目標(biāo)先驗信息用于初始化輪廓的設(shè)置與優(yōu)化LoG 能量項符號的設(shè)計,消除人工干預(yù),增強模型自適應(yīng)能力,而后使用改進的區(qū)域擬合能量模型RBHM 引導(dǎo)演化曲線進一步分割。圖4(a)給出一幅共聚焦顯微鏡圖像和一幅腦血管造影圖像及其初始化,圖4(b)給出僅使用高斯核函數(shù)定義的局部能量項分割圖4(a)中兩幅圖像的結(jié)果,圖4(c)為使用局部灰度測度與高斯核函數(shù)改進的局部能量項的分割結(jié)果,圖4(d)為僅使用自適應(yīng)符號函數(shù)加權(quán)優(yōu)化LoG 全局能量項對應(yīng)的分割結(jié)果,圖4(e)為本文提出的RBHM 模型分割結(jié)果。其中,綠色曲線表示使用ISR 算法獲取的初始化輪廓,紅色曲線表示分割結(jié)果。高斯核標(biāo)準(zhǔn)差=3.3,在分割圖4(b)~(e)第一行圖像時,長度項系數(shù)依次設(shè)置為0.000 5 v、0.000 3 v、0.000 1 v、0.000 3 v,在分割圖4(b)~(e)第二行圖像時,長度項系數(shù)依次設(shè)置為0.000 6 v、0.000 6 v、0.000 1 v、0.000 6 v,其中v=255×255。

    圖4 RBHM 模型的有效性Fig.4 Effectiveness of RBHM model

    改進的局部能量項由于充分考慮了局部灰度域與空域信息的變化,提高了RBHM 模型在弱邊緣處的敏感程度。由圖4(c)第一幅圖像不難看出,改進的局部能量項在圖中弱邊緣(黃色框)處的分割結(jié)果更加連貫,然而僅使用改進的局部能量項在非邊緣目標(biāo)邊界易發(fā)生粘連,如該圖中(藍色框處)所示。且在目標(biāo)分布范圍較廣的弱邊緣圖像中,易陷入局部極值,如圖4(c)中第二幅圖像所示。自適應(yīng)符號函數(shù)加權(quán)的優(yōu)化LoG 能量項由于引入目標(biāo)先驗信息,能夠自動追蹤目標(biāo)邊界方向,且作為一種全局檢測方法,該能量項不易陷入局部極值??紤]到上述兩種能量項的優(yōu)缺點,以線性方式將其結(jié)合起來,一方面保持目標(biāo)輪廓的連貫性和完整性,另一方面降低模型陷入局部極值的可能性,更好地定位目標(biāo)邊界,得到了較為滿意的分割結(jié)果,如圖4(e)所示。

    4 基于ISR算法和RBHM模型的圖像分割算法

    本文首先根據(jù)ISR 算法自動獲取初始輪廓曲線,而后使用本文提出的RBHM 模型引導(dǎo)曲線繼續(xù)演化,直至滿足迭代停止條件。

    本文模型的算法流程圖如圖5 所示,具體實現(xiàn)操作步驟如下:

    圖5 本文模型流程圖Fig.5 Flow chart of proposed model

    (1)演化曲線常規(guī)參數(shù)設(shè)置:=1,=3.0,局部窗口半徑=round(2),=2,===1,=0.8,Δ=0.02,=0.01,=3,長度項系數(shù)需根據(jù)圖像進行調(diào)整。

    (2)使用ISR 算法獲取圖像的初始輪廓,構(gòu)造水平集函數(shù),即:

    其中,=C。

    (3)根據(jù)式(14)~(17)設(shè)計自適應(yīng)符號函數(shù)sgn(,),對優(yōu)化LoG 能量項進行加權(quán),使用式(18)計算自適應(yīng)符號函數(shù)加權(quán)的優(yōu)化LoG 能量項。

    (4)根據(jù)式(20)~(21)計算局部灰度測度,利用式(22)定義改進局部能量項,再根據(jù)式(25)構(gòu)建新的能量泛函。

    (5)根據(jù)式(26)~(27)計算局部擬合均值、。

    (6)固定、,根據(jù)式(28)更新水平集演化方程。

    (7)每迭代6 次,計算演化輪廓曲線所圍面積的變化,當(dāng)滿足|S-S|<10或迭代次數(shù)>500,則停止迭代并輸出分割結(jié)果。

    5 實驗結(jié)果及分析

    本文實驗均在Intel Core i5-4210U 1.70 Hz CPU,4.0 GB RAM 和Windows 7(64 位)操作系統(tǒng)的個人計算機上完成。編程環(huán)境Matlab 2015a。實驗中所使用的示例圖像選自數(shù)據(jù)庫以及四幅真實示例圖像。SPF 模型中的常值速度以及RSF、LGIF、ACML、LPF中的長度項系數(shù)~均需根據(jù)圖像調(diào)整,LPF、ACML 模型中時間步長Δ=0.02 。除特殊說明外,其余參數(shù)均參考原文設(shè)置。

    5.1 ISR 算法有效性分析

    為了評估自動獲取初始化輪廓算法ISR 的有效性,圖6 給出六幅真實示例圖像,前兩幅均為含有弱邊緣且灰度分布不均勻的血管圖像;第三幅為對比度較低的輸卵管造影圖像;第四幅為灰度分布極其不均勻且含有弱邊緣的莖葉圖像;第五幅為一幅含有多目標(biāo)圖像,第六幅為含有極弱邊緣的貝殼圖像。圖6、圖7 分別給出六種模型在人工設(shè)置的初始化輪廓與使用ISR 算法獲取的初始化輪廓下的分割結(jié)果。圖6、圖7 中~表示RSF、LGIF、ACML、LPF 模型中的高斯核標(biāo)準(zhǔn)差,、θ分別表示LGIF模型全局項系數(shù)與ACML 模型全局項系數(shù)。圖6 依次給出人工設(shè)置的初始化輪廓曲線(綠色框)及SPF、RSF、LGIF、ACML、LPF 與RBHM 模型分割結(jié)果(紅色曲線)。表1 為分割圖6、圖7 時模型中所使用的常值速度與長度項系數(shù)值。

    表1 分割圖6、圖7 示例圖像時各模型中v、v1~v5 的取值Table 1 Value of v,v1 to v5 in each model segmenting sample images on Fig.6 and Fig.7

    圖6 人工初始化分割結(jié)果Fig.6 Segmentation results of manual initialization

    圖7 自動初始化分割結(jié)果Fig.7 Segmentation results of automatic initialization

    由圖6 所示分割結(jié)果可以看出,使用人工設(shè)置的矩形框初始化輪廓做引導(dǎo),由于SPF 模型僅考慮圖像全局信息,在處理灰度分布不均勻的弱邊緣圖像中,整體效果欠佳;由于RSF 模型僅考慮局部空間變化信息,對初始化的設(shè)置十分敏感,在圖像4 中出現(xiàn)錯誤分割結(jié)果,且難以將圖像中的多個目標(biāo)分割出來,如該模型在圖像5 中分割結(jié)果所示;LGIF 模型通過引入全局項對RSF 模型進行改進,在圖像5 中取得較優(yōu)的結(jié)果,但仍泄露了部分目標(biāo)弱邊緣,且在分割圖像3~4 時,受背景干擾,該模型仍容易陷入局部極值,這是由于該模型中的全局項在灰度不均勻圖像中處理能力弱,即便增大全局項系數(shù),也難以準(zhǔn)確獲取目標(biāo)邊界;LPF 模型在預(yù)擬合均值函數(shù)的作用下,降低了對初始化輪廓設(shè)置的敏感程度。然而,該模型在分割圖像4 時,由于背景干擾,出現(xiàn)“孤島”區(qū)域;ACML 模型與RBHM 模型均使用優(yōu)化LoG 能量項作為全局項,能夠在灰度不均勻圖像中準(zhǔn)確識別目標(biāo)邊界,如在圖像3~4 中的分割結(jié)果所示。而本文提出的RBHM 模型中的局部能量項充分考慮局部空間和灰度域信息,提高模型在弱邊緣處敏感程度,因此在圖像6 的極弱邊緣處,分割效果明顯優(yōu)于ACML模型。

    圖7 給出使用ISR 算法自動獲取的初始化輪廓(綠色曲線),分別與SPF、RSF、LGIF、ACML、LPF 以及本文提出的RBHM 模型相結(jié)合,依次標(biāo)記為SPF*、RSF*、LGIF*、ACML*、LPF*以及Our(RBHM*)模型。

    由圖7 各模型的分割結(jié)果(紅色曲線)易看出,使用ISR 算法獲取的初始化輪廓,較人工初始輪廓,各模型分割效果大都得到改善。SPF*模型在某些情況下能夠得到部分精確結(jié)果,如在圖像4~6 中分割結(jié)果所示;RSF*模型在灰度不均的圖像3~4 以及含有多目標(biāo)的圖像5 中均取得了較為滿意的分割結(jié)果;LGIF*模型在灰度分布不均的圖像3 與圖像4 上均得到更為準(zhǔn)確的分割結(jié)果;LPF*模型在圖像3~4、6 上的分割效果也得到改善。ACML*模型與本文(RBHM*)模型不僅能夠在圖像1~4 上分割出目標(biāo),在圖像5 中,兩種模型均取得更優(yōu)結(jié)果。表2 給出六種模型在人工設(shè)置的初始化輪廓與使用ISR 算法獲取的初始化輪廓下分割所需迭代時間。

    結(jié)合表2 和圖6、圖7 的分割結(jié)果,不難看出,利用ISR 算法設(shè)置的初始化輪廓曲線,能夠準(zhǔn)確獲取目標(biāo)先驗信息,節(jié)省了整體迭代所需時間,提高了各模型分割效率。

    表2 比較六種模型使用不同初始化輪廓所需的迭代時間Table 2 Iteration time required for six models using different initialized contours s

    5.2 RBHM 模型的有效性分析

    為了進一步驗證本文模型分割結(jié)果的精確度,引入標(biāo)準(zhǔn)數(shù)據(jù)庫。表3 給出模型分割圖8 圖像中所使用的常值速度、長度項系數(shù)。圖8 給出使用ISR 算法自動獲取的初始化輪廓(綠色曲線),再依次使用SPF、RSF、LGIF、ACML、LPF模型及本文提出的RBHM 模型對圖像進一步分割。

    表3 分割圖8 示例圖像時各模型中v、v1~v5 的取值Table 3 Value of v,v1 to v5 in each model segmenting sample images on Fig.8

    由圖8 給出的分割結(jié)果(紅色曲線)可以看出,由于SPF*模型僅考慮全局信息,無法準(zhǔn)確貼合目標(biāo)真實邊緣,如在圖像I 中的分割結(jié)果所示;RSF*模型、LGIF*模型以及ACML*模型均僅利用局部空間信息來定義局部能量項,導(dǎo)致模型在飛機機翼與船桅桿的弱邊緣處發(fā)生泄漏,且無法提取與背景灰度相似的目標(biāo)區(qū)域,如在圖像A~D、I~J 中的分割結(jié)果所示;LPF*模型通過計算三通道(R,G,B)的局部強度均值來定義局部預(yù)擬合函數(shù),在弱邊緣圖像分割中表現(xiàn)較好,但該模型對噪聲較敏感,如在圖像G、J 中分割結(jié)果所示。而本文模型在自適應(yīng)符號函數(shù)加權(quán)的優(yōu)化LoG 能量項的基礎(chǔ)上,加入了改進的局部能量項,兩種能量項的結(jié)合不僅增強了模型的自適應(yīng)性,還能夠有效識別噪聲點,保留目標(biāo)弱邊緣信息,如在圖像E~F、H 中分割結(jié)果所示,在圖8 示例圖像中獲得的分割結(jié)果與真實目標(biāo)邊界更為接近,其中在分割圖像G 時,迭代次數(shù)為1 700。

    圖8 六種模型分割結(jié)果Fig.8 Segmentation results of six models

    采用骰子相似系數(shù)(dice similarity coefficient,DSC)、Jaccard 相似系數(shù)(Jaccard similarity coefficient,JSI)、Hausdorff 距離(Hausdorff distance,HD)三種指標(biāo)來進一步衡量各模型分割結(jié)果的質(zhì)量,定義如式(30)~(32)所示。

    表4 和表5 分別給出六種模型在圖8 示例圖像分割結(jié)果對應(yīng)的DSC、JSI 與HD 的值。結(jié)合表4、表5數(shù)據(jù)及圖8 的實驗結(jié)果可以看出,本文模型能夠在背景復(fù)雜及灰度呈非均勻分布的圖像中有效提取目標(biāo)弱邊緣,且本文模型分割每幅圖片平均迭代時間為16.0 s,與分割效率最高的SPF*模型僅12.1 s 之差,而本文模型取得了更精確的分割結(jié)果??紤]到馬氏距離不受量綱的影響,因此在計算簇與簇之間的差異性時,使用馬氏距離對圖7、圖8 中圖像進行測試,部分圖像無法獲得準(zhǔn)確的目標(biāo)先驗信息,這可能是在計算馬氏距離時,放大了變化微小的變量作用,導(dǎo)致提取的視覺顯著區(qū)域不夠準(zhǔn)確,尋找一種更加有效的距離度量方法來提取視覺顯著區(qū)域也能夠在一定程度上提高本文模型的分割精度。

    表4 六種模型在圖8 圖像上的分割精度Table 4 Segmentation accuracy of six models on Fig.8

    表5 六種模型在圖8 圖像上的分割評價(HD)Table 5 Segmentation evaluation of six models on Fig.8(HD)

    6 結(jié)束語

    本文結(jié)合改進的視覺顯著性檢測算法,為初始化輪廓的設(shè)置提供可靠的依據(jù),并將獲取的目標(biāo)先驗信息用于自適應(yīng)符號函數(shù)的設(shè)計,對優(yōu)化LoG 能量項進行加權(quán),提出自適應(yīng)符號函數(shù)加權(quán)的優(yōu)化LoG能量項,并與改進的局部能量項以線性方式相結(jié)合。一方面,提高模型在弱邊緣的敏感程度,保留目標(biāo)輪廓的連貫性與完整性;另一方面,降低模型陷入局部極值的可能,增強了模型的自適應(yīng)能力。接下來,將進一步優(yōu)化顯著圖提取算法,為主動輪廓模型中的曲線演化提供更準(zhǔn)確的先驗信息。

    猜你喜歡
    測度輪廓灰度
    輪廓錯覺
    采用改進導(dǎo)重法的拓撲結(jié)構(gòu)灰度單元過濾技術(shù)
    三個數(shù)字集生成的自相似測度的乘積譜
    R1上莫朗測度關(guān)于幾何平均誤差的最優(yōu)Vornoi分劃
    基于灰度拉伸的圖像水位識別方法研究
    非等熵Chaplygin氣體測度值解存在性
    Cookie-Cutter集上的Gibbs測度
    基于實時輪廓誤差估算的數(shù)控系統(tǒng)輪廓控制
    基于最大加權(quán)投影求解的彩色圖像灰度化對比度保留算法
    基于灰度線性建模的亞像素圖像抖動量計算
    国产毛片a区久久久久| 亚洲精品成人久久久久久| 看十八女毛片水多多多| 亚洲天堂国产精品一区在线| 亚洲最大成人手机在线| 国产成年人精品一区二区| 男人舔奶头视频| 两性午夜刺激爽爽歪歪视频在线观看| 可以在线观看毛片的网站| 99热这里只有精品一区| 禁无遮挡网站| 欧美性感艳星| 高清日韩中文字幕在线| 日日啪夜夜撸| 寂寞人妻少妇视频99o| 久久精品国产亚洲av涩爱 | 精品久久久久久久末码| 少妇的逼好多水| 婷婷精品国产亚洲av| 高清毛片免费看| 亚洲欧美日韩高清专用| 久久精品综合一区二区三区| 久久精品人妻少妇| 国产免费一级a男人的天堂| 国产午夜精品论理片| av黄色大香蕉| 成人二区视频| 午夜精品在线福利| 综合色丁香网| 精品久久国产蜜桃| 国产精品一及| 亚洲精品日韩av片在线观看| 日本黄大片高清| 亚洲精品一卡2卡三卡4卡5卡| 成人精品一区二区免费| 国内久久婷婷六月综合欲色啪| 亚洲欧美中文字幕日韩二区| 18禁在线无遮挡免费观看视频 | 高清毛片免费看| 亚洲成a人片在线一区二区| 99riav亚洲国产免费| 深爱激情五月婷婷| 国产一区二区亚洲精品在线观看| 国产蜜桃级精品一区二区三区| 美女大奶头视频| 国产av在哪里看| 在线观看av片永久免费下载| 三级男女做爰猛烈吃奶摸视频| 小蜜桃在线观看免费完整版高清| 一本精品99久久精品77| 桃色一区二区三区在线观看| 最近在线观看免费完整版| 性色avwww在线观看| 亚洲高清免费不卡视频| 国产黄色小视频在线观看| 最近视频中文字幕2019在线8| 长腿黑丝高跟| 黄片wwwwww| 亚洲精品日韩av片在线观看| 欧美又色又爽又黄视频| 成人av一区二区三区在线看| 亚洲欧美精品自产自拍| 日韩,欧美,国产一区二区三区 | 国产蜜桃级精品一区二区三区| 3wmmmm亚洲av在线观看| 日本与韩国留学比较| 日本五十路高清| 成人特级黄色片久久久久久久| 少妇高潮的动态图| 亚洲性久久影院| 久久久精品大字幕| 久久精品国产亚洲av天美| 97在线视频观看| 天天一区二区日本电影三级| 亚洲无线观看免费| 91在线观看av| 日韩大尺度精品在线看网址| 久久精品国产鲁丝片午夜精品| 亚洲精品456在线播放app| 久久99热6这里只有精品| 一进一出抽搐动态| 别揉我奶头~嗯~啊~动态视频| 插阴视频在线观看视频| 亚洲高清免费不卡视频| 一夜夜www| 丰满乱子伦码专区| 亚洲18禁久久av| 91麻豆精品激情在线观看国产| 久久精品国产亚洲av香蕉五月| 久久亚洲精品不卡| 国产一区二区在线观看日韩| 日本色播在线视频| 99国产极品粉嫩在线观看| 日韩欧美 国产精品| 日本-黄色视频高清免费观看| 国产精品一区二区三区四区免费观看 | 亚洲美女搞黄在线观看 | 亚洲人成网站在线播放欧美日韩| 国产精品国产高清国产av| 搞女人的毛片| 国产精品国产高清国产av| 国产精品99久久久久久久久| a级一级毛片免费在线观看| 级片在线观看| 国内久久婷婷六月综合欲色啪| 日韩 亚洲 欧美在线| 一本久久中文字幕| a级毛片a级免费在线| 亚洲熟妇熟女久久| 欧美精品国产亚洲| 99热网站在线观看| 一级a爱片免费观看的视频| 日韩av不卡免费在线播放| 91久久精品电影网| 国产精品久久久久久久电影| 国产精品乱码一区二三区的特点| 啦啦啦韩国在线观看视频| 成人鲁丝片一二三区免费| 久99久视频精品免费| 少妇裸体淫交视频免费看高清| 亚洲国产欧洲综合997久久,| 久久99热这里只有精品18| 色在线成人网| 国产午夜精品久久久久久一区二区三区 | 久久久久久久亚洲中文字幕| 老司机福利观看| 欧美极品一区二区三区四区| 女的被弄到高潮叫床怎么办| 色av中文字幕| 亚洲久久久久久中文字幕| 麻豆av噜噜一区二区三区| 欧美xxxx性猛交bbbb| 国产黄色视频一区二区在线观看 | 国内少妇人妻偷人精品xxx网站| 91午夜精品亚洲一区二区三区| 嫩草影视91久久| 久久国产乱子免费精品| 搡女人真爽免费视频火全软件 | 久久6这里有精品| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 男人和女人高潮做爰伦理| 精品人妻熟女av久视频| 97超级碰碰碰精品色视频在线观看| 女人被狂操c到高潮| 中文字幕免费在线视频6| 天堂影院成人在线观看| 久久久精品大字幕| 18+在线观看网站| 中国美白少妇内射xxxbb| 亚洲av五月六月丁香网| 国产熟女欧美一区二区| 国产免费一级a男人的天堂| 亚洲精品日韩在线中文字幕 | 三级国产精品欧美在线观看| 男人舔奶头视频| 97人妻精品一区二区三区麻豆| 亚洲综合色惰| 久久久精品94久久精品| av视频在线观看入口| 毛片女人毛片| 中文字幕精品亚洲无线码一区| 99热这里只有是精品在线观看| 欧美极品一区二区三区四区| av专区在线播放| 少妇的逼好多水| 国产男人的电影天堂91| 国产成人福利小说| 国产精品精品国产色婷婷| 嫩草影院新地址| 午夜福利在线观看吧| 国产精品一区www在线观看| 干丝袜人妻中文字幕| 日本三级黄在线观看| 成人欧美大片| 亚洲,欧美,日韩| 九九在线视频观看精品| 永久网站在线| 色在线成人网| 国产精品免费一区二区三区在线| 亚洲内射少妇av| 黄色视频,在线免费观看| 最近中文字幕高清免费大全6| 欧美zozozo另类| 在现免费观看毛片| 色视频www国产| 精品久久久久久久久久久久久| 午夜精品一区二区三区免费看| 18+在线观看网站| 欧美最新免费一区二区三区| 久久热精品热| 男女下面进入的视频免费午夜| 长腿黑丝高跟| 中国国产av一级| 国产精品久久久久久av不卡| 国产美女午夜福利| 国产v大片淫在线免费观看| 99热这里只有是精品50| 精品久久久久久久人妻蜜臀av| 成人一区二区视频在线观看| 亚洲欧美成人精品一区二区| 久久鲁丝午夜福利片| 熟女电影av网| 欧美bdsm另类| 亚州av有码| 日本免费a在线| 亚洲精品影视一区二区三区av| 乱系列少妇在线播放| 精品久久久久久成人av| 天堂网av新在线| 国产免费男女视频| 精品熟女少妇av免费看| 日韩精品青青久久久久久| 精品久久久久久久久亚洲| 床上黄色一级片| 精品一区二区三区视频在线| 日本撒尿小便嘘嘘汇集6| 少妇人妻精品综合一区二区 | 亚洲精品国产成人久久av| av视频在线观看入口| 久99久视频精品免费| 九九爱精品视频在线观看| 真人做人爱边吃奶动态| av天堂中文字幕网| 日韩av不卡免费在线播放| 一a级毛片在线观看| 国产精品免费一区二区三区在线| 成年女人永久免费观看视频| 欧美xxxx黑人xx丫x性爽| 日韩av在线大香蕉| 亚洲婷婷狠狠爱综合网| 一夜夜www| 午夜老司机福利剧场| 成人鲁丝片一二三区免费| 午夜精品一区二区三区免费看| 久久鲁丝午夜福利片| 男女啪啪激烈高潮av片| 日韩制服骚丝袜av| 成人特级av手机在线观看| 日韩av在线大香蕉| 嫩草影院入口| 成人漫画全彩无遮挡| 熟女电影av网| 一卡2卡三卡四卡精品乱码亚洲| 免费电影在线观看免费观看| 国语自产精品视频在线第100页| 亚洲第一区二区三区不卡| 在线观看av片永久免费下载| 男人舔女人下体高潮全视频| 少妇的逼好多水| 久久天躁狠狠躁夜夜2o2o| 99热全是精品| 联通29元200g的流量卡| 国语自产精品视频在线第100页| 人人妻,人人澡人人爽秒播| 天堂动漫精品| 久久亚洲精品不卡| 日本精品一区二区三区蜜桃| 欧美激情在线99| 99热6这里只有精品| 国产 一区 欧美 日韩| 日本熟妇午夜| 国产乱人视频| 国产精品一区二区三区四区久久| 一级毛片久久久久久久久女| 人妻少妇偷人精品九色| 老师上课跳d突然被开到最大视频| 欧美人与善性xxx| 久久欧美精品欧美久久欧美| 男人舔女人下体高潮全视频| 欧美性感艳星| 国产麻豆成人av免费视频| 国产成人freesex在线 | av卡一久久| 免费看日本二区| 少妇的逼水好多| 午夜福利在线观看免费完整高清在 | 亚洲欧美日韩高清在线视频| 欧美最新免费一区二区三区| 91av网一区二区| 九九热线精品视视频播放| 国产欧美日韩一区二区精品| 国产大屁股一区二区在线视频| 99久久九九国产精品国产免费| eeuss影院久久| 日韩欧美免费精品| 最近的中文字幕免费完整| 人妻丰满熟妇av一区二区三区| 亚洲专区国产一区二区| 在线观看一区二区三区| 观看美女的网站| a级毛色黄片| 99久久精品一区二区三区| 精品熟女少妇av免费看| 国产精品久久久久久精品电影| 亚洲国产精品成人综合色| 91精品国产九色| 国产精品亚洲美女久久久| 国产乱人视频| 久久精品影院6| 国产亚洲av嫩草精品影院| 亚洲天堂国产精品一区在线| 三级经典国产精品| 黄片wwwwww| 亚洲七黄色美女视频| 久久精品国产99精品国产亚洲性色| 欧美高清成人免费视频www| 久久热精品热| 少妇的逼好多水| 欧美色视频一区免费| av在线老鸭窝| 在线天堂最新版资源| 国产成人91sexporn| 成人无遮挡网站| 亚洲第一电影网av| 欧美不卡视频在线免费观看| 深夜a级毛片| 国产高清有码在线观看视频| 国产精品伦人一区二区| 欧美+日韩+精品| 此物有八面人人有两片| 国产男人的电影天堂91| 国产亚洲av嫩草精品影院| 十八禁国产超污无遮挡网站| 成人亚洲欧美一区二区av| 男女做爰动态图高潮gif福利片| 亚洲av熟女| 精品久久久久久久久久久久久| 22中文网久久字幕| 天堂√8在线中文| 俺也久久电影网| 精品久久久久久久久久久久久| 人人妻人人澡人人爽人人夜夜 | 午夜久久久久精精品| 女同久久另类99精品国产91| 黄片wwwwww| 欧美在线一区亚洲| 午夜福利成人在线免费观看| 国产成人一区二区在线| av在线观看视频网站免费| 亚洲精品粉嫩美女一区| 国产真实乱freesex| 99久久精品一区二区三区| 两个人的视频大全免费| 精品久久久久久久久av| 亚洲国产色片| 国产色爽女视频免费观看| 精品福利观看| 深夜a级毛片| 日本a在线网址| 蜜桃亚洲精品一区二区三区| 久久久久久久久久成人| 97碰自拍视频| 中文字幕熟女人妻在线| 亚洲成a人片在线一区二区| 亚洲国产日韩欧美精品在线观看| 日本五十路高清| 超碰av人人做人人爽久久| 精品午夜福利视频在线观看一区| 免费无遮挡裸体视频| 久久亚洲精品不卡| 欧美成人一区二区免费高清观看| 级片在线观看| 欧美三级亚洲精品| www.色视频.com| 国产成人影院久久av| av在线老鸭窝| 永久网站在线| 精品久久国产蜜桃| 少妇人妻精品综合一区二区 | 99久久精品国产国产毛片| 黄色欧美视频在线观看| 国产乱人视频| 欧美zozozo另类| 又爽又黄a免费视频| 我的老师免费观看完整版| 狂野欧美激情性xxxx在线观看| 麻豆乱淫一区二区| 日本熟妇午夜| 日日撸夜夜添| 国产亚洲精品久久久com| 国产精品电影一区二区三区| 久久久午夜欧美精品| 日日干狠狠操夜夜爽| 又爽又黄a免费视频| 午夜福利在线观看免费完整高清在 | 婷婷精品国产亚洲av| 亚洲精品色激情综合| 国产精品美女特级片免费视频播放器| 亚洲欧美成人精品一区二区| 五月伊人婷婷丁香| 午夜老司机福利剧场| 99热精品在线国产| av在线观看视频网站免费| 一个人观看的视频www高清免费观看| 国产又黄又爽又无遮挡在线| 91久久精品电影网| av.在线天堂| 成人特级av手机在线观看| 又爽又黄无遮挡网站| 免费av不卡在线播放| 亚洲精品粉嫩美女一区| 女人被狂操c到高潮| 老熟妇乱子伦视频在线观看| 香蕉av资源在线| 欧美xxxx性猛交bbbb| 欧美绝顶高潮抽搐喷水| 国产精品一及| 亚洲成人久久爱视频| 国产免费男女视频| 免费电影在线观看免费观看| 日韩 亚洲 欧美在线| 免费看光身美女| 99热只有精品国产| 秋霞在线观看毛片| h日本视频在线播放| 久久精品影院6| 一级av片app| 熟女电影av网| 亚洲无线在线观看| 免费黄网站久久成人精品| 成人av一区二区三区在线看| 18禁在线播放成人免费| 国产综合懂色| 日韩欧美 国产精品| 99九九线精品视频在线观看视频| 久久久久久久久久久丰满| 国产高潮美女av| 免费大片18禁| 亚洲天堂国产精品一区在线| 老熟妇仑乱视频hdxx| 一级毛片久久久久久久久女| 又爽又黄a免费视频| 又粗又爽又猛毛片免费看| 亚洲欧美日韩卡通动漫| 青春草视频在线免费观看| av国产免费在线观看| 一本久久中文字幕| 人人妻人人看人人澡| 色噜噜av男人的天堂激情| 亚洲欧美日韩高清专用| 国产精品人妻久久久久久| 草草在线视频免费看| 国产亚洲av嫩草精品影院| 97超碰精品成人国产| 日本 av在线| 亚洲四区av| 久久久久免费精品人妻一区二区| 日韩制服骚丝袜av| 国产人妻一区二区三区在| 国产精品国产高清国产av| 亚洲av美国av| 日本黄色片子视频| .国产精品久久| 丰满乱子伦码专区| 欧美高清性xxxxhd video| 午夜福利高清视频| 亚洲国产精品成人综合色| 中文字幕久久专区| 少妇被粗大猛烈的视频| 简卡轻食公司| 成年版毛片免费区| 天堂动漫精品| 亚洲欧美精品综合久久99| 国产精品一区二区性色av| 岛国在线免费视频观看| 欧美最黄视频在线播放免费| 香蕉av资源在线| 亚洲欧美日韩高清专用| 两性午夜刺激爽爽歪歪视频在线观看| 99久久中文字幕三级久久日本| 麻豆av噜噜一区二区三区| 精品熟女少妇av免费看| 欧美日韩乱码在线| 一本精品99久久精品77| 3wmmmm亚洲av在线观看| 国产成人aa在线观看| 国产老妇女一区| 嫩草影院新地址| 国产激情偷乱视频一区二区| 少妇丰满av| 欧美xxxx性猛交bbbb| 免费人成视频x8x8入口观看| 国产在视频线在精品| 村上凉子中文字幕在线| 亚洲欧美日韩高清专用| 大香蕉久久网| 狂野欧美激情性xxxx在线观看| 丝袜喷水一区| av国产免费在线观看| 三级毛片av免费| 全区人妻精品视频| 日日摸夜夜添夜夜爱| 欧美高清成人免费视频www| 国产精品一区二区三区四区久久| 久久久国产成人精品二区| 悠悠久久av| 久久久久国内视频| 91狼人影院| 国产一区二区在线av高清观看| 级片在线观看| www日本黄色视频网| 韩国av在线不卡| 国产精品一区二区性色av| 亚洲婷婷狠狠爱综合网| 伊人久久精品亚洲午夜| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清在线视频| www.色视频.com| 亚洲精品456在线播放app| 男人舔奶头视频| 久久久国产成人免费| 成人亚洲欧美一区二区av| 天美传媒精品一区二区| 亚洲人成网站在线观看播放| 99热只有精品国产| av黄色大香蕉| 最新中文字幕久久久久| 女人被狂操c到高潮| 91久久精品电影网| 男人的好看免费观看在线视频| 老熟妇乱子伦视频在线观看| 天美传媒精品一区二区| 亚洲人成网站在线观看播放| 久久久久久久久大av| 内地一区二区视频在线| 91av网一区二区| 免费在线观看成人毛片| 亚洲中文字幕日韩| 卡戴珊不雅视频在线播放| 亚洲成av人片在线播放无| 有码 亚洲区| 精品久久久噜噜| 亚洲一级一片aⅴ在线观看| 国产精品久久电影中文字幕| 91av网一区二区| 女人被狂操c到高潮| 久久久午夜欧美精品| 在现免费观看毛片| 小说图片视频综合网站| av视频在线观看入口| 久久久久久久久中文| 国产av不卡久久| 中文字幕人妻熟人妻熟丝袜美| 少妇人妻一区二区三区视频| 国产美女午夜福利| 国产精品一区二区性色av| 男插女下体视频免费在线播放| 亚洲av五月六月丁香网| 成人漫画全彩无遮挡| 午夜免费男女啪啪视频观看 | 午夜日韩欧美国产| 亚洲精品粉嫩美女一区| 精品久久久噜噜| 亚洲av.av天堂| 99久久成人亚洲精品观看| 国内揄拍国产精品人妻在线| 国产精品久久久久久久久免| 日韩av不卡免费在线播放| 久久久色成人| 国产成人精品久久久久久| 一本久久中文字幕| 欧美区成人在线视频| 欧美xxxx黑人xx丫x性爽| 久久久久国产精品人妻aⅴ院| 观看美女的网站| 搡老熟女国产l中国老女人| 床上黄色一级片| 午夜精品一区二区三区免费看| 久久久成人免费电影| 色尼玛亚洲综合影院| 欧美bdsm另类| 欧美极品一区二区三区四区| 美女大奶头视频| 亚洲成a人片在线一区二区| 亚洲精品在线观看二区| 国产片特级美女逼逼视频| 亚洲一区二区三区色噜噜| 午夜福利在线观看吧| 亚洲专区国产一区二区| 男女之事视频高清在线观看| 久久精品国产亚洲av涩爱 | 成年女人永久免费观看视频| 成人综合一区亚洲| 床上黄色一级片| 看片在线看免费视频| 亚洲av成人av| 免费一级毛片在线播放高清视频| 97在线视频观看| 老女人水多毛片| 亚洲成人精品中文字幕电影| 夜夜夜夜夜久久久久| 波多野结衣高清无吗| 日本一二三区视频观看| 人人妻人人澡人人爽人人夜夜 | 热99在线观看视频| 直男gayav资源| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 亚洲人成网站在线观看播放| 欧美性猛交╳xxx乱大交人| 特级一级黄色大片| 亚洲国产精品成人综合色| 欧美区成人在线视频| 禁无遮挡网站| 一区二区三区高清视频在线| 美女 人体艺术 gogo| 深夜a级毛片| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 午夜精品一区二区三区免费看| 99在线人妻在线中文字幕| 97超碰精品成人国产| 日韩欧美免费精品|