• <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)投影求解的彩色圖像灰度化對比度保留算法
    基于灰度線性建模的亞像素圖像抖動量計算
    69精品国产乱码久久久| 国产白丝娇喘喷水9色精品| 啦啦啦啦在线视频资源| 99热6这里只有精品| 狂野欧美激情性xxxx在线观看| 国产成人精品久久久久久| 日韩制服丝袜自拍偷拍| 美国免费a级毛片| 久久久久视频综合| 欧美亚洲 丝袜 人妻 在线| 国产精品 国内视频| 巨乳人妻的诱惑在线观看| 一本久久精品| 99热网站在线观看| 男女免费视频国产| 亚洲成av片中文字幕在线观看 | 边亲边吃奶的免费视频| 久久亚洲国产成人精品v| 在线亚洲精品国产二区图片欧美| 咕卡用的链子| 色吧在线观看| 免费人妻精品一区二区三区视频| 国产av一区二区精品久久| 午夜精品国产一区二区电影| 好男人视频免费观看在线| 亚洲伊人久久精品综合| 午夜福利网站1000一区二区三区| 亚洲欧美色中文字幕在线| 免费看av在线观看网站| 18禁在线无遮挡免费观看视频| 久久这里只有精品19| 老司机影院毛片| 伦精品一区二区三区| 久久久a久久爽久久v久久| 性色av一级| 中文字幕精品免费在线观看视频 | 久久久久久人人人人人| 少妇猛男粗大的猛烈进出视频| 国产欧美另类精品又又久久亚洲欧美| 国产一区二区三区av在线| 日韩熟女老妇一区二区性免费视频| 欧美成人午夜免费资源| av女优亚洲男人天堂| 一本久久精品| 热re99久久国产66热| 丰满乱子伦码专区| 最新中文字幕久久久久| 亚洲色图 男人天堂 中文字幕 | 桃花免费在线播放| 欧美日韩综合久久久久久| a级片在线免费高清观看视频| 亚洲图色成人| 婷婷成人精品国产| 国产乱来视频区| 女性被躁到高潮视频| 18禁在线无遮挡免费观看视频| 一级黄片播放器| 黄片播放在线免费| 午夜久久久在线观看| 国产欧美另类精品又又久久亚洲欧美| 日韩大片免费观看网站| 又大又黄又爽视频免费| 天天躁夜夜躁狠狠久久av| 老司机影院毛片| 黄色配什么色好看| 99re6热这里在线精品视频| 午夜免费观看性视频| 晚上一个人看的免费电影| 久热久热在线精品观看| 熟女人妻精品中文字幕| 黄色毛片三级朝国网站| 欧美精品av麻豆av| 夫妻午夜视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品久久成人aⅴ小说| 飞空精品影院首页| 亚洲精品第二区| 97人妻天天添夜夜摸| www日本在线高清视频| 精品久久久精品久久久| 国产成人a∨麻豆精品| 人人妻人人添人人爽欧美一区卜| 久久人人爽av亚洲精品天堂| 久久久精品94久久精品| 欧美国产精品一级二级三级| 中文字幕最新亚洲高清| 色哟哟·www| 激情五月婷婷亚洲| 亚洲综合色网址| 久久久久精品久久久久真实原创| 国产黄频视频在线观看| 永久免费av网站大全| 波野结衣二区三区在线| 精品少妇黑人巨大在线播放| 亚洲经典国产精华液单| 成年动漫av网址| 国产淫语在线视频| 美女视频免费永久观看网站| 久久久精品区二区三区| 欧美 亚洲 国产 日韩一| 毛片一级片免费看久久久久| 又粗又硬又长又爽又黄的视频| 亚洲精品乱久久久久久| 国产精品国产三级国产av玫瑰| 人人妻人人澡人人看| 这个男人来自地球电影免费观看 | 人人澡人人妻人| 青春草亚洲视频在线观看| 爱豆传媒免费全集在线观看| 日本wwww免费看| 亚洲欧美成人综合另类久久久| 久久这里只有精品19| 国产国语露脸激情在线看| 久久国产亚洲av麻豆专区| 国产爽快片一区二区三区| 国产爽快片一区二区三区| 午夜精品国产一区二区电影| 婷婷成人精品国产| 韩国高清视频一区二区三区| a级毛色黄片| 国产一级毛片在线| 热re99久久精品国产66热6| 少妇熟女欧美另类| 女人被躁到高潮嗷嗷叫费观| 亚洲美女搞黄在线观看| 乱码一卡2卡4卡精品| 国产一区有黄有色的免费视频| 亚洲第一av免费看| 午夜免费观看性视频| 深夜精品福利| 老司机影院毛片| 青春草视频在线免费观看| 丰满乱子伦码专区| a 毛片基地| 国产一区有黄有色的免费视频| 少妇人妻精品综合一区二区| 老司机亚洲免费影院| av国产久精品久网站免费入址| 国产av国产精品国产| av在线播放精品| 成人漫画全彩无遮挡| 成人18禁高潮啪啪吃奶动态图| 国产一区亚洲一区在线观看| 精品国产一区二区久久| 国产精品无大码| 免费观看a级毛片全部| 免费高清在线观看视频在线观看| 满18在线观看网站| 在线看a的网站| 欧美xxⅹ黑人| 亚洲国产av影院在线观看| av天堂久久9| 久久ye,这里只有精品| 青春草亚洲视频在线观看| 欧美精品av麻豆av| 激情五月婷婷亚洲| 久久人人爽人人爽人人片va| av又黄又爽大尺度在线免费看| 日韩av免费高清视频| 黑人猛操日本美女一级片| av又黄又爽大尺度在线免费看| 啦啦啦在线观看免费高清www| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美另类精品又又久久亚洲欧美| 国产亚洲午夜精品一区二区久久| 国产亚洲最大av| 欧美精品av麻豆av| 中文字幕最新亚洲高清| 9色porny在线观看| 国产精品欧美亚洲77777| 少妇的丰满在线观看| 日韩av免费高清视频| 国产精品一区二区在线不卡| 黑人猛操日本美女一级片| 中国国产av一级| 丰满迷人的少妇在线观看| 在线观看国产h片| 国产精品免费大片| 色哟哟·www| 美女xxoo啪啪120秒动态图| 中文字幕最新亚洲高清| 国产精品一区二区在线不卡| av天堂久久9| 久久午夜综合久久蜜桃| xxxhd国产人妻xxx| 国产免费视频播放在线视频| 国产免费一区二区三区四区乱码| 午夜福利影视在线免费观看| 2021少妇久久久久久久久久久| 丝袜脚勾引网站| 久久久久国产精品人妻一区二区| 99热6这里只有精品| 亚洲精品国产av蜜桃| 亚洲av在线观看美女高潮| 日韩制服丝袜自拍偷拍| 乱人伦中国视频| 亚洲第一av免费看| 国产精品麻豆人妻色哟哟久久| 在线观看一区二区三区激情| 精品国产一区二区三区四区第35| 国产成人午夜福利电影在线观看| 久久毛片免费看一区二区三区| 欧美97在线视频| 久久国产精品男人的天堂亚洲 | 人人妻人人澡人人爽人人夜夜| 国产精品.久久久| 欧美最新免费一区二区三区| av有码第一页| 女性生殖器流出的白浆| 成人亚洲欧美一区二区av| 欧美xxⅹ黑人| 亚洲精品aⅴ在线观看| 伊人亚洲综合成人网| 亚洲精品视频女| 老司机影院毛片| 亚洲人成77777在线视频| 熟女电影av网| 男女啪啪激烈高潮av片| 日本av免费视频播放| av免费在线看不卡| 精品久久蜜臀av无| 大码成人一级视频| 国产亚洲午夜精品一区二区久久| 欧美精品av麻豆av| 男男h啪啪无遮挡| kizo精华| 赤兔流量卡办理| 尾随美女入室| 亚洲精品456在线播放app| 国产精品一二三区在线看| 日韩一区二区三区影片| 亚洲av综合色区一区| 在线亚洲精品国产二区图片欧美| 乱人伦中国视频| 两个人免费观看高清视频| 国产在线视频一区二区| 日韩,欧美,国产一区二区三区| 免费高清在线观看日韩| 日韩一本色道免费dvd| 亚洲色图 男人天堂 中文字幕 | 欧美精品国产亚洲| av天堂久久9| 久久影院123| 大片免费播放器 马上看| 国产黄色视频一区二区在线观看| 乱人伦中国视频| 18+在线观看网站| 大陆偷拍与自拍| 国产精品成人在线| 美女国产视频在线观看| 国产成人免费无遮挡视频| 国产精品国产三级国产专区5o| 国产在线一区二区三区精| 97超碰精品成人国产| 中文字幕人妻丝袜制服| 蜜臀久久99精品久久宅男| 最近最新中文字幕免费大全7| 欧美日韩视频高清一区二区三区二| 又黄又粗又硬又大视频| 一级黄片播放器| 精品亚洲成国产av| a级片在线免费高清观看视频| 人体艺术视频欧美日本| 91aial.com中文字幕在线观看| 精品午夜福利在线看| 天美传媒精品一区二区| 老司机影院成人| 免费黄网站久久成人精品| 国产精品国产av在线观看| 成人亚洲精品一区在线观看| freevideosex欧美| 我的女老师完整版在线观看| 国产精品一区二区在线不卡| 一本久久精品| 在线亚洲精品国产二区图片欧美| 免费人成在线观看视频色| 色94色欧美一区二区| 18在线观看网站| 9色porny在线观看| av国产久精品久网站免费入址| 99久久精品国产国产毛片| 狠狠精品人妻久久久久久综合| 考比视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 午夜日本视频在线| 久久久久精品久久久久真实原创| 亚洲内射少妇av| 亚洲婷婷狠狠爱综合网| 丰满饥渴人妻一区二区三| 成人综合一区亚洲| 岛国毛片在线播放| 视频在线观看一区二区三区| 一区二区av电影网| 欧美另类一区| 99九九在线精品视频| 国产精品99久久99久久久不卡 | 久久精品国产亚洲av涩爱| 欧美xxⅹ黑人| 亚洲av男天堂| 午夜免费鲁丝| 亚洲欧美中文字幕日韩二区| 纯流量卡能插随身wifi吗| 国产av一区二区精品久久| 成年美女黄网站色视频大全免费| 欧美亚洲 丝袜 人妻 在线| √禁漫天堂资源中文www| 最后的刺客免费高清国语| 亚洲国产成人一精品久久久| 熟女av电影| 亚洲一区二区三区欧美精品| 精品午夜福利在线看| 嫩草影院入口| 久久久久久人妻| 久久鲁丝午夜福利片| 国产精品国产三级专区第一集| 中文精品一卡2卡3卡4更新| 在线观看免费高清a一片| 激情视频va一区二区三区| 国产不卡av网站在线观看| av网站免费在线观看视频| 日韩三级伦理在线观看| 亚洲欧美清纯卡通| 熟妇人妻不卡中文字幕| 不卡视频在线观看欧美| 欧美精品一区二区免费开放| 久久午夜综合久久蜜桃| 国产乱来视频区| 国产日韩欧美在线精品| 精品熟女少妇av免费看| 中国三级夫妇交换| 大片免费播放器 马上看| 女的被弄到高潮叫床怎么办| 99re6热这里在线精品视频| 汤姆久久久久久久影院中文字幕| 一级片'在线观看视频| 欧美 亚洲 国产 日韩一| 亚洲精品一区蜜桃| 国产成人精品在线电影| 日韩伦理黄色片| 中文字幕av电影在线播放| 女的被弄到高潮叫床怎么办| 日本黄大片高清| 最近最新中文字幕免费大全7| 欧美日韩国产mv在线观看视频| 黑人巨大精品欧美一区二区蜜桃 | 国产国语露脸激情在线看| 人人妻人人添人人爽欧美一区卜| 99热这里只有是精品在线观看| 日韩成人伦理影院| 午夜日本视频在线| 爱豆传媒免费全集在线观看| 丰满少妇做爰视频| 老女人水多毛片| 美女脱内裤让男人舔精品视频| 最近最新中文字幕免费大全7| 夜夜骑夜夜射夜夜干| 久久久久久久亚洲中文字幕| 国产av码专区亚洲av| 精品一区二区免费观看| 精品少妇久久久久久888优播| 蜜桃国产av成人99| 成年人免费黄色播放视频| 日韩av不卡免费在线播放| 少妇的逼好多水| 中文字幕精品免费在线观看视频 | 国产精品成人在线| 国产免费一区二区三区四区乱码| 丝瓜视频免费看黄片| 九九在线视频观看精品| 国产成人91sexporn| 人妻人人澡人人爽人人| 高清在线视频一区二区三区| 黑丝袜美女国产一区| 美女国产视频在线观看| 最近手机中文字幕大全| 日本欧美视频一区| 男的添女的下面高潮视频| 国产男女内射视频| 国产片特级美女逼逼视频| 成人午夜精彩视频在线观看| 丝袜人妻中文字幕| 免费大片18禁| 精品99又大又爽又粗少妇毛片| 国产精品.久久久| 少妇人妻 视频| 亚洲情色 制服丝袜| 国产片特级美女逼逼视频| 国产亚洲一区二区精品| 天堂俺去俺来也www色官网| 又粗又硬又长又爽又黄的视频| 最近中文字幕2019免费版| 欧美 日韩 精品 国产| 80岁老熟妇乱子伦牲交| 男女免费视频国产| 欧美少妇被猛烈插入视频| 三级国产精品片| 久久久久精品人妻al黑| 国产麻豆69| 国产亚洲精品久久久com| 精品人妻在线不人妻| 亚洲精品,欧美精品| 最近中文字幕高清免费大全6| 内地一区二区视频在线| 99视频精品全部免费 在线| 免费播放大片免费观看视频在线观看| 亚洲欧美清纯卡通| 午夜免费观看性视频| 亚洲精品,欧美精品| 日日摸夜夜添夜夜爱| 国产亚洲精品第一综合不卡 | 亚洲高清免费不卡视频| 国产高清三级在线| 日韩av在线免费看完整版不卡| 亚洲国产精品一区三区| 中文字幕制服av| 亚洲第一区二区三区不卡| 中文精品一卡2卡3卡4更新| 亚洲精品第二区| 久久鲁丝午夜福利片| 亚洲国产精品专区欧美| 水蜜桃什么品种好| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 成人综合一区亚洲| 亚洲精品乱码久久久久久按摩| 美女xxoo啪啪120秒动态图| 国产一级毛片在线| 国产高清不卡午夜福利| 国产精品三级大全| 国产成人欧美| 日韩熟女老妇一区二区性免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 99热全是精品| 亚洲国产日韩一区二区| 成人国产av品久久久| 国国产精品蜜臀av免费| 97精品久久久久久久久久精品| 欧美变态另类bdsm刘玥| 久久99蜜桃精品久久| 成人漫画全彩无遮挡| 精品久久国产蜜桃| 丝袜美足系列| 男女无遮挡免费网站观看| 国产精品无大码| 午夜福利影视在线免费观看| 亚洲精品久久成人aⅴ小说| 国产精品无大码| 好男人视频免费观看在线| 午夜影院在线不卡| 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| 亚洲av在线观看美女高潮| 亚洲激情五月婷婷啪啪| 高清不卡的av网站| 女人被躁到高潮嗷嗷叫费观| 女人久久www免费人成看片| 久久久国产一区二区| 国产精品蜜桃在线观看| 狠狠婷婷综合久久久久久88av| 中文字幕最新亚洲高清| 亚洲伊人色综图| 久久久亚洲精品成人影院| av卡一久久| 一区二区日韩欧美中文字幕 | √禁漫天堂资源中文www| 女人被躁到高潮嗷嗷叫费观| 久久人妻熟女aⅴ| 丰满迷人的少妇在线观看| 国产深夜福利视频在线观看| 边亲边吃奶的免费视频| 黄色毛片三级朝国网站| 国产av一区二区精品久久| 美女国产视频在线观看| 国产精品国产三级国产专区5o| 国产免费现黄频在线看| 一边亲一边摸免费视频| 亚洲国产成人一精品久久久| 99国产综合亚洲精品| 在线观看美女被高潮喷水网站| 久久久久久人妻| 日韩av在线免费看完整版不卡| 亚洲成人手机| 亚洲人成网站在线观看播放| 国产精品秋霞免费鲁丝片| 亚洲欧美中文字幕日韩二区| 欧美日韩一区二区视频在线观看视频在线| 国产成人免费无遮挡视频| 国产精品人妻久久久久久| 亚洲综合色惰| 三级国产精品片| 狠狠婷婷综合久久久久久88av| 精品一区在线观看国产| 久久精品国产鲁丝片午夜精品| 日韩av在线免费看完整版不卡| 国产精品久久久久久av不卡| 国产精品不卡视频一区二区| 精品一区二区免费观看| 少妇熟女欧美另类| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| 日本欧美视频一区| 视频区图区小说| 国产精品久久久久久精品电影小说| 成人免费观看视频高清| 男女高潮啪啪啪动态图| 丝瓜视频免费看黄片| 天天操日日干夜夜撸| 久久97久久精品| 免费少妇av软件| 视频在线观看一区二区三区| 日日摸夜夜添夜夜爱| 午夜视频国产福利| 国产成人免费观看mmmm| 国产精品熟女久久久久浪| 国产精品.久久久| 亚洲欧洲精品一区二区精品久久久 | 国产免费一级a男人的天堂| 久久久久久久精品精品| 免费人成在线观看视频色| 欧美日韩精品成人综合77777| 亚洲欧美日韩另类电影网站| 免费观看av网站的网址| 亚洲 欧美一区二区三区| 久久毛片免费看一区二区三区| 99久久中文字幕三级久久日本| 在线观看www视频免费| 中文天堂在线官网| 最后的刺客免费高清国语| www.av在线官网国产| 黑人欧美特级aaaaaa片| 精品一区二区三卡| 日日撸夜夜添| 女的被弄到高潮叫床怎么办| 男女边摸边吃奶| 亚洲av免费高清在线观看| 国产成人精品福利久久| 成人漫画全彩无遮挡| 成年女人在线观看亚洲视频| 色吧在线观看| 国产深夜福利视频在线观看| 午夜av观看不卡| 黄网站色视频无遮挡免费观看| 国产免费一级a男人的天堂| 久久国内精品自在自线图片| 成人免费观看视频高清| 国精品久久久久久国模美| 99国产精品免费福利视频| 国产片内射在线| 国产av一区二区精品久久| 中国三级夫妇交换| 国产爽快片一区二区三区| 美女主播在线视频| 欧美 日韩 精品 国产| 免费高清在线观看视频在线观看| 97人妻天天添夜夜摸| 日韩人妻精品一区2区三区| 夜夜骑夜夜射夜夜干| 亚洲中文av在线| 人人澡人人妻人| 色5月婷婷丁香| 激情五月婷婷亚洲| h视频一区二区三区| 亚洲综合色网址| 国产xxxxx性猛交| 亚洲av福利一区| 久久国产精品男人的天堂亚洲 | 国产精品蜜桃在线观看| 国产精品成人在线| 成人二区视频| 免费观看a级毛片全部| 婷婷色av中文字幕| 国产一区二区三区av在线| a级毛色黄片| 国产精品一国产av| 国产69精品久久久久777片| 一区二区三区精品91| 在线观看www视频免费| 91国产中文字幕| 超碰97精品在线观看| 久久99精品国语久久久| 欧美丝袜亚洲另类| 日本色播在线视频| 婷婷色av中文字幕| av福利片在线| 丝瓜视频免费看黄片| 搡女人真爽免费视频火全软件| 久久精品国产a三级三级三级| 日本vs欧美在线观看视频| tube8黄色片| 亚洲综合色网址| 人妻 亚洲 视频| 黑人猛操日本美女一级片| 亚洲人成网站在线观看播放| 国产亚洲最大av| 国产成人一区二区在线| 视频中文字幕在线观看| 欧美人与善性xxx| 国产av国产精品国产| 亚洲成色77777| 卡戴珊不雅视频在线播放| 欧美人与善性xxx| 亚洲欧美日韩卡通动漫| 人人妻人人澡人人爽人人夜夜| 成人无遮挡网站| 熟女电影av网| 久久久久网色| 一级,二级,三级黄色视频| 性色av一级| 国产极品粉嫩免费观看在线| 国产乱人偷精品视频| 久久久久久伊人网av| 精品一品国产午夜福利视频| 午夜免费鲁丝| 熟女电影av网| 秋霞在线观看毛片| 亚洲美女视频黄频|