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

    基于改進(jìn)三維塊匹配濾波的紅外全息降噪算法

    2023-11-17 07:18:18李亞軒張永安王俊迪郗加民張?jiān)坪?/span>楊永浩
    激光與紅外 2023年10期
    關(guān)鍵詞:維納濾波全息圖散斑

    李亞軒,葉 飛,張永安,王俊迪,郗加民,張?jiān)坪?楊永浩

    (1.昆明理工大學(xué)理學(xué)院,云南 昆明 650500;2.云南省現(xiàn)代信息光學(xué)重點(diǎn)實(shí)驗(yàn)室,云南 昆明 650500;3.昆明榮者光電科技發(fā)展有限公司,云南 昆明 650500)

    1 引 言

    數(shù)字全息技術(shù)是一種將傳統(tǒng)光學(xué)全息技術(shù)與數(shù)字圖像技術(shù)相結(jié)合的技術(shù),因其可以記錄和顯示物體的三維信息并且易于使用,在眾多領(lǐng)域中得到了廣泛的應(yīng)用[1]。隨著數(shù)字全息技術(shù)在更加廣泛領(lǐng)域的應(yīng)用[2],紅外全息檢測(cè)技術(shù)也發(fā)展起來。紅外全息技術(shù)對(duì)比可見光全息技術(shù)有著光源能量高[3]、光源波長(zhǎng)對(duì)抗復(fù)雜環(huán)境能力強(qiáng)、相干距離長(zhǎng)、干涉能力穩(wěn)定的優(yōu)點(diǎn)[4];使得紅外全息有不易受可見光干擾[5]、穿透煙霧火焰能力強(qiáng)[6]、對(duì)實(shí)驗(yàn)穩(wěn)定性要求低[7]、長(zhǎng)距離大目標(biāo)成像能力強(qiáng)[8]的優(yōu)點(diǎn)。使其成為解決無損檢測(cè)[9]與火場(chǎng)救援[10]等現(xiàn)有難以解決領(lǐng)域的最優(yōu)方法。

    紅外全息技術(shù)作為全息技術(shù)與紅外技術(shù)[11]的分支,在長(zhǎng)距離傳播時(shí)會(huì)使光線能量衰減嚴(yán)重,到達(dá)物光系統(tǒng)中的焦平面能量低,導(dǎo)致圖像質(zhì)量較差,圖像噪聲復(fù)雜,影響后續(xù)處理[12]。針對(duì)紅外全息圖進(jìn)行圖像處理的需求也在逐步增長(zhǎng),而對(duì)于紅外全息圖的傳統(tǒng)圖像處理方法主要集中于對(duì)噪聲的抑制;如在2016年王洋等人通過在整幅圖像中尋找與噪聲點(diǎn)具有結(jié)構(gòu)相似的像素區(qū)域作為濾波區(qū)域,將圖像中噪聲點(diǎn)進(jìn)行濾波放到整幅圖像中進(jìn)行考慮,這種非局部均值濾波方法能夠?qū)t外圖像的高斯噪聲進(jìn)行降噪,同時(shí)使得濾波后圖像信息失真程度低,該算法能夠在對(duì)噪聲進(jìn)行良好抑制的同時(shí)具有較好保持邊緣信息的性能[13]。2017年,張勇提出了一種基于動(dòng)靜態(tài)檢測(cè)算法的紅外圖像降噪算法。通過一種動(dòng)靜態(tài)檢測(cè)算法將圖像分成動(dòng)態(tài)圖像和靜態(tài)圖像,用改進(jìn)的自適應(yīng)維納濾波算法處理動(dòng)態(tài)圖像,用改進(jìn)的非局部均值(NL-means)降噪算法處理靜態(tài)圖像,并用現(xiàn)場(chǎng)可編程門陣列實(shí)現(xiàn)紅外圖像降噪系統(tǒng)設(shè)計(jì)。其峰值信噪比(Peak Signal-to-Noise Ratio,PSNR)和細(xì)節(jié)方差-背景方差比(Detail Variance -Background Variance,DV/BV)均高于經(jīng)典降噪算法,該算法能有效減少圖像中的高斯白噪聲,并能很好地保持圖像的邊界細(xì)節(jié)信息[14]。2017年,為了消除紅外全息中顯著的散斑噪聲,劉寧提出了一種基于使用拉普拉斯金字塔技術(shù),并在每層上進(jìn)行雙邊濾波器以區(qū)分噪聲和細(xì)節(jié)信息。從而對(duì)紅外數(shù)字全息圖散斑噪聲進(jìn)行抑制的新技術(shù);由于雙邊濾波器無法區(qū)分所有噪聲,因此設(shè)計(jì)了一個(gè)改進(jìn)的中值濾波器來過濾剩余噪聲。該技術(shù)在噪聲抑制中的紅外數(shù)字全息圖上有效工作,從而去除了紅外全息圖像中大多數(shù)的散斑噪聲[15]。2018年,國(guó)外Bianco,V.使用多種聯(lián)合三維塊匹配濾波(MLDH-BM3D),這是一種特別適合過濾數(shù)字全息(Digital Holography,DH)圖像的方法;它將多個(gè)DH捕獲與三維塊匹配濾波算法(Block-Matching and 3D Filtering,BM3D)相結(jié)合,利用對(duì)象特征的非常稀疏表示而設(shè)計(jì)的一種數(shù)字濾波器;這種方法對(duì)全息圖像中的散斑噪聲與高斯白噪聲有很好的抑制效果,能夠使紅外全息圖像實(shí)現(xiàn)了無噪聲數(shù)字全息重建[16]。而這些對(duì)于紅外圖像的降噪方法對(duì)紅外全息圖像降噪得到結(jié)果的信噪比較低并且細(xì)節(jié)缺失嚴(yán)重[17],并且其運(yùn)算效率與圖像質(zhì)量也存在不足。而BM3D作為較為優(yōu)秀的降噪方法[18],在紅外全息圖像處理中,利用BM3D算法的特點(diǎn)對(duì)圖像中的高斯噪聲與散斑噪聲同時(shí)處理能夠得到對(duì)圖像進(jìn)行更好的降噪處理[19],但是其對(duì)細(xì)節(jié)的保留與對(duì)紅外圖像的適應(yīng)性較差,需要對(duì)BM3D算法的各個(gè)部分進(jìn)行改進(jìn)使其適用于紅外全息圖像的降噪處理。

    本文基于三維塊匹配濾波算法對(duì)紅外全息圖像的降噪原理,提出了一種對(duì)高斯噪聲尤其是加性高斯白噪聲與紅外全息圖散斑噪聲同時(shí)進(jìn)行處理,提升濾波算法的對(duì)細(xì)節(jié)的保留與對(duì)紅外全息圖像的適應(yīng)性,得到針對(duì)紅外全息圖像得到更好降噪效果的改進(jìn)三維塊匹配濾波算法,通過與現(xiàn)有傳統(tǒng)方法進(jìn)行對(duì)比來驗(yàn)證算法的性能。

    2 基本原理

    2.1 紅外全息圖記錄與噪聲分布

    紅外全息技術(shù)是基于數(shù)字全息技術(shù)與紅外技術(shù)結(jié)合的全息技術(shù)。紅外全息技術(shù)的基本原理與可見光波段的激光為光源的全息技術(shù)基本原理相同,都是對(duì)全息基本原理的運(yùn)用,即對(duì)波前信息的采集記錄與重構(gòu)再現(xiàn)。

    紅外激光由于波長(zhǎng)在760 nm以上,不能被肉眼觀察到;為了收集紅外全息圖,基于馬赫-曾德爾(Mach-Zender)干涉光路對(duì)紅外全息光學(xué)裝置進(jìn)行改造,使用綠光作為引導(dǎo)與紅外光進(jìn)行擬合,使紅外光與綠光通過分束鏡所分束的兩束光的光路相同,通過綠光的可見的干涉條紋來判斷紅外條紋的清晰程度,紅外全息光學(xué)裝置如圖1所示。

    圖1 紅外全息光學(xué)裝置

    其中,紅外全息實(shí)驗(yàn)所使用的激光器是中心波長(zhǎng)為1064 nm的在連續(xù)模式下工作的半導(dǎo)體紅外激光器及波長(zhǎng)為632.8 nm的He-Ne激光器,其中He-Ne激光器用于與紅外光擬合后引導(dǎo)光路搭建;紅外激光最大功率為115.3 mW,橫模模式為TEM00,發(fā)出光束發(fā)散角為1.5 mrad,光束直徑為1.329 mm,相干長(zhǎng)度為20 cm。使用的光學(xué)鏡片均為普通透紅外透鏡,分束鏡Ⅰ分光比為1∶2,分束鏡Ⅱ分光比為1∶1,擴(kuò)束鏡放大倍率為25倍,針孔尺寸為10 μm。本文選用CCD分辨率為1920×1080,像素尺寸為2.9 μm,曝光時(shí)間設(shè)置為0.2 ms,再現(xiàn)距離為104.5 mm。樣品為投射式樣品“光”字,材質(zhì)為黑色不透光的亞克力板(如圖2),規(guī)格為50 mm×50 mm,厚度為3 mm,字大小20 mm×15 mm,字筆畫寬度為2.5 mm。

    圖2 透射式樣品“光”圖像

    在數(shù)字全息測(cè)量系統(tǒng)記錄過程中,當(dāng)激光照射到物體表面時(shí),由于激光的高相干性,任意物點(diǎn)的散射光與其他物點(diǎn)產(chǎn)生的散射光發(fā)生干涉,從而形成了無規(guī)則分布的散斑圖樣,這種噪聲被稱為散斑噪聲[20]。

    在數(shù)字全息系統(tǒng)中,散斑噪聲往往被看作是降低全息再現(xiàn)像質(zhì)量的一種乘性光學(xué)噪聲。[21]數(shù)字全息記錄過程中,設(shè)散射面上有N個(gè)獨(dú)立散射點(diǎn)源,散射表面上隨機(jī)散射點(diǎn)q(x,y),此時(shí)接收面P上基元光波的復(fù)振幅為:

    E(P)=α(x,y)exp[jφR(x,y)]exp[jφ(r)]

    (1)

    式中,α(x,y)為散射表面反射光的振幅;φR(x,y)表示散射表面不平整所引入的隨機(jī)相位;φ(r)表示光傳播相位因子。

    接收面P的復(fù)振幅由N個(gè)獨(dú)立散射點(diǎn)源光波傳播到該點(diǎn)的疊加。此時(shí)P點(diǎn)散斑噪聲引起的光強(qiáng)變化表示為:

    I(P)=|E(P)|2=

    (2)

    全息圖上記錄面上光強(qiáng)分布為:

    I(x,y)=U(x,y)U*(x,y)=|O|2+|R|2+O·R*+R·O*+IR(P)

    (3)

    式中,O與R分別代表物光與參考光;IR(P)代表噪聲引起的光強(qiáng)變化。按下式計(jì)算:

    IR(P)=I(P)+IGauss(P)

    (4)

    式中,I(P)為散斑噪聲所引起的光強(qiáng)變化;IGauss(P)為高斯噪聲等其他環(huán)境因素所引起的光強(qiáng)變化。

    2.2 三維塊匹配濾波算法

    三維塊匹配濾波算法(BM3D)借鑒了非局部均值方法的非局部塊匹配思想[22],通過將相似判定的二維圖像塊找出作為參考?jí)K,再將其通過三維變換組合成三維群組并對(duì)群組進(jìn)行協(xié)同濾波處理,最后將處理后的塊聚合到原圖像塊位置得到更好濾波結(jié)果的傳統(tǒng)降噪算法[23]。BM3D濾波算法的過程由兩個(gè)階段所構(gòu)成,即初步估計(jì)階段與最終估計(jì)階段,每一階段由塊匹配、協(xié)同濾波與聚合三部分構(gòu)成,如圖3所示。

    圖3 BM3D流程圖

    如圖3中初步估計(jì)階段為基礎(chǔ)估計(jì)。首先是進(jìn)行按塊估計(jì),對(duì)于噪聲圖像中的每一個(gè)塊,首先利用塊匹配參數(shù)找到所有與圖像相似的相似塊,塊與塊之間的匹配度可以表示為這些系數(shù)的均分誤差,即:

    (5)

    (6)

    (7)

    然后通過相似判定將相似的二維圖像塊三維變換組成三維群組作為塊匹配(Block-matching)。隨后進(jìn)行協(xié)同硬閾值濾波處理,對(duì)三維群組通過硬閾值濾波進(jìn)行降噪,再進(jìn)行三維反變換得到所有分組的塊估計(jì),并將其返回初始位置。BM3D算法中使用的為硬閾值處理,利用硬閾值濾波來對(duì)圖像進(jìn)行初步估計(jì),估計(jì)可用式(8)計(jì)算得到:

    (8)

    (9)

    (10)

    最終估計(jì)階段為使用基礎(chǔ)估計(jì)后進(jìn)行改進(jìn)的分組操作與協(xié)同維納濾波。通過按塊估計(jì),對(duì)于每一個(gè)塊的分組上,使用基礎(chǔ)估計(jì)中得到的塊進(jìn)行塊匹配(BM)找到與當(dāng)前處理塊近似的所有塊的位置,使用這些位置形成源于噪聲圖像與基礎(chǔ)估計(jì)的兩個(gè)數(shù)組,對(duì)上述兩個(gè)數(shù)組進(jìn)行三維變換。隨后利用協(xié)同維納濾波處理,使用基礎(chǔ)估計(jì)的能量譜對(duì)原始噪聲圖像進(jìn)行維納濾波,對(duì)濾波后的三維系數(shù)進(jìn)行三維逆變換來生成所有分組塊的估計(jì),并將其返回初始位置。最后將加權(quán)平均對(duì)所有局部估計(jì)進(jìn)行聚合,計(jì)算出真實(shí)圖像的最終估計(jì)。即為:

    (11)

    (12)

    2.3 改進(jìn)BM3D算法

    為了使BM3D算法更好應(yīng)用于紅外全息圖像的降噪,采用對(duì)硬閾值濾波和維納濾波同時(shí)改進(jìn)的方式來改進(jìn)BM3D算法。

    由于硬閾值濾波需要對(duì)每一個(gè)塊進(jìn)行計(jì)算,而圖像噪聲和細(xì)節(jié)擁有一定的相似度,會(huì)導(dǎo)致降噪的同時(shí)出現(xiàn)細(xì)節(jié)的丟失。在硬閾值濾波前可先利用傅立葉閾值和再濾波對(duì)圖像中的細(xì)節(jié)信息進(jìn)行增強(qiáng),從而將細(xì)節(jié)信息與噪聲區(qū)分開來。然后再濾波,在實(shí)現(xiàn)去噪的同時(shí)對(duì)細(xì)節(jié)信息的有效保留。

    (13)

    (14)

    其中,aΔ為3σ檢驗(yàn),F表示傅里葉變換:

    (15)

    z1=g(xi,yi)*z

    (16)

    (17)

    (18)

    (19)

    (20)

    (21)

    (22)

    通過改進(jìn)BM3D中硬閾值部分與維納濾波部分,使得BM3D算法對(duì)紅外全息圖像的適應(yīng)性與濾波效果提升,并且其保留較多的細(xì)節(jié),能夠更好的對(duì)紅外全息圖像中的噪聲進(jìn)行降噪處理。

    3 實(shí)驗(yàn)驗(yàn)證

    為驗(yàn)證本文提出方法的有效性,選用傳統(tǒng)BM3D濾波、非局部均值濾波、自適應(yīng)維納濾波和傳統(tǒng)BM3D濾波作為對(duì)照組,與本文的改進(jìn)BM3D對(duì)比。實(shí)驗(yàn)中使用實(shí)拍紅外全息圖像(圖4)的濾波結(jié)果與模擬的潔凈環(huán)境下紅外全息圖像(圖5)進(jìn)行了對(duì)比,使用峰值信噪比(PSNR)與結(jié)構(gòu)相似性指標(biāo)(Structural Similarity Index Measure,SSIM)作為濾波效果優(yōu)劣的評(píng)判標(biāo)準(zhǔn)。

    圖4 實(shí)拍紅外全息圖像

    圖5 實(shí)拍圖像再現(xiàn)像

    圖4為圖1所示的紅外光學(xué)裝置對(duì)透射式樣品在大場(chǎng)景下的紅外全息圖像,模擬了在大場(chǎng)景下所呈現(xiàn)的具體圖像,與可見光全息對(duì)比,紅外技術(shù)對(duì)于大視場(chǎng)全息效果更好。圖5為實(shí)拍圖像的再現(xiàn)像,可以看出除了理想環(huán)境下的正負(fù)一級(jí)與零級(jí)像外,還存在其他部分的圖案,這些圖案就是全息成像時(shí)存在的環(huán)境噪聲與乘性散斑噪聲。圖6為添加噪聲方差0.02的乘性噪聲后的紅外全息圖像,模擬了復(fù)雜環(huán)境下噪聲對(duì)紅外全息圖像的影響。圖7為圖6圖像的紅外全息圖像再現(xiàn)像,可以看出噪聲部分明顯變亮并影響了圖中再現(xiàn)像“光”的分辨率,其證明噪聲對(duì)紅外全息圖像無論是原始圖像還是再現(xiàn)像都有所影響。

    圖6 添加噪聲方差0.02的乘性噪聲紅外全息圖像

    圖7 添加噪聲后紅外全息圖像再現(xiàn)像

    為了檢驗(yàn)改進(jìn)算法針對(duì)不同類型噪聲的降噪效果與算法性能,本文在采集的全息圖中人工添加了不同強(qiáng)度的乘性噪聲(包括高斯白噪聲與乘性散斑噪聲),并對(duì)其進(jìn)行了降噪處理。圖8中,由上至下,各行圖像添加了噪聲方差為0,0.02,0.05,0.2,1的乘性噪聲;自左向右,各列圖像所采用的降噪方法分別為3×3鄰域?yàn)V波、自適應(yīng)維納濾波、BM3D濾波和改進(jìn)BM3D濾波。從圖8中可以看出,隨著噪聲方差的增大四種方法的降噪效果均有所下降。其中,3×3鄰域?yàn)V波與BM3D濾波的降噪效果在噪聲方差等于0.05時(shí)明顯下降。自適應(yīng)維納濾波與改進(jìn)BM3D濾波降噪效果則變化較小。當(dāng)噪聲方差達(dá)到0.2時(shí),由于3×3鄰域?yàn)V波與BM3D濾波的降噪性能不足以去除大多數(shù)的噪聲,對(duì)于BM3D濾波方法當(dāng)增加性能則時(shí)也會(huì)導(dǎo)致塊匹配部分過度整合的問題;上述問題使得3×3鄰域?yàn)V波與BM3D濾波降噪后的圖像完全看不到圖中的“光”,表明此時(shí)這兩種濾波已無法實(shí)現(xiàn)有效降噪。自適應(yīng)維納濾波降噪后的圖像,在方差為0.8時(shí),雖能看清“光”字,但圖像清晰度已有明顯下降;在方差為1時(shí),自適應(yīng)維納濾波會(huì)同時(shí)對(duì)圖像與噪聲產(chǎn)生過度抑制,使得圖像清晰度大幅下降,僅部分“光”可隱約看到。表明該方法降噪效果優(yōu)于前兩種方法。改進(jìn)BM3D濾波的降噪圖像隨著噪聲方差的增大,圖像對(duì)比度及“光”字清晰度均有所降低,但圖像質(zhì)量下降小,即使噪聲方差達(dá)到1仍然可以清晰看到圖像中的“光”。表明改進(jìn)BM3D濾波算法可以在降噪的同時(shí)有效保留原圖像中細(xì)節(jié)信息,降噪效果明顯優(yōu)于其它三種降噪算法。

    圖9為采不同濾波方法的PSNR值數(shù)據(jù)圖。自適應(yīng)維納濾波與改進(jìn)BM3D濾波在噪聲方差0.05時(shí)其PSNR值變化較小;而在噪聲方差為0.2時(shí),自適應(yīng)維納濾波的PSNR出現(xiàn)了明顯下降,當(dāng)噪聲方差達(dá)到1時(shí),其PSNR值已經(jīng)下降至25以下,表示此時(shí)圖像質(zhì)量已經(jīng)出現(xiàn)明顯的下降,不能很好的分辨圖中的細(xì)節(jié),但依然優(yōu)于前兩種方法。改進(jìn)BM3D濾波的PSNR值隨著噪聲方差的變大逐漸減小,表示其圖像的分辨率也有所降低,但是其下降的幅度較小并且下降曲線較為平穩(wěn);在0.02噪聲方差之后其PSNR值沒有出現(xiàn)大幅度變化,即使在噪聲方差為1時(shí)PSNR仍然保持在40左右;通過PSNR值得變化,證明改進(jìn)BM3D濾波可以在降噪的同時(shí)能夠有效的保留圖像的細(xì)節(jié)信息,降噪效果明顯優(yōu)于其他三種降噪方法。

    圖9 不同濾波方法的PSNR值數(shù)據(jù)圖

    圖10中(a)為非局部均值濾波對(duì)紅外全息實(shí)拍圖的濾波后的再現(xiàn)像圖像,通過再現(xiàn)像可以觀察到非局部均值濾波算法對(duì)再現(xiàn)像中的噪聲沒有很好的去除效果,其對(duì)紅外全息中的乘性噪聲與環(huán)境噪聲降噪效果不是很理想。(b)為自適應(yīng)維納濾波對(duì)紅外全息實(shí)拍圖的濾波圖像后的再現(xiàn)像圖像,圖像噪聲有部分除去但依然有噪聲的殘留,字體較為清晰可見。(c)與(d)分別為傳統(tǒng)BM3D濾波與改進(jìn)BM3D濾波對(duì)對(duì)紅外全息實(shí)拍圖的濾波后的再現(xiàn)像圖像,圖像中噪聲降噪效果較好,可以表明這兩種濾波已可以對(duì)絕大多數(shù)環(huán)境噪聲與散斑噪聲進(jìn)行有效的降噪。而對(duì)于“光”字體,傳統(tǒng)BM3D濾波的成像效果明出現(xiàn)有部分模糊,其對(duì)再現(xiàn)像的細(xì)節(jié)的保留不夠;而改進(jìn)性BM3D濾波的再現(xiàn)像可以清晰可見的看到圖像中的“光”且噪聲去除程度較大。表明改進(jìn)BM3D濾波算法可以在降噪的同時(shí)有效保留原圖像中細(xì)節(jié)信息,降噪效果明顯優(yōu)于其它三種降噪算法。

    圖10 不同濾波方法對(duì)原始圖像的降噪后的再現(xiàn)像

    4 結(jié) 論

    紅外全息圖像由于紅外全息技術(shù)圖像的環(huán)境噪聲和散斑噪聲較大的缺點(diǎn),而目前的降噪方法無法滿足降噪需求,所以本文提出通過引入全局傅立葉閾值和再濾波技術(shù)與自適應(yīng)維納濾波對(duì)BM3D降噪方法改進(jìn),使其能夠?qū)t外全息圖像進(jìn)行降噪。通過對(duì)比現(xiàn)有方法降噪效果利用PSNR與SSIM進(jìn)行質(zhì)量評(píng)估來驗(yàn)證,并通過全息再現(xiàn)像圖像中噪聲區(qū)域與物體區(qū)域的比對(duì),證明改進(jìn)三維塊匹配技術(shù)能夠?qū)t外全息圖像降噪得到更好的圖像與再現(xiàn)像并且保留更多圖像細(xì)節(jié)。

    猜你喜歡
    維納濾波全息圖散斑
    激光顯示中的彩色散斑測(cè)量研究
    激光投影顯示散斑抑制方法研究
    多級(jí)維納濾波器的快速實(shí)現(xiàn)方法研究
    自適應(yīng)迭代維納濾波算法
    用于檢驗(yàn)散斑協(xié)方差矩陣估計(jì)性能的白化度評(píng)價(jià)方法
    基于多窗譜估計(jì)的改進(jìn)維納濾波語音增強(qiáng)
    基于維納濾波器的去噪研究
    電鏡成像方式對(duì)數(shù)字散斑相關(guān)方法結(jié)果的影響
    能觸摸的全息圖
    能觸摸的全息圖
    精品国产一区二区三区久久久樱花| 三级国产精品欧美在线观看| 精华霜和精华液先用哪个| 又粗又硬又长又爽又黄的视频| 日日爽夜夜爽网站| 日韩一区二区三区影片| 久久毛片免费看一区二区三区| 国产精品久久久久久久电影| 九九久久精品国产亚洲av麻豆| 亚洲精品久久午夜乱码| 韩国高清视频一区二区三区| 欧美人与善性xxx| 22中文网久久字幕| √禁漫天堂资源中文www| 精品少妇内射三级| 欧美 日韩 精品 国产| 一级毛片aaaaaa免费看小| 国产av精品麻豆| 国产精品成人在线| 国产黄色免费在线视频| 中文天堂在线官网| 中文天堂在线官网| av播播在线观看一区| 国产精品免费大片| 日韩亚洲欧美综合| 成人综合一区亚洲| 夜夜看夜夜爽夜夜摸| 国产伦精品一区二区三区视频9| 久久久久久久精品精品| 国产精品蜜桃在线观看| 国产探花极品一区二区| 日本av免费视频播放| 久久久久人妻精品一区果冻| 国产精品一二三区在线看| 亚洲,一卡二卡三卡| 简卡轻食公司| 免费看不卡的av| 成人毛片60女人毛片免费| 我的女老师完整版在线观看| 六月丁香七月| 国产日韩欧美亚洲二区| 中文字幕人妻丝袜制服| 亚洲自偷自拍三级| 夜夜骑夜夜射夜夜干| 日日爽夜夜爽网站| kizo精华| 成人美女网站在线观看视频| 久久国产精品男人的天堂亚洲 | 少妇精品久久久久久久| av黄色大香蕉| 一个人免费看片子| 日韩精品免费视频一区二区三区 | 国产精品秋霞免费鲁丝片| av一本久久久久| 能在线免费看毛片的网站| 国产黄色视频一区二区在线观看| 久久精品国产鲁丝片午夜精品| 少妇被粗大的猛进出69影院 | 嫩草影院入口| 丰满饥渴人妻一区二区三| 国产亚洲av片在线观看秒播厂| 亚洲国产日韩一区二区| a级毛片在线看网站| 亚洲av.av天堂| 午夜福利视频精品| 三级国产精品片| 中文在线观看免费www的网站| 男人添女人高潮全过程视频| 欧美日本中文国产一区发布| 99精国产麻豆久久婷婷| 日韩熟女老妇一区二区性免费视频| 亚洲欧美日韩东京热| h视频一区二区三区| a级一级毛片免费在线观看| 午夜免费男女啪啪视频观看| 美女视频免费永久观看网站| 极品教师在线视频| 日本欧美国产在线视频| 国产亚洲一区二区精品| 伊人久久国产一区二区| h视频一区二区三区| 欧美+日韩+精品| 亚洲成色77777| 午夜福利在线观看免费完整高清在| 伊人亚洲综合成人网| 如日韩欧美国产精品一区二区三区 | 中文精品一卡2卡3卡4更新| 中国国产av一级| 嫩草影院入口| 国产精品久久久久久精品电影小说| 青春草国产在线视频| 精品亚洲乱码少妇综合久久| 成人午夜精彩视频在线观看| 在线播放无遮挡| 欧美xxⅹ黑人| 青春草亚洲视频在线观看| 日韩一区二区三区影片| 日韩成人伦理影院| 全区人妻精品视频| 中文欧美无线码| 精品久久久久久久久亚洲| 少妇丰满av| 五月玫瑰六月丁香| av天堂久久9| 亚洲成人手机| 国产精品国产三级国产专区5o| 高清毛片免费看| 成人综合一区亚洲| 国产伦在线观看视频一区| av福利片在线观看| a级毛片在线看网站| 日本wwww免费看| 亚洲内射少妇av| 一二三四中文在线观看免费高清| 亚洲四区av| 国产精品伦人一区二区| 久久久a久久爽久久v久久| a级片在线免费高清观看视频| 美女脱内裤让男人舔精品视频| 午夜日本视频在线| 日韩强制内射视频| 18禁动态无遮挡网站| 女性生殖器流出的白浆| 国产精品三级大全| 久久久久久久国产电影| 在线观看免费视频网站a站| 国产91av在线免费观看| 亚洲欧美成人综合另类久久久| 欧美精品一区二区免费开放| 九九久久精品国产亚洲av麻豆| 亚洲精品国产色婷婷电影| 美女福利国产在线| av女优亚洲男人天堂| 99久久人妻综合| 丝袜喷水一区| 国产精品99久久久久久久久| 亚洲第一区二区三区不卡| 自拍偷自拍亚洲精品老妇| 啦啦啦在线观看免费高清www| 国产在线男女| 国产白丝娇喘喷水9色精品| 免费人妻精品一区二区三区视频| 国产国拍精品亚洲av在线观看| 91精品国产国语对白视频| 国产69精品久久久久777片| a级一级毛片免费在线观看| 日韩,欧美,国产一区二区三区| 少妇人妻 视频| 日韩成人av中文字幕在线观看| 亚洲欧美中文字幕日韩二区| 精品一区二区三区视频在线| 国产精品久久久久久精品古装| 亚洲性久久影院| 亚洲一区二区三区欧美精品| 亚州av有码| 日日摸夜夜添夜夜爱| 免费久久久久久久精品成人欧美视频 | 国产精品女同一区二区软件| 韩国av在线不卡| 成人美女网站在线观看视频| 国产男女超爽视频在线观看| 精品久久久久久久久亚洲| 一级毛片电影观看| 男男h啪啪无遮挡| 亚洲精品第二区| 赤兔流量卡办理| 在线 av 中文字幕| tube8黄色片| √禁漫天堂资源中文www| 老司机影院成人| 免费大片黄手机在线观看| 国产精品熟女久久久久浪| 最近中文字幕2019免费版| 妹子高潮喷水视频| 欧美另类一区| 男女边摸边吃奶| 18+在线观看网站| 下体分泌物呈黄色| 一边亲一边摸免费视频| 国产精品一区二区性色av| 国产精品久久久久久精品古装| 久久影院123| 国产精品欧美亚洲77777| 韩国高清视频一区二区三区| 欧美bdsm另类| 亚洲国产精品999| 黑人巨大精品欧美一区二区蜜桃 | 最近最新中文字幕免费大全7| 久久鲁丝午夜福利片| 国产成人freesex在线| 国产黄色免费在线视频| 精品一区二区免费观看| 在线看a的网站| 久久精品熟女亚洲av麻豆精品| 曰老女人黄片| 国产精品女同一区二区软件| av播播在线观看一区| 夫妻性生交免费视频一级片| 日本欧美视频一区| 欧美丝袜亚洲另类| av国产久精品久网站免费入址| 国产男人的电影天堂91| 亚洲精品视频女| 麻豆成人av视频| 久久久久久久久久久丰满| 最近的中文字幕免费完整| 丰满少妇做爰视频| 亚洲欧美日韩另类电影网站| 成年人午夜在线观看视频| 久久99热这里只频精品6学生| 日韩强制内射视频| 日本猛色少妇xxxxx猛交久久| 日韩制服骚丝袜av| 少妇丰满av| 深夜a级毛片| 国产黄频视频在线观看| 亚洲欧美精品自产自拍| 熟女电影av网| 一级毛片黄色毛片免费观看视频| .国产精品久久| 久久久久久久久久久免费av| 伦理电影大哥的女人| 国产探花极品一区二区| 午夜视频国产福利| 免费看不卡的av| a级片在线免费高清观看视频| 久久婷婷青草| 久久国产精品男人的天堂亚洲 | 两个人的视频大全免费| 日本免费在线观看一区| 国产av精品麻豆| 国产熟女欧美一区二区| 久久鲁丝午夜福利片| 高清视频免费观看一区二区| 亚洲精品一二三| 亚洲中文av在线| 精品久久久噜噜| 黑人巨大精品欧美一区二区蜜桃 | 又黄又爽又刺激的免费视频.| av女优亚洲男人天堂| 久久久久国产精品人妻一区二区| 国产精品嫩草影院av在线观看| 日韩免费高清中文字幕av| 最近手机中文字幕大全| 亚洲综合精品二区| 一本大道久久a久久精品| 日本av手机在线免费观看| 午夜激情久久久久久久| 在线观看三级黄色| 国产伦精品一区二区三区视频9| 久久精品久久久久久噜噜老黄| 亚洲美女黄色视频免费看| av播播在线观看一区| 国产真实伦视频高清在线观看| 久久久久久久久久久免费av| 欧美97在线视频| av在线app专区| 欧美日韩国产mv在线观看视频| 日韩强制内射视频| 亚洲精品日韩在线中文字幕| 欧美精品一区二区大全| 亚洲精品视频女| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 看非洲黑人一级黄片| 五月伊人婷婷丁香| 亚洲欧洲日产国产| 亚洲国产最新在线播放| 26uuu在线亚洲综合色| 国产国拍精品亚洲av在线观看| 免费观看性生交大片5| 国产精品嫩草影院av在线观看| 纯流量卡能插随身wifi吗| 人人妻人人看人人澡| 熟女电影av网| 一区二区三区精品91| 国产欧美另类精品又又久久亚洲欧美| 黄色毛片三级朝国网站 | 久久女婷五月综合色啪小说| 国产精品福利在线免费观看| 日韩av免费高清视频| 美女福利国产在线| 99热全是精品| 尾随美女入室| 成人国产麻豆网| 一级二级三级毛片免费看| 高清欧美精品videossex| 免费黄网站久久成人精品| 久久久欧美国产精品| 五月开心婷婷网| 这个男人来自地球电影免费观看 | 熟女人妻精品中文字幕| 亚洲av不卡在线观看| 久久久国产欧美日韩av| 国产黄片美女视频| 日本与韩国留学比较| 国产精品人妻久久久影院| 精品卡一卡二卡四卡免费| 高清在线视频一区二区三区| 在线精品无人区一区二区三| 国产黄色免费在线视频| 啦啦啦在线观看免费高清www| 国产av国产精品国产| 国产亚洲一区二区精品| 午夜免费观看性视频| 大话2 男鬼变身卡| 又大又黄又爽视频免费| 成人影院久久| 高清欧美精品videossex| 丝袜在线中文字幕| 国产免费又黄又爽又色| 色网站视频免费| 国产一区二区在线观看日韩| 国产精品免费大片| 国产免费一区二区三区四区乱码| 一级av片app| 国产毛片在线视频| 亚洲精品自拍成人| 97超视频在线观看视频| 中文字幕av电影在线播放| 晚上一个人看的免费电影| 午夜视频国产福利| 成人亚洲欧美一区二区av| 少妇人妻一区二区三区视频| 久久ye,这里只有精品| 亚洲欧美清纯卡通| 纯流量卡能插随身wifi吗| 精品人妻熟女毛片av久久网站| 国产黄色免费在线视频| 激情五月婷婷亚洲| 日韩亚洲欧美综合| 十八禁高潮呻吟视频 | 国产亚洲91精品色在线| 卡戴珊不雅视频在线播放| 三级国产精品片| 精品久久久噜噜| 色视频www国产| 菩萨蛮人人尽说江南好唐韦庄| 午夜激情福利司机影院| 色5月婷婷丁香| 91精品伊人久久大香线蕉| 多毛熟女@视频| 美女xxoo啪啪120秒动态图| 免费黄色在线免费观看| 久久久精品免费免费高清| 99视频精品全部免费 在线| 精品国产一区二区久久| 久久97久久精品| 如何舔出高潮| 国产白丝娇喘喷水9色精品| 久久精品国产自在天天线| 久久国产亚洲av麻豆专区| 亚洲av综合色区一区| 人妻一区二区av| 国产午夜精品久久久久久一区二区三区| 久久午夜综合久久蜜桃| 亚洲人与动物交配视频| 午夜免费观看性视频| 亚洲激情五月婷婷啪啪| 国产亚洲午夜精品一区二区久久| 性色avwww在线观看| videos熟女内射| 日韩电影二区| 国产有黄有色有爽视频| 亚洲真实伦在线观看| 免费不卡的大黄色大毛片视频在线观看| 日韩免费高清中文字幕av| 中文在线观看免费www的网站| 肉色欧美久久久久久久蜜桃| 一本—道久久a久久精品蜜桃钙片| 91精品伊人久久大香线蕉| 欧美成人午夜免费资源| 欧美区成人在线视频| 十分钟在线观看高清视频www | 久久久久久人妻| 亚洲av二区三区四区| 女人精品久久久久毛片| 国内精品宾馆在线| 国产淫语在线视频| 日韩av不卡免费在线播放| 久久影院123| 婷婷色av中文字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 一个人免费看片子| 国产成人精品久久久久久| 国产男人的电影天堂91| 国产 精品1| 午夜福利,免费看| 国产成人精品无人区| 亚洲自偷自拍三级| 在线免费观看不下载黄p国产| 人妻一区二区av| 国产一区二区在线观看日韩| 午夜福利视频精品| 老女人水多毛片| 黄色日韩在线| 美女国产视频在线观看| 国产亚洲精品久久久com| 国产淫语在线视频| 成人黄色视频免费在线看| 成人无遮挡网站| 老女人水多毛片| a级毛片在线看网站| 国产一区亚洲一区在线观看| 五月天丁香电影| 国产伦精品一区二区三区四那| 在线天堂最新版资源| 精品久久久久久久久av| 日本欧美国产在线视频| 久久久国产欧美日韩av| 天天操日日干夜夜撸| av天堂久久9| 午夜精品国产一区二区电影| 三级经典国产精品| 我的女老师完整版在线观看| 国产精品嫩草影院av在线观看| 午夜免费观看性视频| 欧美日韩在线观看h| 国产免费一区二区三区四区乱码| 成人美女网站在线观看视频| 久久热精品热| 99九九线精品视频在线观看视频| 国产精品久久久久久久久免| 免费大片黄手机在线观看| 七月丁香在线播放| 免费黄色在线免费观看| 在线看a的网站| 久久国产精品大桥未久av | 亚洲综合精品二区| 国产淫语在线视频| a级毛片在线看网站| 啦啦啦中文免费视频观看日本| 老熟女久久久| 国产精品一区二区性色av| 国产高清国产精品国产三级| 中文字幕制服av| 国产高清不卡午夜福利| 九色成人免费人妻av| 久久99热6这里只有精品| 内射极品少妇av片p| 色婷婷av一区二区三区视频| 水蜜桃什么品种好| 国产精品99久久久久久久久| 免费播放大片免费观看视频在线观看| 美女cb高潮喷水在线观看| 一个人免费看片子| 国产极品天堂在线| 在线观看免费高清a一片| 成人综合一区亚洲| 汤姆久久久久久久影院中文字幕| 国产免费视频播放在线视频| 一级毛片aaaaaa免费看小| 看非洲黑人一级黄片| 日韩亚洲欧美综合| 亚洲在久久综合| 亚洲精品色激情综合| 免费看av在线观看网站| 中文字幕人妻丝袜制服| 久久ye,这里只有精品| 蜜桃久久精品国产亚洲av| 国产女主播在线喷水免费视频网站| 综合色丁香网| 免费少妇av软件| 亚洲欧美精品专区久久| 乱系列少妇在线播放| 天堂中文最新版在线下载| 亚洲av不卡在线观看| 99热网站在线观看| 五月天丁香电影| 成年人免费黄色播放视频 | 国产综合精华液| 青春草视频在线免费观看| 亚洲精品日韩av片在线观看| 能在线免费看毛片的网站| 两个人免费观看高清视频 | 成人二区视频| 日韩欧美 国产精品| 久久 成人 亚洲| 80岁老熟妇乱子伦牲交| 少妇丰满av| 欧美xxxx性猛交bbbb| 久久av网站| 亚洲av福利一区| 精品99又大又爽又粗少妇毛片| 中文字幕人妻丝袜制服| 久久久久久久精品精品| 日韩精品免费视频一区二区三区 | 国产亚洲最大av| 精品亚洲成a人片在线观看| 欧美精品人与动牲交sv欧美| 国产成人a∨麻豆精品| 国产在线一区二区三区精| 亚洲真实伦在线观看| 亚洲综合精品二区| 成年人免费黄色播放视频 | 国产一区有黄有色的免费视频| 免费观看a级毛片全部| 在线 av 中文字幕| 免费av不卡在线播放| 欧美性感艳星| 亚洲第一av免费看| 亚洲欧美成人综合另类久久久| 少妇 在线观看| 免费观看av网站的网址| 日日爽夜夜爽网站| 少妇被粗大的猛进出69影院 | 亚洲va在线va天堂va国产| 精品少妇黑人巨大在线播放| 一级片'在线观看视频| 老熟女久久久| 国产熟女午夜一区二区三区 | 精品一区二区三区视频在线| 亚洲国产精品一区二区三区在线| 久久午夜福利片| 黄色配什么色好看| 亚洲欧洲精品一区二区精品久久久 | 99九九线精品视频在线观看视频| 久久国产精品男人的天堂亚洲 | 日日啪夜夜爽| 久久99蜜桃精品久久| 人妻少妇偷人精品九色| 大陆偷拍与自拍| 日本爱情动作片www.在线观看| av视频免费观看在线观看| 国产精品人妻久久久影院| 高清在线视频一区二区三区| 少妇丰满av| 欧美日韩亚洲高清精品| 亚洲av中文av极速乱| 国产 精品1| 天天操日日干夜夜撸| 午夜福利影视在线免费观看| 国产探花极品一区二区| 在线观看国产h片| 美女中出高潮动态图| 七月丁香在线播放| 老熟女久久久| 深夜a级毛片| 韩国高清视频一区二区三区| 亚洲精品自拍成人| 亚洲欧美中文字幕日韩二区| 美女脱内裤让男人舔精品视频| 精品少妇内射三级| 九草在线视频观看| 亚洲av电影在线观看一区二区三区| 成人午夜精彩视频在线观看| xxx大片免费视频| 久久久久久久亚洲中文字幕| 日日撸夜夜添| 国产爽快片一区二区三区| 精品久久久久久久久亚洲| 人妻系列 视频| 蜜臀久久99精品久久宅男| 亚洲国产精品成人久久小说| 欧美区成人在线视频| 美女内射精品一级片tv| 国产欧美亚洲国产| 校园人妻丝袜中文字幕| 色5月婷婷丁香| 丰满迷人的少妇在线观看| 国产精品不卡视频一区二区| 日韩强制内射视频| 18+在线观看网站| 免费久久久久久久精品成人欧美视频 | 在现免费观看毛片| 一级毛片黄色毛片免费观看视频| 少妇猛男粗大的猛烈进出视频| 伦精品一区二区三区| 18禁在线播放成人免费| 中文资源天堂在线| 有码 亚洲区| 中文字幕久久专区| 日韩大片免费观看网站| 精品人妻一区二区三区麻豆| 九九久久精品国产亚洲av麻豆| 中文字幕免费在线视频6| 免费观看a级毛片全部| 欧美成人精品欧美一级黄| 国产真实伦视频高清在线观看| 欧美3d第一页| 日日啪夜夜撸| 99热6这里只有精品| 国产永久视频网站| 在现免费观看毛片| 久久久久久久久久久久大奶| 亚洲国产精品成人久久小说| 高清黄色对白视频在线免费看 | 久久鲁丝午夜福利片| 久久青草综合色| www.av在线官网国产| 91精品国产国语对白视频| 岛国毛片在线播放| 一级二级三级毛片免费看| 两个人免费观看高清视频 | 久久精品国产亚洲av天美| 草草在线视频免费看| 老熟女久久久| 亚洲综合色惰| 涩涩av久久男人的天堂| 国产成人免费观看mmmm| 97在线人人人人妻| 成人国产av品久久久| 天堂8中文在线网| 麻豆成人av视频| 最近中文字幕高清免费大全6| 在线精品无人区一区二区三| 午夜免费鲁丝| 老司机影院毛片| 一本一本综合久久| 成人毛片60女人毛片免费| 国产精品熟女久久久久浪| 成人毛片60女人毛片免费| 欧美精品亚洲一区二区| 亚洲精品视频女| 丰满饥渴人妻一区二区三|