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

    一種結(jié)合GVF和CV模型的水平集圖像分割方法

    2016-12-10 09:49:22胡小為劉宏申徐國雄
    關(guān)鍵詞:輪廓全局梯度

    胡小為,劉宏申,徐國雄,阮 越,劉 恒,潘 祥

    (安徽工業(yè)大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,安徽馬鞍山243032)

    一種結(jié)合GVF和CV模型的水平集圖像分割方法

    胡小為,劉宏申,徐國雄,阮 越,劉 恒,潘 祥

    (安徽工業(yè)大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,安徽馬鞍山243032)

    由于CV(Chan-Vese)模型是一個(gè)非凸性泛函,對該泛函求極值只能得到局部最優(yōu)解,運(yùn)用該模型進(jìn)行圖像分割時(shí),很難在全局范圍內(nèi)得到理想的結(jié)果。鑒于此,提出一種結(jié)合梯度矢量流(gradient vector flow,GVF)和CV模型的水平集圖像分割方法。該方法通過GVF將邊緣梯度信息擴(kuò)散至整幅圖像,在保留CV模型基本優(yōu)點(diǎn)的同時(shí),融入GVF的全局性梯度信息,從而引導(dǎo)CV模型在全局范圍內(nèi)演化至準(zhǔn)確的目標(biāo)邊緣。實(shí)驗(yàn)結(jié)果表明,該方法的分割效果和收斂速度均明顯優(yōu)于傳統(tǒng)CV模型。

    圖像分割;水平集方法;CV模型;梯度矢量流

    近20年以來,基于曲線演化理論和偏微分方程圖像分割方法取得了一系列的研究成果[1-4],尤其是結(jié)合Osher等[5]于1988年提出的水平集算法,為此類分割方法提供了有效的數(shù)值解方法。CV模型是一種利用水平集算法進(jìn)行數(shù)值解的圖像分割模型,由于其可以有效分割含噪聲和弱邊緣圖像,受到了國內(nèi)外學(xué)者的廣泛關(guān)注。

    MS模型是Mumford等[6]提出的一種基于區(qū)域信息的分段光滑分割模型。該模型將圖像分割和目標(biāo)平滑結(jié)合在一起,是一種較為理想的分割模型,但是該模型的數(shù)值解計(jì)算復(fù)雜。針對該問題,Chan等[7]提出一種基于精簡MS模型的分段常數(shù)分割模型,即CV模型。CV模型對圖像進(jìn)行二相位分割,分割過程中對目標(biāo)和背景的灰度值進(jìn)行常數(shù)項(xiàng)擬合,由于完全不依賴于圖像的邊緣梯度信息,該模型對含噪聲和弱邊緣圖像具有較好的分割效果,并且在數(shù)值解過程中引入了水平集算法,使得該模型的數(shù)值解計(jì)算方便。但是在實(shí)際應(yīng)用中,CV模型也有其自身的缺陷,主要表現(xiàn)在兩個(gè)方面:(1)該模型把圖像分為目標(biāo)和背景兩個(gè)區(qū)域,并進(jìn)行常數(shù)項(xiàng)擬合,得到二相位的分段常數(shù)圖像,因而不能直接分割多相圖像,該問題已有學(xué)者提出了一系列有效的解決方案[8-11];(2)由于CV模型是一個(gè)非凸性泛函,因此很難得到該泛函的全局最優(yōu)解。針對這類問題,文獻(xiàn)[12]中提出將ROF(Rudin-Osher-Fatemi)去噪模型和CV模型結(jié)合起來,并對能量泛函中的非凸特征函數(shù)進(jìn)行擴(kuò)展,由于ROF是全局凸模型,使得原本的非凸泛函極小化問題轉(zhuǎn)變?yōu)橥狗汉臉O小化問題,該方法的缺點(diǎn)是分割結(jié)果嚴(yán)重依賴初始條件的設(shè)定;文獻(xiàn)[13]中提出利用K-means算法在每一次迭代過程中計(jì)算輪廓的內(nèi)部和外部,利用全局性聚類結(jié)果不斷對CV模型的分割結(jié)果進(jìn)行修正;但是文獻(xiàn)[13]中也指出該方法只是優(yōu)化了CV模型的分割結(jié)果,并不能保證每一次得到全局性的分割結(jié)果。

    針對CV模型的第2個(gè)缺陷,本文提出一種結(jié)合GVF和CV模型的水平集圖像分割模型,首先利用GVF將局部的邊緣梯度擴(kuò)散為全局的梯度矢量場,然后將該全局梯度信息和傳統(tǒng)CV模型結(jié)合起來,并納入能量泛函,從而引導(dǎo)CV模型進(jìn)行全局性的演化。由于添加了全局梯度信息,使得水平集函數(shù)的演化速度大大提高,從而有效減少CV模型的收斂時(shí)間和迭代次數(shù)。

    1 CV模型的原理及分析

    設(shè)Ω為待分割圖像區(qū)域,u0(x,y)為原始圖像(x,y)處的灰度值,C為自定義的活動(dòng)輪廓曲線,Ω1為曲線內(nèi)部區(qū)域,Ω2為曲線外部區(qū)域,c1為曲線內(nèi)部像素灰度的平均值,c2為曲線外部像素灰度的平均值,CV模型定義為如下泛函

    其中:μ為控制長度項(xiàng)的權(quán)重系數(shù);ν為控制面積項(xiàng)的權(quán)重系數(shù);λ1和λ2為控制輪廓內(nèi)部和外部灰度擬合項(xiàng)的權(quán)重系數(shù)。該泛函的第1項(xiàng)為長度項(xiàng),L(C)表示活動(dòng)輪廓的長度,要求活動(dòng)輪廓C的長度盡可能短;第2項(xiàng)為面積項(xiàng),S(C)表示活動(dòng)輪廓的內(nèi)部面積,要求活動(dòng)輪廓C的內(nèi)部面積盡可能小。第1項(xiàng)和第2項(xiàng)的作用是使活動(dòng)輪廓C盡可能收縮且平滑;第3項(xiàng)和第4項(xiàng)為輪廓內(nèi)部和外部的灰度擬合項(xiàng),要求分割結(jié)果為分段常數(shù)圖像且該分段常數(shù)圖像與源圖像盡可能相似。通過求該泛函的極小解,即可得到滿足該泛函條件的最終活動(dòng)輪廓C,再對C的內(nèi)部和外部求像素平均值,就可以得到與原圖像對應(yīng)的分段常數(shù)圖像,即CV模型的分割結(jié)果。

    水平集是一種隱式表達(dá)曲線演化問題的數(shù)值解方法[5],其基本原理是將演化曲線表達(dá)為水平集函數(shù)的零水平集,從而實(shí)現(xiàn)由低維平面的曲線演化問題轉(zhuǎn)化為高維空間的曲面演化問題。與曲線演化相比,水平集演化可以很自然地表達(dá)曲線在演化過程中的拓?fù)浣Y(jié)構(gòu)變化。先定義Heaviside函數(shù)H(z)及其Dirac測度δ(z)

    利用該Heaviside函數(shù)及其Dirac測度可以將能量泛函(1)改寫為以水平集函數(shù)φ表示的形式

    其中:Ω為圖像區(qū)域;?為梯度運(yùn)算符。將Heaviside函數(shù)及其Dirac測度正則化為

    引入時(shí)間變量t并利用變分法可以得到式(3)的梯度下降流為

    其中div為散度運(yùn)算符。利用有限差分法可對式(5)進(jìn)行離散化并求其數(shù)值解。

    針對CV模型分割結(jié)果嚴(yán)重依賴初始活動(dòng)輪廓位置的問題,李俊等在文獻(xiàn)[14]中提出將式(5)中的δε(φ)替換成|?φ(x,y)|,從而其梯度下降流變?yōu)?/p>

    這一替換在一定程度上提高了水平集函數(shù)的演化速度和迭代步長,擴(kuò)大了活動(dòng)輪廓的“捕獲”范圍,但實(shí)驗(yàn)結(jié)果表明,這種優(yōu)化并非是全局性的,如果目標(biāo)邊緣與活動(dòng)輪廓距離過遠(yuǎn),則該方法并不能進(jìn)行有效的全局分割(圖2(c))。

    2 CV模型的改進(jìn)

    從曲線演化的角度看,傳統(tǒng)CV模型不能進(jìn)行全局性分割的原因在于模型本身缺少表示全局信息的變量,每一次演化過程都局限于局部范圍。梯度矢量流[15](gradient vector flow,GVF)是一種對目標(biāo)邊緣梯度進(jìn)行全局?jǐn)U散而產(chǎn)生的靜態(tài)矢量場,它的分布僅依賴于原始圖像目標(biāo)邊緣,而與曲線演化過程中的其他變量無關(guān),可用于解決參數(shù)活動(dòng)輪廓模型無法收斂至深度凹陷邊緣的問題。

    2.1 梯度矢量流

    為了得到GVF,通過式(7)對原始圖像u0(x,y)進(jìn)行邊緣映射,得邊緣映射圖像f(x,y)

    其中:Gσ(x,y)是標(biāo)準(zhǔn)差為σ的二維高斯函數(shù);?為卷積運(yùn)算符。設(shè)梯度矢量流V(x,y)=[m(x,y),n(x,y)],其能量泛函為

    其中η為規(guī)則化參數(shù),η值根據(jù)圖像中噪聲強(qiáng)弱選擇,噪聲強(qiáng)烈增加η值,反之則減少η值??梢钥闯?,在圖像的勻質(zhì)區(qū)域(|?f|值較小),該能量泛函的值主要由第1項(xiàng)決定,產(chǎn)生的擴(kuò)散矢量場變化緩慢。在目標(biāo)邊緣處(|?f|值較大),該能量泛函的值主要由第2項(xiàng)決定,梯度矢量流V(x,y)=?f,利用變分法,求得該泛函的Euler-Lagrange方程組

    其中?2為拉普拉斯算子。最后利用有限差分法對式(9)進(jìn)行離散化,得到數(shù)值解,具體過程見文獻(xiàn)[15]。

    2.2 結(jié)合GVF的改進(jìn)CV模型

    針對傳統(tǒng)CV模型存在的局部最優(yōu)解問題,結(jié)合GVF的全局性信息,本文提出一種改進(jìn)的CV模型

    此式在式(1)的基礎(chǔ)上添加最后一項(xiàng),前4項(xiàng)的作用與式(1)中相同。其中:,V(x,y)為圖像的梯度矢量流;γ為梯度吸引項(xiàng)的權(quán)重系數(shù)。添加項(xiàng)為邊緣梯度吸引項(xiàng),利用GVF將邊緣梯度信息擴(kuò)散至全局范圍,其主要作用是引導(dǎo)活動(dòng)輪廓收斂至邊緣梯度值最大處,即目標(biāo)邊緣。在圖像的非目標(biāo)邊緣區(qū)域,|V(x,y)|值較小,函數(shù)項(xiàng)g值較大;在圖像的目標(biāo)邊緣,|V(x,y)|值較大,函數(shù)項(xiàng)g值較小。運(yùn)用該模型進(jìn)行圖像分割時(shí),收斂過程朝著泛函極小值的方向進(jìn)行,最終停留在g值較小(|V(x,y)|值較大)的目標(biāo)邊緣處,由于該邊緣梯度信息是全局性的,所以厚壁內(nèi)部邊緣、距離活動(dòng)輪廓較遠(yuǎn)的目標(biāo)邊緣均可對活動(dòng)輪廓的演化產(chǎn)生吸引作用,從而引導(dǎo)活動(dòng)輪廓收斂至全局最優(yōu)解。

    同樣對式(10)采用變分水平集方法引入水平集函數(shù)φ和時(shí)間變量t,利用變分法得到梯度下降流

    2.3 改進(jìn)模型的數(shù)值解

    能量泛函的梯度下降流本質(zhì)是一個(gè)偏微分方程,求其解析解十分困難,因此本文利用有限差分法來離散化能量泛函的梯度下降流,求其數(shù)值解。可以用離散網(wǎng)格來表達(dá)水平集函數(shù)φ(x,y,t),設(shè)h為離散網(wǎng)格的空間間隔,Δt為迭代的時(shí)間間隔,則在n時(shí)刻,網(wǎng)格點(diǎn)(i,j)處的水平集函數(shù)。將初始水平集函數(shù)φ定義為活動(dòng)輪廓C的符號距離函數(shù)(signed distance function,SDF),即φ(x,y,t=0)=d(x,y),其中d(x,y)表示圖像中點(diǎn)(x,y)到活動(dòng)輪廓C的符號距離,在C內(nèi)部取正值,在C外部取負(fù)值,在C上取零值。記差分:

    根據(jù)文獻(xiàn)[16-18],為了避免水平集函數(shù)的重新初始化和保持水平集函數(shù)的穩(wěn)定性,選用有限差分半隱式方案進(jìn)行偏微分方程的數(shù)值化計(jì)算,得離散化的迭代式(12)。

    初始化水平集函數(shù)φ和參數(shù) μ,ν,λ1,λ2,γ,Δt及h,對上式進(jìn)行迭代計(jì)算直至分割結(jié)果趨于穩(wěn)定,此時(shí)的活動(dòng)輪廓C={(x,y)|φ(x,y,t)=0}即為最終的目標(biāo)輪廓。

    3 實(shí)驗(yàn)結(jié)果及分析

    為了驗(yàn)證上述改進(jìn)模型的可行性和有效性,通過實(shí)驗(yàn)對本文改進(jìn)模型(記作GVF-CV)與傳統(tǒng)CV模型及李俊改進(jìn)的CV模型(記作Li-CV)在收斂時(shí)間、迭代次數(shù)以及收斂結(jié)果3個(gè)方面進(jìn)行比較。實(shí)驗(yàn)在CPU為Core 2.13 GHz、內(nèi)存為2G的硬件平臺上利用VC++6.0編程實(shí)現(xiàn),參數(shù)取 μ=0.2,ν=0.2,λ1=λ2=1,γ=0.2,Δt=0.1,h=1。以3組具有代表性的實(shí)驗(yàn)來分析3種模型的性能。

    實(shí)驗(yàn)1針對含噪聲厚壁環(huán)狀圖像,分割結(jié)果如圖1。圖1中,初始輪廓在目標(biāo)外圍。使用傳統(tǒng)CV模型無法分割出厚壁內(nèi)部邊緣,Li-CV模型和GVF-CV模型均可實(shí)現(xiàn)有效分割,但是GVF-CV模型在收斂時(shí)間和迭代次數(shù)2個(gè)方面要優(yōu)于Li-CV模型,見表1。

    圖1 含噪聲厚壁環(huán)狀圖像分割結(jié)果Fig.1 Segmentation results of thick walled ring image with noise

    實(shí)驗(yàn)2針對含噪聲多目標(biāo)圖像,分割結(jié)果如圖2。圖2中,初始輪廓包含部分目標(biāo)。使用傳統(tǒng)CV模型和Li-CV模型均只能分割部分目標(biāo),而使用GVF-CV模型可以實(shí)現(xiàn)有效的全局分割,且收斂時(shí)間和迭代次數(shù)均優(yōu)于傳統(tǒng)CV模型和Li-CV模型,見表1。

    圖2 含噪聲多目標(biāo)圖像分割結(jié)果Fig.2 Segmentation results of multi-target image with noise

    實(shí)驗(yàn)3針對無噪聲多目標(biāo)圖像,分割結(jié)果如圖3。圖3中,初始輪廓在兩個(gè)目標(biāo)之間。由于初始輪廓距離內(nèi)部目標(biāo)較遠(yuǎn),傳統(tǒng)CV模型和Li-CV模型均只能分割外部目標(biāo),而使用GVF-CV模型分割的結(jié)果要優(yōu)于傳統(tǒng)CV模型和Li-CV模型,見表1。

    圖3 無噪聲多目標(biāo)圖像分割結(jié)果Fig.3 Segmentation results of multi-target image without noise

    由圖1~3可以看出,模型的改進(jìn)對不同狀況圖像分割結(jié)果均有明顯改善。從表1中可以看出,本文的GVF-CV模型的收斂時(shí)間、迭代次數(shù)以及收斂結(jié)果均優(yōu)于傳統(tǒng)CV模型和Li-CV模型。

    表1 3種模型的性能比較Tab.1 Performance comparison of three models

    4 結(jié) 語

    傳統(tǒng)CV模型是非凸性泛函,其數(shù)值解難以在全局范圍內(nèi)收斂至最優(yōu)。針對這個(gè)缺陷,本文將GVF和傳統(tǒng)CV模型結(jié)合起來,利用GVF將局部邊緣梯度信息擴(kuò)散至全局區(qū)域,并在CV模型中添加邊緣梯度吸引項(xiàng),從而引導(dǎo)活動(dòng)輪廓收斂至正確的目標(biāo)邊緣。實(shí)驗(yàn)結(jié)果表明,本文方法在收斂速度、迭代次數(shù)及分割效果等方面均明顯優(yōu)于傳統(tǒng)CV模型分割方法,其保留了傳統(tǒng)CV模型可以有效分割含噪聲圖像的優(yōu)點(diǎn),且提高了對多目標(biāo)圖像的分割能力。本文改進(jìn)CV模型的不足之處在于:如果圖像噪聲過于強(qiáng)烈且目標(biāo)邊緣過于虛弱,則會影響GVF的分布,從而降低該模型的分割效果。這是后續(xù)工作將要研究解決的一個(gè)問題。

    [1]WEICKERT J.Efficient image segmentation using partial differential equations and morphology[J].Pattern Recognition,2001, 34(9)∶1813-1824.

    [2]劉利雄,馬忠梅,趙恒博.一種基于主動(dòng)輪廓模型的心臟核磁共振圖像分割方法[J].計(jì)算機(jī)學(xué)報(bào),2012,35(1)∶146-153.

    [3]KAUR J,AGRAWAL S,VIG R.Integration of clustering,optimization and partial differential equation method for improved image segmentation[J].International Journal of Image,Graphics and Signal Processing(IJIGSP),2012,4(11)∶26-33.

    [4]NIANG O,THIOUNE A,GUEIREA M C E.Partial differential equation-based approach for empirical mode decomposition∶application on image analysis[J].Image Processing,IEEE Transactions on,2012,21(9)∶3991-4001.

    [5]OSHER S,SETHIAN JA.Fronts propagating with curvature dependent speed∶algorithms based on Hamilton-Jacobi formulations [J].Journal of Computational Physics,1988,79(1)∶12-49.

    [6]MUMFORD D,SHAH J.Optimal approximations by piecewise smooth functions and associated variational problems[J].Communications on Pure andApplied Mathematics,1989,42(5)∶577-685.

    [7]CHAN T F,VESE LA.Active contours without edges[J].Image Processing,IEEE Transactions on,2001,10(2)∶266-277.

    [8]VESE L A,CHAN T F.A multiphase level set framework for image segmentation using the Mumford and Shah model[J].International Journal of Computer Vision,2002,50(3)∶271-293.

    [9]GAO S,BUI T D.Image segmentation and selective smoothing by using Mumford-Shah model[J].Image Processing,IEEE Transactions on,2005,14(10)∶1537-1549.

    [10]CHUNG G,VESE L A.Image segmentation using a multilayer level-set approach[J].Computing and Visualization in Science, 2009,12(6)∶267-285.

    [11]LIU S,PENGY.Alocal region-based Chan-Vese model for image segmentation[J].Pattern Recognition,2012,45(7)∶2769-2779.

    [12]CHAN T F,ESEDOGLU S,NIKOLOVA M.Algorithms for finding global minimizers of image segmentation and denoising models[J].SIAM Journal onApplied Mathematics,2006,66(5)∶1632-1648.

    [13]BROWN E S,CHAN T F,BRESSON X.Completely convex formulation of the Chan-Vese image segmentation model[J].International Journal of Computer Vision,2012,98(1)∶103-121.

    [14]李俊,楊新.基于Mumford-Shah模型的快速水平集圖像分割方法[J].計(jì)算機(jī)學(xué)報(bào),2002,25(11)∶1175-1183.

    [15]XU C,PRINCE J L.Snakes,shapes,and gradient vector flow[J].Image Processing,IEEE Transactions on,1998,7(3)∶359-369.

    [16]LI C,XU C,GUI C.Level set evolution without re-initialization∶a new variational formulation[C]//2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition.[s.l.]∶IEEE,2005,1∶430-436.

    [17]SMEREKA P.Semi-implicit level set methods for curvature and surface diffusion motion[J].Journal of Scientific Computing, 2003,19(1)∶439-456.

    [18]KOLAHDOUZ E M,SALAC D.A semi-implicit gradient augmented level set method[J].SIAM Journal on Scientific Computing,2013,35(1)∶A231-A254.

    責(zé)任編輯:丁吉海

    ALevel Set Image Segmentation Method Combined with GVF and CV Model

    HU Xiaowei,LIU Hongshen,XU Guoxiong,RUAN Yue,LIU Heng,PAN Xiang
    (School of Computer Science and Technology,Anhui University of Technology,Ma'anshan 243032,China)

    Owning to the non-convex functional of with the Chan-Vese(CV)model,one can only obtain a local optimal solution.It is difficult to achieve an ideal result for image segmentation in the global range.Therefore a new level set based image segmentation method that combining CV model and gradient vector flow(GVF)was proposed.The edge gradient information is spreaded to the entire image with GVF,which guides the evolution of CV model to the correcttarget edge in the global range and retains the basic advantages of CV model.The experimental results indicate that the present method are obviously better than the traditional CV model.

    image segmentation;level set method;CV model;gradient vector flow

    TP 391.41

    A

    10.3969/j.issn.1671-7872.2016.03.017

    2015-11-10

    安徽省自然科學(xué)基金項(xiàng)目(1308085QF113)

    胡小為(1989-),男,安徽安慶人,碩士生,研究方向?yàn)閳D像處理。

    劉宏申(1962-),男,安徽宣城人,教授,研究方向?yàn)閳D像處理與模式識別。

    1671-7872(2016)03-0289-06

    猜你喜歡
    輪廓全局梯度
    Cahn-Hilliard-Brinkman系統(tǒng)的全局吸引子
    量子Navier-Stokes方程弱解的全局存在性
    一個(gè)改進(jìn)的WYL型三項(xiàng)共軛梯度法
    OPENCV輪廓識別研究與實(shí)踐
    一種自適應(yīng)Dai-Liao共軛梯度法
    基于實(shí)時(shí)輪廓誤差估算的數(shù)控系統(tǒng)輪廓控制
    一類扭積形式的梯度近Ricci孤立子
    落子山東,意在全局
    金橋(2018年4期)2018-09-26 02:24:54
    在線學(xué)習(xí)機(jī)制下的Snake輪廓跟蹤
    新思路:牽一發(fā)動(dòng)全局
    午夜激情福利司机影院| 国产精品野战在线观看| 亚洲中文字幕日韩| 在线视频色国产色| 国产精品日韩av在线免费观看| 国产高清有码在线观看视频 | 可以在线观看毛片的网站| 亚洲欧美日韩东京热| 18禁观看日本| 国产私拍福利视频在线观看| 午夜a级毛片| 亚洲中文av在线| cao死你这个sao货| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久久久久| 欧美绝顶高潮抽搐喷水| 亚洲自拍偷在线| 床上黄色一级片| 99riav亚洲国产免费| 成人国语在线视频| 成人永久免费在线观看视频| 亚洲真实伦在线观看| 美女高潮喷水抽搐中文字幕| 十八禁人妻一区二区| 欧美色视频一区免费| 法律面前人人平等表现在哪些方面| 国产av麻豆久久久久久久| 757午夜福利合集在线观看| 全区人妻精品视频| 欧美午夜高清在线| 久久久久国内视频| 亚洲欧洲精品一区二区精品久久久| 搞女人的毛片| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄色小视频在线观看| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 国产精品久久久久久久电影 | 人人妻人人澡欧美一区二区| 亚洲成人国产一区在线观看| 日韩高清综合在线| 中文字幕最新亚洲高清| 欧美日韩精品网址| 欧美黑人精品巨大| a在线观看视频网站| 欧美中文综合在线视频| 成人永久免费在线观看视频| 日本 av在线| 亚洲全国av大片| 男女床上黄色一级片免费看| 在线观看日韩欧美| 啦啦啦观看免费观看视频高清| 韩国av一区二区三区四区| av天堂在线播放| 亚洲精品国产精品久久久不卡| 欧美在线一区亚洲| 欧美另类亚洲清纯唯美| xxx96com| 久热爱精品视频在线9| 亚洲精品久久成人aⅴ小说| 两个人视频免费观看高清| 国产69精品久久久久777片 | av在线播放免费不卡| 91av网站免费观看| 婷婷丁香在线五月| av免费在线观看网站| 亚洲精品久久成人aⅴ小说| 亚洲aⅴ乱码一区二区在线播放 | 在线免费观看的www视频| 叶爱在线成人免费视频播放| 久久精品91无色码中文字幕| 欧美+亚洲+日韩+国产| 搡老岳熟女国产| 国产精品精品国产色婷婷| 在线观看免费午夜福利视频| 亚洲欧洲精品一区二区精品久久久| 国产伦一二天堂av在线观看| 亚洲 国产 在线| 国产不卡一卡二| 亚洲熟妇熟女久久| or卡值多少钱| 国产精品爽爽va在线观看网站| 叶爱在线成人免费视频播放| 欧美日韩福利视频一区二区| 高清毛片免费观看视频网站| 美女 人体艺术 gogo| 美女扒开内裤让男人捅视频| 精品久久久久久久久久免费视频| 成人精品一区二区免费| 亚洲午夜理论影院| 国产99白浆流出| 999久久久精品免费观看国产| 身体一侧抽搐| 男女午夜视频在线观看| 亚洲片人在线观看| 久久精品国产清高在天天线| 亚洲无线在线观看| 免费看十八禁软件| 亚洲成a人片在线一区二区| 亚洲成a人片在线一区二区| 国产欧美日韩一区二区三| 欧美黑人欧美精品刺激| 亚洲人成77777在线视频| 精品久久久久久久末码| www.www免费av| 精品一区二区三区四区五区乱码| 精品久久久久久成人av| 免费在线观看日本一区| 国产免费男女视频| 国产亚洲精品久久久久久毛片| 国产精品综合久久久久久久免费| 亚洲精品在线观看二区| 国产精品久久久av美女十八| 日韩欧美三级三区| 日本五十路高清| 久久精品影院6| 亚洲欧美一区二区三区黑人| 久久这里只有精品19| 中文字幕av在线有码专区| 此物有八面人人有两片| 91成年电影在线观看| xxxwww97欧美| 老司机在亚洲福利影院| 九色成人免费人妻av| 国产成人aa在线观看| 亚洲午夜理论影院| 亚洲激情在线av| 男女视频在线观看网站免费 | 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 精品久久久久久久久久久久久| 欧美日韩精品网址| 999久久久国产精品视频| 正在播放国产对白刺激| 精品久久久久久久人妻蜜臀av| 亚洲自偷自拍图片 自拍| 精品久久久久久久久久免费视频| 欧美精品亚洲一区二区| 亚洲片人在线观看| 露出奶头的视频| 1024视频免费在线观看| 久99久视频精品免费| 99热这里只有精品一区 | АⅤ资源中文在线天堂| 精品国产亚洲在线| 国产伦在线观看视频一区| 悠悠久久av| 亚洲精品美女久久av网站| 亚洲黑人精品在线| 首页视频小说图片口味搜索| 精品国内亚洲2022精品成人| 国产亚洲精品久久久久久毛片| av中文乱码字幕在线| 久久久久久久午夜电影| 1024视频免费在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 狂野欧美激情性xxxx| 看黄色毛片网站| 变态另类成人亚洲欧美熟女| 欧美日本亚洲视频在线播放| 亚洲av美国av| 中文字幕久久专区| 人妻丰满熟妇av一区二区三区| 日韩国内少妇激情av| 黄色 视频免费看| 亚洲美女黄片视频| av有码第一页| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品久久男人天堂| 好看av亚洲va欧美ⅴa在| 国产三级在线视频| 午夜免费激情av| 两个人看的免费小视频| 无遮挡黄片免费观看| 欧美成人一区二区免费高清观看 | 亚洲专区字幕在线| 男女床上黄色一级片免费看| 欧美3d第一页| 老汉色∧v一级毛片| or卡值多少钱| 亚洲在线自拍视频| 国产精品久久久av美女十八| 麻豆国产av国片精品| 午夜精品一区二区三区免费看| 亚洲熟女毛片儿| 视频区欧美日本亚洲| 精品久久久久久成人av| 精品人妻1区二区| 亚洲国产精品成人综合色| 亚洲人成网站高清观看| 巨乳人妻的诱惑在线观看| 一本久久中文字幕| 亚洲欧美日韩东京热| 亚洲国产欧洲综合997久久,| 中出人妻视频一区二区| 亚洲电影在线观看av| 亚洲狠狠婷婷综合久久图片| 正在播放国产对白刺激| 美女免费视频网站| 日韩有码中文字幕| 国语自产精品视频在线第100页| 亚洲激情在线av| 91国产中文字幕| 精品一区二区三区四区五区乱码| av有码第一页| 亚洲一码二码三码区别大吗| 免费在线观看亚洲国产| 黑人欧美特级aaaaaa片| 久久久久久久久中文| 亚洲精品在线美女| 国产免费av片在线观看野外av| 成在线人永久免费视频| 嫁个100分男人电影在线观看| 99热6这里只有精品| 精品国内亚洲2022精品成人| 亚洲午夜精品一区,二区,三区| 亚洲精品一区av在线观看| 床上黄色一级片| 久9热在线精品视频| 国产精品永久免费网站| 手机成人av网站| 伊人久久大香线蕉亚洲五| 久久久久久久午夜电影| 首页视频小说图片口味搜索| 亚洲最大成人中文| 国产精品日韩av在线免费观看| 人妻久久中文字幕网| 欧美av亚洲av综合av国产av| 精品电影一区二区在线| 毛片女人毛片| 超碰成人久久| 国产亚洲欧美在线一区二区| 亚洲天堂国产精品一区在线| 国产成人系列免费观看| 在线观看日韩欧美| 熟女电影av网| 亚洲一码二码三码区别大吗| 亚洲专区字幕在线| 国产黄a三级三级三级人| 欧美极品一区二区三区四区| 老熟妇仑乱视频hdxx| 亚洲国产欧美一区二区综合| 成人高潮视频无遮挡免费网站| 国产欧美日韩精品亚洲av| 精品日产1卡2卡| 亚洲一区二区三区色噜噜| 搡老妇女老女人老熟妇| 男男h啪啪无遮挡| 宅男免费午夜| 午夜成年电影在线免费观看| 神马国产精品三级电影在线观看 | 免费观看精品视频网站| 草草在线视频免费看| 高清毛片免费观看视频网站| 亚洲国产日韩欧美精品在线观看 | 午夜影院日韩av| a在线观看视频网站| 在线十欧美十亚洲十日本专区| 天堂√8在线中文| svipshipincom国产片| 久久久久亚洲av毛片大全| 毛片女人毛片| 国产黄片美女视频| 国产精品av久久久久免费| 两个人看的免费小视频| 欧美午夜高清在线| 日本在线视频免费播放| 久久精品夜夜夜夜夜久久蜜豆 | 国产精品久久久久久久电影 | 国内毛片毛片毛片毛片毛片| 50天的宝宝边吃奶边哭怎么回事| 婷婷六月久久综合丁香| 免费在线观看亚洲国产| 嫩草影视91久久| 国模一区二区三区四区视频 | 亚洲片人在线观看| 人妻久久中文字幕网| 午夜免费成人在线视频| 亚洲人成电影免费在线| 久久香蕉精品热| 久久中文字幕一级| 真人做人爱边吃奶动态| 五月玫瑰六月丁香| 久久伊人香网站| 999久久久国产精品视频| cao死你这个sao货| 特大巨黑吊av在线直播| 99久久无色码亚洲精品果冻| 国产主播在线观看一区二区| 老司机在亚洲福利影院| 村上凉子中文字幕在线| 精品国产乱码久久久久久男人| 老司机福利观看| 亚洲av片天天在线观看| 无人区码免费观看不卡| 欧美日本视频| 久久国产乱子伦精品免费另类| 国产三级黄色录像| 国产日本99.免费观看| 久久精品亚洲精品国产色婷小说| 97碰自拍视频| 99国产综合亚洲精品| 99久久精品国产亚洲精品| 9191精品国产免费久久| 欧美av亚洲av综合av国产av| 在线观看美女被高潮喷水网站 | 国产成人aa在线观看| 窝窝影院91人妻| 人人妻人人澡欧美一区二区| 人成视频在线观看免费观看| 高潮久久久久久久久久久不卡| 午夜a级毛片| 桃色一区二区三区在线观看| 男女午夜视频在线观看| 最新在线观看一区二区三区| 人人妻人人澡欧美一区二区| 亚洲成a人片在线一区二区| 91字幕亚洲| 后天国语完整版免费观看| 久热爱精品视频在线9| 欧美在线一区亚洲| 国产乱人伦免费视频| 久久精品国产亚洲av高清一级| 亚洲av成人一区二区三| 淫秽高清视频在线观看| 波多野结衣巨乳人妻| 亚洲18禁久久av| 成人18禁在线播放| 亚洲一区二区三区色噜噜| av有码第一页| 亚洲第一电影网av| 99国产极品粉嫩在线观看| 欧美色视频一区免费| 在线观看www视频免费| 18禁裸乳无遮挡免费网站照片| 国产精品国产高清国产av| 精品国产乱码久久久久久男人| 亚洲中文日韩欧美视频| av在线播放免费不卡| 欧美色视频一区免费| 岛国在线免费视频观看| 亚洲人与动物交配视频| 特级一级黄色大片| 国产主播在线观看一区二区| 丝袜人妻中文字幕| 黑人操中国人逼视频| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 免费在线观看日本一区| 亚洲成人精品中文字幕电影| 欧美成人一区二区免费高清观看 | 欧美最黄视频在线播放免费| 一级a爱片免费观看的视频| 91国产中文字幕| 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 首页视频小说图片口味搜索| 久久精品91蜜桃| 国产精品av久久久久免费| 亚洲天堂国产精品一区在线| 亚洲精品一区av在线观看| 日本一二三区视频观看| 久久精品亚洲精品国产色婷小说| 亚洲精品美女久久久久99蜜臀| 久久草成人影院| 动漫黄色视频在线观看| 欧美色视频一区免费| 人妻夜夜爽99麻豆av| 青草久久国产| 亚洲国产高清在线一区二区三| 桃红色精品国产亚洲av| 岛国在线观看网站| 久久久精品欧美日韩精品| 免费在线观看日本一区| 天堂√8在线中文| 欧美极品一区二区三区四区| 无遮挡黄片免费观看| 可以免费在线观看a视频的电影网站| 亚洲国产精品sss在线观看| 成年女人毛片免费观看观看9| 搡老岳熟女国产| 成人亚洲精品av一区二区| 男人舔女人下体高潮全视频| 天堂影院成人在线观看| 久久这里只有精品中国| 亚洲成a人片在线一区二区| 999精品在线视频| 丝袜人妻中文字幕| 少妇熟女aⅴ在线视频| 久久久久久国产a免费观看| 欧美国产日韩亚洲一区| 在线免费观看的www视频| 99热这里只有是精品50| 久久精品影院6| 国产精品野战在线观看| 叶爱在线成人免费视频播放| 日韩精品中文字幕看吧| 国模一区二区三区四区视频 | 亚洲精品色激情综合| 国产aⅴ精品一区二区三区波| 亚洲真实伦在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲av五月六月丁香网| av有码第一页| 久久久久免费精品人妻一区二区| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| 亚洲精品久久国产高清桃花| 一区二区三区激情视频| 后天国语完整版免费观看| 精品电影一区二区在线| 桃色一区二区三区在线观看| 日韩精品免费视频一区二区三区| 亚洲免费av在线视频| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 成人高潮视频无遮挡免费网站| 99国产极品粉嫩在线观看| 男女视频在线观看网站免费 | 一进一出抽搐gif免费好疼| 一夜夜www| 免费在线观看日本一区| 久久精品国产清高在天天线| 精品日产1卡2卡| 亚洲精品在线观看二区| 1024手机看黄色片| 又黄又粗又硬又大视频| 波多野结衣高清作品| 欧美大码av| 国产一区二区在线观看日韩 | 国产亚洲精品久久久久久毛片| 亚洲欧美一区二区三区黑人| 亚洲aⅴ乱码一区二区在线播放 | 亚洲成av人片在线播放无| 国产黄a三级三级三级人| 精品一区二区三区av网在线观看| 美女黄网站色视频| 亚洲天堂国产精品一区在线| 国产真人三级小视频在线观看| 国产成人aa在线观看| 一区二区三区高清视频在线| 青草久久国产| 国产精品久久久久久久电影 | 在线观看一区二区三区| av福利片在线| 久久久久久九九精品二区国产 | 欧美在线一区亚洲| 欧美乱码精品一区二区三区| 90打野战视频偷拍视频| 波多野结衣巨乳人妻| 久久精品影院6| 午夜精品一区二区三区免费看| 老汉色av国产亚洲站长工具| 成人18禁高潮啪啪吃奶动态图| 18禁国产床啪视频网站| 丰满人妻熟妇乱又伦精品不卡| 美女午夜性视频免费| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频| 十八禁网站免费在线| 50天的宝宝边吃奶边哭怎么回事| 国产成年人精品一区二区| 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| e午夜精品久久久久久久| 不卡一级毛片| 久久香蕉激情| 国产高清视频在线播放一区| 国产欧美日韩一区二区精品| 国产精品一区二区三区四区免费观看 | tocl精华| 动漫黄色视频在线观看| 日本一二三区视频观看| 日本成人三级电影网站| 黄频高清免费视频| netflix在线观看网站| 久久香蕉国产精品| 午夜精品一区二区三区免费看| 久久天堂一区二区三区四区| 欧美中文综合在线视频| 日本 av在线| √禁漫天堂资源中文www| 久久国产乱子伦精品免费另类| 亚洲一区高清亚洲精品| 精品电影一区二区在线| 久久香蕉激情| 一区二区三区高清视频在线| 99在线视频只有这里精品首页| 久久精品91蜜桃| 一本久久中文字幕| 在线观看免费视频日本深夜| av视频在线观看入口| 熟妇人妻久久中文字幕3abv| videosex国产| 亚洲专区字幕在线| 免费看美女性在线毛片视频| 老司机深夜福利视频在线观看| 人成视频在线观看免费观看| 女同久久另类99精品国产91| 久久久久九九精品影院| 麻豆成人av在线观看| 制服诱惑二区| 精品电影一区二区在线| 国产成人aa在线观看| 亚洲专区中文字幕在线| 欧美日韩一级在线毛片| 欧美一区二区精品小视频在线| 欧美日韩一级在线毛片| 男人舔女人的私密视频| aaaaa片日本免费| 国产野战对白在线观看| 欧美丝袜亚洲另类 | 国产1区2区3区精品| 19禁男女啪啪无遮挡网站| 成人一区二区视频在线观看| 久久久国产精品麻豆| 午夜日韩欧美国产| 国产精品自产拍在线观看55亚洲| 欧美乱码精品一区二区三区| 成在线人永久免费视频| 9191精品国产免费久久| 黄片小视频在线播放| 欧美久久黑人一区二区| 91成年电影在线观看| av免费在线观看网站| 国产亚洲av嫩草精品影院| 亚洲欧美精品综合一区二区三区| 69av精品久久久久久| 动漫黄色视频在线观看| x7x7x7水蜜桃| 一进一出好大好爽视频| 欧美乱码精品一区二区三区| 在线观看免费日韩欧美大片| 此物有八面人人有两片| 精品久久久久久成人av| 免费在线观看完整版高清| 中文字幕高清在线视频| 老鸭窝网址在线观看| 久久久久久九九精品二区国产 | 亚洲精品美女久久久久99蜜臀| 欧美 亚洲 国产 日韩一| tocl精华| 男人舔女人下体高潮全视频| 天堂av国产一区二区熟女人妻 | 亚洲自偷自拍图片 自拍| 亚洲成av人片免费观看| 精华霜和精华液先用哪个| 精品久久久久久久久久免费视频| 日本a在线网址| 欧美日韩黄片免| 国产成人欧美在线观看| 深夜精品福利| www日本在线高清视频| 国产真实乱freesex| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲美女黄片视频| 国产成人av激情在线播放| 国产高清有码在线观看视频 | 毛片女人毛片| 亚洲成人中文字幕在线播放| 成人欧美大片| 日本 欧美在线| 久久欧美精品欧美久久欧美| 最近最新中文字幕大全电影3| 听说在线观看完整版免费高清| 成人特级黄色片久久久久久久| 国产亚洲欧美98| 在线免费观看的www视频| 老司机深夜福利视频在线观看| 久久久国产成人免费| 久久亚洲精品不卡| 中文字幕高清在线视频| 久久久久九九精品影院| 一级a爱片免费观看的视频| 日日爽夜夜爽网站| 一个人观看的视频www高清免费观看 | 国产精品1区2区在线观看.| 国内精品久久久久久久电影| 亚洲激情在线av| 熟妇人妻久久中文字幕3abv| 午夜精品久久久久久毛片777| 国产不卡一卡二| 色综合婷婷激情| 久久性视频一级片| 国产不卡一卡二| 真人做人爱边吃奶动态| a在线观看视频网站| 美女免费视频网站| 麻豆一二三区av精品| 亚洲国产精品成人综合色| 人人妻,人人澡人人爽秒播| 国产又黄又爽又无遮挡在线| 久久精品国产亚洲av高清一级| 亚洲全国av大片| 最新美女视频免费是黄的| 制服丝袜大香蕉在线| 国产爱豆传媒在线观看 | 哪里可以看免费的av片| 国内精品一区二区在线观看| tocl精华| av在线播放免费不卡| 亚洲成av人片在线播放无| 后天国语完整版免费观看| 久久这里只有精品中国| 国产伦一二天堂av在线观看| 欧美日韩一级在线毛片| 91麻豆精品激情在线观看国产| 国产伦一二天堂av在线观看| av福利片在线| 久久中文看片网| 欧美性长视频在线观看| 99久久国产精品久久久| 88av欧美| 国产精品1区2区在线观看.| 国产成人欧美在线观看|