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

    應(yīng)用植被變化與火燒跡地對(duì)森林火災(zāi)反演1)

    2021-05-26 07:56:58王金鑫桑學(xué)鋒劉鑫常家軒鄭陽(yáng)李子恒
    關(guān)鍵詞:跡地火燒植被指數(shù)

    王金鑫 桑學(xué)鋒 劉鑫 常家軒 鄭陽(yáng) 李子恒

    (中國(guó)水利水電科學(xué)研究院水資源所,北京,100038)

    準(zhǔn)確評(píng)估森林火災(zāi)后的損失情況對(duì)于預(yù)判林火蔓延發(fā)展態(tài)勢(shì)、組織林火撲救、安全滅火、災(zāi)后重建、維護(hù)生態(tài)平衡以及保護(hù)森林資源具有現(xiàn)實(shí)意義[1]。隨著遙感技術(shù)的不斷發(fā)展,衛(wèi)星的時(shí)空分辨率的不斷提高,可以通過(guò)衛(wèi)星影像實(shí)時(shí)觀測(cè)火災(zāi)過(guò)程并且提取有效信息,實(shí)現(xiàn)對(duì)森林火災(zāi)的全過(guò)程進(jìn)行動(dòng)態(tài)反演以及災(zāi)害評(píng)估。災(zāi)區(qū)的植被變化與火燒跡地面積是火災(zāi)反演的重要評(píng)估指標(biāo)[2]。近年來(lái),國(guó)內(nèi)外許多學(xué)者提出了基于遙感技術(shù)的火災(zāi)反演評(píng)估方法,其中常用的是利用遙感影像中的光譜構(gòu)造植被指數(shù)與燃燒指數(shù)來(lái)突出植被損傷區(qū)域與火災(zāi)的火燒跡地面積,其原理是基于火災(zāi)后植被燃燒減少和燃燒殘留物沉積這兩個(gè)變化特征。Allen et al.[3]通過(guò)Landsat TM/ETM+數(shù)據(jù)計(jì)算了火災(zāi)前后的歸一化燃燒比(NBR)與差分歸一化燃燒比(dNBR)對(duì)阿拉斯加北部森林火災(zāi)進(jìn)行了評(píng)估; Veraverbeke et al.[4]通過(guò)Landsat TM/ETM+與MODIS數(shù)據(jù),利用差分歸一化燃燒比對(duì)希臘2007年伯羅奔尼撒半島的森林大火進(jìn)行不同時(shí)序的火災(zāi)反演以及火災(zāi)嚴(yán)重程度評(píng)估;Chang et al.[5]通過(guò)構(gòu)建歸一化燃燒比指數(shù)對(duì)黑龍江森林大火進(jìn)行了火災(zāi)嚴(yán)重程度的評(píng)估;Escuin et al.[6]從LANDSAT TM/ETM圖像中分析提取的歸一化燃燒比和歸一化植被指數(shù),對(duì)西班牙南部發(fā)生的3起森林火災(zāi)進(jìn)行了火災(zāi)嚴(yán)重程度評(píng)估;張春桂等[7]利用分辨率250 m的 MODIS可見(jiàn)光數(shù)據(jù),根據(jù)MODIS的歸一化植被指數(shù),在火災(zāi)發(fā)生前后的變化特征估算了林火面積。傳統(tǒng)方法只依據(jù)單一指數(shù)指標(biāo)對(duì)火災(zāi)損失進(jìn)行評(píng)價(jià),往往忽視了災(zāi)區(qū)的植被損失與火燒跡地的區(qū)別;實(shí)際上,在植被損失估算中往往只計(jì)算出災(zāi)后植被指數(shù),通過(guò)目視解譯而估算植被損失,實(shí)際操作易產(chǎn)生誤差。在火燒跡地面積評(píng)估中,通過(guò)計(jì)算火燒跡地的燃燒指數(shù)來(lái)區(qū)別過(guò)火區(qū)與未過(guò)火區(qū),對(duì)如何能夠準(zhǔn)確快速追蹤識(shí)別與提取火燒跡地面積的方法很少有考慮。

    本文基于Sentinel-2影像,以木里森林火災(zāi)為例,采用基于火災(zāi)前后歸一化差分植被指數(shù)、差分歸一化燃燒比(dNBR)、迭代加權(quán)多元變化檢測(cè)算法(IR-MAD)、追蹤識(shí)別提取算法等分析火災(zāi)前后的植被變化和火燒跡地面積,為森林火災(zāi)快速反演與評(píng)估提供數(shù)據(jù)支持。

    1 研究區(qū)概況

    2020年3月28日四川省涼山州木里藏族自治縣項(xiàng)腳鄉(xiāng)發(fā)生森林火災(zāi),火場(chǎng)位于木里縣東南部,西側(cè)距離木里縣城約5 km。其地理位置為101°18′~101°31′E,27°51′~28°2′N。災(zāi)區(qū)位于高山峽谷地帶,森林覆蓋率高,植被多為松樹(shù)、桉樹(shù)等易燃植被,其枝葉中均含有大量油脂,且地表可燃物復(fù)雜(包括倒木、腐木、高山腐殖層、泥炭等),一旦燃燒,會(huì)分解產(chǎn)生大量揮發(fā)性易燃?xì)怏w,導(dǎo)致燃燒劇烈。木里森林火災(zāi)和平原、盆地、丘陵等地區(qū)不同,涼山州地處偏僻,森林植被處于的高山峽谷,給滅火工作帶來(lái)極大的挑戰(zhàn),地面囿于視野限制,難以判斷火勢(shì)走向,衛(wèi)星“站得高看得遠(yuǎn)”,可以快速發(fā)現(xiàn)火情,判斷火情趨勢(shì),從宏觀的角度把握火災(zāi)全貌,有利于應(yīng)急救援指揮與決策。

    根據(jù)涼山州氣象信息中心的數(shù)據(jù)顯示,火災(zāi)期間,正值干風(fēng)季節(jié),晴熱天氣異常,連續(xù)20 d無(wú)任何降水,空氣濕度為5%~10%;特別是3月下旬已連續(xù)4 d高溫天氣,最高氣溫達(dá)到31.2 ℃,風(fēng)力7~8級(jí),火災(zāi)期間全州平均溫度較歷史同期偏高2.0 ℃;全州平均降水量12 mm,木里比歷史同期減少64%的降水量。在高溫、大風(fēng)和干燥氣候加持下,引起了木里特大森林火災(zāi)。

    2 研究方法

    2.1 影像數(shù)據(jù)

    Sentinel-2影像具有 13 個(gè)光譜波段,空間分辨率為10 m[8]。采用SNAP軟件對(duì)數(shù)據(jù)進(jìn)行輻射定標(biāo)、幾何校正、波段融合和影像裁剪。Sentinel-2在重訪周期、空間分辨率以及光譜分辨率等方面均優(yōu)于大多數(shù)遙感影像,能夠?yàn)樯只馂?zāi)反演評(píng)估提供更高精度的數(shù)據(jù)。本文使用Sentinel-2災(zāi)前(3月25日)、災(zāi)中(3月30日)、災(zāi)后(4月9日)無(wú)云影像數(shù)據(jù)各2幅,Landsat-8災(zāi)后(4月14日)影像一幅?;饒?chǎng)影像如圖1所示。

    圖1 不同時(shí)間的火災(zāi)區(qū)的影像

    2.2 迭代加權(quán)多元檢測(cè)算法的植被變化檢測(cè)

    通過(guò)計(jì)算火災(zāi)前后的歸一化植被指數(shù)的遙感影像后,利用迭代加權(quán)多元檢測(cè)算法(IR-MAD)進(jìn)行火災(zāi)前后的植被變化檢測(cè)。通過(guò)變化檢測(cè)火災(zāi)前后植被的變化像元,分離出健康植被與受災(zāi)植被,制作出植被受損專題圖。

    歸一化植被指數(shù):反映土地覆蓋植被狀況的遙感指標(biāo),也是反映植被長(zhǎng)勢(shì)和營(yíng)養(yǎng)信息的重要參數(shù)[9]。通過(guò)計(jì)算火災(zāi)前后的歸一化植被指數(shù),能夠?yàn)橛行гu(píng)估植被受損區(qū)域提供支持。遙感影像中,歸一化植被指數(shù)的計(jì)算公式為:INDVI=(NR-R)/(NR+R)。式中,NR為近紅外波段的反射值,R為紅光波段的反射值。

    迭代加權(quán)多元檢測(cè)算法(IR-MAD):多元變化檢測(cè)算法(MAD)[10]其數(shù)學(xué)本質(zhì),主要是多元統(tǒng)計(jì)分析中的典型相關(guān)分析(CCA)以及波段差值運(yùn)算,但該算法仍然不能完全改善目前多元遙感影像處理中的局限性。Canty et al.[11]在MAD算法的基礎(chǔ)上,并提出了迭代加權(quán)多元檢測(cè)算法。IR-MAD是一種檢測(cè)多元影像變化的方法,具有準(zhǔn)確的獲取變化信息以及受外界因素影響較小的特點(diǎn),因此,在多元影像檢測(cè)變化中被廣泛應(yīng)用[12]。IR-MAD核心思想是每個(gè)像元初始權(quán)重為1,每一次迭代均賦予2幅影像中每個(gè)像元新的權(quán)重,通過(guò)計(jì)算,未發(fā)生變化的像元具有較大的權(quán)重,最終得到的權(quán)重是決定各個(gè)像元是否發(fā)生變化的依據(jù)。經(jīng)過(guò)若干次迭代后,每個(gè)像元的權(quán)重會(huì)趨于穩(wěn)定,直到變化小于設(shè)定的閾值或不再變化則停止迭代。研究證明,IR-MAD通過(guò)迭代更新權(quán)值,得到的變化區(qū)域更精確,噪聲更少,能夠更好地將變化信息從非變化區(qū)域中提取出來(lái)[13]。

    2.3 追蹤識(shí)別提取算法

    采用差分歸一化燃燒比和追蹤識(shí)別提取算法對(duì)研究區(qū)進(jìn)行火燒跡地面積的提取。

    采用差分歸一化燃燒比:采用差分歸一化燃燒比是在歸一化燃燒比方法的基礎(chǔ)上建立的[14]。通常運(yùn)用差分歸一化燃燒比和歸一化燃燒比生成森林大火燃燒后的圖像,可以獲得對(duì)燃燒嚴(yán)重性的初步評(píng)估數(shù)據(jù),并支持現(xiàn)場(chǎng)救災(zāi)工作。本文采用差分歸一化燃燒比對(duì)受災(zāi)區(qū)域制作火災(zāi)等級(jí)分布專題圖[15]。差分歸一化燃燒比和歸一化燃燒比的計(jì)算公式為:INBR=(NR-SR)/(NR+SR);IdNBR=IF-IP。式中,INBR為歸一化燃燒比,IdNBR 為差分歸一化燃燒比,NR為近紅外波段的反射值,SR為短波紅外波段的反射值,IF為災(zāi)區(qū)火災(zāi)前的歸一化燃燒比,IP為火災(zāi)后的歸一化燃燒比。

    運(yùn)用目標(biāo)追蹤識(shí)別提取算法能夠進(jìn)行追蹤影像中所有連通區(qū),并根據(jù)閾值提取目標(biāo)的面積信息,與兩邊掃描法、區(qū)域增長(zhǎng)法相比,掃描重復(fù)率低、精度高[16]。差分歸一化燃燒比和火燒跡地面積追蹤識(shí)別提取影像信息,是在差分歸一化燃燒比的突出燃燒區(qū)域信息基礎(chǔ)上,繼續(xù)火燒跡地信息的追蹤識(shí)別。算法原理是在差分歸一化燃燒比提取后的影像中,自主搜尋種子點(diǎn),遍歷整幅影像;根據(jù)每個(gè)點(diǎn)的性質(zhì)判斷該點(diǎn)是否符合要求,決定種子點(diǎn);再根據(jù)種子點(diǎn)的4-鄰域追蹤連通區(qū)域,將種子點(diǎn)的鄰域判斷結(jié)束后,確定該點(diǎn)和4-鄰域的連通性。根據(jù)火燒跡地的特有光譜信息范圍作為閾值,采用雙閾值法來(lái)分類過(guò)火區(qū),對(duì)過(guò)火區(qū)域進(jìn)行標(biāo)記,分別計(jì)算不同連通區(qū)域面積對(duì)應(yīng)的影像中的火燒跡地面積,從而實(shí)現(xiàn)整幅影像中的火燒跡地面積一次性識(shí)別提取。主要方法流程如下:

    (1)從影像的第二行第二列開(kāi)始掃描,若掃描到A點(diǎn),先對(duì)A點(diǎn)的像素性質(zhì)進(jìn)行判斷,若符合種子點(diǎn)要求則對(duì)其鄰域進(jìn)行判別;

    (2)若鄰域內(nèi)有與A點(diǎn)具有相同性質(zhì)的點(diǎn)則進(jìn)行標(biāo)記、儲(chǔ)存;

    (3)4-鄰域內(nèi)像素判斷完后,將儲(chǔ)存的具有相同性質(zhì)的點(diǎn)提取出來(lái),一個(gè)連通區(qū)域就被提取出來(lái),表明該區(qū)域追蹤完成;

    (4)接著對(duì)下一個(gè)像素進(jìn)行步驟(2)和步驟(3),若遇到的點(diǎn)已經(jīng)被標(biāo)記則跳過(guò);

    (5)所有圖像掃描后,將所有的連通區(qū)域排列,根據(jù)火燒跡地閾值提取各個(gè)火燒跡地連通區(qū)并計(jì)算像元點(diǎn)數(shù),將其自動(dòng)與影像分辨率相乘得出火燒跡地面積。

    3 結(jié)果與分析

    3.1 火災(zāi)前后植被變化檢測(cè)

    如圖2所示,通過(guò)對(duì)災(zāi)前災(zāi)后的影像進(jìn)行歸一化植被指數(shù)計(jì)算得到火災(zāi)前后的兩幅影像,觀察研究區(qū)災(zāi)前(3月25日)與災(zāi)后(4月9日)的無(wú)云兩幅影像,這一時(shí)間區(qū)域跨越了整個(gè)火災(zāi)的過(guò)程,將災(zāi)前災(zāi)后的植被變化顯現(xiàn)出來(lái)。

    如圖3所示,通過(guò)基于IDL語(yǔ)言編寫集成的IR-MAD功能模塊在ENVI5.4軟件中進(jìn)行加載使用,對(duì)火災(zāi)前后歸一化植被指數(shù)的影像變化進(jìn)行檢測(cè),變化結(jié)果被認(rèn)定為受損植被;將結(jié)果影像用閾值分割進(jìn)行色彩增強(qiáng)制成專題圖,突出了健康植被與受損植被的區(qū)別,經(jīng)過(guò)統(tǒng)計(jì)受損植被像元數(shù)量為1 464 408個(gè)像素。

    圖2 火災(zāi)前后NDVI影像

    圖3 火災(zāi)前后植被變化影像

    3.2 火燒跡地面積追蹤識(shí)別提取

    采用dNBR+火燒跡地面積追蹤識(shí)別提取算法是在dNBR結(jié)果基礎(chǔ)上對(duì)火燒跡地進(jìn)行二次處理分析,并對(duì)火燒跡地面積進(jìn)行一次性快速準(zhǔn)確提取與計(jì)算。使用火災(zāi)的前中后3個(gè)時(shí)段的影像,對(duì)災(zāi)中(3月30日)、災(zāi)后(4月9日)火燒跡地的面積進(jìn)行了快速識(shí)別提取。流程為:(1)計(jì)算3個(gè)時(shí)段火災(zāi)區(qū)域影像的NBR;(2)通過(guò)計(jì)算災(zāi)前與災(zāi)中、災(zāi)后的NBR差值得到災(zāi)中(3月30日)、災(zāi)后(4月9日)dNBR影像(見(jiàn)圖4);(3)再通過(guò)基于(2)中dNBR影像進(jìn)行火燒跡地面積追蹤識(shí)別提取,將災(zāi)中(3月30日)、災(zāi)后(4月9日)兩個(gè)時(shí)序的災(zāi)區(qū)火燒跡地面積一次性提出并計(jì)算,結(jié)果見(jiàn)圖5、圖6,圖中白色區(qū)域?yàn)榛馃E地。

    圖4 災(zāi)中災(zāi)后dNBR影像

    圖5 火災(zāi)等級(jí)分布專題圖

    圖6 dNBR+火燒跡地面積識(shí)別提取成果影像

    通過(guò)dNBR+火燒跡地面積追蹤識(shí)別提取算法進(jìn)行實(shí)驗(yàn),提取火災(zāi)中(3月30日)火燒跡地斑塊18個(gè),火燒跡地面積6.01 km2;提取災(zāi)后(4月9日)火燒跡地斑塊24個(gè),識(shí)別出火燒跡地面積171.73 km2。

    3.3 精度評(píng)定

    采用災(zāi)后(4月14日)的Landsat-8的重采樣為10 m分辨率影像作為精度評(píng)定數(shù)據(jù)。通過(guò)計(jì)算災(zāi)后歸一化差分植被指數(shù),并通過(guò)人工目視解譯獲取受損植被面積,以此作為實(shí)際結(jié)果進(jìn)行植被變化檢測(cè)的精度評(píng)定;同時(shí)通過(guò)多次目視解譯獲得火燒跡地典型斑塊,并與dNBR+火燒跡地面積追蹤識(shí)別提取方法結(jié)果進(jìn)行對(duì)比評(píng)價(jià),在研究區(qū)域范圍隨機(jī)生成驗(yàn)證像元點(diǎn),結(jié)合Landsat-8影像數(shù)據(jù)確定驗(yàn)證像元點(diǎn)的屬性值構(gòu)建混淆矩陣,計(jì)算火燒跡地面積提取精度的kappa系數(shù)。因此,采用的NDVI+IR-MAD檢測(cè)植被變化方法,NDVI+IR-MAD檢測(cè)植被變化像元數(shù)1 464 408個(gè),實(shí)際像元數(shù)1 592 631個(gè),檢測(cè)到受損植被精度達(dá)到91.95%;dNBR+火燒跡地追蹤識(shí)別提取火燒跡地面積171.73 km2,火燒基地實(shí)際面積為189.19 km2,提取精度90.77%,kappa系數(shù)為0.87。評(píng)價(jià)結(jié)果表明該方法具有較好的效果與精度,驗(yàn)證了方法的可靠性。

    4 結(jié)論

    由于西南林區(qū)的高森林覆蓋率以及冬春季節(jié)干燥風(fēng)大,每年都會(huì)發(fā)生如木里這樣的森林火災(zāi),木里火災(zāi)是四川乃至西南林區(qū)近年來(lái)經(jīng)常發(fā)生的典型森林火災(zāi)。本文采用Sentinel-2遙感影像,以2020年3月28日四川涼山州木里縣發(fā)生的典型西南林區(qū)森林大火為例,進(jìn)行植被變化與火燒跡地提取的森林火災(zāi)反演實(shí)驗(yàn)。得出以下結(jié)論:該方法是利用所有森林火災(zāi)后都會(huì)產(chǎn)生的植被變化與火燒跡地信息來(lái)反演火災(zāi)情況,具有普適性。通過(guò)對(duì)四川木里火災(zāi)的反演,較好地證明了該方法的可行性。該方法流程能夠準(zhǔn)確分離出受損植被的區(qū)域與快速準(zhǔn)確提取計(jì)算火燒跡地像元面積,有效排除其他噪聲的干擾,且精度較高;能夠?qū)λ谢馂?zāi)區(qū)域進(jìn)行災(zāi)后快速評(píng)估,提供植被變化與火燒跡地的參考數(shù)據(jù),反演整個(gè)森林火災(zāi)的時(shí)空發(fā)展過(guò)程。根據(jù)木里地區(qū)的良好的遙感反演情況,也證明了該方法能夠判斷出火災(zāi)的發(fā)展勢(shì)態(tài),為森林火災(zāi)的安全撲救提供現(xiàn)實(shí)參考。

    猜你喜歡
    跡地火燒植被指數(shù)
    封二、封三:寒蔥嶺戰(zhàn)跡地 (伏擊戰(zhàn)) 新考
    火燒毒品
    南方周末(2019-06-27)2019-06-27 19:47:54
    AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
    我國(guó)科學(xué)家在國(guó)際上首發(fā)30 m分辨率全球火燒跡地?cái)?shù)據(jù)產(chǎn)品
    驢肉火燒
    火燒跡地造林技術(shù)
    基于Landsat影像提取大興安嶺區(qū)域火燒跡地★
    山西建筑(2017年15期)2017-07-01 22:03:22
    不怕火燒的手帕
    莫愁(2017年15期)2017-05-11 01:21:55
    “小火燒”成為大產(chǎn)業(yè)
    麻豆精品久久久久久蜜桃| 人妻 亚洲 视频| 免费大片黄手机在线观看| 亚洲av中文av极速乱| 国产成人免费观看mmmm| 一级毛片我不卡| 美女视频免费永久观看网站| 亚洲国产精品999| a级毛片在线看网站| 国产免费现黄频在线看| 久久精品国产亚洲av涩爱| 国产精品人妻久久久久久| 国产精品国产av在线观看| 人妻 亚洲 视频| 日韩精品免费视频一区二区三区 | 国产av一区二区精品久久| 亚洲丝袜综合中文字幕| 国产精品人妻久久久久久| 久久人妻熟女aⅴ| 九草在线视频观看| 色94色欧美一区二区| 久久久国产精品麻豆| 老司机影院毛片| 亚洲成人av在线免费| 精品国产乱码久久久久久小说| 天堂8中文在线网| 久久精品夜色国产| 日韩成人伦理影院| 免费人妻精品一区二区三区视频| av一本久久久久| 狂野欧美激情性xxxx在线观看| 波多野结衣一区麻豆| 国产亚洲欧美精品永久| 黄色视频在线播放观看不卡| 午夜免费男女啪啪视频观看| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 三级国产精品片| 赤兔流量卡办理| 一级毛片黄色毛片免费观看视频| 天堂中文最新版在线下载| 国产福利在线免费观看视频| 中国三级夫妇交换| 久久久久精品性色| 狠狠精品人妻久久久久久综合| 亚洲成av片中文字幕在线观看 | 国产精品成人在线| 久久热在线av| 中文字幕制服av| 久久鲁丝午夜福利片| 中文欧美无线码| 全区人妻精品视频| 国产成人精品婷婷| 少妇 在线观看| 菩萨蛮人人尽说江南好唐韦庄| 十分钟在线观看高清视频www| 日韩制服骚丝袜av| 桃花免费在线播放| 爱豆传媒免费全集在线观看| 国产精品欧美亚洲77777| 成人18禁高潮啪啪吃奶动态图| av天堂久久9| 久久午夜福利片| 免费日韩欧美在线观看| 美女国产高潮福利片在线看| 欧美成人精品欧美一级黄| 国产成人精品久久久久久| 一级,二级,三级黄色视频| 免费人妻精品一区二区三区视频| 亚洲国产av新网站| 美国免费a级毛片| 国产乱人偷精品视频| 大香蕉久久网| 草草在线视频免费看| 91精品伊人久久大香线蕉| 九九爱精品视频在线观看| 夫妻性生交免费视频一级片| 大码成人一级视频| 在现免费观看毛片| 久久这里有精品视频免费| 成人18禁高潮啪啪吃奶动态图| 国产成人av激情在线播放| 日本黄色日本黄色录像| 青春草视频在线免费观看| 九草在线视频观看| 亚洲成人av在线免费| 亚洲av日韩在线播放| 欧美日韩成人在线一区二区| 久久精品久久久久久噜噜老黄| 国产精品人妻久久久影院| 亚洲少妇的诱惑av| 一级爰片在线观看| 97超碰精品成人国产| 90打野战视频偷拍视频| 亚洲欧美中文字幕日韩二区| 国产深夜福利视频在线观看| 国产av码专区亚洲av| 成人国产av品久久久| 日韩不卡一区二区三区视频在线| 欧美日本中文国产一区发布| 观看av在线不卡| 色94色欧美一区二区| 搡女人真爽免费视频火全软件| 国产成人av激情在线播放| 婷婷色av中文字幕| 插逼视频在线观看| 精品久久国产蜜桃| 在线观看www视频免费| 日韩,欧美,国产一区二区三区| 国产一区二区在线观看日韩| 久久国产亚洲av麻豆专区| 久久久久久久国产电影| 精品国产露脸久久av麻豆| 波多野结衣一区麻豆| 看非洲黑人一级黄片| 国产欧美另类精品又又久久亚洲欧美| 欧美性感艳星| 中文字幕制服av| 精品国产一区二区三区久久久樱花| 日韩制服骚丝袜av| 亚洲色图综合在线观看| 国产精品一区www在线观看| 男女午夜视频在线观看 | 成人无遮挡网站| videos熟女内射| 肉色欧美久久久久久久蜜桃| 久久影院123| 午夜精品国产一区二区电影| 高清黄色对白视频在线免费看| 国产片内射在线| 丰满饥渴人妻一区二区三| 精品人妻熟女毛片av久久网站| 国产一区二区在线观看av| 中文字幕精品免费在线观看视频 | 日韩成人av中文字幕在线观看| 久久人人97超碰香蕉20202| 肉色欧美久久久久久久蜜桃| 国产成人av激情在线播放| 日韩制服丝袜自拍偷拍| 老司机影院毛片| 久久人人爽人人爽人人片va| 激情视频va一区二区三区| 九草在线视频观看| 国产精品欧美亚洲77777| 国产无遮挡羞羞视频在线观看| 久久久国产精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 九色亚洲精品在线播放| 性高湖久久久久久久久免费观看| 色吧在线观看| 少妇熟女欧美另类| 九九在线视频观看精品| 久久热在线av| 国产又爽黄色视频| 中文字幕人妻熟女乱码| 国产精品久久久久久精品电影小说| 18禁国产床啪视频网站| 我的女老师完整版在线观看| 超碰97精品在线观看| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图| 熟妇人妻不卡中文字幕| 亚洲精品日韩在线中文字幕| 美女福利国产在线| 成人午夜精彩视频在线观看| 草草在线视频免费看| 男人爽女人下面视频在线观看| 两个人免费观看高清视频| 国产一区二区三区av在线| 黄色毛片三级朝国网站| 又黄又爽又刺激的免费视频.| 伊人亚洲综合成人网| 中国国产av一级| 18禁动态无遮挡网站| 精品国产露脸久久av麻豆| 黑人猛操日本美女一级片| 午夜精品国产一区二区电影| freevideosex欧美| av卡一久久| 亚洲欧美色中文字幕在线| 免费观看av网站的网址| 亚洲av成人精品一二三区| av在线播放精品| 国产日韩欧美在线精品| 国产一区二区在线观看日韩| a级毛片黄视频| 亚洲久久久国产精品| 亚洲综合色惰| 日本黄大片高清| 青春草亚洲视频在线观看| 九九爱精品视频在线观看| 交换朋友夫妻互换小说| 热re99久久国产66热| 国产高清三级在线| 看非洲黑人一级黄片| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看 | 蜜桃在线观看..| 精品一区二区三区四区五区乱码 | 寂寞人妻少妇视频99o| 久久综合国产亚洲精品| 成人综合一区亚洲| 少妇人妻 视频| 97人妻天天添夜夜摸| 免费大片黄手机在线观看| 午夜激情久久久久久久| 18禁观看日本| 亚洲av中文av极速乱| 欧美日韩视频高清一区二区三区二| 亚洲欧美中文字幕日韩二区| 毛片一级片免费看久久久久| 黄色视频在线播放观看不卡| 十八禁高潮呻吟视频| 国产在线视频一区二区| 黄片无遮挡物在线观看| 亚洲欧洲日产国产| 亚洲精品日本国产第一区| 日韩视频在线欧美| 中文字幕另类日韩欧美亚洲嫩草| 亚洲美女视频黄频| 免费观看性生交大片5| 少妇熟女欧美另类| 久久狼人影院| 午夜福利在线观看免费完整高清在| 99热网站在线观看| 青春草国产在线视频| 中文字幕另类日韩欧美亚洲嫩草| 卡戴珊不雅视频在线播放| 亚洲国产毛片av蜜桃av| 亚洲一码二码三码区别大吗| 高清av免费在线| 最近2019中文字幕mv第一页| 国产色爽女视频免费观看| 老司机影院毛片| 日韩av免费高清视频| 婷婷色av中文字幕| 最新的欧美精品一区二区| 咕卡用的链子| 免费看光身美女| 免费久久久久久久精品成人欧美视频 | 成人二区视频| 汤姆久久久久久久影院中文字幕| 国产亚洲精品久久久com| 黄色配什么色好看| 在线观看美女被高潮喷水网站| 日韩视频在线欧美| 七月丁香在线播放| 久久人人97超碰香蕉20202| 最新的欧美精品一区二区| 久久国内精品自在自线图片| 90打野战视频偷拍视频| 久久综合国产亚洲精品| 国产午夜精品一二区理论片| 亚洲内射少妇av| 蜜桃国产av成人99| 国产黄色免费在线视频| 黄色视频在线播放观看不卡| 寂寞人妻少妇视频99o| 在线观看免费高清a一片| 亚洲欧美成人精品一区二区| 午夜久久久在线观看| 老熟女久久久| 日韩中字成人| 侵犯人妻中文字幕一二三四区| 大陆偷拍与自拍| 日本爱情动作片www.在线观看| 精品视频人人做人人爽| 五月玫瑰六月丁香| 男女午夜视频在线观看 | 亚洲国产毛片av蜜桃av| 日日撸夜夜添| 久久精品国产亚洲av天美| av播播在线观看一区| 下体分泌物呈黄色| 国产成人精品无人区| 美女福利国产在线| 少妇高潮的动态图| 精品第一国产精品| 宅男免费午夜| 国产亚洲午夜精品一区二区久久| 夜夜爽夜夜爽视频| 大码成人一级视频| 国产精品久久久av美女十八| av免费在线看不卡| 国产亚洲精品第一综合不卡 | 国产一区有黄有色的免费视频| 26uuu在线亚洲综合色| a级毛色黄片| 母亲3免费完整高清在线观看 | 国产亚洲一区二区精品| 亚洲av中文av极速乱| 老司机亚洲免费影院| av有码第一页| 国精品久久久久久国模美| 尾随美女入室| 久久午夜综合久久蜜桃| 男女无遮挡免费网站观看| 桃花免费在线播放| 丁香六月天网| 观看av在线不卡| 日本爱情动作片www.在线观看| 国产成人精品婷婷| 欧美精品亚洲一区二区| 高清视频免费观看一区二区| 黄片无遮挡物在线观看| 最后的刺客免费高清国语| a级毛色黄片| 久久婷婷青草| 成人手机av| 午夜av观看不卡| 日韩大片免费观看网站| 亚洲欧洲精品一区二区精品久久久 | 成人手机av| 男女边吃奶边做爰视频| 一级黄片播放器| 国产一区二区激情短视频 | av国产久精品久网站免费入址| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 80岁老熟妇乱子伦牲交| 91久久精品国产一区二区三区| 日韩电影二区| 色哟哟·www| 国产色婷婷99| 最近最新中文字幕免费大全7| 女人被躁到高潮嗷嗷叫费观| 成人黄色视频免费在线看| 亚洲欧美日韩另类电影网站| av网站免费在线观看视频| 在现免费观看毛片| 亚洲av日韩在线播放| 美国免费a级毛片| 国产探花极品一区二区| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品电影小说| 午夜免费男女啪啪视频观看| 91午夜精品亚洲一区二区三区| 我要看黄色一级片免费的| 少妇高潮的动态图| 亚洲婷婷狠狠爱综合网| www.av在线官网国产| 亚洲精品久久午夜乱码| 国产精品人妻久久久久久| 一级片'在线观看视频| 国产色爽女视频免费观看| 国产一区二区在线观看av| 国产精品久久久久久久久免| 欧美人与性动交α欧美精品济南到 | 人妻人人澡人人爽人人| 精品国产露脸久久av麻豆| 黄色毛片三级朝国网站| 在线天堂中文资源库| 精品第一国产精品| 亚洲性久久影院| 中国国产av一级| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| a级片在线免费高清观看视频| 欧美国产精品一级二级三级| 午夜视频国产福利| 国产极品粉嫩免费观看在线| 亚洲精品av麻豆狂野| 中国三级夫妇交换| 夜夜骑夜夜射夜夜干| 国产黄色免费在线视频| a级毛色黄片| 免费人妻精品一区二区三区视频| 最近2019中文字幕mv第一页| 亚洲av中文av极速乱| 51国产日韩欧美| 伊人亚洲综合成人网| 曰老女人黄片| 永久网站在线| 国产精品99久久99久久久不卡 | 国产精品人妻久久久久久| 欧美另类一区| 久久久欧美国产精品| 97在线视频观看| 久热这里只有精品99| √禁漫天堂资源中文www| 久久狼人影院| 中文字幕最新亚洲高清| 高清视频免费观看一区二区| xxxhd国产人妻xxx| 亚洲成色77777| 综合色丁香网| 国产成人精品福利久久| 精品卡一卡二卡四卡免费| 伦理电影大哥的女人| 国产1区2区3区精品| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 黄色一级大片看看| 丰满少妇做爰视频| 国产精品久久久久久久电影| 大香蕉97超碰在线| 亚洲国产精品一区二区三区在线| 国产片特级美女逼逼视频| 成人毛片60女人毛片免费| 春色校园在线视频观看| 永久网站在线| 国产黄色视频一区二区在线观看| a级毛片在线看网站| 韩国av在线不卡| 婷婷色麻豆天堂久久| 麻豆精品久久久久久蜜桃| 如日韩欧美国产精品一区二区三区| 超碰97精品在线观看| 看十八女毛片水多多多| 国产免费视频播放在线视频| 日本免费在线观看一区| 老司机亚洲免费影院| 交换朋友夫妻互换小说| 99热6这里只有精品| 亚洲精品视频女| 51国产日韩欧美| 18在线观看网站| 欧美成人午夜精品| 哪个播放器可以免费观看大片| 精品99又大又爽又粗少妇毛片| av免费观看日本| 久久久久久久国产电影| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频在线观看免费| 国产激情久久老熟女| 一区二区日韩欧美中文字幕 | 久久99热6这里只有精品| 国产av码专区亚洲av| 久久人人爽av亚洲精品天堂| 日本91视频免费播放| 日本av免费视频播放| 美女国产高潮福利片在线看| 天天躁夜夜躁狠狠久久av| 国产精品一二三区在线看| 黄片无遮挡物在线观看| 内地一区二区视频在线| 国产男人的电影天堂91| 精品一品国产午夜福利视频| 国产日韩一区二区三区精品不卡| 捣出白浆h1v1| 亚洲高清免费不卡视频| 边亲边吃奶的免费视频| 九色亚洲精品在线播放| 免费人妻精品一区二区三区视频| 内地一区二区视频在线| 成人国语在线视频| h视频一区二区三区| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 91在线精品国自产拍蜜月| 国产高清不卡午夜福利| 99久久精品国产国产毛片| 纵有疾风起免费观看全集完整版| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 内地一区二区视频在线| 我的女老师完整版在线观看| 亚洲精品久久成人aⅴ小说| 久久久精品94久久精品| 欧美人与善性xxx| 大片免费播放器 马上看| 狠狠婷婷综合久久久久久88av| 久久久久久人人人人人| 日本午夜av视频| 韩国精品一区二区三区 | 欧美变态另类bdsm刘玥| 国产成人av激情在线播放| 女性生殖器流出的白浆| av不卡在线播放| 边亲边吃奶的免费视频| 婷婷色av中文字幕| 精品一区二区免费观看| 精品人妻在线不人妻| 天堂8中文在线网| 久久精品人人爽人人爽视色| 老司机影院成人| 亚洲av国产av综合av卡| 欧美亚洲日本最大视频资源| 美女国产视频在线观看| 日韩一区二区视频免费看| 美女xxoo啪啪120秒动态图| 亚洲中文av在线| 2021少妇久久久久久久久久久| 亚洲av欧美aⅴ国产| 激情视频va一区二区三区| 亚洲美女视频黄频| 男女无遮挡免费网站观看| 黄片播放在线免费| 黑丝袜美女国产一区| 欧美精品亚洲一区二区| 成人黄色视频免费在线看| 精品亚洲成国产av| 夫妻午夜视频| 五月开心婷婷网| av不卡在线播放| 国产精品偷伦视频观看了| 亚洲美女搞黄在线观看| 国产精品三级大全| 亚洲国产日韩一区二区| 国产精品嫩草影院av在线观看| 2022亚洲国产成人精品| 亚洲精品一区蜜桃| 欧美国产精品一级二级三级| 免费大片黄手机在线观看| 婷婷色综合www| 国产精品国产三级专区第一集| 久久久久久久国产电影| 日本猛色少妇xxxxx猛交久久| 欧美激情国产日韩精品一区| 性色avwww在线观看| 日韩av不卡免费在线播放| 精品99又大又爽又粗少妇毛片| 国产一区亚洲一区在线观看| 尾随美女入室| 丰满少妇做爰视频| 汤姆久久久久久久影院中文字幕| 韩国精品一区二区三区 | 久久久久国产网址| 春色校园在线视频观看| 久久99精品国语久久久| 高清在线视频一区二区三区| 一级a做视频免费观看| 大香蕉久久网| a级毛片黄视频| 99香蕉大伊视频| 亚洲国产看品久久| 看免费av毛片| 纵有疾风起免费观看全集完整版| 久热久热在线精品观看| 日韩中文字幕视频在线看片| 高清视频免费观看一区二区| 婷婷色综合www| 欧美日韩成人在线一区二区| 在线观看人妻少妇| 天堂8中文在线网| 国产有黄有色有爽视频| 黄网站色视频无遮挡免费观看| 菩萨蛮人人尽说江南好唐韦庄| 一区二区日韩欧美中文字幕 | 久久精品国产亚洲av涩爱| 最近手机中文字幕大全| 国产精品一国产av| 亚洲第一av免费看| 一区二区三区四区激情视频| 天堂中文最新版在线下载| 国产精品国产三级国产专区5o| 亚洲av男天堂| 中国美白少妇内射xxxbb| 丰满少妇做爰视频| 视频在线观看一区二区三区| 亚洲少妇的诱惑av| 精品一区在线观看国产| 精品国产露脸久久av麻豆| 麻豆精品久久久久久蜜桃| 亚洲av男天堂| 伦精品一区二区三区| 国产精品久久久久成人av| 国产精品成人在线| videos熟女内射| 综合色丁香网| 街头女战士在线观看网站| 韩国av在线不卡| 亚洲精品日韩在线中文字幕| 亚洲精品日本国产第一区| 天天躁夜夜躁狠狠躁躁| 精品熟女少妇av免费看| 丝袜美足系列| 久久精品人人爽人人爽视色| 欧美精品亚洲一区二区| 亚洲国产精品专区欧美| 久久久欧美国产精品| 最近手机中文字幕大全| 观看av在线不卡| 成年女人在线观看亚洲视频| 久久 成人 亚洲| a级毛片在线看网站| 十八禁高潮呻吟视频| 色哟哟·www| 久久国产精品大桥未久av| 人妻系列 视频| 99久久中文字幕三级久久日本| 午夜福利,免费看| 亚洲四区av| 一区在线观看完整版| 午夜免费观看性视频| 卡戴珊不雅视频在线播放| 亚洲精品久久久久久婷婷小说| 91成人精品电影| 日本爱情动作片www.在线观看| 纵有疾风起免费观看全集完整版| av线在线观看网站| 亚洲国产精品一区三区| 国产成人精品久久久久久| 纯流量卡能插随身wifi吗| 丰满乱子伦码专区| 大陆偷拍与自拍| 一边摸一边做爽爽视频免费| 国产成人91sexporn| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国语对白做爰xxxⅹ性视频网站| 成人手机av| 青春草视频在线免费观看| 18禁在线无遮挡免费观看视频| 欧美日韩成人在线一区二区| 精品国产一区二区三区四区第35| 久久97久久精品| 亚洲精品自拍成人| 免费黄频网站在线观看国产| 色视频在线一区二区三区| 欧美精品一区二区大全| 国产在线一区二区三区精| 亚洲情色 制服丝袜| 国产免费一级a男人的天堂| 日韩制服丝袜自拍偷拍| 青春草亚洲视频在线观看|