• <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)全局
    在线a可以看的网站| 色尼玛亚洲综合影院| 亚洲国产精品成人综合色| 欧美性感艳星| 精品免费久久久久久久清纯| 欧美日韩乱码在线| 精品久久久久久久毛片微露脸| 偷拍熟女少妇极品色| 最近视频中文字幕2019在线8| 少妇的逼好多水| 制服丝袜大香蕉在线| а√天堂www在线а√下载| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧洲综合997久久,| 69人妻影院| 嫁个100分男人电影在线观看| 最近在线观看免费完整版| av女优亚洲男人天堂| 中出人妻视频一区二区| 欧美黄色淫秽网站| 在线观看午夜福利视频| 性色av乱码一区二区三区2| 国产欧美日韩一区二区三| 久久人人精品亚洲av| 久久精品国产亚洲av香蕉五月| 热99在线观看视频| 黄片小视频在线播放| 国产高清视频在线观看网站| 国产国拍精品亚洲av在线观看 | 亚洲18禁久久av| 国产成人福利小说| 在线免费观看不下载黄p国产 | 内射极品少妇av片p| 精品国产超薄肉色丝袜足j| 日韩欧美一区二区三区在线观看| 亚洲无线在线观看| 在线观看美女被高潮喷水网站 | 国产不卡一卡二| 久久婷婷人人爽人人干人人爱| 日韩欧美三级三区| 中文字幕熟女人妻在线| 淫秽高清视频在线观看| a级毛片a级免费在线| 亚洲精品国产精品久久久不卡| 成年女人永久免费观看视频| 亚洲国产欧洲综合997久久,| 噜噜噜噜噜久久久久久91| 久久天躁狠狠躁夜夜2o2o| 给我免费播放毛片高清在线观看| 日韩亚洲欧美综合| 首页视频小说图片口味搜索| 内地一区二区视频在线| 99热这里只有是精品50| 久久精品国产综合久久久| 又黄又粗又硬又大视频| 欧美日韩精品网址| 人妻夜夜爽99麻豆av| 国产精品影院久久| 国产伦精品一区二区三区视频9 | 天美传媒精品一区二区| 男女视频在线观看网站免费| 99久久精品国产亚洲精品| 日韩欧美国产在线观看| 看片在线看免费视频| 成人国产综合亚洲| 国产亚洲精品久久久久久毛片| 又粗又爽又猛毛片免费看| 免费看a级黄色片| 国产高清激情床上av| 国产亚洲欧美在线一区二区| 亚洲av第一区精品v没综合| 男插女下体视频免费在线播放| 国产视频内射| 国产精品久久电影中文字幕| 亚洲第一电影网av| 日日摸夜夜添夜夜添小说| 精品福利观看| 不卡一级毛片| 久久人人精品亚洲av| 亚洲av免费在线观看| 老汉色av国产亚洲站长工具| a级毛片a级免费在线| 好男人在线观看高清免费视频| 老司机午夜福利在线观看视频| 99热这里只有是精品50| 91九色精品人成在线观看| 动漫黄色视频在线观看| 亚洲美女黄片视频| 少妇丰满av| 琪琪午夜伦伦电影理论片6080| 在线观看免费午夜福利视频| 久久精品国产亚洲av涩爱 | or卡值多少钱| 国产麻豆成人av免费视频| 成人亚洲精品av一区二区| 国产伦精品一区二区三区四那| 看黄色毛片网站| 夜夜看夜夜爽夜夜摸| www.色视频.com| 日韩av在线大香蕉| aaaaa片日本免费| 黑人欧美特级aaaaaa片| 搡老熟女国产l中国老女人| 伊人久久大香线蕉亚洲五| 亚洲av日韩精品久久久久久密| 国产淫片久久久久久久久 | 黄片大片在线免费观看| 国产成人aa在线观看| 国产欧美日韩一区二区精品| 黄色丝袜av网址大全| 国产视频一区二区在线看| 长腿黑丝高跟| 欧美丝袜亚洲另类 | 亚洲av中文字字幕乱码综合| 一区二区三区高清视频在线| 国产91精品成人一区二区三区| 露出奶头的视频| 亚洲国产精品成人综合色| 日本黄色片子视频| 亚洲中文字幕日韩| 色综合站精品国产| 成人鲁丝片一二三区免费| 无限看片的www在线观看| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 午夜激情欧美在线| 高清在线国产一区| 黄片小视频在线播放| 亚洲av成人不卡在线观看播放网| 一a级毛片在线观看| 欧美成人一区二区免费高清观看| 看片在线看免费视频| 日韩欧美精品v在线| 母亲3免费完整高清在线观看| 最近最新中文字幕大全电影3| 两个人的视频大全免费| 国产精品久久久久久人妻精品电影| 中文在线观看免费www的网站| 亚洲成人久久性| 欧美乱色亚洲激情| 男女做爰动态图高潮gif福利片| 给我免费播放毛片高清在线观看| 一个人看视频在线观看www免费 | 51国产日韩欧美| 国产亚洲精品av在线| 中文字幕人成人乱码亚洲影| 精品欧美国产一区二区三| 亚洲国产欧美网| 国产不卡一卡二| 亚洲最大成人手机在线| 在线观看舔阴道视频| 午夜a级毛片| 黄片大片在线免费观看| 国产免费男女视频| 一二三四社区在线视频社区8| 很黄的视频免费| 欧美不卡视频在线免费观看| 久久久国产精品麻豆| 精品日产1卡2卡| 国产精品久久久久久久电影 | 婷婷亚洲欧美| 国产三级黄色录像| 日韩欧美在线二视频| 亚洲一区二区三区色噜噜| 露出奶头的视频| 法律面前人人平等表现在哪些方面| 亚洲精品在线美女| а√天堂www在线а√下载| xxx96com| 亚洲欧美日韩无卡精品| 久久久色成人| 欧美激情久久久久久爽电影| 19禁男女啪啪无遮挡网站| 国产精品亚洲av一区麻豆| 三级男女做爰猛烈吃奶摸视频| 床上黄色一级片| www国产在线视频色| 亚洲电影在线观看av| a级毛片a级免费在线| 欧洲精品卡2卡3卡4卡5卡区| 神马国产精品三级电影在线观看| 最后的刺客免费高清国语| 久久久久久久久中文| 最新中文字幕久久久久| 黄色片一级片一级黄色片| 少妇高潮的动态图| 亚洲av电影在线进入| 法律面前人人平等表现在哪些方面| 久久久久久人人人人人| 18禁裸乳无遮挡免费网站照片| 真人做人爱边吃奶动态| 97超视频在线观看视频| 欧美中文日本在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 听说在线观看完整版免费高清| 日本在线视频免费播放| 美女大奶头视频| 一本一本综合久久| 少妇高潮的动态图| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 欧美午夜高清在线| 日本黄色片子视频| 国产乱人视频| 男插女下体视频免费在线播放| 亚洲无线观看免费| 国产一区二区激情短视频| 老司机福利观看| 人妻夜夜爽99麻豆av| 狠狠狠狠99中文字幕| 淫妇啪啪啪对白视频| 午夜福利18| 亚洲最大成人手机在线| av天堂中文字幕网| 欧美zozozo另类| 午夜视频国产福利| 天堂动漫精品| 人妻久久中文字幕网| 一区二区三区高清视频在线| 嫁个100分男人电影在线观看| 成人永久免费在线观看视频| 九色成人免费人妻av| 91麻豆av在线| 人人妻人人澡欧美一区二区| 亚洲成人久久爱视频| 手机成人av网站| 国产一区二区激情短视频| 亚洲国产精品久久男人天堂| 神马国产精品三级电影在线观看| 99在线视频只有这里精品首页| 亚洲最大成人中文| 欧美一区二区亚洲| 天堂av国产一区二区熟女人妻| 黄色女人牲交| 国产亚洲欧美98| 日本黄色视频三级网站网址| 啪啪无遮挡十八禁网站| 99国产精品一区二区三区| 国产国拍精品亚洲av在线观看 | 日韩欧美精品v在线| 久久久色成人| 亚洲国产高清在线一区二区三| 51午夜福利影视在线观看| 99热精品在线国产| 深夜精品福利| 精品久久久久久成人av| 亚洲成人中文字幕在线播放| 国产日本99.免费观看| 熟女电影av网| 午夜福利在线观看吧| 韩国av一区二区三区四区| 日韩欧美三级三区| 日韩欧美国产在线观看| 国产伦精品一区二区三区视频9 | 村上凉子中文字幕在线| 免费在线观看成人毛片| 亚洲精品亚洲一区二区| 给我免费播放毛片高清在线观看| 在线观看免费视频日本深夜| 一个人观看的视频www高清免费观看| 色综合婷婷激情| 黄片小视频在线播放| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 欧美成狂野欧美在线观看| 激情在线观看视频在线高清| а√天堂www在线а√下载| 真人做人爱边吃奶动态| 少妇的逼好多水| 亚洲天堂国产精品一区在线| 九九久久精品国产亚洲av麻豆| 最新在线观看一区二区三区| 此物有八面人人有两片| 欧美日韩福利视频一区二区| 波野结衣二区三区在线 | 一区二区三区免费毛片| 老司机午夜十八禁免费视频| netflix在线观看网站| 此物有八面人人有两片| 在线观看美女被高潮喷水网站 | 日本免费a在线| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| 噜噜噜噜噜久久久久久91| 制服人妻中文乱码| 国产免费av片在线观看野外av| 蜜桃亚洲精品一区二区三区| 亚洲美女黄片视频| 最新在线观看一区二区三区| 国产三级在线视频| e午夜精品久久久久久久| 悠悠久久av| 国产精品久久视频播放| 白带黄色成豆腐渣| 岛国在线免费视频观看| 国产伦在线观看视频一区| 国产av不卡久久| av黄色大香蕉| 国产91精品成人一区二区三区| 成人亚洲精品av一区二区| 亚洲久久久久久中文字幕| 精品人妻1区二区| 日韩有码中文字幕| ponron亚洲| 国产精华一区二区三区| 精品国产亚洲在线| 亚洲成av人片免费观看| 精品一区二区三区视频在线观看免费| 国产成人影院久久av| 又黄又爽又免费观看的视频| 久久久久久久精品吃奶| 白带黄色成豆腐渣| 嫩草影院入口| 一个人看的www免费观看视频| 69人妻影院| 久久精品亚洲精品国产色婷小说| 免费观看人在逋| netflix在线观看网站| 久久久久久国产a免费观看| 国内精品久久久久精免费| 美女黄网站色视频| 国产成人av激情在线播放| 久久99热这里只有精品18| 国产精品久久电影中文字幕| 男女做爰动态图高潮gif福利片| 婷婷丁香在线五月| 久久精品国产自在天天线| 成年女人毛片免费观看观看9| 久久欧美精品欧美久久欧美| 午夜激情欧美在线| 女人高潮潮喷娇喘18禁视频| 欧美丝袜亚洲另类 | 女同久久另类99精品国产91| 日韩高清综合在线| 成人特级黄色片久久久久久久| 一个人观看的视频www高清免费观看| 在线十欧美十亚洲十日本专区| 我要搜黄色片| 好看av亚洲va欧美ⅴa在| 九九热线精品视视频播放| 一本精品99久久精品77| 亚洲av免费在线观看| 久久国产精品人妻蜜桃| 日本精品一区二区三区蜜桃| 啦啦啦韩国在线观看视频| 国产91精品成人一区二区三区| 黄色丝袜av网址大全| 99久久精品热视频| 最近最新中文字幕大全电影3| 免费av毛片视频| 亚洲国产中文字幕在线视频| 国产精品久久久久久精品电影| 日本在线视频免费播放| 色吧在线观看| 精品久久久久久久久久免费视频| 亚洲18禁久久av| 午夜精品一区二区三区免费看| 黄片小视频在线播放| 最近最新中文字幕大全电影3| 国产91精品成人一区二区三区| 一本精品99久久精品77| 日本精品一区二区三区蜜桃| 高潮久久久久久久久久久不卡| 亚洲熟妇中文字幕五十中出| 欧美最黄视频在线播放免费| 九色成人免费人妻av| av专区在线播放| 麻豆国产97在线/欧美| 免费在线观看日本一区| 色综合站精品国产| 好看av亚洲va欧美ⅴa在| 国内毛片毛片毛片毛片毛片| 内射极品少妇av片p| 69人妻影院| 又黄又爽又免费观看的视频| 18美女黄网站色大片免费观看| 在线a可以看的网站| 十八禁人妻一区二区| 国产乱人伦免费视频| 天天躁日日操中文字幕| 小说图片视频综合网站| 久久久色成人| 日韩精品青青久久久久久| 欧美色欧美亚洲另类二区| 色在线成人网| 欧美中文日本在线观看视频| 一级a爱片免费观看的视频| 亚洲第一欧美日韩一区二区三区| avwww免费| 亚洲男人的天堂狠狠| 一区福利在线观看| 国产综合懂色| 19禁男女啪啪无遮挡网站| 波多野结衣巨乳人妻| 乱人视频在线观看| 岛国在线观看网站| 一区二区三区免费毛片| 国产一区二区三区视频了| 亚洲国产色片| 日韩欧美国产一区二区入口| xxxwww97欧美| 啦啦啦免费观看视频1| 亚洲男人的天堂狠狠| 他把我摸到了高潮在线观看| 国产高清三级在线| 欧美一区二区国产精品久久精品| 欧美日本视频| 狠狠狠狠99中文字幕| 国产乱人伦免费视频| 性色avwww在线观看| 国产一区二区三区视频了| 在线观看免费午夜福利视频| 人人妻人人澡欧美一区二区| 亚洲人成网站在线播| 黄色女人牲交| 啦啦啦观看免费观看视频高清| 亚洲国产欧美人成| 国产免费男女视频| 国模一区二区三区四区视频| 国产一级毛片七仙女欲春2| 国产淫片久久久久久久久 | 热99在线观看视频| 亚洲在线观看片| 欧美一区二区国产精品久久精品| 久久精品国产亚洲av涩爱 | 国产aⅴ精品一区二区三区波| 精品乱码久久久久久99久播| 99精品在免费线老司机午夜| or卡值多少钱| 99久久九九国产精品国产免费| 亚洲精品影视一区二区三区av| 日本黄色视频三级网站网址| 51国产日韩欧美| 日韩欧美精品v在线| 亚洲久久久久久中文字幕| 国产私拍福利视频在线观看| 国产探花在线观看一区二区| 偷拍熟女少妇极品色| 激情在线观看视频在线高清| 欧美乱妇无乱码| 日韩人妻高清精品专区| 久久精品国产99精品国产亚洲性色| av国产免费在线观看| 精品久久久久久久久久久久久| 在线免费观看的www视频| 啦啦啦韩国在线观看视频| 乱人视频在线观看| 欧美日韩精品网址| 最近在线观看免费完整版| 婷婷精品国产亚洲av在线| 亚洲第一欧美日韩一区二区三区| 91麻豆精品激情在线观看国产| 夜夜躁狠狠躁天天躁| 亚洲av电影在线进入| 成人高潮视频无遮挡免费网站| 欧美日韩精品网址| 亚洲五月天丁香| 国产探花极品一区二区| 嫩草影院精品99| 黑人欧美特级aaaaaa片| 草草在线视频免费看| 色视频www国产| 人妻夜夜爽99麻豆av| 国产成人影院久久av| 免费看光身美女| 天天添夜夜摸| 99热精品在线国产| 中文字幕人妻丝袜一区二区| 亚洲第一电影网av| 色噜噜av男人的天堂激情| x7x7x7水蜜桃| 1000部很黄的大片| 国产精品嫩草影院av在线观看 | 欧美成人性av电影在线观看| 免费av观看视频| 亚洲无线在线观看| 91九色精品人成在线观看| 母亲3免费完整高清在线观看| 日本黄大片高清| 日本撒尿小便嘘嘘汇集6| 国产一区二区在线观看日韩 | 亚洲美女视频黄频| 精品久久久久久久久久免费视频| 极品教师在线免费播放| 麻豆成人av在线观看| 国产精品美女特级片免费视频播放器| www日本黄色视频网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99久久九九国产精品国产免费| 给我免费播放毛片高清在线观看| 欧美日本亚洲视频在线播放| 日韩欧美国产在线观看| 国产精品爽爽va在线观看网站| 国产成+人综合+亚洲专区| 久久精品91无色码中文字幕| ponron亚洲| 色综合欧美亚洲国产小说| 亚洲精品色激情综合| 人人妻人人看人人澡| 一个人观看的视频www高清免费观看| 日本黄大片高清| 欧美丝袜亚洲另类 | 日韩 欧美 亚洲 中文字幕| 国产又黄又爽又无遮挡在线| 亚洲性夜色夜夜综合| 三级男女做爰猛烈吃奶摸视频| 国产精品美女特级片免费视频播放器| 国产不卡一卡二| 中文字幕久久专区| av视频在线观看入口| 国产午夜福利久久久久久| 免费在线观看成人毛片| www.色视频.com| 小蜜桃在线观看免费完整版高清| 91久久精品国产一区二区成人 | a在线观看视频网站| 精品电影一区二区在线| 午夜福利成人在线免费观看| 久久精品影院6| 欧美精品啪啪一区二区三区| 麻豆国产av国片精品| 亚洲欧美日韩无卡精品| 淫妇啪啪啪对白视频| 天美传媒精品一区二区| 色哟哟哟哟哟哟| 99久国产av精品| 免费高清视频大片| 国产毛片a区久久久久| 亚洲欧美日韩高清专用| 亚洲欧美日韩无卡精品| 国产色婷婷99| 日本a在线网址| 日韩中文字幕欧美一区二区| e午夜精品久久久久久久| 色吧在线观看| 欧美+日韩+精品| 亚洲av免费高清在线观看| 日韩欧美精品v在线| 久久久精品大字幕| 男女之事视频高清在线观看| 1000部很黄的大片| 99热这里只有精品一区| 女生性感内裤真人,穿戴方法视频| 成人国产综合亚洲| 狂野欧美白嫩少妇大欣赏| 国产主播在线观看一区二区| a级一级毛片免费在线观看| 母亲3免费完整高清在线观看| 久久精品人妻少妇| 日本三级黄在线观看| 精品久久久久久久毛片微露脸| 欧美区成人在线视频| 美女高潮的动态| 国产精品一区二区三区四区久久| 亚洲国产精品合色在线| 一进一出抽搐动态| 亚洲美女视频黄频| 在线播放无遮挡| 日本一二三区视频观看| 18禁黄网站禁片免费观看直播| 亚洲在线观看片| 久久久久久久精品吃奶| 男人舔奶头视频| 一本精品99久久精品77| 一个人免费在线观看电影| 国产99白浆流出| 亚洲av电影不卡..在线观看| 成年女人永久免费观看视频| 麻豆成人av在线观看| 色视频www国产| 亚洲熟妇中文字幕五十中出| 亚洲av五月六月丁香网| 日韩高清综合在线| 一区二区三区高清视频在线| 色播亚洲综合网| 麻豆成人午夜福利视频| 久久精品91蜜桃| 女警被强在线播放| 亚洲乱码一区二区免费版| www国产在线视频色| 国产一区在线观看成人免费| 99热这里只有精品一区| 在线免费观看的www视频| 欧美大码av| 亚洲熟妇熟女久久| 久久久久久久久中文| 叶爱在线成人免费视频播放| 亚洲专区国产一区二区| 日本 欧美在线| 99久久久亚洲精品蜜臀av| 97超级碰碰碰精品色视频在线观看| 男女午夜视频在线观看| 色av中文字幕| 草草在线视频免费看| 亚洲欧美日韩高清在线视频| 色综合站精品国产| 一个人看视频在线观看www免费 | 精品国内亚洲2022精品成人| 中文字幕熟女人妻在线| 国产精品一区二区三区四区免费观看 | 成年女人永久免费观看视频| 国产午夜精品论理片| 久久6这里有精品| 夜夜看夜夜爽夜夜摸| 少妇人妻精品综合一区二区 | 国产精品电影一区二区三区| 免费大片18禁| h日本视频在线播放| 国产中年淑女户外野战色| 国产精品一区二区三区四区久久| 网址你懂的国产日韩在线| 成人鲁丝片一二三区免费| 国产三级黄色录像| 亚洲 欧美 日韩 在线 免费| 色av中文字幕|