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

    改進(jìn)多尺度的Retinex紅外圖像增強(qiáng)

    2021-04-10 05:55:08魏然然詹偉達(dá)朱德鵬
    液晶與顯示 2021年3期
    關(guān)鍵詞:光暈濾波器尺度

    魏然然, 詹偉達(dá), 朱德鵬, 田 永

    (長(zhǎng)春理工大學(xué) 電子信息工程學(xué)院,吉林 長(zhǎng)春 130022)

    1 引 言

    近年來,紅外成像系統(tǒng)由于具有良好的隱蔽性、較強(qiáng)的抗干擾能力以及可晝夜工作等特點(diǎn)而被廣泛應(yīng)用在軍事和民用領(lǐng)域[1-3]。紅外圖像在獲取時(shí)由于受到紅外成像系統(tǒng)器件性能及波長(zhǎng)等因素的影響,常常會(huì)表現(xiàn)出對(duì)比度較低、邊緣細(xì)節(jié)不突出、視覺效果模糊等缺點(diǎn)[4]。這些缺點(diǎn)將嚴(yán)重影響到后續(xù)紅外圖像的觀察、目標(biāo)識(shí)別和目標(biāo)探測(cè)等工作。因此,利用紅外圖像增強(qiáng)算法增強(qiáng)圖像對(duì)比度、細(xì)節(jié)等已成為一種必不可少的圖像處理操作。

    常用的紅外圖像處理算法包括空域算法和頻域算法??沼蛩惴òㄖ狈綀D均衡化算法和灰度變換等;頻域算法包括傅里葉變換、小波變換、Culvelet變換和Contourlet變換等[5]。雖然這些算法都可以或多或少地對(duì)紅外圖像進(jìn)行增強(qiáng),但也存在一些缺點(diǎn),如直方圖均衡化算法和灰度變換在對(duì)紅外圖像增強(qiáng)的同時(shí)會(huì)放大噪聲。Land等人首次提出Retinex理論,Retinex理論是一種非線性算法,其處理紅外圖像時(shí)可以在保留細(xì)節(jié)信息的同時(shí)降低圖像噪聲[6]。在Land的理論基礎(chǔ)上,Jobson等人提出了一種單尺度的Retinex算法(SSR),因?yàn)樵摲椒ㄖ挥幸粋€(gè)參數(shù),故其不能很好地兼顧動(dòng)態(tài)范圍壓縮和色感一致性[7]。為了改進(jìn)SSR算法,Rahman等人隨后提出了多尺度Retinex算法(MSR)[8]。在該算法的基礎(chǔ)上研究者們提出了一系列的改進(jìn)算法。劉輝等人改進(jìn)傳統(tǒng)MSR算法,將算法中的指數(shù)用灰度線性拉伸代替,該算法雖然可以提高圖像的對(duì)比度,改善圖像的視覺效果,但其在操作過程中會(huì)擴(kuò)大圖像噪聲,同時(shí)產(chǎn)生“光暈”現(xiàn)象[9]。張承泓等人使用自適應(yīng)引導(dǎo)濾波器代替MSR算法中的高斯濾波器,該改進(jìn)算法可以較好地去除圖像中的“光暈”和噪聲[10]。本文在前人研究的基礎(chǔ)上,將改進(jìn)引導(dǎo)濾波器模型使用到傳統(tǒng)MSR算法中,以提出一種新的紅外圖像增強(qiáng)算法。

    2 基礎(chǔ)理論

    2.1 多尺度Retinex理論

    根據(jù)多尺度Retinex理論,一幅紅外圖像I(x,y)可以由入射分量L(x,y)和反射分量R(x,y)乘積得到,其數(shù)學(xué)表達(dá)式如下所示:

    I(x,y)=L(x,y)R(x,y)

    (1)

    對(duì)式(1)兩邊同時(shí)取對(duì)數(shù)可得:

    lnI(x,y)=ln(L(x,y)R(x,y))=
    lnL(x,y)+lnR(x,y)

    (2)

    當(dāng)L(x,y)=I(x,y)*G(x,y)時(shí),反射分量可表示為:

    lnR(x,y)=lnI(x,y)-
    ln(I(x,y)*G(x,y))

    (3)

    式中,*表示卷積運(yùn)算,G(x,y)為高斯濾波函數(shù),其定義為:

    (4)

    式中,λ為滿足時(shí)的歸一化常數(shù),δ為環(huán)繞參數(shù)。

    式(3)又被稱為單尺度Retinex算法(SSR)。由式(3)和式(4)可以看出,SSR只有單一的參數(shù)δ,該算法不能同時(shí)兼顧動(dòng)態(tài)范圍壓縮和色感一致性。為了解決這一問題,研究者們提出了多尺度Retinex算法(MSR),其主要思想是對(duì)每個(gè)尺度的SSR結(jié)果進(jìn)行不同權(quán)重的加權(quán)平均,其數(shù)學(xué)表達(dá)式為:

    (5)

    2.2 引導(dǎo)濾波

    引導(dǎo)濾波器可以對(duì)圖像進(jìn)行線性濾波,其具有一定的邊緣保持特性[11]。其數(shù)學(xué)模型如下:

    qi=akIi+bk,?i∈ωk

    (6)

    式中,qi為輸出圖像;Ii為引導(dǎo)圖像;i、k表示像素索引;ωk為半徑為r的局部方形窗口;ak和bk為ωk對(duì)應(yīng)的線性系數(shù)。

    對(duì)式(6)兩邊同時(shí)求梯度可得:

    qi=aIi

    (7)

    式中,a為常數(shù)。

    由式(7)可以看出qi與Ii具有相似的邊緣梯度,而Ii的邊緣平滑程度取決于ak。為了確定其值,需要使qi與pi的差異最小,則其代價(jià)函數(shù)如下所示[10]:

    (8)

    式中,ε為防止ak和bk過大的正則化因子。

    利用線性回歸法求解式(8)取最小值時(shí),ak和bk的值為:

    (9)

    (10)

    式中,|ω|為ωk中像素的數(shù)量;δk和μk為ωk中Ii的標(biāo)準(zhǔn)差和均值。

    假設(shè)Ii=pi,則式(9)和式(10)可化為:

    (11)

    bk=μk(1-ak)

    (12)

    對(duì)式(6)、式(11)、式(12)進(jìn)行分析可得,當(dāng)ε=0時(shí),ak=1,bk=0,故此時(shí)圖像保持不變;當(dāng)ε?δk時(shí),ak≈0,bk≈μk,qi≈μk,故當(dāng)ωk處于低方差區(qū)域時(shí),原始圖像得到了平滑處理;當(dāng)ε?δk時(shí),ak≈1,bk≈0,qi≈Ii,故當(dāng)ωk處于高方差區(qū)域時(shí),圖像可以進(jìn)行邊緣保持操作。由以上操作可以看出,ε大小可以決定圖像被平滑的程度,故ε又被稱為平滑因子,對(duì)比高斯濾波的特點(diǎn),引導(dǎo)濾波器也有決定其大小的尺度因子s。

    3 本文算法

    紅外圖像在獲取時(shí)由于受到紅外成像系統(tǒng)器件性能及波長(zhǎng)等因素的影響,常常會(huì)表現(xiàn)出對(duì)比度較低、邊緣細(xì)節(jié)不突出、視覺效果模糊等缺點(diǎn)。傳統(tǒng)的MSR在解決這些問題時(shí),會(huì)出現(xiàn)光暈現(xiàn)象和放大噪聲等問題。為了解決這些問題,本文提出一種改進(jìn)的MSR算法對(duì)紅外圖像進(jìn)行圖像增強(qiáng)。設(shè)待增強(qiáng)的紅外圖像為F(x,y),增強(qiáng)后的紅外圖像為f(x,y),改進(jìn)引導(dǎo)濾波器輸出的入射分量為L(zhǎng)(x,y),經(jīng)過MSR計(jì)算得到的反射分量為R(x,y),增強(qiáng)后的入射和反射分量分別為l(x,y)和r(x,y),其具體算法實(shí)現(xiàn)框圖如圖1所示。

    圖1 本文算法實(shí)現(xiàn)框圖Fig.1 Algorithm block diagram of this paper

    3.1 入射分量估計(jì)

    傳統(tǒng)多尺度Retinex算法在圖像邊緣發(fā)生躍變處,不能精確估計(jì)光照分量即反射分量,進(jìn)而產(chǎn)生“光暈”現(xiàn)象。為了解決這一問題,本文采用具有保邊平滑的改進(jìn)引導(dǎo)濾波器代替?zhèn)鹘y(tǒng)的高斯濾波器對(duì)光照層進(jìn)行精確估計(jì)。

    引導(dǎo)濾波器與高斯濾波器相比,雖然可以改善圖像中的“光暈”現(xiàn)象,但其每個(gè)窗口對(duì)權(quán)重的取值均相同,沒有考慮每個(gè)窗口的灰度差異。針對(duì)這一不足,對(duì)引導(dǎo)濾波器進(jìn)行改進(jìn),具體方法為:首先計(jì)算窗內(nèi)每個(gè)像素點(diǎn)的均值、標(biāo)準(zhǔn)差、平均梯度來作為該點(diǎn)特征,然后根據(jù)第k次的該點(diǎn)特征來求得權(quán)重值,最后對(duì)輸出的引導(dǎo)圖像加權(quán)平均。

    圖像像素點(diǎn)的均值、標(biāo)準(zhǔn)差、平均梯度定義如下:

    (13)

    (14)

    (15)

    則權(quán)重值為:

    (16)

    (17)

    式中,Uk(i,j)為質(zhì)量測(cè)度,其由上述均值、標(biāo)準(zhǔn)差、平均梯度的一個(gè)或多個(gè)組成。對(duì)不同圖像進(jìn)行處理時(shí),其側(cè)重點(diǎn)不同,α,β,γ取值不同。α,β,γ其值可以為0或1。當(dāng)α,β,γ中一個(gè)或兩個(gè)為0時(shí),說明圖像的質(zhì)量測(cè)度不采用此特性;當(dāng)α,β,γ其中一個(gè)或全部為1時(shí),表明此特性占有很大權(quán)重,圖像的質(zhì)量測(cè)度采取此特性。ε為一個(gè)很小的正數(shù)。

    改進(jìn)引導(dǎo)濾波輸出圖像為:

    (18)

    式中,Lk(i,j)為3次引導(dǎo)濾波對(duì)應(yīng)的輸出圖像。

    將原始圖像F和引導(dǎo)圖像I通過改進(jìn)引導(dǎo)濾波器處理后得到輸出圖像即入射圖像L為:

    L=L(F,I)

    (19)

    式中,L(·)為改進(jìn)引導(dǎo)濾波器操作。

    為了驗(yàn)證本文改進(jìn)引導(dǎo)濾波器的優(yōu)越性,使用本文改進(jìn)引導(dǎo)濾波器與傳統(tǒng)高斯濾波器處理估計(jì)得到的入射分量圖像,如圖2所示(高斯濾波器參數(shù)選擇:尺度因子取值為80;引導(dǎo)濾波參數(shù)選擇:尺度因子為80,平滑因子取值為0.002)。根據(jù)圖2可以看出,本文改進(jìn)的引導(dǎo)濾波器可以更準(zhǔn)確地對(duì)入射分量進(jìn)行估計(jì)。

    (a)原圖(a)Origin image

    (b) 傳統(tǒng)高斯濾波器估計(jì)(b) Traditional Gaussian filter estimation

    (c)改進(jìn)的引導(dǎo)濾波器估計(jì)(c)Improved guide filter estimation圖2 入射分量估計(jì)結(jié)果對(duì)比圖Fig.2 Comparison of estimated results of incident components

    3.2 入射分量自適應(yīng)灰度拉伸

    為了更好地提高入射分量的全局對(duì)比度,本文使用全局直方圖均衡化來實(shí)現(xiàn)對(duì)入射分量的自適應(yīng)灰度拉伸。對(duì)入射圖像進(jìn)行歸一化處理后的分布函數(shù)為:

    (20)

    式中,Lk為第k級(jí)像素的像素值,nk為第k級(jí)像素的個(gè)數(shù),n為像素總數(shù),N為最大像素級(jí)。

    若原入射分量的像素值為L(zhǎng),直方圖均衡化后的像素值為l,則其灰度變換公式為:

    (21)

    式中,Li表示第i級(jí)像素的原像素值,p(Li)表示Li的概率。

    根據(jù)以上理論,用MATLAB對(duì)圖2(c)進(jìn)行入射分量自適應(yīng)灰度拉伸,結(jié)果如圖3所示。

    圖3 圖2(c)入射分量自適應(yīng)灰度拉伸結(jié)果圖Fig.3 Result of adaptive grayscale stretching of incident component of Fig.2(c)

    3.3 反射分量非線性校正

    反射分量中包含豐富的細(xì)節(jié)信息,但在其輸出時(shí),總體的亮度偏暗,對(duì)比度較低。為了調(diào)節(jié)反射分量的對(duì)比度,增強(qiáng)人眼視覺效果,大多數(shù)學(xué)者選擇使用gamma曲線進(jìn)行校正,但該校正方法容易產(chǎn)生過增強(qiáng)的現(xiàn)象。為了解決這一問題,本文使用以下非線性公式對(duì)其進(jìn)行校正,具體公式如下:

    (22)

    (23)

    式中,r為校正后的反射分量,RN為歸一化后的反射分量,R為原始反射分量,Rmin和Rmax分別為R中像素最小值和最大值,t為調(diào)整參數(shù)。

    在本次實(shí)驗(yàn)中,當(dāng)t取值為0.1,gamma曲線參數(shù)取值為0.3。如圖4所示,本文非線性校正曲線與gamma曲線校正相比,gamma曲線在低像素區(qū)域會(huì)出現(xiàn)拉伸過大的情況,進(jìn)而導(dǎo)致細(xì)節(jié)信息出現(xiàn)過增強(qiáng)。本文校正曲線相對(duì)比較緩和,既可以很好地對(duì)圖像中的細(xì)節(jié)信息進(jìn)行圖像增強(qiáng),又可以避免圖像失真。

    圖4 本文曲線與gamma曲線對(duì)比圖Fig.4 Comparison between the curve in this paper and gamma curve

    4 實(shí)驗(yàn)結(jié)果與分析

    為了驗(yàn)證本文算法,選取4個(gè)紅外圖像場(chǎng)景進(jìn)行對(duì)比實(shí)驗(yàn)(紅外圖像取自O(shè)TCBVS數(shù)據(jù)集)。首先,分別對(duì)這3個(gè)場(chǎng)景(如圖5所示)使用直方圖均衡化算法(HE)、單尺度Retinex算法(SSR)、多尺度Retinex算法(MSR)、基于引導(dǎo)濾波器的多尺度Retinex算法(未改進(jìn)引導(dǎo)濾波+MSR)、文獻(xiàn)[10]、本文算法來進(jìn)行紅外圖像的增強(qiáng);然后,通過對(duì)增強(qiáng)后的紅外圖像進(jìn)行客觀和主觀評(píng)價(jià)來驗(yàn)證本文算法的優(yōu)越性和可行性。本文所有算法均使用Matlab R2016b編譯環(huán)境進(jìn)行編程實(shí)現(xiàn)(DELL系列電腦,操作系統(tǒng)為Windows7,CPU為3.1 GHz,內(nèi)存為4 GB)。本文算法的改進(jìn)引導(dǎo)濾波器選用3個(gè)尺度,其相應(yīng)窗口大小分別為20、70、115;SSR選擇的參數(shù)σ=80;MSR選取3個(gè)尺度,每個(gè)尺度的高斯濾波參數(shù)為σ1=20,σ2=80,σ3=110,且本文中所有的多尺度MSR均采用相同的3個(gè)尺度。

    圖5 用于測(cè)試的4個(gè)不同場(chǎng)景的紅外圖像Fig.5 Infrared images of 4 different scenes for testing

    4.1 主觀評(píng)價(jià)

    本文對(duì)選取的3幅場(chǎng)景進(jìn)行不同算法比較,結(jié)果如圖6~8所示。

    由圖6~8可以看出,HE算法對(duì)紅外圖像的對(duì)比度、亮度、邊緣都有一定的增強(qiáng)。但圖6、7中后面的背景細(xì)節(jié)信息嚴(yán)重丟失,圖8中的船和人細(xì)節(jié)信息丟失。同時(shí),圖6~8中的噪聲也很嚴(yán)重。SSR算法對(duì)紅外圖像中的亮度的到明顯增強(qiáng),但其增強(qiáng)圖像的對(duì)比度較差,且有些地方產(chǎn)生了“光暈”現(xiàn)象,人眼視覺效果較差。圖6中的左邊飛機(jī)場(chǎng)窗口出現(xiàn)“光暈”。MSR算法增強(qiáng)的紅外圖像亮度有一定提高,“光暈”現(xiàn)象有所改善,但其對(duì)比度較低,部分背景信息丟失,紅外圖像整體視覺較差。圖6和圖7背景模糊,部分細(xì)節(jié)信息丟失。未改進(jìn)的引導(dǎo)濾波器+MSR算法增強(qiáng)的紅外圖像邊緣清晰度高,亮度有所增強(qiáng),但其整體圖像效果比較模糊,如圖6和圖8所示。文獻(xiàn)[10]所得圖像效果相對(duì)前面幾種算法亮度有一定提升,有效改善了圖像中的“光暈”現(xiàn)象,但其主要細(xì)節(jié)信息有所丟失。本文算法在有效去除“光暈”現(xiàn)象的同時(shí),很好地提升了圖像的對(duì)比度、細(xì)節(jié)和清晰度等,人眼整體視覺效果較好。

    圖6 場(chǎng)景1的各個(gè)算法比較實(shí)驗(yàn)結(jié)果Fig.6 Comparison experimental results of various algorithms in scene 1.

    圖7 場(chǎng)景2的各個(gè)算法比較實(shí)驗(yàn)結(jié)果Fig.7 Comparison experimental results of various algorithms in scene 2.

    圖8 場(chǎng)景3的各個(gè)算法比較實(shí)驗(yàn)結(jié)果Fig.8 Comparison experimental results of various algorithms in scene 3.

    4.2 客觀評(píng)價(jià)

    為了進(jìn)一步驗(yàn)證本文算法的有效性,本文選用圖像標(biāo)準(zhǔn)差(SD)[12]、平均梯度(AVG)[13]、峰值信噪比(PSNR)[14]、信息熵(IE)[15]客觀指標(biāo)對(duì)增強(qiáng)的紅外圖像進(jìn)行定量評(píng)價(jià)。

    對(duì)以上3個(gè)場(chǎng)景分別進(jìn)行客觀指標(biāo)評(píng)價(jià),結(jié)果如表1~3所示。

    表1 場(chǎng)景1不同算法所得紅外增強(qiáng)圖像的客觀評(píng)價(jià)指標(biāo)Tab.1 Objective evaluation indicators of infrared enhanced images obtained by different algorithms in scene 1

    表2 場(chǎng)景2不同算法所得紅外增強(qiáng)圖像的客觀評(píng)價(jià)指標(biāo)Tab.2 Objective evaluation indicators of infrared enhanced images obtained by different algorithms in scene 2

    表3 場(chǎng)景3不同算法所得紅外增強(qiáng)圖像的客觀評(píng)價(jià)指標(biāo)Tab.3 Objective evaluation indicators of infrared enhanced images obtained by different algorithms in scene 3

    由表1~3可以看出:本文算法并不能同時(shí)兼顧SD、AVG、PSNR、IE,但本文算法相對(duì)其他算法相對(duì)占優(yōu)。HE算法所得紅外增強(qiáng)圖像的標(biāo)準(zhǔn)差較高,但其得到的增強(qiáng)圖像細(xì)節(jié)信息丟失,邊緣模糊。本文算法處理得到的紅外圖像與SSR、MSR、未改進(jìn)引導(dǎo)濾波器+MSR算法、文獻(xiàn)[10]相比平均梯度(AVG)、信息熵(IE)較高,表明本文算法增強(qiáng)后的紅外圖像有更加豐富的信息量,圖像更清晰。同時(shí)本文算法處理后的圖像峰值信噪比(PSNR)較高,表明增強(qiáng)后的圖像去除噪聲更徹底,整體處理效果更好,對(duì)比度更強(qiáng),人眼整體視覺效果更好。綜上可以看出本文算法紅外圖像增強(qiáng)效果較好,且具有可行性和優(yōu)化性。

    5 結(jié) 論

    本文提出一種基于改進(jìn)多尺度Retinex的紅外圖像細(xì)節(jié)增強(qiáng)算法。該算法實(shí)驗(yàn)結(jié)果表明:利用改進(jìn)引導(dǎo)濾波器代替高斯濾波器對(duì)入射分量進(jìn)行精確估計(jì),解決了原始算法放大噪聲及產(chǎn)生“光暈”等問題;對(duì)反射分量進(jìn)行校正時(shí),使用一種平緩的非線性曲線代替常用的gamma曲線進(jìn)行校正,解決了圖像過增強(qiáng)問題,進(jìn)一步提高了反射圖像細(xì)節(jié)信息。同時(shí),本文算法與HE、SSR、MSR、引導(dǎo)濾波+MSR算法相比,對(duì)紅外圖像的對(duì)比度、細(xì)節(jié)信息、人眼視覺效果都有很好的增強(qiáng)效果。

    猜你喜歡
    光暈濾波器尺度
    基于無擾濾波器和AED-ADT的無擾切換控制
    盧美涵:夕陽
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    開關(guān)電源EMI濾波器的應(yīng)用方法探討
    電子制作(2018年16期)2018-09-26 03:26:50
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    基于TMS320C6678的SAR方位向預(yù)濾波器的并行實(shí)現(xiàn)
    微軟大作《光暈》首現(xiàn)iOS射擊視角變了
    9
    室外雕塑的尺度
    雕塑(2000年2期)2000-06-22 16:13:30
    伊人久久大香线蕉亚洲五| 亚洲精品久久久久久婷婷小说| 国产麻豆69| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品第二区| 日韩 欧美 亚洲 中文字幕| 久久久久久久久久久免费av| 精品一区二区三区四区五区乱码 | 午夜福利视频精品| 久久人人97超碰香蕉20202| 欧美黑人欧美精品刺激| 婷婷色综合大香蕉| 啦啦啦视频在线资源免费观看| 亚洲三区欧美一区| 亚洲视频免费观看视频| 亚洲精品,欧美精品| 国产精品欧美亚洲77777| 五月天丁香电影| av在线app专区| 看十八女毛片水多多多| 精品久久久久久电影网| av电影中文网址| 熟妇人妻不卡中文字幕| 欧美最新免费一区二区三区| 视频区图区小说| 啦啦啦中文免费视频观看日本| 亚洲熟女毛片儿| 欧美黑人精品巨大| 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 久久久久国产精品人妻一区二区| 免费看不卡的av| av网站免费在线观看视频| 一级片'在线观看视频| 婷婷色综合www| 精品亚洲乱码少妇综合久久| 51午夜福利影视在线观看| 日日撸夜夜添| 夜夜骑夜夜射夜夜干| 精品一区二区三区四区五区乱码 | 一本久久精品| 日本欧美视频一区| 母亲3免费完整高清在线观看| av又黄又爽大尺度在线免费看| 91精品伊人久久大香线蕉| 在线亚洲精品国产二区图片欧美| 久久久久精品人妻al黑| 亚洲欧美成人综合另类久久久| 啦啦啦在线免费观看视频4| 丝袜在线中文字幕| 午夜福利免费观看在线| 亚洲av成人精品一二三区| 黄片播放在线免费| 亚洲av综合色区一区| 黄片小视频在线播放| 日韩一本色道免费dvd| 国产激情久久老熟女| 超碰成人久久| 精品一区二区三卡| 亚洲七黄色美女视频| 亚洲成国产人片在线观看| 99九九在线精品视频| 日日啪夜夜爽| 男女下面插进去视频免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲熟妇少妇任你| 欧美黑人精品巨大| 色精品久久人妻99蜜桃| 国产人伦9x9x在线观看| 中文字幕制服av| 亚洲久久久国产精品| 久久99精品国语久久久| 久久久久久久精品精品| 老熟女久久久| 国产有黄有色有爽视频| 国产成人91sexporn| 中文字幕高清在线视频| 成年av动漫网址| 亚洲中文av在线| 亚洲美女视频黄频| 国产精品熟女久久久久浪| 国产一区有黄有色的免费视频| 亚洲第一区二区三区不卡| 色综合欧美亚洲国产小说| 男人爽女人下面视频在线观看| 精品亚洲成a人片在线观看| 99精国产麻豆久久婷婷| 2018国产大陆天天弄谢| 韩国高清视频一区二区三区| 亚洲国产日韩一区二区| 一级毛片 在线播放| 亚洲自偷自拍图片 自拍| 日本欧美视频一区| 日韩人妻精品一区2区三区| 国产成人啪精品午夜网站| 国产 精品1| 1024香蕉在线观看| 99热国产这里只有精品6| 卡戴珊不雅视频在线播放| 久久精品国产a三级三级三级| 天天添夜夜摸| 久久久精品94久久精品| 国产一级毛片在线| 你懂的网址亚洲精品在线观看| 狂野欧美激情性xxxx| 在现免费观看毛片| 久久av网站| 观看av在线不卡| 久久韩国三级中文字幕| 日本欧美视频一区| 欧美日韩视频精品一区| 国产免费视频播放在线视频| 国产免费又黄又爽又色| 中文字幕人妻丝袜制服| 一区二区av电影网| 美国免费a级毛片| 久久人人爽人人片av| 国产av一区二区精品久久| 91国产中文字幕| 久久精品久久精品一区二区三区| 日日啪夜夜爽| 色精品久久人妻99蜜桃| 久久97久久精品| 国产爽快片一区二区三区| 欧美最新免费一区二区三区| www.精华液| 国产亚洲av高清不卡| 亚洲人成网站在线观看播放| 99精品久久久久人妻精品| av国产久精品久网站免费入址| 国产精品麻豆人妻色哟哟久久| 成人毛片60女人毛片免费| 久久久久精品国产欧美久久久 | 久久热在线av| 女人高潮潮喷娇喘18禁视频| 1024视频免费在线观看| 丝袜美足系列| 欧美黑人精品巨大| 日韩不卡一区二区三区视频在线| 日本欧美国产在线视频| 人人妻人人爽人人添夜夜欢视频| 日日撸夜夜添| 久久精品久久久久久久性| 大香蕉久久网| 午夜激情久久久久久久| 999精品在线视频| 国产高清不卡午夜福利| av天堂久久9| 飞空精品影院首页| 交换朋友夫妻互换小说| 黄色怎么调成土黄色| 高清在线视频一区二区三区| 1024视频免费在线观看| 亚洲欧美一区二区三区久久| 两个人免费观看高清视频| 亚洲美女黄色视频免费看| av国产精品久久久久影院| www.自偷自拍.com| 午夜福利网站1000一区二区三区| 久久影院123| 国产爽快片一区二区三区| 成年美女黄网站色视频大全免费| 国产精品 欧美亚洲| 热re99久久精品国产66热6| 亚洲成色77777| xxx大片免费视频| 建设人人有责人人尽责人人享有的| 91精品三级在线观看| 国产精品秋霞免费鲁丝片| 成年人午夜在线观看视频| 久久婷婷青草| 国产精品秋霞免费鲁丝片| a级毛片黄视频| 亚洲久久久国产精品| 女人精品久久久久毛片| 在线观看www视频免费| 免费高清在线观看日韩| 黑人猛操日本美女一级片| 啦啦啦在线观看免费高清www| 亚洲欧美精品综合一区二区三区| 日韩熟女老妇一区二区性免费视频| 777米奇影视久久| a级毛片黄视频| 蜜桃在线观看..| 日韩精品有码人妻一区| 午夜激情av网站| 国产毛片在线视频| 欧美日韩一级在线毛片| 777久久人妻少妇嫩草av网站| 一本一本久久a久久精品综合妖精| 亚洲美女视频黄频| 亚洲精品在线美女| 成人毛片60女人毛片免费| 亚洲精品久久午夜乱码| 一区二区日韩欧美中文字幕| 成人国产麻豆网| av女优亚洲男人天堂| 精品人妻熟女毛片av久久网站| 七月丁香在线播放| 999精品在线视频| 久久人妻熟女aⅴ| 久久人人爽av亚洲精品天堂| 精品一区二区三卡| bbb黄色大片| 下体分泌物呈黄色| 国产av精品麻豆| 欧美精品av麻豆av| 在线观看一区二区三区激情| 在线观看免费午夜福利视频| 亚洲精品成人av观看孕妇| 日韩熟女老妇一区二区性免费视频| 亚洲国产中文字幕在线视频| 人人妻,人人澡人人爽秒播 | 亚洲美女视频黄频| 亚洲欧美日韩另类电影网站| 免费女性裸体啪啪无遮挡网站| 亚洲欧洲精品一区二区精品久久久 | 国产一区二区激情短视频 | 最新在线观看一区二区三区 | 国产乱人偷精品视频| 亚洲国产欧美一区二区综合| 黄色毛片三级朝国网站| 80岁老熟妇乱子伦牲交| 日韩大片免费观看网站| 成年av动漫网址| 久久久亚洲精品成人影院| 国产免费视频播放在线视频| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 老汉色∧v一级毛片| 香蕉丝袜av| 国产高清国产精品国产三级| 日韩 欧美 亚洲 中文字幕| 看非洲黑人一级黄片| 777久久人妻少妇嫩草av网站| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 欧美在线黄色| 日韩制服骚丝袜av| 欧美日韩亚洲高清精品| 久久99精品国语久久久| 国产精品一区二区精品视频观看| 亚洲欧美精品自产自拍| 色播在线永久视频| 国产探花极品一区二区| av天堂久久9| 国产成人系列免费观看| 日本wwww免费看| 青春草国产在线视频| 免费高清在线观看日韩| 一二三四中文在线观看免费高清| 国产精品国产av在线观看| 亚洲av男天堂| 中文字幕亚洲精品专区| 一区二区三区精品91| 超碰成人久久| 国产精品av久久久久免费| 精品一区二区三区av网在线观看 | 国产麻豆69| 日韩中文字幕欧美一区二区 | 日韩中文字幕视频在线看片| 99久久99久久久精品蜜桃| 九九爱精品视频在线观看| 精品国产乱码久久久久久小说| 少妇人妻 视频| 熟女少妇亚洲综合色aaa.| 99久久99久久久精品蜜桃| 男人舔女人的私密视频| 老汉色av国产亚洲站长工具| 99精品久久久久人妻精品| 午夜日韩欧美国产| 新久久久久国产一级毛片| 成年美女黄网站色视频大全免费| 美女扒开内裤让男人捅视频| 狂野欧美激情性bbbbbb| 国产精品久久久人人做人人爽| 一区二区三区精品91| 久久精品久久精品一区二区三区| 九九爱精品视频在线观看| 90打野战视频偷拍视频| 日韩中文字幕视频在线看片| 9热在线视频观看99| 国产毛片在线视频| 极品人妻少妇av视频| 啦啦啦中文免费视频观看日本| 不卡av一区二区三区| 亚洲精品国产av蜜桃| 在线观看免费视频网站a站| 亚洲精品视频女| 精品第一国产精品| 如何舔出高潮| 99精国产麻豆久久婷婷| 女人精品久久久久毛片| 老司机在亚洲福利影院| 亚洲,欧美,日韩| 少妇人妻久久综合中文| 五月天丁香电影| 免费在线观看完整版高清| 日韩不卡一区二区三区视频在线| 99久国产av精品国产电影| 免费看不卡的av| 精品国产乱码久久久久久小说| 毛片一级片免费看久久久久| 看免费成人av毛片| 日韩成人av中文字幕在线观看| 热99久久久久精品小说推荐| 激情五月婷婷亚洲| 巨乳人妻的诱惑在线观看| 免费不卡黄色视频| 亚洲美女搞黄在线观看| 国产欧美日韩综合在线一区二区| 欧美 日韩 精品 国产| 超色免费av| 9热在线视频观看99| 欧美xxⅹ黑人| 欧美日韩国产mv在线观看视频| 亚洲国产日韩一区二区| 亚洲精品av麻豆狂野| 国产毛片在线视频| 欧美日韩视频高清一区二区三区二| 国产精品国产av在线观看| 2018国产大陆天天弄谢| 午夜福利乱码中文字幕| 国产成人啪精品午夜网站| 免费少妇av软件| 国产成人91sexporn| 国产爽快片一区二区三区| 波多野结衣一区麻豆| 各种免费的搞黄视频| 9热在线视频观看99| 国产av精品麻豆| 午夜日韩欧美国产| 色网站视频免费| 丁香六月欧美| 人人妻人人澡人人爽人人夜夜| 亚洲情色 制服丝袜| 国产又色又爽无遮挡免| 日本猛色少妇xxxxx猛交久久| 国产又爽黄色视频| 午夜激情久久久久久久| 午夜福利乱码中文字幕| 男女国产视频网站| 亚洲欧美成人精品一区二区| 成人漫画全彩无遮挡| 国产不卡av网站在线观看| 国产精品欧美亚洲77777| 日韩精品免费视频一区二区三区| 久久久国产一区二区| 精品亚洲成a人片在线观看| 亚洲国产av影院在线观看| 女的被弄到高潮叫床怎么办| 日韩伦理黄色片| 国产精品一二三区在线看| 最新的欧美精品一区二区| 水蜜桃什么品种好| 电影成人av| 亚洲欧美中文字幕日韩二区| 一级片免费观看大全| 久久久久精品性色| 亚洲激情五月婷婷啪啪| 青春草国产在线视频| 十八禁人妻一区二区| 一边摸一边做爽爽视频免费| 中文欧美无线码| 热re99久久精品国产66热6| 街头女战士在线观看网站| 黄网站色视频无遮挡免费观看| 亚洲精品视频女| 久热爱精品视频在线9| 久久天躁狠狠躁夜夜2o2o | 亚洲美女搞黄在线观看| 99久久精品国产亚洲精品| www.自偷自拍.com| 色吧在线观看| 国产精品久久久av美女十八| 啦啦啦在线观看免费高清www| 美女午夜性视频免费| 国产亚洲欧美精品永久| 午夜免费观看性视频| 欧美国产精品一级二级三级| 宅男免费午夜| av在线播放精品| 美女高潮到喷水免费观看| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频 | 精品一区二区三区av网在线观看 | 亚洲av电影在线观看一区二区三区| 国产又色又爽无遮挡免| 人人妻人人添人人爽欧美一区卜| www.av在线官网国产| 亚洲 欧美一区二区三区| 精品人妻熟女毛片av久久网站| 亚洲 欧美一区二区三区| 欧美激情极品国产一区二区三区| 99九九在线精品视频| 国产熟女欧美一区二区| 成年人午夜在线观看视频| 美女大奶头黄色视频| 老司机在亚洲福利影院| 大香蕉久久成人网| 亚洲国产成人一精品久久久| av片东京热男人的天堂| 国产淫语在线视频| 伦理电影免费视频| 亚洲在久久综合| av在线播放精品| 亚洲成色77777| 午夜福利影视在线免费观看| 天天操日日干夜夜撸| 国产黄色视频一区二区在线观看| 国产精品偷伦视频观看了| 黄色 视频免费看| 日韩不卡一区二区三区视频在线| 波多野结衣av一区二区av| 国产精品人妻久久久影院| 久久婷婷青草| 黄片播放在线免费| 看十八女毛片水多多多| av不卡在线播放| 满18在线观看网站| 麻豆乱淫一区二区| 99香蕉大伊视频| 亚洲av成人不卡在线观看播放网 | 日本午夜av视频| www.自偷自拍.com| 亚洲av福利一区| 精品人妻在线不人妻| 久久久久人妻精品一区果冻| 日韩,欧美,国产一区二区三区| 精品少妇内射三级| 亚洲欧美一区二区三区黑人| 九色亚洲精品在线播放| 亚洲精品aⅴ在线观看| 国产免费一区二区三区四区乱码| 欧美日韩综合久久久久久| 国产精品熟女久久久久浪| 免费高清在线观看日韩| 久久国产亚洲av麻豆专区| 看十八女毛片水多多多| www日本在线高清视频| 欧美精品人与动牲交sv欧美| 亚洲成人一二三区av| 91国产中文字幕| 成年av动漫网址| 一级a爱视频在线免费观看| 天天操日日干夜夜撸| 国产乱来视频区| 婷婷色综合大香蕉| 波多野结衣av一区二区av| 国产成人啪精品午夜网站| 国产毛片在线视频| 国产免费视频播放在线视频| 亚洲欧美一区二区三区黑人| 亚洲熟女毛片儿| 免费在线观看黄色视频的| 丝袜喷水一区| 亚洲国产中文字幕在线视频| 十八禁网站网址无遮挡| 国产在视频线精品| 亚洲精品日韩在线中文字幕| 在线亚洲精品国产二区图片欧美| 爱豆传媒免费全集在线观看| 亚洲成人国产一区在线观看 | 亚洲图色成人| 久久久久国产一级毛片高清牌| 久久久欧美国产精品| 我的亚洲天堂| 亚洲成人av在线免费| 亚洲国产精品成人久久小说| 香蕉国产在线看| 亚洲国产精品国产精品| 2018国产大陆天天弄谢| 中文字幕人妻丝袜一区二区 | a级毛片在线看网站| 97在线人人人人妻| 色综合欧美亚洲国产小说| 国产不卡av网站在线观看| 啦啦啦视频在线资源免费观看| 99精国产麻豆久久婷婷| 久久久久国产精品人妻一区二区| 国产欧美日韩综合在线一区二区| 少妇精品久久久久久久| 亚洲av电影在线观看一区二区三区| 女人久久www免费人成看片| 热99国产精品久久久久久7| 亚洲精品自拍成人| 亚洲专区中文字幕在线 | 久久av网站| 好男人视频免费观看在线| 亚洲精品aⅴ在线观看| 婷婷色av中文字幕| 桃花免费在线播放| 亚洲精品久久成人aⅴ小说| 欧美日韩视频高清一区二区三区二| 天天添夜夜摸| 亚洲精品成人av观看孕妇| 视频区图区小说| 久久狼人影院| 日韩成人av中文字幕在线观看| 最近中文字幕2019免费版| 日韩 亚洲 欧美在线| 一区在线观看完整版| 亚洲国产毛片av蜜桃av| av一本久久久久| 岛国毛片在线播放| av女优亚洲男人天堂| 国产熟女午夜一区二区三区| 久久国产精品男人的天堂亚洲| 两个人看的免费小视频| 黑丝袜美女国产一区| 国产黄色视频一区二区在线观看| 人人妻人人澡人人看| 国产精品香港三级国产av潘金莲 | 女的被弄到高潮叫床怎么办| 777米奇影视久久| 男女床上黄色一级片免费看| 高清欧美精品videossex| 亚洲人成电影观看| videos熟女内射| 看免费av毛片| 狠狠婷婷综合久久久久久88av| 精品少妇一区二区三区视频日本电影 | 黄片播放在线免费| 免费观看av网站的网址| 国产精品久久久久久精品古装| e午夜精品久久久久久久| 考比视频在线观看| 好男人视频免费观看在线| 日韩av在线免费看完整版不卡| 男人舔女人的私密视频| 黑人猛操日本美女一级片| 成人亚洲精品一区在线观看| 一级a爱视频在线免费观看| 午夜老司机福利片| 赤兔流量卡办理| 王馨瑶露胸无遮挡在线观看| 久久久亚洲精品成人影院| 日韩制服丝袜自拍偷拍| 777米奇影视久久| 欧美日本中文国产一区发布| 亚洲国产精品成人久久小说| 国产极品粉嫩免费观看在线| 午夜福利一区二区在线看| 大片免费播放器 马上看| 男女边摸边吃奶| 亚洲美女搞黄在线观看| 国产视频首页在线观看| 叶爱在线成人免费视频播放| 久久久久久免费高清国产稀缺| 可以免费在线观看a视频的电影网站 | 午夜激情久久久久久久| 不卡视频在线观看欧美| 亚洲精品国产av成人精品| 亚洲av国产av综合av卡| 国产野战对白在线观看| 我的亚洲天堂| 亚洲精品aⅴ在线观看| 18禁观看日本| 男女国产视频网站| 叶爱在线成人免费视频播放| 久久久久久免费高清国产稀缺| 大码成人一级视频| www.av在线官网国产| 国产精品一区二区精品视频观看| 精品一区二区免费观看| 亚洲男人天堂网一区| 国产精品嫩草影院av在线观看| 最黄视频免费看| 日韩欧美一区视频在线观看| netflix在线观看网站| 国产精品成人在线| 免费少妇av软件| 黄色毛片三级朝国网站| 亚洲av国产av综合av卡| 一个人免费看片子| 亚洲色图综合在线观看| 伊人亚洲综合成人网| 自线自在国产av| 久久午夜综合久久蜜桃| 制服丝袜香蕉在线| 成年av动漫网址| 久久久久人妻精品一区果冻| 欧美xxⅹ黑人| 国产爽快片一区二区三区| 久久ye,这里只有精品| 免费观看人在逋| 中文字幕亚洲精品专区| 亚洲精品国产色婷婷电影| 中文字幕人妻丝袜制服| 在线亚洲精品国产二区图片欧美| 国产片内射在线| 男女之事视频高清在线观看 | 精品一区二区三卡| 嫩草影院入口| 亚洲美女黄色视频免费看| 久久狼人影院| av国产久精品久网站免费入址| 两性夫妻黄色片| 亚洲精品久久成人aⅴ小说| 久久精品国产综合久久久| 国产精品无大码| 久久精品aⅴ一区二区三区四区| 亚洲国产最新在线播放| 日本午夜av视频| 午夜免费男女啪啪视频观看| 2021少妇久久久久久久久久久| 精品国产乱码久久久久久小说| 国产一卡二卡三卡精品 | 夫妻性生交免费视频一级片| 免费看av在线观看网站| 国产精品 欧美亚洲| 狠狠婷婷综合久久久久久88av| 999久久久国产精品视频| 久久影院123| 天天操日日干夜夜撸|