• <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)投影求解的彩色圖像灰度化對比度保留算法
    基于灰度線性建模的亞像素圖像抖動量計算
    精品久久久久久久久久久久久| 亚洲,欧美,日韩| 一区二区三区激情视频| 国产精品,欧美在线| 久久精品夜夜夜夜夜久久蜜豆| 好男人在线观看高清免费视频| 波多野结衣巨乳人妻| 男女之事视频高清在线观看| 欧美潮喷喷水| 国产高清激情床上av| 精品人妻熟女av久视频| 亚洲精华国产精华精| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 久久精品人妻少妇| 国产aⅴ精品一区二区三区波| 国产精品一及| 在线观看舔阴道视频| 亚洲人与动物交配视频| 国产精品av视频在线免费观看| 日韩欧美国产一区二区入口| av在线老鸭窝| 久久亚洲真实| 此物有八面人人有两片| 国产欧美日韩精品一区二区| 啪啪无遮挡十八禁网站| 欧美+亚洲+日韩+国产| 美女 人体艺术 gogo| 日本黄大片高清| 国产精品三级大全| 久久精品国产亚洲网站| 成人三级黄色视频| 天堂av国产一区二区熟女人妻| 精品一区二区三区视频在线| 国产精品免费一区二区三区在线| 国产精品99久久久久久久久| 少妇人妻精品综合一区二区 | 熟女电影av网| 99久久精品热视频| 啦啦啦韩国在线观看视频| 精品久久国产蜜桃| 国产伦精品一区二区三区四那| 少妇的逼好多水| 国产探花极品一区二区| 天堂√8在线中文| 在线国产一区二区在线| 国产精品久久久久久av不卡| 亚洲第一区二区三区不卡| 夜夜夜夜夜久久久久| 人人妻人人澡欧美一区二区| 日本在线视频免费播放| 日本爱情动作片www.在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 亚洲精华国产精华液的使用体验 | 国产精品不卡视频一区二区| 日韩精品中文字幕看吧| 国产精品人妻久久久影院| 中国美女看黄片| 国产主播在线观看一区二区| 又紧又爽又黄一区二区| 免费在线观看日本一区| 欧美3d第一页| 国产91精品成人一区二区三区| 国产亚洲精品久久久久久毛片| 亚洲av成人精品一区久久| 两个人视频免费观看高清| 亚洲成人精品中文字幕电影| 国产一区二区亚洲精品在线观看| 国产私拍福利视频在线观看| 亚洲色图av天堂| 国产成人aa在线观看| 床上黄色一级片| 欧美日韩综合久久久久久 | 最新中文字幕久久久久| 亚洲av熟女| 国产精品乱码一区二三区的特点| 国产精品一及| 亚洲欧美激情综合另类| 嫩草影院新地址| 久久精品夜夜夜夜夜久久蜜豆| 在线免费十八禁| 日韩欧美精品免费久久| 国产精品一及| 午夜影院日韩av| 伦理电影大哥的女人| 国产精华一区二区三区| 99视频精品全部免费 在线| 国产精品女同一区二区软件 | 国产乱人视频| 日韩欧美 国产精品| 久久久久久大精品| 精品无人区乱码1区二区| 国产欧美日韩精品一区二区| 午夜福利在线观看吧| avwww免费| 精品不卡国产一区二区三区| 久久精品国产亚洲网站| 亚洲专区中文字幕在线| 国产视频一区二区在线看| 成熟少妇高潮喷水视频| 欧美一级a爱片免费观看看| 亚洲熟妇熟女久久| 男人和女人高潮做爰伦理| 国内精品美女久久久久久| а√天堂www在线а√下载| 中出人妻视频一区二区| 99久久成人亚洲精品观看| 国产 一区精品| 毛片女人毛片| 亚洲中文字幕一区二区三区有码在线看| 久久久精品大字幕| 毛片一级片免费看久久久久 | 男人狂女人下面高潮的视频| 久久久久久国产a免费观看| 欧美日韩瑟瑟在线播放| 成年女人毛片免费观看观看9| 日韩欧美三级三区| or卡值多少钱| 亚洲美女黄片视频| 性色avwww在线观看| 精品国产三级普通话版| 成人性生交大片免费视频hd| 午夜福利18| 毛片女人毛片| 日本五十路高清| 成人av一区二区三区在线看| 久久精品国产亚洲av涩爱 | 精品久久久久久久久久久久久| 69人妻影院| 少妇丰满av| 亚洲专区中文字幕在线| .国产精品久久| 琪琪午夜伦伦电影理论片6080| 大又大粗又爽又黄少妇毛片口| 可以在线观看毛片的网站| 色综合婷婷激情| 深爱激情五月婷婷| 日日摸夜夜添夜夜添小说| 日韩强制内射视频| 国内少妇人妻偷人精品xxx网站| 国产色婷婷99| 午夜老司机福利剧场| 久久精品国产自在天天线| 在线播放国产精品三级| 五月玫瑰六月丁香| 午夜免费男女啪啪视频观看 | 欧美最新免费一区二区三区| 又爽又黄无遮挡网站| 国产精品永久免费网站| 亚洲七黄色美女视频| 2021天堂中文幕一二区在线观| 亚洲av中文字字幕乱码综合| 亚洲男人的天堂狠狠| 久久香蕉精品热| 在线观看午夜福利视频| 天天一区二区日本电影三级| 亚洲综合色惰| av国产免费在线观看| 91狼人影院| 日本欧美国产在线视频| 少妇被粗大猛烈的视频| 久久精品人妻少妇| 一个人看的www免费观看视频| 亚洲久久久久久中文字幕| 国产大屁股一区二区在线视频| 麻豆国产av国片精品| 亚洲aⅴ乱码一区二区在线播放| 欧美又色又爽又黄视频| 国产日本99.免费观看| 性插视频无遮挡在线免费观看| 少妇人妻一区二区三区视频| 少妇猛男粗大的猛烈进出视频 | 婷婷色综合大香蕉| 综合色av麻豆| 久99久视频精品免费| av国产免费在线观看| 亚洲av成人精品一区久久| 18禁在线播放成人免费| 欧美bdsm另类| 久久久久久九九精品二区国产| 12—13女人毛片做爰片一| 日本-黄色视频高清免费观看| 极品教师在线视频| 亚洲av.av天堂| 99热这里只有精品一区| 熟女电影av网| 91av网一区二区| 看片在线看免费视频| 亚洲性夜色夜夜综合| 亚洲欧美日韩无卡精品| 欧美日本视频| 亚洲成人久久爱视频| 男女视频在线观看网站免费| 国产免费av片在线观看野外av| 看片在线看免费视频| 久久久久久大精品| 在线观看一区二区三区| 老司机深夜福利视频在线观看| 亚洲四区av| 在线观看av片永久免费下载| av黄色大香蕉| 国产精品自产拍在线观看55亚洲| 99久久九九国产精品国产免费| 无人区码免费观看不卡| 国产精品人妻久久久影院| 成人精品一区二区免费| 午夜福利在线观看吧| 一进一出抽搐gif免费好疼| 亚洲av成人精品一区久久| 色视频www国产| 不卡一级毛片| 动漫黄色视频在线观看| 黄色日韩在线| 18+在线观看网站| 美女黄网站色视频| 狠狠狠狠99中文字幕| 少妇人妻一区二区三区视频| 国产精品福利在线免费观看| 国语自产精品视频在线第100页| 日韩欧美在线乱码| 久久精品影院6| 综合色av麻豆| 高清日韩中文字幕在线| 国产爱豆传媒在线观看| 国产探花极品一区二区| 国产亚洲精品综合一区在线观看| 久久草成人影院| 欧美性猛交╳xxx乱大交人| 亚洲美女视频黄频| 成人国产综合亚洲| 韩国av在线不卡| 日韩 亚洲 欧美在线| 免费在线观看日本一区| 99久久成人亚洲精品观看| 国产精品永久免费网站| 日日夜夜操网爽| 亚洲精华国产精华液的使用体验 | 人妻制服诱惑在线中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 夜夜看夜夜爽夜夜摸| avwww免费| 长腿黑丝高跟| 久久精品国产亚洲av天美| 欧美黑人巨大hd| 制服丝袜大香蕉在线| 99久久中文字幕三级久久日本| 国产伦精品一区二区三区四那| 精品久久久久久成人av| 精品久久久久久久久av| 国产三级在线视频| www日本黄色视频网| 免费一级毛片在线播放高清视频| 午夜激情欧美在线| 免费大片18禁| 亚洲av中文av极速乱 | 男女啪啪激烈高潮av片| 少妇人妻一区二区三区视频| 国产一区二区三区视频了| 搞女人的毛片| av在线老鸭窝| 韩国av在线不卡| 99在线视频只有这里精品首页| 女人十人毛片免费观看3o分钟| 国产欧美日韩精品亚洲av| 国产伦精品一区二区三区视频9| 12—13女人毛片做爰片一| 精品久久久久久久久久免费视频| 又粗又爽又猛毛片免费看| 人人妻人人看人人澡| 中文字幕av在线有码专区| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 免费搜索国产男女视频| 亚洲avbb在线观看| 成人精品一区二区免费| 国产视频一区二区在线看| 在线观看av片永久免费下载| 婷婷色综合大香蕉| 欧美精品国产亚洲| 男女之事视频高清在线观看| 18禁裸乳无遮挡免费网站照片| 91精品国产九色| 久久久久久久亚洲中文字幕| 小说图片视频综合网站| 又黄又爽又免费观看的视频| 干丝袜人妻中文字幕| 赤兔流量卡办理| 国产探花极品一区二区| 久久精品人妻少妇| 热99在线观看视频| 91麻豆精品激情在线观看国产| 黄色视频,在线免费观看| 国产伦人伦偷精品视频| 欧美区成人在线视频| 淫秽高清视频在线观看| 国产精品人妻久久久影院| 国产私拍福利视频在线观看| 精品人妻偷拍中文字幕| 欧美+亚洲+日韩+国产| 搡女人真爽免费视频火全软件 | 久久人妻av系列| 国产伦一二天堂av在线观看| 夜夜看夜夜爽夜夜摸| 久久国产乱子免费精品| 十八禁网站免费在线| 女人十人毛片免费观看3o分钟| 无人区码免费观看不卡| 在线播放国产精品三级| 变态另类成人亚洲欧美熟女| 精品人妻1区二区| 精品国产三级普通话版| 日韩欧美三级三区| 观看美女的网站| 亚洲国产欧美人成| 一本久久中文字幕| 亚洲人成伊人成综合网2020| 国产高清三级在线| 成人高潮视频无遮挡免费网站| 亚洲一区二区三区色噜噜| av天堂在线播放| 亚洲自偷自拍三级| 69人妻影院| 能在线免费观看的黄片| 国产精品伦人一区二区| 日韩欧美国产在线观看| 观看免费一级毛片| 午夜日韩欧美国产| 91狼人影院| 国产极品精品免费视频能看的| 成人毛片a级毛片在线播放| 欧美+日韩+精品| 国产欧美日韩一区二区精品| 亚洲美女搞黄在线观看 | 亚洲最大成人av| 国产精品自产拍在线观看55亚洲| 草草在线视频免费看| 亚洲精华国产精华液的使用体验 | 最新中文字幕久久久久| 精品人妻偷拍中文字幕| 国产 一区精品| 色噜噜av男人的天堂激情| 国内少妇人妻偷人精品xxx网站| 在线观看免费视频日本深夜| 白带黄色成豆腐渣| 亚洲va日本ⅴa欧美va伊人久久| 国产视频一区二区在线看| 美女高潮的动态| 一本精品99久久精品77| 在线天堂最新版资源| 男女边吃奶边做爰视频| 国产一区二区在线观看日韩| 日韩一本色道免费dvd| 亚洲一级一片aⅴ在线观看| 色播亚洲综合网| 狂野欧美白嫩少妇大欣赏| 午夜精品一区二区三区免费看| 国产成人一区二区在线| av福利片在线观看| 精品人妻一区二区三区麻豆 | av在线天堂中文字幕| 久久久国产成人免费| 国产伦精品一区二区三区四那| 永久网站在线| 在线免费观看的www视频| 露出奶头的视频| 99riav亚洲国产免费| 女人被狂操c到高潮| 别揉我奶头 嗯啊视频| 久久久久久久亚洲中文字幕| 国产精品99久久久久久久久| 国产精品电影一区二区三区| 免费人成在线观看视频色| 国产精品人妻久久久久久| 特级一级黄色大片| 亚洲三级黄色毛片| 国产亚洲精品av在线| 国产av在哪里看| 哪里可以看免费的av片| 99在线视频只有这里精品首页| 欧美潮喷喷水| 日韩人妻高清精品专区| 亚洲18禁久久av| 桃红色精品国产亚洲av| 欧美潮喷喷水| 国产精品无大码| 亚洲国产精品成人综合色| 欧美成人免费av一区二区三区| 99久久九九国产精品国产免费| 国产蜜桃级精品一区二区三区| 偷拍熟女少妇极品色| 国产探花极品一区二区| 久久久久久久久久黄片| 久久久久性生活片| 一级黄色大片毛片| 黄色视频,在线免费观看| 午夜久久久久精精品| 欧美又色又爽又黄视频| 久久久久国内视频| 赤兔流量卡办理| 国产精品久久久久久精品电影| 久久99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 嫩草影院新地址| 国产欧美日韩精品一区二区| 九九久久精品国产亚洲av麻豆| 免费在线观看成人毛片| 成年女人看的毛片在线观看| 少妇裸体淫交视频免费看高清| 久久久久久久亚洲中文字幕| 村上凉子中文字幕在线| 99riav亚洲国产免费| 国产午夜精品论理片| 亚洲av免费在线观看| 国产69精品久久久久777片| 中文亚洲av片在线观看爽| 亚洲五月天丁香| 亚洲av五月六月丁香网| 校园人妻丝袜中文字幕| 国产毛片a区久久久久| 一区二区三区免费毛片| 又紧又爽又黄一区二区| 午夜福利视频1000在线观看| av在线蜜桃| 亚洲av免费高清在线观看| 亚洲一区二区三区色噜噜| 男人和女人高潮做爰伦理| 身体一侧抽搐| 麻豆成人午夜福利视频| av天堂在线播放| 亚洲成a人片在线一区二区| 亚洲成人精品中文字幕电影| 老司机福利观看| 男人舔女人下体高潮全视频| 欧美高清性xxxxhd video| 一进一出好大好爽视频| 亚洲av成人精品一区久久| 精品乱码久久久久久99久播| 亚洲在线自拍视频| 国产黄a三级三级三级人| 午夜免费男女啪啪视频观看 | 久久精品综合一区二区三区| 色综合婷婷激情| 嫩草影院精品99| 国产伦一二天堂av在线观看| 女人十人毛片免费观看3o分钟| av福利片在线观看| 成人av在线播放网站| 欧美日本亚洲视频在线播放| 亚洲人成网站在线播| 精品久久久噜噜| 天堂网av新在线| 国产国拍精品亚洲av在线观看| 老师上课跳d突然被开到最大视频| 免费av毛片视频| 乱码一卡2卡4卡精品| 97超视频在线观看视频| 国产色婷婷99| 在线观看美女被高潮喷水网站| h日本视频在线播放| 夜夜爽天天搞| 日韩欧美 国产精品| 成年免费大片在线观看| 热99在线观看视频| 欧美潮喷喷水| 狠狠狠狠99中文字幕| 国产精品亚洲一级av第二区| 国产真实乱freesex| 亚洲成人精品中文字幕电影| 午夜久久久久精精品| 亚洲国产欧美人成| 老司机深夜福利视频在线观看| 国模一区二区三区四区视频| 国产中年淑女户外野战色| 亚洲欧美日韩无卡精品| 一本精品99久久精品77| 亚洲精品国产成人久久av| 欧美绝顶高潮抽搐喷水| 最后的刺客免费高清国语| 日韩人妻高清精品专区| 欧美日韩亚洲国产一区二区在线观看| 国产男人的电影天堂91| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 国产极品精品免费视频能看的| 一级黄片播放器| 国产主播在线观看一区二区| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放| 国产午夜福利久久久久久| 热99re8久久精品国产| 97人妻精品一区二区三区麻豆| 国产aⅴ精品一区二区三区波| 最近最新中文字幕大全电影3| 精品午夜福利在线看| 99久久精品一区二区三区| 亚洲七黄色美女视频| 日本黄色视频三级网站网址| 非洲黑人性xxxx精品又粗又长| 亚洲精品成人久久久久久| 观看免费一级毛片| 深爱激情五月婷婷| 欧美区成人在线视频| 99热6这里只有精品| 精品久久久久久久久亚洲 | 久久精品综合一区二区三区| 国产精品一区二区三区四区免费观看 | av专区在线播放| 美女免费视频网站| 国产精品久久久久久久久免| 久久精品影院6| 69av精品久久久久久| 国产免费男女视频| 亚洲av成人av| 直男gayav资源| 国产精品综合久久久久久久免费| 亚洲精品久久国产高清桃花| 精品国内亚洲2022精品成人| 我的女老师完整版在线观看| 99热这里只有是精品在线观看| 日韩欧美 国产精品| 俄罗斯特黄特色一大片| 免费av观看视频| 久久精品综合一区二区三区| 最近最新中文字幕大全电影3| 久久久久久伊人网av| 午夜精品久久久久久毛片777| 麻豆久久精品国产亚洲av| 免费观看的影片在线观看| 国产精品亚洲美女久久久| 偷拍熟女少妇极品色| 大又大粗又爽又黄少妇毛片口| 亚洲av美国av| 国产精品久久电影中文字幕| 久久亚洲精品不卡| 熟女人妻精品中文字幕| 免费无遮挡裸体视频| 啪啪无遮挡十八禁网站| 国产精品野战在线观看| 精品久久久久久久久av| 午夜久久久久精精品| 中文字幕熟女人妻在线| av福利片在线观看| 最近在线观看免费完整版| 日韩中字成人| 成人精品一区二区免费| 国产高清视频在线播放一区| 精品日产1卡2卡| www日本黄色视频网| 窝窝影院91人妻| 婷婷六月久久综合丁香| 一区二区三区四区激情视频 | 91久久精品国产一区二区三区| 少妇熟女aⅴ在线视频| 久久精品夜夜夜夜夜久久蜜豆| 午夜视频国产福利| 99热网站在线观看| netflix在线观看网站| 中文字幕熟女人妻在线| 亚洲国产高清在线一区二区三| 一本一本综合久久| 久久精品久久久久久噜噜老黄 | 午夜福利在线观看免费完整高清在 | 国产激情偷乱视频一区二区| 亚洲av熟女| 老熟妇乱子伦视频在线观看| 国产淫片久久久久久久久| 欧美日韩黄片免| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩综合久久久久久 | 中国美女看黄片| 日韩一本色道免费dvd| 97超级碰碰碰精品色视频在线观看| 香蕉av资源在线| 国产 一区 欧美 日韩| 一级毛片久久久久久久久女| 国产午夜精品论理片| 好男人在线观看高清免费视频| 网址你懂的国产日韩在线| 在线看三级毛片| 国产成年人精品一区二区| 内射极品少妇av片p| 精品欧美国产一区二区三| 日韩精品有码人妻一区| 成人欧美大片| 国产一区二区三区在线臀色熟女| 日本-黄色视频高清免费观看| 美女xxoo啪啪120秒动态图| 亚洲欧美激情综合另类| 色吧在线观看| 日日夜夜操网爽| 欧美成人a在线观看| 免费av不卡在线播放| 啦啦啦啦在线视频资源| 1000部很黄的大片| 在线播放国产精品三级| 男女之事视频高清在线观看| 国产精品一区二区三区四区免费观看 | 熟妇人妻久久中文字幕3abv| 一本久久中文字幕| 国产精品一区二区性色av| 一卡2卡三卡四卡精品乱码亚洲| 搡女人真爽免费视频火全软件 | 尤物成人国产欧美一区二区三区| 免费av毛片视频| 伊人久久精品亚洲午夜| 男女之事视频高清在线观看| 婷婷亚洲欧美| 一本一本综合久久| 欧美又色又爽又黄视频| 两人在一起打扑克的视频| 一a级毛片在线观看| 成年女人毛片免费观看观看9| 一进一出好大好爽视频| 成人高潮视频无遮挡免费网站| 国内精品久久久久久久电影|