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

    基于薄云厚度分布評估的遙感影像高保真薄云去除方法

    2019-01-07 12:16:40汪月云
    計(jì)算機(jī)應(yīng)用 2018年12期
    關(guān)鍵詞:薄云波段藍(lán)色

    汪月云,黃 微,王 睿

    (上海大學(xué) 通信與信息工程學(xué)院,上海 200444)(*通信作者電子郵箱yueyunW@163.com)

    0 引言

    光學(xué)遙感影像往往受到薄云污染的影響,造成影像地物細(xì)節(jié)對比度下降,這給影像的應(yīng)用[1-3]帶來了很多困難。因此,如何有效去除降質(zhì)影像的薄云影響從而提高影像的利用率是遙感影像應(yīng)用的一項(xiàng)重要研究。

    目前,光學(xué)遙感影像薄云去除主要有頻域法和空域法兩類。頻域法去云主要根據(jù)薄云信息處于影像低頻區(qū)域的原理,所以可以設(shè)計(jì)合適的高通濾波器以濾除低頻成分從而有效去云。比如:閻慶等[4]結(jié)合小波變換和同態(tài)濾波、馮春等[5]利用同態(tài)濾波都達(dá)到了有效去云效果。

    空域法去云將薄云污染分為乘性污染和加性污染兩種模型。乘性污染經(jīng)典的去云方法有暗通道先驗(yàn)法[6-7]和薄云優(yōu)化變換(Haze Optimized Transformation, HOT)法[8-9]。暗通道先驗(yàn)法通過對清晰影像庫統(tǒng)計(jì)得到暗通道分布規(guī)律,即清晰影像至少在某一個(gè)顏色通道中有灰度值接近于0的像素,利用該先驗(yàn)規(guī)律得出云模型的大氣光和透射率從而進(jìn)行薄云去除。Pan等[10]研究發(fā)現(xiàn)無云遙感影像的暗通道數(shù)值很低,但不接近于0,通過改進(jìn)遙感影像暗通道分布從而有效去除薄云。HOT法是利用遙感影像紅光波段和藍(lán)光波段相關(guān)性高的特點(diǎn)尋找合適的晴空線,用像元到晴空線的距離表征該點(diǎn)云量進(jìn)而有效去云。He等[11]提出了AHOT(Advanced HOT)算法,先利用空間信息剔除易混淆地物,提高云檢測精度,再考慮大氣散射乘性效應(yīng),提出基于AHOT的虛擬云點(diǎn)(Virtual Cloud Point, VCP)技術(shù)擴(kuò)大由薄云污染降低的數(shù)字量化值(Digital Number, DN)方差從而有效去云。Makarau等[12]利用加性薄云污染模型,在傳統(tǒng)的暗目標(biāo)減法(Dark Object Subtraction, DOS)基礎(chǔ)上[13],通過局部搜索整個(gè)影像的暗目標(biāo)構(gòu)建影像的薄云厚度分布(Haze Thickness Map, HTM),最后用原降質(zhì)影像減去薄云厚度分布得到無云影像。這些去云算法大多關(guān)注去云后影像細(xì)節(jié)的銳化、地物間對比度的增強(qiáng),如何在去云的同時(shí)保證數(shù)據(jù)保真是一項(xiàng)具有挑戰(zhàn)性的問題。

    近年來,已有學(xué)者提出了一些考慮保真的薄云去除算法[14-16]。Shen等[14]基于同態(tài)濾波法半自動化地確定不同波段的截止頻率,對薄云像素和無云像素分別處理從而獲得高保真的無云影像。Singh等[15]基于暗通道先驗(yàn)法提出改進(jìn)的無云影像恢復(fù)算法,通過改進(jìn)聯(lián)合三邊濾波器以提高大氣光的精細(xì)估計(jì)效果,并重新定義透射率從而降低無云影像的顏色失真程度。Liu等[16]基于加性薄云污染模型,將原始薄云厚度分布先經(jīng)過紋理去除,再進(jìn)行大面積地物反射抑制,然后用導(dǎo)向?yàn)V波器進(jìn)行平滑,最后用線性回歸求取不同波段的薄云厚度從而得到無云影像。文獻(xiàn)[16]算法的去云結(jié)果整體具有保真效果,但是依舊存在無云區(qū)域過度校正,以及影像中藍(lán)色地物色彩嚴(yán)重失真的問題。

    因此,本文在文獻(xiàn)[16]算法的基礎(chǔ)上提出了一種改進(jìn)的薄云去除方法,即基于薄云厚度分布(HTM)評估的高保真薄云去除方法。該方法首先利用基于暗通道先驗(yàn)薄云去除算法求得影像的透射率分布圖[10],根據(jù)透射率分布圖中像素?cái)?shù)值接近于1來自動確定無云區(qū)域的范圍。然后,在HTM中,求取無云區(qū)域的平均值,將整個(gè)HTM減去無云區(qū)域平均值,使得HTM滿足無云區(qū)域的薄云厚度接近于零,以確保無云區(qū)域的保真度。針對藍(lán)色地物區(qū)域云量被過估計(jì),造成去云結(jié)果中色彩失真的問題,本文利用藍(lán)色地物RGB波段的光譜特性,結(jié)合閾值法檢測出藍(lán)色地物,用藍(lán)色地物附近的HTM平均值作為藍(lán)色地物區(qū)域的云量估計(jì)。最后,用降質(zhì)影像減去不同波段優(yōu)化的HTM得到最終去云結(jié)果。

    通過多幅影像的實(shí)驗(yàn)發(fā)現(xiàn),與文獻(xiàn)[14]算法去云結(jié)果、文獻(xiàn)[15]算法去云結(jié)果,以及文獻(xiàn)[16]算法去云結(jié)果相比,本文方法的結(jié)果影像去云效果良好,同時(shí)具有數(shù)據(jù)高保真的優(yōu)點(diǎn)。

    1 加性模型薄云去除

    1.1 云成像模型

    薄云污染可以認(rèn)為是加性污染[12,16],并且薄云污染對不同波段的影響不同[16-17],所以云成像模型可以描述為:

    Ic(x)=Jc(x)+Hc(x);c∈{R,G,B}

    (1)

    其中:x代表像素的位置;I代表降質(zhì)影像;J代表無云影像;H代表薄云厚度分布;c是{R,G,B}中的任一波段。

    1.2 薄云估計(jì)

    在Makarau等[12]的研究中,薄云厚度分布(HTM)估計(jì)如下:

    (2)

    其中:Ib(x)、Ig(x)分別是降質(zhì)影像的藍(lán)光波段和綠光波段;W(x)是以x為中心、大小為w×w像素的窗口。本文實(shí)驗(yàn)中該窗口大小取3×3像素。

    Liu等[16]發(fā)現(xiàn)由式(2)得到的Hraw(x)存在將地物反射光作為薄云組成成分,造成降質(zhì)影像被過度矯正的問題。為了減少地物反射光的影響,文獻(xiàn)[16]對Hraw(x)進(jìn)一步改進(jìn):首先為了去除Hraw(x)中的地物紋理,文獻(xiàn)[16]采用了可以去除無用紋理和細(xì)節(jié)的全變差(Total Variation, TV)正則化圖像修復(fù)模型[18-19]得到平滑的薄云厚度分布Hs(x);接著,為解決大面積高亮地物的薄云厚度過估計(jì)問題,文獻(xiàn)[16]用Cannny算子檢測出高亮地物邊緣,用局部窗口對Hs(x)中高亮地物區(qū)域作數(shù)值抑制得到Hd(x);然后,應(yīng)用導(dǎo)向?yàn)V波器[20]去除Hd(x)中的塊狀效應(yīng)得到Hg(x);最后,求取不同波段的薄云厚度分布。

    Himpr(x)=Hg(x)-minHg

    (3)

    Hc(x)=kcHimpr(x);c∈{R,G,B}

    (4)

    為了去除氣溶膠的影響,式(3)用Hg減去其最小值minHg得到最終的薄云厚度估計(jì)Himpr(x)。式(4)中:kc是不同波段的云估計(jì)強(qiáng)度,它是在Hg低于平均值的像素集與降質(zhì)影像I的不同波段中低于平均值的像素集的交集中,以Hg的像素為自變量,以影像的不同波段的像素為因變量,通過線性回歸得到;Hc(x)是不同波段的薄云厚度估計(jì)。得到Hc(x)后,通過式(1)即可得到無云影像J。

    圖1是根據(jù)文獻(xiàn)[16]算法得到的實(shí)驗(yàn)結(jié)果。其中:圖1(a)是大小為300×400像素的GF-2降質(zhì)影像;圖1(b)是最終的云厚度分布估計(jì)Himpr(x);圖1(c)是去云結(jié)果。由圖1可以看出,Himpr(x)中盡管云厚度分布范圍估計(jì)比較準(zhǔn)確,但是無云區(qū)域的云厚度值并不為0或接近于0,同時(shí)無云區(qū)域的藍(lán)色和白色物體估計(jì)的云厚度分布值也很高,這就造成了去云影像的無云區(qū)域存在過度矯正和藍(lán)色地物區(qū)域嚴(yán)重失真的問題。

    圖1 文獻(xiàn)[16]算法的薄云估計(jì)和去云結(jié)果Fig.1 Haze estimation and haze removal results of algorithm in literature [16]

    2 薄云去除改進(jìn)

    為解決圖1去云結(jié)果中反映出來的無云區(qū)域過度矯正和藍(lán)色地物失真嚴(yán)重問題,基于加性薄云污染模型,本文在文獻(xiàn)[16]得到的最終薄云厚度分布Himpr(x)基礎(chǔ)上進(jìn)行改進(jìn),在獲得良好去云結(jié)果的同時(shí)確保影像的高保真。本文方法的流程如圖2所示。

    圖2 本文方法流程Fig. 2 Flow chart of the proposed method

    2.1 無云區(qū)域保真

    首先針對無云區(qū)域薄云厚度被過估計(jì),造成無云區(qū)域數(shù)據(jù)失真的問題,本文對Himpr(x)作了如下改進(jìn):

    Hnew(x)=max(Himpr(x)-K,0)

    (5)

    其中:

    (6)

    式中:Dclean是Himpr(x)中的無云區(qū)域;M是Dclean中的像素個(gè)數(shù);K是無云區(qū)域的平均值。理論上,Hnew(x)的無云區(qū)域DN值接近或等于0,有云區(qū)域DN值大于0,為了避免由于Himpr(x)-K導(dǎo)致負(fù)數(shù)的出現(xiàn),這里使用了最大值算子max。

    無云區(qū)域Dclean的確定影響著最終無云影像無云區(qū)域的保真度。為了有效確定無云區(qū)域Dclean的范圍,本文利用基于暗通道先驗(yàn)薄云去除算法[6-7]中的透射率來估計(jì)該區(qū)域。其中透射率為:

    t(x)∈(0,1]

    球形工頻電場檢測傳感器采用球形電容檢測感應(yīng)電壓的方法來檢測電場,一維球形傳感器的結(jié)構(gòu)如圖1所示,將一個(gè)空心金屬球殼切割成兩部分,通過絕緣物質(zhì)將兩個(gè)半球連接在一起,并在兩個(gè)半球中間并聯(lián)一個(gè)測量電容,兩個(gè)半球分別構(gòu)成了傳感器的兩個(gè)電極[6-7]。

    (7)

    2.2 藍(lán)色地物云量估計(jì)

    對于降質(zhì)影像中的藍(lán)色地物,由于像素在藍(lán)光波段的亮度值遠(yuǎn)大于綠光波段,因此根據(jù)式(1)得到的薄云厚度分布Hraw(x)中容易造成藍(lán)色地物薄云厚度被過估計(jì),最終導(dǎo)致去云結(jié)果中藍(lán)色地物色彩失真。

    通過大量實(shí)驗(yàn)數(shù)據(jù)發(fā)現(xiàn),在優(yōu)化后的薄云厚度分布Hnew(x)中,白色地物和藍(lán)色地物區(qū)域的DN值都很高。因此,本文先利用閾值選定白色區(qū)域和藍(lán)色區(qū)域,再根據(jù)原始影像中藍(lán)色地物藍(lán)光波段像素值遠(yuǎn)大于紅光波段像素值的特點(diǎn)[21],剔除白色地物區(qū)域最終得到藍(lán)色地物的像素集:

    P1={x|Hnew(x)≥T}

    (8)

    P2={x|Ib(x)-Ir(x)>0}

    (9)

    Pf=P1∩P2

    (10)

    其中:x代表像素坐標(biāo)位置;P1是包含藍(lán)色地物和白色地物的像素集;P2是包含偏藍(lán)地物的像素集;Pf是藍(lán)色地物像素集。通過閾值T可以選定白色地物區(qū)域和藍(lán)色地物區(qū)域,閾值T的選定不能過大,否則影像中的藍(lán)色地物區(qū)域不能全部選中,閾值T也不能過小,否則容易包含非藍(lán)色地物。本文中閾值T取Hnew(x)中最大像素值的70%,即T=max(Hnew(x))*0.7。Ib(x)、Ir(x)分別是原始影像的藍(lán)光波段和紅光波段。由于薄云的分布一般是連續(xù)和平緩的,因此本文將藍(lán)色地物區(qū)域的薄云厚度用周圍區(qū)域薄云厚度的均值代替,其表達(dá)式為:

    (11)

    其中:Hblue是藍(lán)色地物薄云厚度;Hnew(x)是優(yōu)化的薄云厚度分布;Dnear是在Hnew(x)中藍(lán)色地物附近人工選取出的薄云區(qū)域;n是Dnear中像素點(diǎn)個(gè)數(shù)。

    經(jīng)過以上的優(yōu)化,可得到完整的薄云厚度分布Htotal(x):

    (12)

    為了去除Htotal(x)中的塊狀效應(yīng),本文應(yīng)用導(dǎo)向?yàn)V波器[20]將其平滑得到最終的薄云厚度分布Hfinal(x)。

    本文方法的每一個(gè)優(yōu)化步驟結(jié)果如圖3所示。其中:圖3(b)是根據(jù)透射率得到的無云區(qū)域二值圖像,白色表示無云區(qū)域。圖3(c)是優(yōu)化后的薄云厚度分布Hnew(x)??梢钥闯?,Hnew(x)在無云區(qū)域的DN值接近或等于0,這能保證恢復(fù)的無云影像在無云區(qū)域的高保真度。圖3(d)是藍(lán)色地物的檢測結(jié)果, 可以發(fā)現(xiàn),圖中的藍(lán)色地物位置檢測比較準(zhǔn)確。圖3(e)是最終的薄云厚度分布Hfinal(x)。與圖3(c)相比,圖3(e)的Hfinal(x)在藍(lán)色地物區(qū)域的DN值相比Hnew(x)小很多,可以避免藍(lán)色地物云量被過估計(jì)。

    圖3 本文方法各步驟結(jié)果Fig. 3 Results of each step of the proposed method

    2.3 高保真無云影像

    通過逐步優(yōu)化薄云分布,單獨(dú)估計(jì)藍(lán)色地物云量,最終得到的無云影像J表示如下:

    Jc(x)=Ic(x)-kcHfinal(x)

    (13)

    其中:I是降質(zhì)影像;kc是根據(jù)文獻(xiàn)[16]算法中通過線性回歸得到的不同波段云估計(jì)強(qiáng)度;Hfinal(x)是最終優(yōu)化的薄云厚度分布。圖3(f)是本文改進(jìn)方法得到的無云影像,對比圖3(a)的降質(zhì)影像可以發(fā)現(xiàn),圖3(f)的無云影像不但具有良好的去云效果,而且無云區(qū)域的地物色彩基本不變,藍(lán)色地物色彩自然,即本文方法得到的無云云影像具有高保真的特點(diǎn)。

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

    本文采用三幅大小均是512×512遙感影像作為實(shí)驗(yàn)對象。圖4是分辨率為30 m的Landsat7遙感影像;圖5是分辨率為5.8 m的ZY3遙感影像;圖6是分辨率為2.0 m的WV-2遙感影像。數(shù)據(jù)包括了城市、農(nóng)田、草地等多種地物類型。本文采用文獻(xiàn)[14]方法、文獻(xiàn)[15]方法以及文獻(xiàn)[16]方法與本文方法進(jìn)行比較。

    圖4 不同方法去云結(jié)果比較(Landsat7遙感影像)Fig. 4 Haze removal result comparison of different methods (Landsat7 remote sensing image)

    圖5 不同方法去云結(jié)果比較(ZY3遙感影像)Fig. 5 Haze removal result comparison of different methods (ZY3 remote sensing image)

    從圖4~6可看出,4種算法都具有良好的去云效果。文獻(xiàn)[14]去云方法會使降質(zhì)影像的云區(qū)變得有些模糊,使得地物細(xì)節(jié)不夠清晰,且去云影像的無云區(qū)邊界色彩過渡不自然。文獻(xiàn)[15]去云方法對影像地物數(shù)據(jù)的保真度較差,偏綠地物太過鮮艷,色彩不自然。文獻(xiàn)[16]去云方法對非藍(lán)色地物整體有一定的保真效果,然而在藍(lán)色地物區(qū)域(圖5(d)和圖6(d)中方塊)去云后色彩與原圖中的藍(lán)色地物(圖5(a)和圖6(a)中方塊)相比過于鮮艷,失真嚴(yán)重。本文方法得到的去云結(jié)果不但在無云區(qū)域有較高的保真度,而且影像中的藍(lán)色地物(圖5(e)和圖6(e)中方塊)去云后色彩自然,達(dá)到高保真去云目的。

    為了定量評估本文方法的去云效果,采用廣義奇異值分解(Generalized Singular Value Decomposition, GSVD)作為評價(jià)因子[22]。GSVD是用于評估圖像受到模糊或噪聲影響的程度,薄云覆蓋是光學(xué)遙感影像被模糊的主要原因。GSVD值越小,則影像的去云效果越好。

    不同方法的GSVD統(tǒng)計(jì)值如表1所示。從表1可看出,四種去云方法中,本文方法對應(yīng)的三波段GSVD平均值總體較小,表明本文方法具有理想的去云效果。

    圖6 不同方法去云結(jié)果比較(WV-2遙感影像)Fig. 6 Haze removal result comparison of different methods (WV-2 remote sensing image)表1 不同方法的GSVD統(tǒng)計(jì)值對比Tab. 1 Comparison of GSVD statistical values of different methods

    圖像波段降質(zhì)影像文獻(xiàn)[14]方法文獻(xiàn)[15]方法文獻(xiàn)[16]方法本文方法圖4藍(lán)波段25.8524.6225.1824.7925.14綠波段24.5824.2323.9824.0424.08紅波段23.3523.5623.3123.2123.20平均值24.6024.1424.1624.0124.14圖5藍(lán)波段22.5120.8118.6719.1812.15綠波段17.3514.8615.4716.9317.30紅波段19.6819.1518.5019.6919.44平均值19.8518.2717.5518.6016.30圖6藍(lán)波段15.1713.6313.7613.9514.35綠波段13.2913.7212.8113.0912.37紅波段10.029.459.709.9810.01平均值12.8312.2712.0912.3412.24

    為了定量評價(jià)降質(zhì)影像無云區(qū)域的保真度,本文采用平均相對誤差(Mean Relative Error, MRE)作為評價(jià)指標(biāo),結(jié)果記作δ:

    (14)

    其中:Dclea是降質(zhì)影像I(x)和恢復(fù)的無云影像J(x)的同一片無云區(qū)域;m是無云區(qū)內(nèi)像素點(diǎn)個(gè)數(shù)。δ值越小保真效果越好,該評價(jià)指標(biāo)MRE的統(tǒng)計(jì)值如表2所示。

    表2 不同方法MRE統(tǒng)計(jì)值對比Tab. 2 Comparison of MRE statistical values of different methods

    由表2可看出,四種去云方法都具有高度的保真特性,文獻(xiàn)[15]方法去云結(jié)果的數(shù)據(jù)保真性能相比其他方法略差一些。本文方法的無云區(qū)域保真性能最好。綜合兩個(gè)評價(jià)指標(biāo)可以看出,本文方法能得到高保真的最優(yōu)去云結(jié)果。

    通過對不同數(shù)據(jù)范圍、不同分辨率大小、多種地物類型的薄云污染影像數(shù)據(jù)實(shí)驗(yàn)得出,本文方法得到的無云影像都具有數(shù)據(jù)高保真效果,表明本文方法對不同薄云污染影像去云具有高度適用性。

    4 結(jié)語

    本文提出了一種基于薄云厚度評估的高保真去薄云方法,不僅能夠滿足無云影像中無云區(qū)域數(shù)據(jù)的高保真度,還能夠解決傳統(tǒng)算法藍(lán)色地物色彩失真問題。針對影響算法的兩大因素:無云區(qū)域Dclean的確定和閾值T的選取,本文通過透射率分布圖中像素?cái)?shù)值接近于1實(shí)現(xiàn)了自動化確定無云區(qū)域的范圍,然而在確定藍(lán)色區(qū)域的過程中,本文是以經(jīng)驗(yàn)閾值T選取藍(lán)色地物和白色地物區(qū)域,自適應(yīng)求取最優(yōu)閾值T將是我們后續(xù)的改進(jìn)工作。

    猜你喜歡
    薄云波段藍(lán)色
    春日暖陽
    藍(lán)色的祝愿
    青年歌聲(2021年2期)2021-03-05 09:02:14
    藍(lán)色的祝愿
    青年歌聲(2020年11期)2020-11-24 06:57:20
    藍(lán)色使命
    商周刊(2018年13期)2018-07-11 03:34:10
    藍(lán)色的五月
    海峽姐妹(2017年5期)2017-06-05 08:53:18
    M87的多波段輻射過程及其能譜擬合
    日常維護(hù)對L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    L波段雷達(dá)磁控管的使用與維護(hù)
    河南科技(2014年18期)2014-02-27 14:14:53
    利用小波變換去除遙感影像薄云研究
    蘇木山觀日出
    老年世界(2013年22期)2013-12-16 08:09:38
    亚洲在线自拍视频| 毛片女人毛片| 99久久精品国产国产毛片| 99在线视频只有这里精品首页| 久久人人爽人人片av| 一区二区三区乱码不卡18| 1024手机看黄色片| a级毛色黄片| 亚洲精品一区蜜桃| 日日摸夜夜添夜夜添av毛片| 久久鲁丝午夜福利片| 欧美激情久久久久久爽电影| 91久久精品电影网| 国内揄拍国产精品人妻在线| 99久久中文字幕三级久久日本| 久久婷婷人人爽人人干人人爱| 亚州av有码| av免费观看日本| 色综合亚洲欧美另类图片| 深爱激情五月婷婷| 亚洲四区av| 精品99又大又爽又粗少妇毛片| 久久久欧美国产精品| 小说图片视频综合网站| 成人三级黄色视频| 亚洲国产欧美人成| 99久久中文字幕三级久久日本| 99九九线精品视频在线观看视频| 久久久久久大精品| 一级毛片久久久久久久久女| 国产精品久久久久久av不卡| 日韩高清综合在线| 超碰97精品在线观看| 最近最新中文字幕大全电影3| 免费av观看视频| 亚洲,欧美,日韩| 日韩欧美 国产精品| 国产69精品久久久久777片| 国产成人a区在线观看| 欧美不卡视频在线免费观看| 只有这里有精品99| 亚洲av免费在线观看| 九九久久精品国产亚洲av麻豆| 久久精品影院6| 欧美性猛交╳xxx乱大交人| 听说在线观看完整版免费高清| 好男人在线观看高清免费视频| 干丝袜人妻中文字幕| 国产69精品久久久久777片| 18禁在线无遮挡免费观看视频| 亚洲精品国产av成人精品| 在线免费十八禁| 99热全是精品| 国产69精品久久久久777片| 色尼玛亚洲综合影院| 99热精品在线国产| 亚洲欧美成人精品一区二区| 岛国毛片在线播放| 免费观看在线日韩| 国产极品天堂在线| 久久99热这里只有精品18| 久久精品久久久久久噜噜老黄 | 色视频www国产| 午夜免费男女啪啪视频观看| 免费看a级黄色片| 久久久久久久午夜电影| 亚洲一区高清亚洲精品| 亚洲,欧美,日韩| 国产真实伦视频高清在线观看| 最近视频中文字幕2019在线8| 三级毛片av免费| 精品久久久久久久人妻蜜臀av| 亚洲精品乱码久久久久久按摩| 国产精品三级大全| 99热这里只有是精品在线观看| 尾随美女入室| 久久99蜜桃精品久久| 一个人免费在线观看电影| 亚洲一区高清亚洲精品| 夫妻性生交免费视频一级片| 亚洲18禁久久av| 能在线免费看毛片的网站| 日韩欧美国产在线观看| 我的老师免费观看完整版| 变态另类丝袜制服| 成人av在线播放网站| 精品国内亚洲2022精品成人| 午夜福利在线在线| 久久韩国三级中文字幕| 十八禁国产超污无遮挡网站| 天天躁夜夜躁狠狠久久av| 国产视频首页在线观看| 极品教师在线视频| 欧美区成人在线视频| 亚洲精品色激情综合| 精品欧美国产一区二区三| 青春草视频在线免费观看| 国产女主播在线喷水免费视频网站 | 少妇人妻一区二区三区视频| av卡一久久| 午夜福利成人在线免费观看| 欧美区成人在线视频| 亚洲综合色惰| 免费黄色在线免费观看| 日韩,欧美,国产一区二区三区 | 亚洲在线自拍视频| 男女啪啪激烈高潮av片| 国产免费视频播放在线视频 | 在线观看66精品国产| 国产精品,欧美在线| 亚洲av免费高清在线观看| 秋霞伦理黄片| 91av网一区二区| 岛国在线免费视频观看| 国产乱人偷精品视频| 你懂的网址亚洲精品在线观看 | 赤兔流量卡办理| 国产极品天堂在线| 成人无遮挡网站| 国产精品麻豆人妻色哟哟久久 | 国国产精品蜜臀av免费| 久久99热6这里只有精品| 不卡视频在线观看欧美| 在线免费观看的www视频| 国产精品熟女久久久久浪| 一级毛片aaaaaa免费看小| 神马国产精品三级电影在线观看| 一个人观看的视频www高清免费观看| 高清在线视频一区二区三区 | 亚洲国产最新在线播放| 啦啦啦韩国在线观看视频| 久久欧美精品欧美久久欧美| 六月丁香七月| 麻豆一二三区av精品| 日本与韩国留学比较| av在线播放精品| 99久久成人亚洲精品观看| 精品熟女少妇av免费看| 丰满乱子伦码专区| 蜜桃久久精品国产亚洲av| 精品国产露脸久久av麻豆 | 91在线精品国自产拍蜜月| 国产av一区在线观看免费| 国产极品精品免费视频能看的| 精品久久久久久久久av| 色视频www国产| 国语自产精品视频在线第100页| 最近最新中文字幕免费大全7| 欧美色视频一区免费| 亚洲国产日韩欧美精品在线观看| 精品一区二区免费观看| 午夜激情欧美在线| 国语自产精品视频在线第100页| 亚洲18禁久久av| 麻豆精品久久久久久蜜桃| 三级经典国产精品| 亚洲国产欧洲综合997久久,| 国产精品久久久久久久久免| 国产免费一级a男人的天堂| 我要看日韩黄色一级片| av黄色大香蕉| 亚洲真实伦在线观看| 亚洲av熟女| 免费黄网站久久成人精品| 人人妻人人澡人人爽人人夜夜 | 青春草亚洲视频在线观看| 亚洲av电影不卡..在线观看| 欧美三级亚洲精品| 天美传媒精品一区二区| 免费人成在线观看视频色| 精品人妻视频免费看| 亚洲高清免费不卡视频| av专区在线播放| 久久国产乱子免费精品| 欧美色视频一区免费| 国产精品99久久久久久久久| 亚洲国产欧洲综合997久久,| 欧美3d第一页| 人妻夜夜爽99麻豆av| 视频中文字幕在线观看| 国产成人精品一,二区| 日韩精品青青久久久久久| av视频在线观看入口| 亚洲人成网站高清观看| 久久久久精品久久久久真实原创| 少妇人妻精品综合一区二区| 精品国产露脸久久av麻豆 | 天天躁日日操中文字幕| 日韩欧美在线乱码| 在线免费观看的www视频| a级毛色黄片| 色网站视频免费| 亚洲精品国产成人久久av| 内地一区二区视频在线| av.在线天堂| 国产高清国产精品国产三级 | 日本猛色少妇xxxxx猛交久久| 一级av片app| 在线天堂最新版资源| 久久精品综合一区二区三区| av黄色大香蕉| 国产一区二区在线观看日韩| 国产黄色视频一区二区在线观看 | 日韩亚洲欧美综合| 麻豆成人av视频| 99久久精品热视频| 韩国高清视频一区二区三区| 亚洲电影在线观看av| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久久免| 国产三级中文精品| 精品欧美国产一区二区三| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产欧美人成| av在线老鸭窝| 搡女人真爽免费视频火全软件| 最近的中文字幕免费完整| 国产精品一区二区在线观看99 | 亚洲aⅴ乱码一区二区在线播放| 亚洲天堂国产精品一区在线| 国产一区二区在线观看日韩| 三级国产精品片| 在线观看66精品国产| 国产高清国产精品国产三级 | 乱人视频在线观看| 成人一区二区视频在线观看| 国产精品久久久久久精品电影| 狂野欧美白嫩少妇大欣赏| 精品国内亚洲2022精品成人| 欧美日韩精品成人综合77777| 国产精品熟女久久久久浪| 免费看光身美女| 国产亚洲91精品色在线| 亚洲国产精品国产精品| 欧美一区二区精品小视频在线| 岛国在线免费视频观看| 小蜜桃在线观看免费完整版高清| 久久久久久久久久黄片| 女的被弄到高潮叫床怎么办| 亚洲精品日韩在线中文字幕| 99久国产av精品国产电影| 人妻制服诱惑在线中文字幕| 深爱激情五月婷婷| 国产色婷婷99| 久久婷婷人人爽人人干人人爱| 三级国产精品片| 成年女人看的毛片在线观看| 99久久无色码亚洲精品果冻| 一个人看视频在线观看www免费| 九色成人免费人妻av| 久热久热在线精品观看| 欧美又色又爽又黄视频| 午夜爱爱视频在线播放| 亚洲中文字幕一区二区三区有码在线看| 成人午夜高清在线视频| 久久久久免费精品人妻一区二区| 久久精品影院6| av线在线观看网站| 日韩中字成人| 成年免费大片在线观看| 精品人妻偷拍中文字幕| 久久久久九九精品影院| 人人妻人人澡人人爽人人夜夜 | 一级毛片我不卡| 搞女人的毛片| av在线天堂中文字幕| 99热精品在线国产| 国产高清视频在线观看网站| 丝袜喷水一区| 国产私拍福利视频在线观看| 国产精品永久免费网站| 精品久久久久久久久亚洲| 国产精品美女特级片免费视频播放器| 我的老师免费观看完整版| 好男人视频免费观看在线| 日韩成人伦理影院| 日韩强制内射视频| 国产激情偷乱视频一区二区| av又黄又爽大尺度在线免费看 | 2021天堂中文幕一二区在线观| 精品欧美国产一区二区三| av福利片在线观看| 亚洲av电影在线观看一区二区三区 | 99久国产av精品国产电影| 看十八女毛片水多多多| 国产黄片视频在线免费观看| 欧美日本视频| 嫩草影院精品99| 国产探花在线观看一区二区| 爱豆传媒免费全集在线观看| 中文在线观看免费www的网站| 国产不卡一卡二| 国产黄色视频一区二区在线观看 | 国产精品一二三区在线看| av.在线天堂| 99国产精品一区二区蜜桃av| 天天一区二区日本电影三级| 我要搜黄色片| 麻豆成人av视频| 丰满乱子伦码专区| 赤兔流量卡办理| 欧美精品国产亚洲| 欧美97在线视频| 99久久精品热视频| 亚洲国产精品国产精品| www.色视频.com| 嘟嘟电影网在线观看| 男女啪啪激烈高潮av片| 狂野欧美白嫩少妇大欣赏| 大话2 男鬼变身卡| 国产精品熟女久久久久浪| 欧美日韩综合久久久久久| 日韩av在线大香蕉| 亚洲色图av天堂| 91在线精品国自产拍蜜月| 日韩成人伦理影院| 久久久亚洲精品成人影院| 中文字幕熟女人妻在线| 日韩欧美在线乱码| 中文字幕熟女人妻在线| 国产亚洲精品av在线| 亚州av有码| 级片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 桃色一区二区三区在线观看| 午夜福利在线在线| 一级黄片播放器| 欧美精品一区二区大全| 国产精品一区www在线观看| 高清午夜精品一区二区三区| 国产精华一区二区三区| 亚洲欧美精品自产自拍| 日韩欧美精品免费久久| 国产老妇女一区| 18禁在线无遮挡免费观看视频| 国产精品一区二区三区四区免费观看| 狠狠狠狠99中文字幕| 国产免费男女视频| 中国国产av一级| 少妇的逼水好多| 国产亚洲91精品色在线| 美女被艹到高潮喷水动态| av又黄又爽大尺度在线免费看 | 91午夜精品亚洲一区二区三区| 国产淫片久久久久久久久| 你懂的网址亚洲精品在线观看 | 男人舔奶头视频| 亚洲国产色片| 亚洲激情五月婷婷啪啪| 国产中年淑女户外野战色| 在线观看av片永久免费下载| 国产成人aa在线观看| 身体一侧抽搐| 午夜精品在线福利| 亚洲电影在线观看av| 男插女下体视频免费在线播放| 人人妻人人澡人人爽人人夜夜 | 久99久视频精品免费| 边亲边吃奶的免费视频| 亚洲成人久久爱视频| 身体一侧抽搐| 亚洲欧美一区二区三区国产| 午夜日本视频在线| 中文字幕精品亚洲无线码一区| 日韩av不卡免费在线播放| 水蜜桃什么品种好| 日产精品乱码卡一卡2卡三| 麻豆乱淫一区二区| 成人三级黄色视频| 国产单亲对白刺激| 亚洲欧美日韩无卡精品| 国产精品久久久久久久久免| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| 九九久久精品国产亚洲av麻豆| 国产久久久一区二区三区| 在线观看一区二区三区| 又爽又黄a免费视频| 精品久久久久久电影网 | 国国产精品蜜臀av免费| 欧美日本视频| 91在线精品国自产拍蜜月| 久久久久国产网址| 亚洲国产精品成人综合色| 国产精品99久久久久久久久| 日日干狠狠操夜夜爽| 国产毛片a区久久久久| 少妇猛男粗大的猛烈进出视频 | 一个人看视频在线观看www免费| 国产黄色小视频在线观看| 久久鲁丝午夜福利片| 亚洲av熟女| 三级经典国产精品| 美女cb高潮喷水在线观看| 亚洲第一区二区三区不卡| 久久久色成人| 国产女主播在线喷水免费视频网站 | 超碰97精品在线观看| 亚洲精品日韩av片在线观看| 亚洲久久久久久中文字幕| 久久99热这里只有精品18| 乱码一卡2卡4卡精品| 亚洲av.av天堂| 久久久久久久亚洲中文字幕| 高清av免费在线| 久久精品熟女亚洲av麻豆精品 | 一夜夜www| 黄色一级大片看看| 久久韩国三级中文字幕| 久久久久久久久久久免费av| 日韩高清综合在线| 中文亚洲av片在线观看爽| 日本一二三区视频观看| 久久热精品热| 国产单亲对白刺激| 久久久成人免费电影| 国产高清三级在线| 久久精品综合一区二区三区| 精品少妇黑人巨大在线播放 | 亚洲中文字幕一区二区三区有码在线看| 老司机影院成人| 精品不卡国产一区二区三区| 日本欧美国产在线视频| 精品熟女少妇av免费看| 午夜激情欧美在线| 五月伊人婷婷丁香| 在线播放国产精品三级| АⅤ资源中文在线天堂| 国产精品人妻久久久久久| 精品久久久久久久久av| 伦精品一区二区三区| 嫩草影院新地址| 欧美不卡视频在线免费观看| 国产亚洲精品久久久com| 22中文网久久字幕| 成人鲁丝片一二三区免费| 在线观看一区二区三区| 亚洲人成网站高清观看| 国产精品爽爽va在线观看网站| 国产精品99久久久久久久久| videossex国产| 国内精品美女久久久久久| 热99re8久久精品国产| 欧美日韩在线观看h| 内地一区二区视频在线| 亚洲av福利一区| 看黄色毛片网站| 国产精品女同一区二区软件| 少妇高潮的动态图| 亚洲av免费高清在线观看| 日韩精品青青久久久久久| 女的被弄到高潮叫床怎么办| www日本黄色视频网| 尤物成人国产欧美一区二区三区| 亚洲熟妇中文字幕五十中出| 日日摸夜夜添夜夜添av毛片| 一个人观看的视频www高清免费观看| 国产精品久久久久久久电影| 最近2019中文字幕mv第一页| 国产真实伦视频高清在线观看| 亚洲精品乱码久久久v下载方式| 在线免费观看的www视频| 联通29元200g的流量卡| 天堂√8在线中文| 中文资源天堂在线| 久久韩国三级中文字幕| 我的老师免费观看完整版| 国产精品国产三级国产专区5o | 国产成人精品一,二区| 综合色av麻豆| www.av在线官网国产| 男插女下体视频免费在线播放| 成人毛片60女人毛片免费| 建设人人有责人人尽责人人享有的 | 国产一区二区在线av高清观看| 亚洲欧美成人精品一区二区| 天天躁日日操中文字幕| 欧美精品一区二区大全| 人妻少妇偷人精品九色| 寂寞人妻少妇视频99o| 亚洲欧美清纯卡通| 成人高潮视频无遮挡免费网站| 亚洲国产欧美人成| 中文在线观看免费www的网站| 菩萨蛮人人尽说江南好唐韦庄 | 看片在线看免费视频| 国产不卡一卡二| 99热这里只有精品一区| 久久这里只有精品中国| 亚洲图色成人| 搡女人真爽免费视频火全软件| av线在线观看网站| 热99在线观看视频| 最近的中文字幕免费完整| 我的老师免费观看完整版| 最近手机中文字幕大全| 国产精品一区二区三区四区免费观看| 淫秽高清视频在线观看| 亚洲国产精品久久男人天堂| www日本黄色视频网| 非洲黑人性xxxx精品又粗又长| 全区人妻精品视频| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久电影| 国产精品国产三级国产专区5o | 老师上课跳d突然被开到最大视频| 亚洲国产欧美人成| 欧美高清成人免费视频www| 久久综合国产亚洲精品| 亚洲av二区三区四区| 欧美色视频一区免费| 深夜a级毛片| 中文在线观看免费www的网站| 91久久精品国产一区二区三区| 久久久久久久久久久免费av| 水蜜桃什么品种好| 搡女人真爽免费视频火全软件| 国产欧美日韩精品一区二区| 一二三四中文在线观看免费高清| 我的老师免费观看完整版| 国产一级毛片在线| 中文精品一卡2卡3卡4更新| 国产精品人妻久久久久久| 好男人视频免费观看在线| 日韩av不卡免费在线播放| 欧美成人精品欧美一级黄| 久久久久久久久久黄片| a级一级毛片免费在线观看| 一个人免费在线观看电影| 成年女人看的毛片在线观看| 桃色一区二区三区在线观看| 亚洲一级一片aⅴ在线观看| 色5月婷婷丁香| 欧美成人一区二区免费高清观看| 国产极品天堂在线| 国产综合懂色| 国产一区二区亚洲精品在线观看| 国产精品国产高清国产av| 国产人妻一区二区三区在| 亚洲精品日韩在线中文字幕| 欧美另类亚洲清纯唯美| 亚洲18禁久久av| 国产成人freesex在线| 1024手机看黄色片| 男人舔女人下体高潮全视频| 午夜视频国产福利| 97热精品久久久久久| 日本五十路高清| 高清在线视频一区二区三区 | 国产高清三级在线| 精品免费久久久久久久清纯| 国产一级毛片七仙女欲春2| 极品教师在线视频| 国产精品三级大全| av福利片在线观看| 久久99热6这里只有精品| 亚洲精品乱久久久久久| 欧美日韩国产亚洲二区| 午夜福利视频1000在线观看| av免费观看日本| 一级毛片aaaaaa免费看小| 99久久无色码亚洲精品果冻| 亚洲av中文av极速乱| 乱码一卡2卡4卡精品| 免费观看人在逋| 午夜福利网站1000一区二区三区| 国产精品久久电影中文字幕| 亚洲最大成人av| 欧美日韩一区二区视频在线观看视频在线 | 国产一区二区在线观看日韩| a级一级毛片免费在线观看| 亚洲婷婷狠狠爱综合网| 人人妻人人澡欧美一区二区| 十八禁国产超污无遮挡网站| av在线播放精品| 男女啪啪激烈高潮av片| 国产综合懂色| 日日干狠狠操夜夜爽| 男女啪啪激烈高潮av片| 国产又黄又爽又无遮挡在线| 久久99热这里只频精品6学生 | 男插女下体视频免费在线播放| 美女大奶头视频| 国产激情偷乱视频一区二区| 一级毛片我不卡| 亚洲精品日韩av片在线观看| 免费一级毛片在线播放高清视频| 一边摸一边抽搐一进一小说| 欧美日韩在线观看h| 一级毛片aaaaaa免费看小| 午夜福利在线观看吧| 国产亚洲5aaaaa淫片| 国产成人a区在线观看| 丝袜美腿在线中文| 神马国产精品三级电影在线观看| 国产一区有黄有色的免费视频 | 精品一区二区三区视频在线| 九九久久精品国产亚洲av麻豆| 麻豆成人av视频| 成人性生交大片免费视频hd| 久久久久久伊人网av| 国产精品乱码一区二三区的特点| 简卡轻食公司| 久久99热这里只有精品18| 欧美+日韩+精品| 国产精品一区www在线观看| 日日撸夜夜添| 国内精品一区二区在线观看| 少妇熟女欧美另类| 亚洲国产欧洲综合997久久,| 黄色日韩在线| 麻豆av噜噜一区二区三区| 国产成人福利小说| 精品国产一区二区三区久久久樱花 | 国产免费男女视频|