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

    小波結(jié)合雙邊濾波抑制全息相干噪聲

    2016-10-14 02:13:07王慶伍
    光電工程 2016年8期
    關(guān)鍵詞:方法

    肖 文,王慶伍,潘 鋒

    ?

    小波結(jié)合雙邊濾波抑制全息相干噪聲

    肖 文,王慶伍,潘 鋒

    ( 北京航空航天大學(xué) 儀器科學(xué)與光電工程學(xué)院,北京 100191 )

    抑制數(shù)字全息中相干噪聲可以獲得更高質(zhì)量的數(shù)字全息再現(xiàn)圖像。本文提出一種結(jié)合雙邊濾波和小波鄰域閾值的數(shù)字全息相干噪聲抑制方法。對再現(xiàn)像進行小波變換,利用小波鄰域方法對高頻分量小波閾值去噪,對低頻分量做雙邊濾波處理,最后做小波逆變換得到去噪再現(xiàn)像。通過仿真分析將該方法與中值、均值、維納濾波等方法進行了比較,去噪后所得圖像的峰值信噪比提高了0.4 dB~2 dB,與原始圖像的相關(guān)性更高。利用全息實驗數(shù)據(jù)驗證了方法,結(jié)果表明本文方法能夠有效抑制散斑噪聲,同時更好的保留再現(xiàn)像的細節(jié)信息。

    數(shù)字全息;相干噪聲;散斑抑制;小波變換;雙邊濾波

    0 引 言

    數(shù)字全息技術(shù)[1]是基于傳統(tǒng)全息原理,借助光電圖像探測和計算機圖像處理技術(shù),利用CCD或CMOS傳感器以數(shù)字形式記錄全息圖,然后通過計算機模擬光學(xué)衍射過程來實現(xiàn)全息再現(xiàn),能夠定量獲得物體的振幅和相位像。數(shù)字全息作為一種相干成像技術(shù),通常采用激光器作為光源。由于激光具有較好的時間和空間相干性,當(dāng)物體具有粗糙表面或折射率不均勻等特性,以及光路中存在多層反射面、劃痕、灰塵顆粒等因素,將產(chǎn)生寄生條紋、散斑分布等相干噪聲[2],并與物光場同時記錄在全息圖中。因此數(shù)字再現(xiàn)時,相干噪聲也被重構(gòu)再現(xiàn),這將嚴重影響再現(xiàn)像質(zhì)量,使得再現(xiàn)圖像分辨率降低,圖像細節(jié)結(jié)構(gòu)模糊,降低測量精度。因此,在數(shù)字全息中抑制相干噪聲是非常必要的。

    目前全息相干噪聲抑制方法主要有兩類。一類是在全息圖獲取過程中運用光學(xué)方法[3-7],例如采用弱相干或非相干光源,但是該方法對光路的結(jié)構(gòu)及其穩(wěn)定性有嚴格要求。另外通過多波長、多偏振態(tài)、多角度照明等方式記錄多幅全息圖。這些方法的物理基礎(chǔ)是各全息圖記錄了同一物的光場,但具有不同相干噪聲分布。進而通過平均疊加處理實現(xiàn)相干噪聲抑制。這些方法的局限性在于增加了光路結(jié)構(gòu)和記錄過程的復(fù)雜性。另一類方法是運用圖像處理技術(shù)對再現(xiàn)像進行去噪處理。例如Garcia-Sucerquia[8]等人采用中值濾波和均值濾波減小激光相干散斑顆粒。Maycock[9]等人采用傅里葉濾波法濾除相干散斑噪聲。Sharma[10]等人采用小波變換濾除相干散斑噪聲。蔡曉鷗[11]等用帶孔徑函數(shù)的維納濾波法濾除散斑噪聲。賈勤[12]等人利用邊緣檢測結(jié)合小波變換的方法有效抑制了全息再現(xiàn)像的激光散斑。Uzan等人利用NLM濾波算法[13-14]對數(shù)字全息散相干噪聲進行抑制。Montresor等人將均值、維納、小波等濾波方法對全息圖降噪效果進行了對比。與第一類方法相比,這類方法的優(yōu)勢在于不需要任何硬件支持,僅通過對全息再現(xiàn)像的后處理來進行噪聲抑制,其代價最小。

    然而基于圖像處理技術(shù)來抑制相干噪聲方法面臨的共同問題是抑制噪聲的同時會損失部分再現(xiàn)像的細節(jié)信息,即影響圖像的分辨率。針對這一矛盾,本文提出一種結(jié)合雙邊濾波的小波鄰域閾值相干噪聲抑制方法。針對全息再現(xiàn)像的小波變換高頻分量利用鄰域閾值進行閾值萎縮,而對其低頻分量利用雙邊濾波進行濾波,這樣將兩種濾波方法的優(yōu)勢都利用起來,同時彌補了其各自不足,實現(xiàn)全息再現(xiàn)像去噪的同時更好的保留圖像的細節(jié)信息。

    1 結(jié)合雙邊濾波的小波鄰域閾值抑制全息相干噪聲算法

    基于Goodman完全發(fā)育的相干散斑噪聲特性,數(shù)字全息相干散斑噪聲可近似為一種均勻分布乘性噪聲模型。其模型為

    式中:表示受噪聲影響的再現(xiàn)像,表示物體再現(xiàn)像,表示散斑噪聲,、分別代表再現(xiàn)像橫縱像素。對式(1)作對數(shù)運算,將再現(xiàn)像乘性噪聲轉(zhuǎn)換為加性噪聲進行處理,得到:

    此時散斑噪聲轉(zhuǎn)化為加性噪聲ln,是近似服從零均值高斯分布的。

    基于小波變換具有多分辨率、去相關(guān)性、選基靈活性等特點,本文將采用小波閾值濾波結(jié)合雙邊濾波來進行相干噪聲抑制。首先,對全息再現(xiàn)像進行小波變換,得到一個低頻分量及若干尺度的三方向高頻分量。然后,對高頻分量利用鄰域閾值進行閾值萎縮,保留高頻細節(jié)的同時抑制噪聲。一般小波濾波方法中低頻分量是不進行濾波處理的。由于相干散斑場與物光場疊加是隨機的,使得散斑頻譜與物光頻譜產(chǎn)生混淆,導(dǎo)致再現(xiàn)像進行小波變換后,其低頻分量中也含有部分顆粒狀相干噪聲?;谶@一特點,在利用小波變換抑制全息相干噪聲過程中,有必要對低頻分量進行去噪處理。雙邊濾波[15]在保留圖像的基本特性,如結(jié)構(gòu)、紋理、平均強度等信息是有優(yōu)勢的,因此選用雙邊濾波來對再現(xiàn)像的低頻分量進行處理。通過綜合利用小波閾值方法以及雙邊濾波方法,在一定程度上緩解了圖像處理抑制相干噪聲的共同矛盾,抑制散斑的同時較好地保留了數(shù)字全息再現(xiàn)像的細節(jié)信息,從而得到更高質(zhì)量的數(shù)字全息再現(xiàn)像。

    本文綜合濾波算法的流程如圖1所示。

    具體步驟如下:

    1) 再現(xiàn)像對數(shù)化 通過對數(shù)化操作,將再現(xiàn)像中的乘性散斑噪聲轉(zhuǎn)化為近似高斯分布的加性噪聲。

    圖1 濾波流程圖

    2) 再現(xiàn)像小波分解為了獲得原圖像大小一致的高頻分量,采用平穩(wěn)小波變換,對對數(shù)化再現(xiàn)像進行小波分解。得到3+1個子帶(為分解尺度),其中低頻子帶LL一個,高頻子帶共3個,分別為水平方向高頻分量HL,垂直方向高頻分量VH,對角方向高頻分量HH(=1,2,…,)。

    3) 低頻分量LL雙邊濾波去噪運用平移不變高斯濾波器對小波分解的低頻分量進行處理,在實現(xiàn)濾除噪聲的同時保護了低頻分量中的圖像特征。利用低頻分量小波系數(shù)區(qū)域位置相似性和絕對值相似性構(gòu)造了雙邊濾波器中區(qū)域濾波器和值域濾波器。

    其中:為待濾波低頻分量像素,為鄰域內(nèi)像素,鄰域大小為(2+1)×(2+1);及分別為低頻分量空間臨近度函數(shù)及小波系數(shù)值相似函數(shù),均取為參數(shù)間歐幾里德距離的高斯函數(shù);參數(shù)為鄰域空間方差,決定了鄰域權(quán)重因子衰減系數(shù);為灰度方差,控制著值域權(quán)重因子衰減系數(shù)。

    4) 高頻分量小波鄰域閾值去噪小波去噪中閾值的選擇對去噪結(jié)果影響明顯。閾值過大則容易造成過度平滑,閾值太小則噪聲去除效果不佳。考慮到再現(xiàn)像的鄰域小波系數(shù)相互之間具有強相關(guān)性,而噪聲是獨立分布的,其鄰域內(nèi)小波系數(shù)幾乎沒有相關(guān)性這一特點。本文選用小波NeighShrink閾值[16]對再現(xiàn)像小波分解的高頻分量進行閾值處理,實現(xiàn)噪聲與圖像高頻信息的分離。

    NeighShrink閾值是根據(jù)鄰域內(nèi)小波系數(shù)相關(guān)性來確定閾值,與系數(shù)越近,相關(guān)性越強。綜合考慮相關(guān)性以及運算效率,鄰域窗口選擇為3×3或5×5。以鄰域窗口內(nèi)小波系數(shù)的平方總和與通用閾值進行比較,為圖像總像素數(shù)。當(dāng)平方總和小于,則認為窗口中心的小波系數(shù)由噪聲產(chǎn)生,直接置零;反之,則對窗口中心的小波系數(shù)作收縮運算。小波系數(shù)的平方和和收縮因子及收縮運算的計算表達式分別為

    5) 小波逆變換利用處理后的低頻分量及高頻分量作小波逆變換,得到抑制噪聲后的再現(xiàn)像。

    2 仿真分析

    2.1 仿真數(shù)據(jù)

    本文方法將與文獻[8]中均值濾波、中值濾波算法、文獻[11]中所用維納濾波算法、文獻[16]中提出的雙邊濾波、文獻[12]中所用結(jié)合邊緣檢測的小波濾波等算法進行比較。圖像選擇cameraman、lena、mandrill、barbara四幅圖片,通過將均值為0,方差為0.1的均勻分布的隨機Speckle乘性噪聲模擬相干散斑噪聲,再對含噪圖像進行對數(shù)化處理,將乘性噪聲轉(zhuǎn)化為加性噪聲,此時仿真圖像中噪聲模型實際全息圖可用式(2)來表示。為保持統(tǒng)一,用于比較的均值、中值、維納小波鄰域濾波算法中鄰域窗口一律設(shè)定為5×5。雙邊濾波鄰近閾值及系數(shù)閾值分別設(shè)定為2及0.2。采用定量評價參數(shù)均方差MSE(用表示)、峰值信噪比PSNR(用表示)和相似度來衡量去噪結(jié)果,其表達式分別為:

    2.2 仿真結(jié)果及分析

    表1為加噪圖像以及經(jīng)六種不同方法去噪后的圖像與原始圖像進行比較,計算出的各自峰值信噪比以及結(jié)構(gòu)相似度統(tǒng)計結(jié)果。圖2為紋理、細節(jié)信息最豐富的mandrill圖像胡須附近區(qū)域去噪結(jié)果對比圖,圖3是lena圖像去噪結(jié)果對比圖。圖2及圖3中(a)是原始圖像,(b)是加噪后圖像,(c)是加噪圖像經(jīng)均值濾波結(jié)果,(d)是加噪圖像經(jīng)中值濾波結(jié)果,(e)是加噪圖像經(jīng)維納濾波結(jié)果,(f)是加噪圖像經(jīng)雙邊濾波結(jié)果,(g)是加噪圖像經(jīng)小波結(jié)合邊緣檢測方法濾波結(jié)果,(h)為加噪圖像經(jīng)本文方法濾波所得結(jié)果。

    表1 濾波結(jié)果對比表Table 1 Comparison table of filtering results

    圖2 mandrill區(qū)域截圖

    圖3 Lena濾波組圖

    通過表1,可以看出運用本文方法去噪后圖像的PSNR值比其它方法都要高0.5 dB至4 dB,同時與原始圖像相比其結(jié)構(gòu)相似度也最高,細節(jié)信息保持最好。從lena去噪組圖也可以看出,本文方法相干噪聲去除效果好,同時從lena的帽子、頭發(fā)和眼睫毛等也可看出,圖像細節(jié)信息保持較好。而在mandrill細節(jié)信息最為豐富的胡須區(qū)域截圖也可以看出,本文方法處理后狒狒(mandrill)的胡須邊緣輪廓清晰、細節(jié)結(jié)構(gòu)保持較好。由此可見本文方法與相比較的其它傳統(tǒng)去噪方法不僅去噪效果好,在保持細節(jié)信息上具有優(yōu)勢。

    3 實驗驗證

    3.1 實驗系統(tǒng)

    實驗對象選定色子、硬幣和分辨率板,三者材質(zhì)分別為塑料、金屬、玻璃,表面粗糙程度差異較大。實驗光路圖如圖4。激光器發(fā)出的激光束進入偏振分光棱鏡(PBS)分為兩束。兩束光分別經(jīng)空間濾波器擴束并經(jīng)透鏡L1和L2準直,其中一束光經(jīng)被測物反射后攜帶物面信息,即作為物光與另一束參考光沿入射方向被合光棱鏡(BS)接收。兩光束以小角度在CCD相機光敏面上相干疊加,產(chǎn)生離軸全息圖,進而傳送到計算機進行再現(xiàn)計算。其中光源為Nd:YAG固體激光器,波長為532 nm,輸出功率為15 mW,相干長度約為50 m。CCD相機像素數(shù)1 024 pixel×1 024 pixel,像素尺寸6.7 μm×6.7 μm。

    圖4 實驗光路結(jié)構(gòu)示意圖

    3.2 實驗結(jié)果與分析

    圖5及圖6分別是硬幣以及色子的再現(xiàn)圖像濾波結(jié)果對比圖。圖5及圖6中(a)是原始再現(xiàn)像,(b)是經(jīng)均值濾波結(jié)果,(c)是經(jīng)中值濾波結(jié)果,(d)是經(jīng)維納濾波結(jié)果,(e)是經(jīng)雙邊濾波結(jié)果,(f)是經(jīng)小波濾波結(jié)果,(g)為經(jīng)本文方法濾波所得結(jié)果。

    由圖5和6可以看出,相比于其它方法,本文方法去噪后硬幣和色子表面顆粒狀相干散斑噪聲消除得最好。同時,為了定量比較,用相干散斑指數(shù)來評價濾波算法對散斑的抑制能力。散斑指數(shù)(Speckle index,用表示)[17]定義為

    圖5 硬幣濾波對比圖

    圖6 色子濾波對比圖

    表2 散斑指數(shù)統(tǒng)計表Table 2 Speckle index statistics

    進一步利用分辨率板再現(xiàn)像進行比較。圖7中(a)、(b)、(c)、(d)分別是分辨率板再現(xiàn)像原始圖像,利用雙邊濾波去噪后的再現(xiàn)像,利用小波鄰域閾值去噪后的再現(xiàn)像,結(jié)合小波鄰域閾值及雙邊濾波去噪后的再現(xiàn)像。圖7(a)中白線位置屬于分辨率板中平坦區(qū)域,在不受相干噪聲影響的理想狀態(tài)下,其再現(xiàn)像上灰度值理應(yīng)趨于一個定值。圖8是針對圖7(a)中白線位置,對圖7(b)、7(c)、7(d)同樣位置直線上像元灰度值分析圖。圖8橫坐標為像素位置,縱坐標為像素歸一化灰度值,實線為小波鄰域去噪結(jié)果、短虛線為雙邊濾波去噪結(jié)果、雙點間斷虛線為本文所用方法去噪結(jié)果,單點間斷虛線為白線上灰度值平均值??煽闯霰疚姆椒ㄈピ牒?,藍色曲線像元灰度值變化減緩,平坦區(qū)域基本處于同一水平,更逼近于灰度平均值,去噪效果最好。

    圖7 分辨率板濾波對比圖

    圖8 灰度變化對比圖

    4 結(jié) 論

    本文提出了一種結(jié)合小波閾值去噪及雙邊濾波來抑制數(shù)字全息相干散斑噪聲的方法。針對此方法,首先進行仿真分析,利用均方差MSE、峰值信噪比PSNR和相似度三項指標對仿真結(jié)果進行了定量分析,結(jié)果表明本文方法均方差更小,峰值信噪比更大,去噪后的圖像在邊緣、紋理、結(jié)構(gòu)與原始圖像更接近。然后,開展數(shù)字全息實驗,對本方法的有效性和優(yōu)越性進行了實驗驗證。基于表面較粗糙的色子、硬幣和表面光滑的分辨率板的再現(xiàn)像,將本方法與其它方法的去噪結(jié)果進行定性和定量的比較分析。結(jié)果表明,相比于其它方法,本文結(jié)合雙邊濾波的鄰域小波閾值去噪方法,不僅抑制相干散斑噪聲效果好,同時去噪后圖像與原圖像相關(guān)度最高,對再現(xiàn)像的細節(jié)信息保留更完整。

    [1] Goodman J W,Awrence R W. Digital image fortnulation form electronically detected holograms [J]. Applied Physics Letters (S0003-6951),1967,11(2):77-79.

    [2] 張亦卓. 生物樣品的數(shù)字全息顯微相襯成像技術(shù)研究 [D]. 北京:北京工業(yè)大學(xué),2012:40-41.

    ZHANG Yizhuo. Digital holographic microscopy study of biological samples phase contrast imaging technique [D]. Beijing:Beijing University of Technology,2012:40-41.

    [3] XIAO Wen,ZHANG Jing,RONG Lu,. Improvement of speckle noise suppression in digital holography by rotating linear polarization state [J]. Chinese Optics Letters(S1671-7694),2011,9(6):0609011-0609013.

    [4] PAN Feng,XIAO Wen,LIU Shuo,. Coherent noise reduction in digital holographic phase contrast microscopy by slitly shifting object [J]. Optics Express(S1094-4087),2011,19(5):3862-3869.

    [5] 王大勇,王云新,郭莎,等. 基于多角度無透鏡傅里葉變換數(shù)字全息的散斑噪聲抑制成像研究 [J]. 物理學(xué)報,2014,63(15):175-180.

    WANG Dayong,WANG Yunxin,GUO Sha,. Research on speckle denoising by lensless Fourier transform holographic imaging with angular diversity [J]. Acta Physica Sinica,2014,63(15):175-180.

    [6] LI Dayan,Damien P Kelly,John T Sheridan. Speckle suppression by doubly scattering systems [J]. Applied Optics (S1539-4522),2013,52(35):8617-8626.

    [7] 張雯,周皓,顧濟華,等. 多光束數(shù)字全息的研究 [J]. 光子學(xué)報,2010,39(3):533-536.

    ZHANG Wen,ZHOU Hao,GU Jihua,. Multi-beamdigital holography [J]. Acta Photonics Sinica,2010,39(3):533-536.

    [8] Jorge Garcia-Sucerquia,Jorge Alexis Herrera Ramírez,Daniel Velásquez Prieto. Reduction of speckle noise in digital holography by using digital image processing [J]. Optik-International Journal for Light and Electron Optics(S0030-4026),2005,116(1):44-48.

    [9] Jonathan Maycock,Bryan M Hennelly,John B Mcdonald,. Reduction of speckle in digital holography by discrete Fourier filtering [J]. Journal of the Optical Society of America A(S1084-7529),2007,24(6):1617-1622.

    [10] Akshay Sharma,Gyanendra Sheoran,Jaffery Z A,. Improvement of signal-to-noise ratio in digital holography using wavelet transform [J]. Optics & Lasers in Engineering(S0143-8166),2008,46(1):42-47.

    [11] CAI Xiao¢ou. Reduction of speckle noise in the reconstructed image of digital holography [J]. Optik-International Journal for Light and Electron Optics(S0030-4026),2010,121(4):394-399.

    [12] 賈勤,李志全. 數(shù)字全息再現(xiàn)圖像散斑噪聲消除新方法 [J]. 計算機工程與應(yīng)用,2012,48(9):205-207.

    JIA Qin,LI Zhiquan. Speckle denoising new method for reconstructed image of digital holography [J]. Computer Engineering and Applications,2012,48(9):205-207.

    [13] Uzan A,Rivenson Y,Stern A. Speckle denoising in digital holography by nonlocal means filtering [J]. Applied Optics (S1539-4522),2013,52(1):195-200.

    [14] Jaehong Aum,Ji-hyun Kim,Jeong Jichai. Effective speckle noise suppression in optical coherence tomography images using nonlocal means denoising filter with double Gaussian anisotropic kernels [J]. Applied Optics(S1539-4522),2015,54(13):D43-D50.

    [15] Tomasi C,Manduchi R. Bilateral filtering for gray and color images [C]// Proceedings of the Sixth International Conference on Computer Vision,New Delhi,India,Jan 4-7,1998:839-846.

    [16] Dongwook C,Tien D B,CHEN Guangyi. Image Denoising based on Wavelet Shrinkage Using Neighbor and Level Dependency [J]. International Journal of Wavelets Multiresolution and Information Processing(S0219-6913),2009,7(3):299-311.

    [17] Dewaele P,Wambacq P,Oosterlinck A,. Comparison of some speckle reduction techniques for SAR images [Z]. Geoscience and Remote Sensing Symposium. Maryland,USA:the University of Maryland college Park,1990:2417-2422.

    Suppression of Coherent Noise by Wavelet Combined with Bilateral Filtering in Digital Holography

    XIAO Wen,WANG Qingwu,PAN Feng

    ( School of Instrumentation Science and Opto-electronics Engineering, Beihang University, Beijing 100191, China )

    Suppression of coherent noise in digital holography can obtain higher quality digital holographic reconstructed image. A suppression of holographic coherent noise method which combines bilateral filter and wavelet filter is proposed. Firstly, the holographic reconstructed image is decomposed by wavelet transformation. Then the high-frequency component in wavelet domain is suppressed by the wavelet shrinkage algorithm. The low-pass approximation component is filtered by the bilateral filter. Finally, the inverse wavelet transform is processed. Through the simulation analysis, this method improves the image peak signal-to-noise ratio 0.4 dB~2 dB, and got a higher correlation coefficientcloser with the original image. The method is verified by holographic experiment data,and the results show that this method can effectively suppress the speckle noise, and at the same time, retain more details information of the reconstructed image.

    digital holographic; coherent noise; speckle suppression; wavelet transform; bilateral filtering

    1003-501X(2016)08-0039-08

    TN26

    A

    10.3969/j.issn.1003-501X.2016.08.007

    2015-12-04;

    2016-01-27

    國家自然科學(xué)基金資助項目(61177006)

    肖文(1962-),男(漢族),陜西西安人。教授,博士,主要研究方向為數(shù)字全息三維成像及相干檢測。E-mail: xiaow@buaa.edu.cn。

    王慶伍(1985-),男(漢族),安徽蕪湖人。碩士研究生,主要研究方向為數(shù)字全息成像優(yōu)化。E-mail: wangqingwu8510@163.com。

    猜你喜歡
    方法
    中醫(yī)特有的急救方法
    中老年保健(2021年9期)2021-08-24 03:52:04
    高中數(shù)學(xué)教學(xué)改革的方法
    河北畫報(2021年2期)2021-05-25 02:07:46
    化學(xué)反應(yīng)多變幻 “虛擬”方法幫大忙
    變快的方法
    兒童繪本(2020年5期)2020-04-07 17:46:30
    學(xué)習(xí)方法
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    最有效的簡單方法
    山東青年(2016年1期)2016-02-28 14:25:23
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    韩国精品一区二区三区 | 美女福利国产在线| 777米奇影视久久| 欧美精品一区二区免费开放| 亚洲精品aⅴ在线观看| 国产有黄有色有爽视频| 十八禁网站网址无遮挡| 久久国产精品男人的天堂亚洲 | 亚洲av电影在线进入| 精品一品国产午夜福利视频| 久久国产精品男人的天堂亚洲 | 这个男人来自地球电影免费观看 | 国产成人a∨麻豆精品| 亚洲欧洲日产国产| 十分钟在线观看高清视频www| 少妇被粗大猛烈的视频| 亚洲色图 男人天堂 中文字幕 | 亚洲国产精品专区欧美| 欧美日韩视频精品一区| 日韩av在线免费看完整版不卡| 中文字幕人妻丝袜制服| 黄色配什么色好看| 欧美人与善性xxx| 宅男免费午夜| 午夜视频国产福利| 日韩视频在线欧美| av在线老鸭窝| 久久精品国产亚洲av天美| 国产深夜福利视频在线观看| 国产精品国产三级国产专区5o| 秋霞在线观看毛片| 在线 av 中文字幕| 国产精品熟女久久久久浪| 精品少妇内射三级| 嫩草影院入口| 国语对白做爰xxxⅹ性视频网站| 国产综合精华液| 色视频在线一区二区三区| 国产成人av激情在线播放| 精品第一国产精品| 亚洲av电影在线观看一区二区三区| 欧美少妇被猛烈插入视频| 久久ye,这里只有精品| 国产精品一区二区在线不卡| 丰满迷人的少妇在线观看| 韩国高清视频一区二区三区| 丝袜美足系列| videossex国产| 欧美人与性动交α欧美软件 | 亚洲欧美精品自产自拍| 人体艺术视频欧美日本| 大陆偷拍与自拍| 在线天堂最新版资源| 一级毛片黄色毛片免费观看视频| 亚洲第一区二区三区不卡| 蜜臀久久99精品久久宅男| 满18在线观看网站| 成人毛片60女人毛片免费| 亚洲激情五月婷婷啪啪| 精品人妻偷拍中文字幕| 男的添女的下面高潮视频| a级毛片在线看网站| 国产永久视频网站| 亚洲av中文av极速乱| 观看av在线不卡| 亚洲美女视频黄频| www.色视频.com| 免费不卡的大黄色大毛片视频在线观看| 少妇人妻久久综合中文| 亚洲国产av新网站| 成人国产av品久久久| 美女主播在线视频| 国产熟女午夜一区二区三区| 日韩制服骚丝袜av| 日韩不卡一区二区三区视频在线| 久久久久久伊人网av| 国产精品一区二区在线观看99| 在线观看www视频免费| 丝袜在线中文字幕| 亚洲av在线观看美女高潮| 少妇精品久久久久久久| 成人午夜精彩视频在线观看| 少妇被粗大的猛进出69影院 | 精品人妻偷拍中文字幕| 欧美日韩成人在线一区二区| 两性夫妻黄色片 | 午夜视频国产福利| 久热这里只有精品99| 少妇熟女欧美另类| 一二三四在线观看免费中文在 | 国产精品一二三区在线看| av不卡在线播放| 国产精品三级大全| 久久人人爽人人爽人人片va| 日本猛色少妇xxxxx猛交久久| 免费av中文字幕在线| 亚洲精品国产色婷婷电影| 在线观看美女被高潮喷水网站| 亚洲欧美日韩卡通动漫| 亚洲国产最新在线播放| 免费日韩欧美在线观看| 下体分泌物呈黄色| 国产精品国产三级国产专区5o| 亚洲精品456在线播放app| 女的被弄到高潮叫床怎么办| 欧美日韩亚洲高清精品| 涩涩av久久男人的天堂| 国产xxxxx性猛交| 精品视频人人做人人爽| 国产不卡av网站在线观看| 亚洲经典国产精华液单| 99国产精品免费福利视频| 免费看光身美女| 99国产综合亚洲精品| 亚洲欧美精品自产自拍| 亚洲精品久久久久久婷婷小说| 亚洲天堂av无毛| 精品一区二区三区四区五区乱码 | 成人毛片60女人毛片免费| 男女啪啪激烈高潮av片| 欧美日韩视频精品一区| 色吧在线观看| 成人18禁高潮啪啪吃奶动态图| 观看美女的网站| 国语对白做爰xxxⅹ性视频网站| 国产精品熟女久久久久浪| 人妻系列 视频| 飞空精品影院首页| 日本91视频免费播放| 美女中出高潮动态图| 欧美激情 高清一区二区三区| 九草在线视频观看| 免费观看av网站的网址| 精品人妻偷拍中文字幕| 一区在线观看完整版| 日韩不卡一区二区三区视频在线| 欧美精品国产亚洲| 王馨瑶露胸无遮挡在线观看| 大香蕉97超碰在线| 亚洲精品第二区| 成年av动漫网址| 免费人成在线观看视频色| 99国产精品免费福利视频| 香蕉丝袜av| 黄片无遮挡物在线观看| 一本色道久久久久久精品综合| 日韩大片免费观看网站| av一本久久久久| av在线播放精品| 又黄又粗又硬又大视频| 免费黄色在线免费观看| 天天躁夜夜躁狠狠躁躁| 精品国产一区二区三区四区第35| 国产精品国产av在线观看| 国产白丝娇喘喷水9色精品| 满18在线观看网站| 69精品国产乱码久久久| 中国三级夫妇交换| av天堂久久9| 我要看黄色一级片免费的| 妹子高潮喷水视频| 久久久久久久亚洲中文字幕| 中国美白少妇内射xxxbb| 人体艺术视频欧美日本| 日韩一本色道免费dvd| 精品久久久久久电影网| 中国国产av一级| 国产精品久久久久成人av| 亚洲国产欧美日韩在线播放| 久久国产精品男人的天堂亚洲 | 激情五月婷婷亚洲| 美女福利国产在线| 欧美成人午夜精品| 亚洲欧美色中文字幕在线| 91精品三级在线观看| 成人毛片60女人毛片免费| 亚洲欧洲国产日韩| 欧美精品av麻豆av| 国产在视频线精品| 久久99热6这里只有精品| 亚洲精品一区蜜桃| 秋霞伦理黄片| 伦精品一区二区三区| 十八禁网站网址无遮挡| 亚洲av中文av极速乱| 亚洲精品乱久久久久久| 久久99热这里只频精品6学生| 丝袜脚勾引网站| 久久 成人 亚洲| 亚洲精品第二区| 一个人免费看片子| 色5月婷婷丁香| 永久网站在线| 人人妻人人澡人人看| 国产视频首页在线观看| 国内精品宾馆在线| 97在线人人人人妻| 欧美日本中文国产一区发布| 大话2 男鬼变身卡| 亚洲av男天堂| 高清不卡的av网站| 麻豆精品久久久久久蜜桃| 天天影视国产精品| 国产白丝娇喘喷水9色精品| 国产一区二区三区av在线| 男女高潮啪啪啪动态图| 又黄又粗又硬又大视频| 欧美日韩亚洲高清精品| 免费看不卡的av| 纵有疾风起免费观看全集完整版| 亚洲精品日韩在线中文字幕| 纵有疾风起免费观看全集完整版| 亚洲av在线观看美女高潮| 在线观看三级黄色| 久久久久久久久久久久大奶| 久久久国产欧美日韩av| 老司机影院成人| 欧美日韩av久久| 黑人巨大精品欧美一区二区蜜桃 | 午夜久久久在线观看| 99久久综合免费| 国产69精品久久久久777片| 欧美变态另类bdsm刘玥| 欧美xxⅹ黑人| 久久这里只有精品19| 一二三四中文在线观看免费高清| 久久99热6这里只有精品| 赤兔流量卡办理| 国产精品久久久久久久电影| 男人操女人黄网站| 国产高清国产精品国产三级| 国产精品久久久久成人av| 男人舔女人的私密视频| 熟女av电影| 日韩av不卡免费在线播放| 亚洲国产欧美在线一区| videossex国产| 在线亚洲精品国产二区图片欧美| 97超碰精品成人国产| 国产免费福利视频在线观看| 日本vs欧美在线观看视频| 热re99久久国产66热| 婷婷色av中文字幕| 国产不卡av网站在线观看| 亚洲一码二码三码区别大吗| 日韩制服骚丝袜av| 国产亚洲最大av| a 毛片基地| 夜夜爽夜夜爽视频| 精品熟女少妇av免费看| 日韩中文字幕视频在线看片| 嫩草影院入口| 香蕉精品网在线| 宅男免费午夜| 亚洲精品视频女| 在线观看三级黄色| av电影中文网址| 伊人久久国产一区二区| 国产免费现黄频在线看| av免费观看日本| 人人妻人人添人人爽欧美一区卜| 免费黄色在线免费观看| 成人综合一区亚洲| 国产精品久久久久久久电影| 美女国产视频在线观看| 91精品国产国语对白视频| 婷婷成人精品国产| 午夜福利,免费看| 日产精品乱码卡一卡2卡三| 亚洲内射少妇av| 99久国产av精品国产电影| 久久国产精品男人的天堂亚洲 | 亚洲,一卡二卡三卡| 建设人人有责人人尽责人人享有的| 97超碰精品成人国产| 亚洲精品乱码久久久久久按摩| 2018国产大陆天天弄谢| 欧美 亚洲 国产 日韩一| 亚洲伊人色综图| 亚洲经典国产精华液单| 亚洲国产精品专区欧美| 久久精品国产综合久久久 | 男人操女人黄网站| 黄色怎么调成土黄色| 国产精品久久久久久久久免| 国产探花极品一区二区| 99香蕉大伊视频| 国产免费视频播放在线视频| 久久久久久人人人人人| www.av在线官网国产| 26uuu在线亚洲综合色| 最近中文字幕高清免费大全6| 男人舔女人的私密视频| 性高湖久久久久久久久免费观看| 日本午夜av视频| 黄色毛片三级朝国网站| 18禁动态无遮挡网站| 女人久久www免费人成看片| 亚洲国产av新网站| 国产1区2区3区精品| 天天躁夜夜躁狠狠躁躁| videosex国产| 国产成人91sexporn| 亚洲综合精品二区| 一个人免费看片子| 狠狠婷婷综合久久久久久88av| 18禁裸乳无遮挡动漫免费视频| 热99国产精品久久久久久7| 久久人人爽人人爽人人片va| 嫩草影院入口| 精品亚洲乱码少妇综合久久| 熟女电影av网| xxx大片免费视频| 丝袜脚勾引网站| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 黑人欧美特级aaaaaa片| 欧美3d第一页| 性高湖久久久久久久久免费观看| av免费在线看不卡| 在现免费观看毛片| 一区二区日韩欧美中文字幕 | 国产亚洲欧美精品永久| 中文乱码字字幕精品一区二区三区| 亚洲三级黄色毛片| 最新的欧美精品一区二区| 看免费av毛片| 亚洲伊人色综图| 满18在线观看网站| 欧美少妇被猛烈插入视频| 婷婷成人精品国产| 精品国产一区二区三区久久久樱花| 黑人高潮一二区| 久久精品夜色国产| 国产精品不卡视频一区二区| 欧美人与性动交α欧美精品济南到 | 99热这里只有是精品在线观看| 伦精品一区二区三区| 亚洲欧洲国产日韩| 国产在视频线精品| 精品人妻在线不人妻| av在线观看视频网站免费| 午夜福利视频精品| 日韩视频在线欧美| 18在线观看网站| av卡一久久| 成人综合一区亚洲| 国产黄频视频在线观看| 国产1区2区3区精品| 日日啪夜夜爽| 精品人妻偷拍中文字幕| 天美传媒精品一区二区| 亚洲少妇的诱惑av| 一级片免费观看大全| 在线观看美女被高潮喷水网站| 亚洲图色成人| 久久久久久久久久久久大奶| 青青草视频在线视频观看| 亚洲国产精品国产精品| av网站免费在线观看视频| 久久青草综合色| 免费黄频网站在线观看国产| 22中文网久久字幕| 人人澡人人妻人| av在线app专区| 伊人亚洲综合成人网| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久久久免| 全区人妻精品视频| 伦理电影免费视频| 国产欧美日韩一区二区三区在线| 欧美日韩一区二区视频在线观看视频在线| 午夜激情久久久久久久| 99热6这里只有精品| 久久国产精品大桥未久av| 一级毛片电影观看| 天堂俺去俺来也www色官网| videosex国产| 亚洲欧洲日产国产| 国产欧美日韩一区二区三区在线| 国产成人精品久久久久久| av一本久久久久| 欧美成人午夜免费资源| 日本黄大片高清| 91精品国产国语对白视频| 国产亚洲av片在线观看秒播厂| 成人毛片60女人毛片免费| 男女高潮啪啪啪动态图| 久久精品国产亚洲av天美| 1024视频免费在线观看| 久久97久久精品| 一区二区日韩欧美中文字幕 | 天美传媒精品一区二区| 久久久久久久亚洲中文字幕| 国产精品国产三级专区第一集| 18禁在线无遮挡免费观看视频| 成年女人在线观看亚洲视频| 免费在线观看黄色视频的| 欧美国产精品va在线观看不卡| av视频免费观看在线观看| 免费播放大片免费观看视频在线观看| 男人爽女人下面视频在线观看| 大香蕉久久网| 国产男女内射视频| 18禁观看日本| 永久免费av网站大全| 日本爱情动作片www.在线观看| 欧美97在线视频| 两性夫妻黄色片 | 天天影视国产精品| 最黄视频免费看| 熟女人妻精品中文字幕| 色94色欧美一区二区| 久久久精品94久久精品| av黄色大香蕉| 婷婷色综合大香蕉| 一级黄片播放器| 日本-黄色视频高清免费观看| 看免费av毛片| 国产av一区二区精品久久| 亚洲欧洲国产日韩| 建设人人有责人人尽责人人享有的| 香蕉丝袜av| 日本爱情动作片www.在线观看| 亚洲三级黄色毛片| 99热国产这里只有精品6| 51国产日韩欧美| 99国产综合亚洲精品| 久久综合国产亚洲精品| tube8黄色片| 99香蕉大伊视频| 一边摸一边做爽爽视频免费| 国产av码专区亚洲av| 黄片播放在线免费| 男人添女人高潮全过程视频| 视频区图区小说| 天堂俺去俺来也www色官网| 久久人妻熟女aⅴ| 精品卡一卡二卡四卡免费| 亚洲精品一二三| 亚洲av福利一区| 九九爱精品视频在线观看| 欧美日韩成人在线一区二区| 五月天丁香电影| 老司机影院毛片| 国产综合精华液| 赤兔流量卡办理| 国产精品人妻久久久影院| 亚洲美女搞黄在线观看| 免费人妻精品一区二区三区视频| 欧美国产精品va在线观看不卡| 免费av不卡在线播放| 黄色怎么调成土黄色| 考比视频在线观看| 爱豆传媒免费全集在线观看| 亚洲国产成人一精品久久久| 国产黄色免费在线视频| 日韩电影二区| 精品国产乱码久久久久久小说| 日韩一区二区三区影片| 亚洲精品aⅴ在线观看| av电影中文网址| 人人妻人人澡人人爽人人夜夜| 天美传媒精品一区二区| 久久精品国产a三级三级三级| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久午夜乱码| 在线观看一区二区三区激情| 亚洲欧洲国产日韩| 男人爽女人下面视频在线观看| 一个人免费看片子| 99久久人妻综合| 丝袜人妻中文字幕| 国产老妇伦熟女老妇高清| 成人亚洲欧美一区二区av| 观看美女的网站| 国产精品人妻久久久久久| 九色成人免费人妻av| 亚洲国产精品国产精品| 69精品国产乱码久久久| 亚洲精品成人av观看孕妇| 亚洲欧美精品自产自拍| 丝袜美足系列| 99re6热这里在线精品视频| 欧美精品av麻豆av| 成年女人在线观看亚洲视频| 成人国产av品久久久| 如何舔出高潮| 亚洲精品,欧美精品| 水蜜桃什么品种好| 欧美xxⅹ黑人| 精品久久久精品久久久| 亚洲,欧美,日韩| 精品一区在线观看国产| 日韩视频在线欧美| 亚洲一级一片aⅴ在线观看| 久久精品国产鲁丝片午夜精品| av免费观看日本| 日韩一区二区三区影片| 中文天堂在线官网| 最近手机中文字幕大全| 亚洲经典国产精华液单| 看免费av毛片| 欧美人与善性xxx| 欧美性感艳星| 久久这里有精品视频免费| 精品久久蜜臀av无| 国产色爽女视频免费观看| 久久这里只有精品19| 久久久久久久国产电影| 日本av免费视频播放| 久久久精品区二区三区| 日韩人妻精品一区2区三区| 女性被躁到高潮视频| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| av视频免费观看在线观看| 国产国拍精品亚洲av在线观看| av在线老鸭窝| 日本vs欧美在线观看视频| 国产欧美日韩一区二区三区在线| 久久99蜜桃精品久久| 国产日韩欧美视频二区| 久久免费观看电影| 国产 精品1| 欧美精品一区二区大全| 热99国产精品久久久久久7| 亚洲一码二码三码区别大吗| 亚洲高清免费不卡视频| 99国产精品免费福利视频| 国产高清国产精品国产三级| 国产不卡av网站在线观看| 免费av中文字幕在线| 亚洲欧美色中文字幕在线| 各种免费的搞黄视频| 新久久久久国产一级毛片| 精品一区二区三区四区五区乱码 | 麻豆乱淫一区二区| 亚洲精品日韩在线中文字幕| 人妻少妇偷人精品九色| 欧美xxⅹ黑人| 2021少妇久久久久久久久久久| 日本-黄色视频高清免费观看| 一本大道久久a久久精品| 插逼视频在线观看| 国产免费一区二区三区四区乱码| 性高湖久久久久久久久免费观看| 久久97久久精品| 精品酒店卫生间| 国产成人a∨麻豆精品| 久久久国产精品麻豆| 精品卡一卡二卡四卡免费| 久久久久久久久久人人人人人人| 国产av国产精品国产| 国产白丝娇喘喷水9色精品| 丝袜美足系列| 成人免费观看视频高清| 好男人视频免费观看在线| 亚洲av电影在线进入| videossex国产| 欧美日韩综合久久久久久| 美国免费a级毛片| 天堂中文最新版在线下载| 国产老妇伦熟女老妇高清| 人成视频在线观看免费观看| 国产精品熟女久久久久浪| 交换朋友夫妻互换小说| 国产精品99久久99久久久不卡 | 激情视频va一区二区三区| 两性夫妻黄色片 | www.av在线官网国产| 美女福利国产在线| 99国产综合亚洲精品| 大香蕉97超碰在线| 制服丝袜香蕉在线| 9热在线视频观看99| 久久久久久久久久久久大奶| 狂野欧美激情性bbbbbb| 超色免费av| 国产精品蜜桃在线观看| av福利片在线| 久久精品国产亚洲av天美| 七月丁香在线播放| 巨乳人妻的诱惑在线观看| 国产精品免费大片| a级片在线免费高清观看视频| 精品第一国产精品| 欧美日本中文国产一区发布| 中文字幕av电影在线播放| 波野结衣二区三区在线| 久久国内精品自在自线图片| 在线看a的网站| 久久精品夜色国产| 欧美日韩视频精品一区| 午夜精品国产一区二区电影| 视频区图区小说| 免费黄频网站在线观看国产| 久久久国产一区二区| 国产成人午夜福利电影在线观看| 国产成人av激情在线播放| 免费大片黄手机在线观看| 久久 成人 亚洲| 免费高清在线观看日韩| 久久97久久精品| 欧美+日韩+精品| 男男h啪啪无遮挡| av.在线天堂| 少妇 在线观看| 欧美日韩精品成人综合77777| 老熟女久久久| 亚洲成国产人片在线观看| 亚洲三级黄色毛片| 久久这里只有精品19| 精品99又大又爽又粗少妇毛片| 成人国语在线视频|