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

    基于QUEST決策樹的火燒跡地自動(dòng)提取

    2018-09-13 10:48:30岳彩榮
    四川林業(yè)科技 2018年4期
    關(guān)鍵詞:跡地林火火燒

    霍 鵬,岳彩榮

    (西南林業(yè)大學(xué),云南 昆明 650224)

    森林火災(zāi)是影響森林生態(tài)系統(tǒng)的一個(gè)重要的擾動(dòng)因素[1]。林火的發(fā)生不僅會(huì)直接導(dǎo)致森林面積減少、森林結(jié)構(gòu)和生態(tài)系統(tǒng)平衡被破壞,而且樹木燃燒的過程中向大氣中釋放大量的CO2以及其他有害氣體對(duì)氣候產(chǎn)生了深遠(yuǎn)的影響,也給當(dāng)?shù)厣鐣?huì)經(jīng)濟(jì)帶來了嚴(yán)重的損失。因此,森林火災(zāi)已成為森林資源和自然環(huán)境破壞的重大威脅,其發(fā)生原因復(fù)雜且難以預(yù)防。如何對(duì)森林火災(zāi)及其遭受的損失進(jìn)行及時(shí)、準(zhǔn)確和有效地監(jiān)測(cè)和事后損失評(píng)估,對(duì)預(yù)防火災(zāi)發(fā)生、減少災(zāi)害損失、維護(hù)生態(tài)平衡和保護(hù)森林資源都具有非常重要的現(xiàn)實(shí)意義[2]。

    近年來,隨著遙感技術(shù)的逐步成熟,利用其空間分辨率高、監(jiān)測(cè)范圍廣、時(shí)效性強(qiáng)、成本較低的特點(diǎn),將森林火災(zāi)受損狀況監(jiān)測(cè)與遙感技術(shù)結(jié)合起來,已成為我國(guó)森林火災(zāi)監(jiān)測(cè)的一個(gè)重要手段。決策樹算法是一種數(shù)據(jù)挖掘監(jiān)督分類方法,通過對(duì)訓(xùn)練樣本集進(jìn)行數(shù)據(jù)挖掘和知識(shí)發(fā)現(xiàn)建立決策樹,利用規(guī)則對(duì)遙感影像進(jìn)行分類[3]。不僅直觀簡(jiǎn)單、運(yùn)算速度快、分類精度高,而且能自動(dòng)選取特征并融入影像以外的各種知識(shí)、能高效地處理大量數(shù)據(jù)和高維數(shù)據(jù)。傳統(tǒng)的決策樹分類在決策樹的建立、特征指數(shù)的選擇和閾值的設(shè)定上都較為煩瑣,即使有諸如CART決策樹能自動(dòng)分類的決策樹,但在變量選擇時(shí)都是有偏的[4]。

    QUEST算法是Loh和Shih于1997年提出的,該算法將分支變量和分割點(diǎn)的選擇分開進(jìn)行,使其具有兼顧適用于連續(xù)型的變量和離散型的變量的特點(diǎn)。此外,它還克服了其他決策樹算法更傾向于選擇那些具有更多的潛在分割點(diǎn)的預(yù)測(cè)變量[5~8]。

    CRUISE 2D算法是Kim和Loh于2001年提出的,它由FACT 算法改進(jìn)而得來,它融合了FACT、QUEST 和CART 算法的思想。CRUISE 2D算法通過多分割減少樹的深度,并采用單因子或雙因子交互方法選擇預(yù)測(cè)變量,同時(shí)提高了樹的可讀性[9-10]。

    這兩種決策樹分類方法皆可以快速、無偏、高效提取具有復(fù)雜地形的遙感影像的不同土地類型,與普通決策樹相比具有方便、快捷、分類精度高的特點(diǎn)[11]。在特征指數(shù)的適用性上也更趨近于較小的研究范圍,一般僅在火災(zāi)發(fā)生地的局部區(qū)域提取火燒跡地,對(duì)大尺度范圍內(nèi)的火燒跡地提取較少。根據(jù)指數(shù)適用的尺度范圍、林火發(fā)生后植被減少、燃燒殘留物沉積以及地形、紋理特征等,在已提出的修正植被指數(shù)、火燒跡地指數(shù)等指數(shù)中選取最佳指數(shù),為QUEST決策樹的建立提供堅(jiān)實(shí)的基礎(chǔ)。

    針對(duì)2006年云南省昆明市的“3·29”林火事件,由于該燃燒區(qū)從安寧市溫泉鎮(zhèn)西北部的古朗大箐下段森林蔓延至西山區(qū),因此本文采用QUEST算法、CRUISE 2D算法的決策樹以及最大似然法在大尺度范圍內(nèi)提取Landsat影像中的火燒跡地,并估測(cè)每種分類方法的精度,并選擇出最適合大尺度范圍內(nèi)火燒跡地提取的方法。結(jié)合2015年云南省昆明市的“3·02”林火事件驗(yàn)證該方法的普適性。這為在中低分辨率的遙感影像中林火損失面積的自動(dòng)提取、長(zhǎng)時(shí)間大尺度范圍內(nèi)的森林火災(zāi)監(jiān)測(cè)、林火發(fā)生的時(shí)空規(guī)律、林火的預(yù)防等提供數(shù)據(jù)準(zhǔn)備和參考。

    1 研究區(qū)和研究數(shù)據(jù)

    1.1 研究區(qū)概況

    本文實(shí)驗(yàn)研究區(qū)為云南省昆明市,它地處我國(guó)云貴高原中部,地理位置介于東經(jīng)102°10′~103°40′,北緯24°23′~26°32′;平均海拔在 1 500 m~2 800 m之間,總體地勢(shì)北高南低,中部隆起,東西兩側(cè)較低(見圖1)。

    昆明市屬北緯低緯度亞熱帶高原山地季風(fēng)氣候,境內(nèi)大多數(shù)地區(qū)夏無酷暑,冬無嚴(yán)寒,四季如春,氣候宜人,素以“春城”之稱而享譽(yù)中外。按全年降水量分布分為干、濕兩季,每年的5月至10月為雨季,11月至次年的4月為干季。因整個(gè)干季降水量極少僅占全年降水量的3%~5%,每逢3月、4月份蒸發(fā)旺盛致使春季為林火多發(fā)季節(jié),昆明市曾多次發(fā)生森林火災(zāi)事件,因此選該地區(qū)為本實(shí)驗(yàn)的研究區(qū)。

    1.2 研究數(shù)據(jù)

    本文的研究數(shù)據(jù)主要為美國(guó)陸地衛(wèi)星(Landsat)影像和DEM數(shù)據(jù),所用數(shù)據(jù)皆取自于地理空間數(shù)據(jù)云平臺(tái)(http://www.gscloud.cn),其中Landsat數(shù)據(jù)產(chǎn)品為L(zhǎng)evel1T的地形校正數(shù)據(jù),其投影坐標(biāo)為UTM/WGS84,空間分辨率為30 m;DEM數(shù)據(jù)是2009年采集的同為UTM/WGS84投影坐標(biāo)的數(shù)字高程模型,其空間分辨率為30 m。

    2 研究方法

    2.1 數(shù)據(jù)預(yù)處理

    實(shí)驗(yàn)數(shù)據(jù)選用是2005年和2006年云南省昆明市的Landsat TM數(shù)據(jù),驗(yàn)證數(shù)據(jù)采用的是2015年1月和3月的Landsat8 OLI數(shù)據(jù)。由于從地理空間數(shù)據(jù)云獲取的數(shù)據(jù)為L(zhǎng)evel1T級(jí)地形校正數(shù)據(jù),已經(jīng)過系統(tǒng)輻射校正和地面控制點(diǎn)幾何校正,因此只需對(duì)研究區(qū)域做大氣校正、影像鑲嵌和裁切。

    2.2 Cart決策樹火燒跡地信息提取

    2.2.1 植被指數(shù)

    針對(duì)大范圍內(nèi)的火燒跡地提取,對(duì)其他地物的分離度也要求較高??紤]到電磁波輻射、大氣、植被覆蓋土壤背景等因素的影響,本文主要對(duì)以下幾種植被指數(shù)進(jìn)行研究,包括歸一化差值植被指數(shù)(Normalized Difference Vegetation,NDVI)、多時(shí)相差值植被指數(shù)(Multi-temporal Normalized Difference Vegetation,MNDVI)、土壤調(diào)節(jié)植被指數(shù)(Soil Adjusted Vegetation Index,SAVI)、修改型土壤調(diào)節(jié)植被指數(shù)(Modified Soil Adjusted Vegetation Index,MSAVI)、全球環(huán)境監(jiān)測(cè)指數(shù)(Global Environment Monitoring Index,GEMI)以及亮度植被指數(shù)(Bright Vegetation Index,BVI)、綠度植被指數(shù)(Green Vegetation Index,GVI)、濕度植被指數(shù)(Wetness Vegetation Index,WVI)。

    經(jīng)過對(duì)不同指數(shù)進(jìn)行比較發(fā)現(xiàn),NDVI對(duì)綠色植被表現(xiàn)最為敏感,即使在低密度的植被區(qū)域仍能很好地提取出來[12]。MNDVI是將同一區(qū)域的兩個(gè)不同時(shí)期的NDVI進(jìn)行波段相減。對(duì)于小范圍內(nèi)的發(fā)生過火災(zāi)的區(qū)域,能很好的將植被和其他地類區(qū)分開來;但對(duì)于大范圍的地物提取,即使整景圖像中地類變化較小,但也可能會(huì)受到影像拍攝時(shí)間不同以及陰影等其他因素的影響,降低了整個(gè)區(qū)域中植被的提取精度。幾種植被指數(shù)中僅有MSAVI不受云層影響。MSAVI對(duì)火燒跡地區(qū)域的區(qū)分較為明顯,但是植被和城鎮(zhèn)會(huì)難以區(qū)分,且受陰影影響較大,而SAVI能將植被和城鎮(zhèn)較好的區(qū)分開,總體受陰影的影響較小。GEMI與NDVI相比,既保存了比NDVI指數(shù)相對(duì)低密度至濃密度覆蓋更大的動(dòng)態(tài)范圍,又在不用改變植被信息的前提下減小了大氣影響,但仍受到裸土的亮度和顏色相當(dāng)大的影響,對(duì)于稀疏或中密度植被覆蓋的區(qū)域?qū)χ脖坏奶崛⌒Ч^差[13]。而根據(jù)經(jīng)驗(yàn)方法得到土壤亮度指數(shù)(BVI)、綠度植被指數(shù)(GVI)、黃度植被指數(shù)(YVI),雖然忽略了大氣、土壤、植被間相互作用,但就不同地類區(qū)分度來說,BVI和WVI植被和水體嚴(yán)重混分,GVI與GEMI區(qū)分效果相似,但受傳感器的影響因素較大。

    2.2.2 水體指數(shù)

    水體指數(shù)是同植被指數(shù)相似的在提取水體領(lǐng)域的一種指數(shù)。常見的水體指數(shù)主要包括歸一化差值水體指數(shù)(Normalized Difference Water Index,NDWI)、修改型歸一化差異水體指數(shù)(Modified Normalized Difference Water Index,MNDWI)、增強(qiáng)型水體指數(shù)(Enhanced Water Index,EWI)、新型水體指數(shù)(New Water Index,NWI)。

    這四種水體指數(shù)中NDWI、EWI受云的影響較大。本研究采用影像大約為4月、5月左右,滇池的湖水邊緣生長(zhǎng)了許多漂浮的水生植物,而EWI、NWI對(duì)水體中的植被比較敏感,在提取水體時(shí)對(duì)水質(zhì)要求較高。MNDWI在幾個(gè)植被指數(shù)中提取水體的效果最好,受云、城市影響較小,不受水質(zhì)的限制,但對(duì)陰影以及覆有大棚的農(nóng)田影響仍難以剔除[14]。除以上的水體指數(shù)以外,纓帽變換中的濕度植被指數(shù)對(duì)面積較大的水體提取較好,但對(duì)于大范圍中的水體提取,效果并不理想易與城市等其他地物發(fā)生混淆。

    2.2.3 火燒跡地指數(shù)

    用于林火跡地提取的遙感指數(shù)主要有歸一化燃燒率(Normalized Burn Ratio,NBR)、差值歸一化燃燒率(Difference Normalized Burn Ratio,dNBR)、燃燒面積指數(shù)(Burn Area Index,BAI)、碳土壤指數(shù)(The Char Soil Index,CSI)。對(duì)于顏色較深的水體,這幾種火燒跡地指數(shù)都難以將之與火燒跡地區(qū)分開來。dNBR比NBR受其他地物影響較少,區(qū)分度較高[15]。

    2.2.4 特征融合

    將全球環(huán)境監(jiān)測(cè)指數(shù)(GEMI)、修改型歸一化差異水體指數(shù)(MNDWI)差值歸一化燃燒率(dNBR)、灰度共生矩陣的對(duì)比度(contrast)以及數(shù)字高程數(shù)據(jù)(DEM)同Landsat 的多光譜影像做波段融合。

    2.3 基于QUEST算法的分類

    本實(shí)驗(yàn)將研究區(qū)分為6種主要土地類型,包括林地、農(nóng)地、水體、城市、火燒跡地和裸地,并采用目視解譯的方法,在特征融合后的影像上選取訓(xùn)練樣本(550個(gè)感興趣區(qū),總計(jì) 82 836個(gè)像元)。目視解譯是在解譯者具備影像解譯專業(yè)知識(shí)、區(qū)域背景知識(shí)及遙感系統(tǒng)知識(shí)的基礎(chǔ)上,根據(jù)圖像基本解譯要素和解譯標(biāo)志來識(shí)別目標(biāo)地物。訓(xùn)練樣本的選取是構(gòu)建決策樹的關(guān)鍵部分,它的好壞直接影響后期分類精度。因此,在選取訓(xùn)練樣本時(shí),考慮了其代表性、典型性及完備性,也就是說訓(xùn)練樣本代表了整個(gè)區(qū)域內(nèi)不同地類間的光譜特征差異。

    3 研究結(jié)果與分析

    3.1 火燒跡地提取精度評(píng)價(jià)

    為了驗(yàn)證火燒跡地自動(dòng)識(shí)別的精度,根據(jù)野外采集數(shù)據(jù)和目視解譯隨機(jī)選取550個(gè)驗(yàn)證樣本點(diǎn),分別計(jì)算QUEST決策樹、CRUISE決策樹和極大似然分類方法的分類精度,建立混淆矩陣進(jìn)行精度評(píng)價(jià),如表1所示。結(jié)果表明,QUEST決策樹的分類總體精度為84.5%,Kappa系數(shù)為0.8161;CRUISE 2D決策樹的分類總體精度為83.5%,Kappa系數(shù)為0.8041;最大似然分類的分類總體精度為80.0%,Kappa系數(shù)為0.7624。分類總體精度是分類器正確分類的像元總數(shù)與總像元數(shù)之比,若總體精度越大,分類器分類精度越高;Kappa系數(shù)是綜合了生產(chǎn)精度和用戶精度的指標(biāo),若Kappa系數(shù)越大,分類結(jié)果與真實(shí)參考之間差異越小,則分類器分類精度越高。相比其他兩種分類方法,火燒跡地的生產(chǎn)者精度和用戶精度也是3種方法中最高的。因此,QUEST決策樹在大尺度范圍內(nèi)提取火燒跡地的效果最為理想(圖2和表1)。

    表1 QUEST決策樹、CRUISE決策樹和最大似然分類精度評(píng)價(jià)表

    分類方法地類用戶精度生產(chǎn)者精度Kappa系數(shù)總體精度QUEST決策樹森林61%99%0.816184.5%水體99%91%城市86%79%火燒跡地99%98%裸地86%93%農(nóng)地88%82%CRUISE決策樹森林62%99%0.804183.5%水體95%93%城市90%84%火燒跡地98%96%裸地87%95%農(nóng)地96%80%最大似然分類森林63%99%0.762480.0%水體99%88%城市97%62%火燒跡地86%98%裸地69%96%農(nóng)地90%83%

    3.2 火燒跡地提取精度分析

    3種分類方法中QUEST和CRUISE 2D分類精度相差較小,但在提取火燒跡地的面積上存在很大的差異。僅針對(duì)“3·29”事件發(fā)生區(qū)域,利用火燒跡地提取精度驗(yàn)證公式(如公式(1)和公式(2))計(jì)算3種方法的火燒跡地提取精度。其中,QUEST決策樹的火燒跡地提取精度為97.9%;CRUISE 2D決策樹的火燒跡地提取精度為80.5%;最大似然分類的火燒跡地提取精度為64.8%。QUEST算法在提取局部火燒跡地與實(shí)際面積最為吻合,如圖3(a)所示。

    圖2 分類結(jié)果

    識(shí)別出的火燒跡地面積<實(shí)際面積時(shí):

    識(shí)別精度=識(shí)別面積/實(shí)際面積×100%

    (1)

    識(shí)別出的火燒跡地面積>實(shí)際面積時(shí):

    (2)

    但是QUEST決策樹在大尺度范圍內(nèi)對(duì)火燒跡地提取的效果并不理想(如圖4),在分類時(shí)會(huì)產(chǎn)生與其他地類混分的現(xiàn)象,如火燒跡地與裸地混分較為嚴(yán)重,特別是因水土流失后的裸土,與林火發(fā)生后裸露的土壤區(qū)域在目視解譯上顏色較為相近;春季農(nóng)地焚燒的區(qū)域也會(huì)被提取出來,判斷是否是林火跡地還需有相關(guān)地類調(diào)查信息,如圖4中的(a)和(b)的紅框?yàn)橐陨蟽煞N情況誤分的區(qū)域。此外,這3種方法在提取火燒跡地的過程中皆會(huì)產(chǎn)生通過分類后處理無法去除的大量的碎斑,這些碎斑需根據(jù)國(guó)家規(guī)定的森林火災(zāi)等級(jí)的面積,采用手動(dòng)的方式對(duì)碎斑進(jìn)行剔除。與其相比,最大似然分類雖在總體分類精度上與QUEST決策樹相當(dāng),且混分的地類較少,但提取局部區(qū)域的火燒跡地面積與實(shí)際面積差距巨大。

    圖3 火燒跡地提取結(jié)果

    由于自然修復(fù)或人為種植會(huì)使得火燒跡地在一定時(shí)間范圍內(nèi)緩慢恢復(fù),雖無法將生物量恢復(fù)為原來的數(shù)值,但超過一定的時(shí)間限制會(huì)使得在目視解譯上難以辨認(rèn),這就造成了火燒跡地提取的時(shí)效性。因此,本實(shí)驗(yàn)方法僅可以識(shí)別出較短時(shí)間范圍內(nèi)火燒跡地的狀況,時(shí)間范圍的選取需根據(jù)研究區(qū)域的氣候環(huán)境以及林火災(zāi)后修復(fù)手段加以判斷。

    3.3 普適性驗(yàn)證

    本文選取2015年云南省昆明市的為驗(yàn)證區(qū)域采用同樣的方法對(duì)研究區(qū)進(jìn)行分類(如圖4),研究發(fā)現(xiàn)2015年的總體分類精度也達(dá)到97%,Kappa系數(shù)為0.957,提取火燒跡地的面積占實(shí)際精度的89%??梢婋m采用不同傳感器的同系列衛(wèi)星的產(chǎn)品,與Landsat TM影像相比的Landsat8 OLI影像多了一個(gè)波段,不僅能保證分類的精度,而且Landsat8 OLI數(shù)據(jù)的多種假彩色波段的組合更有利于通過目視解譯區(qū)分不同的地類。

    兩個(gè)研究區(qū)在時(shí)間選取上的差異性對(duì)整個(gè)影像的分類具有一定的影響。Landsat5 TM影像采用的是2006年4月初的數(shù)據(jù),而Landsat8 OLI采用的是2015年3月初的數(shù)據(jù)。雖僅存在一月之差,但對(duì)于中國(guó)西南部的氣候來說,進(jìn)入四月氣候轉(zhuǎn)暖快滇池等一些區(qū)域的水質(zhì)明顯下降,水生生物增多造成水體邊界模糊;又正值農(nóng)耕播種的季節(jié),3月初明顯未更耕種的農(nóng)地會(huì)被分成裸地,使得與火燒跡地間隙顏色接近的裸地被誤分的可能性增加,這使得在同一區(qū)域的不同時(shí)期的土地分類產(chǎn)生了差異。又因?yàn)闆Q策樹的精度對(duì)感興趣區(qū)的選取有較高的依賴性,因此這對(duì)目視解譯人員有較高的要求。

    圖4 大范圍火燒跡地提取圖

    4 討論

    本文采用QUEST算法、CRUISE算法的決策樹以及最大似然分類法對(duì)大范圍內(nèi)的火燒跡地進(jìn)行提取,經(jīng)實(shí)驗(yàn)驗(yàn)證QUEST算法的決策樹在三種分類方法中對(duì)分類精度最高、火燒跡地的提取的效果上最為理想。在實(shí)現(xiàn)大尺度范圍內(nèi)對(duì)火燒跡地自動(dòng)提取方面,QUEST決策樹雖在這三種方法中效果最好,但仍有許多待改進(jìn)之處。

    但本實(shí)驗(yàn)中的不足之處包括以下幾點(diǎn):根據(jù)森林火災(zāi)中火種的蔓延特點(diǎn),在提取火燒跡地的結(jié)構(gòu)中會(huì)有許多零散的誤分區(qū)域,還需要在分類后人工剔除,尚未實(shí)現(xiàn)火燒跡地提取的完全自動(dòng)化;僅利用地物自身的生物特性和整體的生長(zhǎng)特征,提取指標(biāo)單一,檢測(cè)結(jié)果易受干旱、植被季節(jié)變化等影響;能否可以通過增加其他特征開降低某些顏色的裸土與火燒跡地混分的狀況;為了降低決策樹的深度,本文并未采用過多指數(shù)的波段融合,僅選取最佳指數(shù),過少的特征是否會(huì)影響分類提取的精度;在決策樹分到多少級(jí)時(shí)分類效果最佳等問題還需要以后的研究中逐步驗(yàn)證。

    猜你喜歡
    跡地林火火燒
    無錫惠山區(qū)堅(jiān)持“六抓六強(qiáng)” 構(gòu)建林火防治銅墻鐵壁
    封二、封三:寒蔥嶺戰(zhàn)跡地 (伏擊戰(zhàn)) 新考
    林火監(jiān)測(cè)系統(tǒng)在森林防火中的應(yīng)用與發(fā)展
    半邊天
    火燒毒品
    南方周末(2019-06-27)2019-06-27 19:47:54
    我國(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
    黄色一级大片看看| 三级毛片av免费| 亚洲欧美激情综合另类| 亚洲欧美日韩高清在线视频| 熟女电影av网| 欧美bdsm另类| 国产精品女同一区二区软件 | 在线国产一区二区在线| 成人特级黄色片久久久久久久| 欧美bdsm另类| 天天躁日日操中文字幕| 内地一区二区视频在线| 亚洲最大成人中文| eeuss影院久久| 国产欧美日韩精品亚洲av| 97热精品久久久久久| 成熟少妇高潮喷水视频| h日本视频在线播放| 美女免费视频网站| 少妇猛男粗大的猛烈进出视频 | 久久九九热精品免费| 亚洲美女搞黄在线观看 | 色噜噜av男人的天堂激情| 深爱激情五月婷婷| 久久久国产成人精品二区| 在线国产一区二区在线| 深夜a级毛片| 乱人视频在线观看| 又爽又黄无遮挡网站| 欧美日本视频| 亚洲精品日韩av片在线观看| 狠狠狠狠99中文字幕| 成人欧美大片| 不卡一级毛片| 国产精品三级大全| aaaaa片日本免费| 日本一本二区三区精品| 悠悠久久av| 深夜精品福利| 麻豆国产av国片精品| 高清在线国产一区| 国产精品1区2区在线观看.| 亚洲精品成人久久久久久| 国产精品精品国产色婷婷| aaaaa片日本免费| 欧美高清成人免费视频www| 亚洲 国产 在线| 九色成人免费人妻av| 久久国内精品自在自线图片| 淫妇啪啪啪对白视频| 日本黄色视频三级网站网址| 天堂√8在线中文| 日本欧美国产在线视频| 国产在线男女| 亚洲精品一区av在线观看| 91久久精品国产一区二区三区| 国产淫片久久久久久久久| 日韩大尺度精品在线看网址| 欧美+亚洲+日韩+国产| 亚洲精品国产成人久久av| 国产成人aa在线观看| 啦啦啦韩国在线观看视频| 日本黄大片高清| 99热只有精品国产| 亚洲成人久久爱视频| 日日夜夜操网爽| 国产女主播在线喷水免费视频网站 | 午夜福利在线在线| 精品免费久久久久久久清纯| 欧美极品一区二区三区四区| 亚洲专区国产一区二区| 亚洲国产欧洲综合997久久,| 亚洲人与动物交配视频| 免费av观看视频| 国产精品久久久久久亚洲av鲁大| 国产精品日韩av在线免费观看| 乱系列少妇在线播放| 99热精品在线国产| 2021天堂中文幕一二区在线观| 中文字幕熟女人妻在线| 国产精品伦人一区二区| 亚洲av二区三区四区| 一进一出好大好爽视频| 99精品久久久久人妻精品| 免费人成在线观看视频色| 国产高潮美女av| 少妇的逼好多水| 午夜爱爱视频在线播放| 韩国av一区二区三区四区| 免费看a级黄色片| 舔av片在线| 亚洲成a人片在线一区二区| 亚洲,欧美,日韩| 亚洲精品日韩av片在线观看| 午夜老司机福利剧场| 亚洲天堂国产精品一区在线| 一进一出抽搐gif免费好疼| 国产精品人妻久久久影院| 国产亚洲精品久久久com| 淫妇啪啪啪对白视频| 国产爱豆传媒在线观看| 嫩草影院精品99| 两人在一起打扑克的视频| 亚洲va日本ⅴa欧美va伊人久久| eeuss影院久久| 亚洲精华国产精华精| 久久久久精品国产欧美久久久| 国产精品久久久久久久电影| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 国产白丝娇喘喷水9色精品| 黄片wwwwww| videossex国产| 国产又黄又爽又无遮挡在线| 麻豆久久精品国产亚洲av| 国产一区二区激情短视频| 国产中年淑女户外野战色| 韩国av一区二区三区四区| 精品久久久久久久久亚洲 | 国产高潮美女av| or卡值多少钱| 精品午夜福利视频在线观看一区| 欧美日韩黄片免| 亚洲精品粉嫩美女一区| 99热这里只有是精品在线观看| 成人永久免费在线观看视频| 免费av毛片视频| 男人的好看免费观看在线视频| 99精品在免费线老司机午夜| eeuss影院久久| www.www免费av| 成人美女网站在线观看视频| 国产精品无大码| 99久久无色码亚洲精品果冻| 国产激情偷乱视频一区二区| 久久精品国产自在天天线| 夜夜爽天天搞| 免费一级毛片在线播放高清视频| 亚洲美女黄片视频| 最后的刺客免费高清国语| 久久久精品大字幕| 亚洲精品日韩av片在线观看| 99riav亚洲国产免费| 免费观看精品视频网站| 97人妻精品一区二区三区麻豆| 少妇猛男粗大的猛烈进出视频 | 黄色丝袜av网址大全| 午夜福利在线观看免费完整高清在 | 中文字幕av成人在线电影| 日韩欧美免费精品| 久久国产精品人妻蜜桃| 亚洲自拍偷在线| 亚洲精品影视一区二区三区av| 国产在视频线在精品| 国产爱豆传媒在线观看| 亚洲五月天丁香| 99热精品在线国产| 日本 av在线| 日本一二三区视频观看| 国产成人一区二区在线| 在线免费观看的www视频| 可以在线观看的亚洲视频| 国产成人一区二区在线| 国产高清视频在线观看网站| 色av中文字幕| 午夜福利高清视频| 国产午夜福利久久久久久| 99在线人妻在线中文字幕| 国产精品久久久久久av不卡| 免费观看精品视频网站| 欧美又色又爽又黄视频| 少妇的逼好多水| 久久久久性生活片| 99视频精品全部免费 在线| 麻豆国产97在线/欧美| 久久久国产成人精品二区| 在线观看66精品国产| 哪里可以看免费的av片| 别揉我奶头 嗯啊视频| 一进一出抽搐gif免费好疼| 男女之事视频高清在线观看| 中文字幕久久专区| 亚洲精品日韩av片在线观看| 简卡轻食公司| 女同久久另类99精品国产91| 日韩欧美免费精品| 99热精品在线国产| 丰满的人妻完整版| 国产成人aa在线观看| 国产毛片a区久久久久| 五月玫瑰六月丁香| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 国产成人影院久久av| 小蜜桃在线观看免费完整版高清| 婷婷六月久久综合丁香| 别揉我奶头~嗯~啊~动态视频| 伦精品一区二区三区| 他把我摸到了高潮在线观看| 亚洲成a人片在线一区二区| 欧美区成人在线视频| 综合色av麻豆| 国产综合懂色| 一级a爱片免费观看的视频| 毛片女人毛片| 婷婷亚洲欧美| 亚洲熟妇熟女久久| 国产av一区在线观看免费| 男女啪啪激烈高潮av片| 桃红色精品国产亚洲av| 亚洲av五月六月丁香网| 日韩欧美免费精品| 99热精品在线国产| 国产精品,欧美在线| 最新在线观看一区二区三区| 1000部很黄的大片| 一夜夜www| 人妻丰满熟妇av一区二区三区| 精华霜和精华液先用哪个| 久久久色成人| 在线观看免费视频日本深夜| 日日干狠狠操夜夜爽| 欧美最黄视频在线播放免费| 亚洲av.av天堂| 久久99热6这里只有精品| 成人特级黄色片久久久久久久| 国内毛片毛片毛片毛片毛片| 免费人成在线观看视频色| 久久精品久久久久久噜噜老黄 | 欧美bdsm另类| 日本黄色片子视频| 日韩欧美 国产精品| av视频在线观看入口| 日日摸夜夜添夜夜添小说| 黄色女人牲交| 精品日产1卡2卡| 国产高潮美女av| 国产精品久久久久久av不卡| 亚洲无线观看免费| 97超视频在线观看视频| 国产精品永久免费网站| 18禁黄网站禁片免费观看直播| 国产真实乱freesex| 国语自产精品视频在线第100页| 淫秽高清视频在线观看| 亚洲av成人精品一区久久| 婷婷色综合大香蕉| 狂野欧美白嫩少妇大欣赏| 国产精品98久久久久久宅男小说| 狠狠狠狠99中文字幕| 国产激情偷乱视频一区二区| 亚洲18禁久久av| 精品日产1卡2卡| 欧美+亚洲+日韩+国产| a在线观看视频网站| 国产精品亚洲美女久久久| 十八禁网站免费在线| 特大巨黑吊av在线直播| 无遮挡黄片免费观看| 精品一区二区三区视频在线| 久久精品91蜜桃| 嫩草影院新地址| 天堂av国产一区二区熟女人妻| 久久久成人免费电影| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区 | 日本五十路高清| 久久久色成人| 小说图片视频综合网站| 啪啪无遮挡十八禁网站| 亚洲最大成人中文| 免费看美女性在线毛片视频| 赤兔流量卡办理| 97超视频在线观看视频| 天堂影院成人在线观看| 国内揄拍国产精品人妻在线| 成熟少妇高潮喷水视频| 成人三级黄色视频| 久久精品国产亚洲网站| 精品人妻偷拍中文字幕| 国产精品嫩草影院av在线观看 | 亚洲av第一区精品v没综合| 午夜视频国产福利| 88av欧美| 级片在线观看| 99热网站在线观看| 亚洲美女搞黄在线观看 | 岛国在线免费视频观看| 国产精品亚洲美女久久久| 久久久久久国产a免费观看| 在线播放国产精品三级| 少妇丰满av| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| 18禁黄网站禁片午夜丰满| 国产午夜福利久久久久久| 亚洲五月天丁香| 欧美精品啪啪一区二区三区| 69人妻影院| 午夜免费男女啪啪视频观看 | 久99久视频精品免费| 精品不卡国产一区二区三区| videossex国产| 亚洲性久久影院| 他把我摸到了高潮在线观看| 好男人在线观看高清免费视频| 亚洲成a人片在线一区二区| 久久久久久久精品吃奶| 精品日产1卡2卡| 国产精品98久久久久久宅男小说| 十八禁网站免费在线| 久9热在线精品视频| 小蜜桃在线观看免费完整版高清| 十八禁网站免费在线| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 亚洲人成网站在线播| 久久人人精品亚洲av| 三级毛片av免费| 91久久精品电影网| 亚洲最大成人手机在线| 亚洲国产高清在线一区二区三| 毛片一级片免费看久久久久 | 亚洲精品成人久久久久久| 成年免费大片在线观看| ponron亚洲| 日韩亚洲欧美综合| 欧美日韩黄片免| 国产精品人妻久久久久久| АⅤ资源中文在线天堂| 免费人成在线观看视频色| 亚洲av美国av| 国产精品一区二区免费欧美| 97超视频在线观看视频| 黄色丝袜av网址大全| 精品国产三级普通话版| 国产精品永久免费网站| 亚洲综合色惰| 亚洲天堂国产精品一区在线| 免费不卡的大黄色大毛片视频在线观看 | 成人美女网站在线观看视频| 波多野结衣高清作品| 亚洲人成网站高清观看| 久久人人爽人人爽人人片va| 精品国内亚洲2022精品成人| 丰满人妻一区二区三区视频av| 又爽又黄a免费视频| avwww免费| 蜜桃亚洲精品一区二区三区| 欧美bdsm另类| 日本在线视频免费播放| 国产私拍福利视频在线观看| 久久国产精品人妻蜜桃| 欧美精品国产亚洲| 久久国产精品人妻蜜桃| 国产亚洲91精品色在线| 看片在线看免费视频| 亚州av有码| 国产又黄又爽又无遮挡在线| 天堂动漫精品| 中文字幕久久专区| 欧美一区二区国产精品久久精品| 看片在线看免费视频| 国产精品日韩av在线免费观看| 免费观看精品视频网站| 有码 亚洲区| 免费人成在线观看视频色| 少妇的逼好多水| 欧美潮喷喷水| 中文字幕熟女人妻在线| 老司机深夜福利视频在线观看| 校园人妻丝袜中文字幕| 久久久久国产精品人妻aⅴ院| 岛国在线免费视频观看| 国产精品1区2区在线观看.| 日本 欧美在线| 成年版毛片免费区| 99riav亚洲国产免费| 日韩中文字幕欧美一区二区| 麻豆成人午夜福利视频| 欧美又色又爽又黄视频| 免费看日本二区| 精品久久久久久久久亚洲 | 长腿黑丝高跟| 中文字幕av在线有码专区| 国产一区二区三区在线臀色熟女| 欧美激情国产日韩精品一区| 91在线精品国自产拍蜜月| 欧美激情久久久久久爽电影| 欧美成人一区二区免费高清观看| 久久香蕉精品热| 一本久久中文字幕| 色哟哟哟哟哟哟| 人妻少妇偷人精品九色| 此物有八面人人有两片| 热99re8久久精品国产| 99国产极品粉嫩在线观看| 啦啦啦啦在线视频资源| 男女下面进入的视频免费午夜| 国产久久久一区二区三区| 国产精品av视频在线免费观看| 亚洲美女搞黄在线观看 | 99久久精品国产国产毛片| 国内精品美女久久久久久| 在现免费观看毛片| 亚洲av中文av极速乱 | 99久久精品热视频| 天堂网av新在线| 国产一区二区激情短视频| 少妇的逼水好多| 18禁在线播放成人免费| 少妇高潮的动态图| 男人的好看免费观看在线视频| 欧美日韩黄片免| 亚洲成人久久爱视频| 99国产极品粉嫩在线观看| 性欧美人与动物交配| 亚洲成人中文字幕在线播放| 精品久久久久久久久久久久久| 色吧在线观看| 国产精品电影一区二区三区| 国产精品女同一区二区软件 | 国产精品嫩草影院av在线观看 | 午夜精品在线福利| 国产精品久久久久久亚洲av鲁大| 亚洲欧美精品综合久久99| 婷婷六月久久综合丁香| 黄色欧美视频在线观看| 综合色av麻豆| 欧美激情久久久久久爽电影| 午夜老司机福利剧场| 国内精品美女久久久久久| 国产精品美女特级片免费视频播放器| 亚洲美女黄片视频| 日本三级黄在线观看| 亚洲av免费在线观看| 国产av麻豆久久久久久久| 听说在线观看完整版免费高清| 欧美人与善性xxx| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 亚洲国产色片| 国产真实伦视频高清在线观看 | 国产一区二区在线av高清观看| netflix在线观看网站| 欧美成人性av电影在线观看| 夜夜看夜夜爽夜夜摸| 成人av在线播放网站| 神马国产精品三级电影在线观看| 亚洲精品456在线播放app | 麻豆一二三区av精品| 偷拍熟女少妇极品色| 欧美不卡视频在线免费观看| 午夜精品在线福利| 2021天堂中文幕一二区在线观| 亚洲av电影不卡..在线观看| 五月玫瑰六月丁香| 色5月婷婷丁香| 毛片一级片免费看久久久久 | 岛国在线免费视频观看| 性欧美人与动物交配| 国产真实伦视频高清在线观看 | 亚洲美女视频黄频| 99热精品在线国产| 丰满乱子伦码专区| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久末码| 午夜福利高清视频| 老司机福利观看| 日韩欧美三级三区| 国产亚洲精品av在线| 啪啪无遮挡十八禁网站| 1000部很黄的大片| 国产单亲对白刺激| 99热只有精品国产| 搡老熟女国产l中国老女人| 99久久无色码亚洲精品果冻| 麻豆久久精品国产亚洲av| 99热这里只有是精品在线观看| 国内揄拍国产精品人妻在线| 国产探花在线观看一区二区| 国产黄色小视频在线观看| 成年版毛片免费区| 成人性生交大片免费视频hd| 亚洲精品成人久久久久久| 久久国产精品人妻蜜桃| 国产精品自产拍在线观看55亚洲| 99久久九九国产精品国产免费| 午夜免费激情av| 精品人妻偷拍中文字幕| 精品99又大又爽又粗少妇毛片 | 特级一级黄色大片| 美女cb高潮喷水在线观看| 亚洲精品在线观看二区| 亚洲美女搞黄在线观看 | 国产麻豆成人av免费视频| 少妇人妻一区二区三区视频| 又黄又爽又刺激的免费视频.| 国产一区二区三区在线臀色熟女| 黄色日韩在线| 99精品久久久久人妻精品| 无遮挡黄片免费观看| 乱码一卡2卡4卡精品| 国产精品爽爽va在线观看网站| 欧美精品啪啪一区二区三区| 在现免费观看毛片| 亚洲自拍偷在线| 99久国产av精品| av国产免费在线观看| 国产单亲对白刺激| 精品人妻视频免费看| 一级av片app| ponron亚洲| 麻豆国产av国片精品| 在线观看舔阴道视频| 18禁在线播放成人免费| 国产精品亚洲美女久久久| 免费看美女性在线毛片视频| av在线亚洲专区| 一个人免费在线观看电影| 国产精品乱码一区二三区的特点| 搡老妇女老女人老熟妇| 尾随美女入室| 国产伦人伦偷精品视频| 国产毛片a区久久久久| 69av精品久久久久久| 国产精品野战在线观看| 国产一区二区在线av高清观看| 赤兔流量卡办理| 欧美黑人巨大hd| ponron亚洲| 亚洲人成网站在线播| 精品久久国产蜜桃| 国产黄片美女视频| 国产中年淑女户外野战色| 久久久成人免费电影| 免费av观看视频| 久久久久国内视频| 亚洲乱码一区二区免费版| 久久久久精品国产欧美久久久| 噜噜噜噜噜久久久久久91| 国产精品98久久久久久宅男小说| 非洲黑人性xxxx精品又粗又长| 精品国产三级普通话版| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看吧| 亚洲精品亚洲一区二区| 欧美国产日韩亚洲一区| 嫩草影院精品99| 日韩中文字幕欧美一区二区| 国产单亲对白刺激| 哪里可以看免费的av片| 亚洲色图av天堂| 热99re8久久精品国产| 国产精品久久久久久精品电影| 中文字幕高清在线视频| 性色avwww在线观看| 国产精品福利在线免费观看| 成人美女网站在线观看视频| 一区二区三区四区激情视频 | 日韩国内少妇激情av| 麻豆国产av国片精品| 日韩精品有码人妻一区| 中文字幕高清在线视频| a在线观看视频网站| 亚洲欧美日韩无卡精品| 午夜精品一区二区三区免费看| 岛国在线免费视频观看| 男女啪啪激烈高潮av片| 亚洲天堂国产精品一区在线| 欧美日韩乱码在线| 一区二区三区高清视频在线| 精品99又大又爽又粗少妇毛片 | 国产乱人伦免费视频| 亚洲人成伊人成综合网2020| 日韩欧美三级三区| 久久久久国产精品人妻aⅴ院| 熟女人妻精品中文字幕| 一进一出好大好爽视频| 熟妇人妻久久中文字幕3abv| 久久午夜亚洲精品久久| 国内揄拍国产精品人妻在线| 亚洲成av人片在线播放无| 国产亚洲精品久久久com| 精品久久久久久久久亚洲 | 亚洲黑人精品在线| 久久精品人妻少妇| 国产 一区 欧美 日韩| 在线观看美女被高潮喷水网站| 免费在线观看影片大全网站| 亚洲avbb在线观看| 久久久久久大精品| 精品一区二区三区视频在线观看免费| 午夜福利成人在线免费观看| 日韩强制内射视频| 不卡一级毛片| 日本三级黄在线观看| 偷拍熟女少妇极品色| 1000部很黄的大片| 亚洲第一电影网av| 欧美一区二区亚洲| 三级国产精品欧美在线观看| 久久精品国产亚洲av天美| 超碰av人人做人人爽久久| 一进一出抽搐gif免费好疼| 精品久久久久久久久久久久久| 国产伦人伦偷精品视频| 国产乱人视频| 老师上课跳d突然被开到最大视频| 日本a在线网址| 亚洲av中文av极速乱 | 一个人免费在线观看电影| 精品人妻偷拍中文字幕| 久久国产精品人妻蜜桃| 91久久精品国产一区二区成人|