• <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高清不卡| 色94色欧美一区二区| 最黄视频免费看| 欧美乱码精品一区二区三区| 国产日韩欧美亚洲二区| 久久免费观看电影| 丝袜美足系列| 中文字幕人妻丝袜一区二区| 久久精品国产99精品国产亚洲性色 | 三上悠亚av全集在线观看| 黄频高清免费视频| 精品一区二区三区av网在线观看 | 国产欧美日韩一区二区三区在线| 精品国产乱码久久久久久男人| 久久中文字幕一级| 亚洲avbb在线观看| 欧美+亚洲+日韩+国产| 亚洲国产欧美一区二区综合| 亚洲 欧美一区二区三区| 成年动漫av网址| 亚洲人成伊人成综合网2020| 色老头精品视频在线观看| 日韩三级视频一区二区三区| 国产精品免费大片| 亚洲伊人色综图| 久久久国产成人免费| 亚洲欧美激情在线| 99国产精品一区二区蜜桃av | 丝袜美腿诱惑在线| 久久热在线av| 国产精品香港三级国产av潘金莲| 国产成人一区二区三区免费视频网站| 久久久精品区二区三区| 夜夜爽天天搞| 岛国在线观看网站| 91成年电影在线观看| 怎么达到女性高潮| 午夜免费成人在线视频| 国产极品粉嫩免费观看在线| 不卡一级毛片| 国产精品偷伦视频观看了| 五月天丁香电影| 精品国产乱码久久久久久男人| 高清在线国产一区| 美女午夜性视频免费| 欧美黄色淫秽网站| 国产男靠女视频免费网站| 中文字幕人妻熟女乱码| 亚洲精品中文字幕在线视频| 啪啪无遮挡十八禁网站| av不卡在线播放| 久久精品人人爽人人爽视色| 大片免费播放器 马上看| 精品国产乱子伦一区二区三区| 一区在线观看完整版| 国产一区二区激情短视频| 999久久久国产精品视频| 99精品在免费线老司机午夜| 国产三级黄色录像| 亚洲成a人片在线一区二区| 日韩制服丝袜自拍偷拍| 免费在线观看黄色视频的| 每晚都被弄得嗷嗷叫到高潮| 亚洲一卡2卡3卡4卡5卡精品中文| 夜夜骑夜夜射夜夜干| 丝袜美腿诱惑在线| 亚洲av成人不卡在线观看播放网| 欧美精品啪啪一区二区三区| 久久久国产一区二区| 精品福利永久在线观看| 考比视频在线观看| 悠悠久久av| 国产精品秋霞免费鲁丝片| 精品亚洲成a人片在线观看| 亚洲国产欧美在线一区| 久久精品aⅴ一区二区三区四区| 美女扒开内裤让男人捅视频| 久久精品国产亚洲av香蕉五月 | 在线永久观看黄色视频| 国产精品av久久久久免费| 午夜福利视频精品| 亚洲九九香蕉| 高清欧美精品videossex| 黄色视频,在线免费观看| 久久精品亚洲精品国产色婷小说| 精品免费久久久久久久清纯 | 18在线观看网站| 黄网站色视频无遮挡免费观看| 亚洲人成电影观看| 免费在线观看影片大全网站| 久久这里只有精品19| 搡老岳熟女国产| 国产精品欧美亚洲77777| 午夜福利视频在线观看免费| 亚洲成国产人片在线观看| 精品一品国产午夜福利视频| 制服人妻中文乱码| 久久久久精品人妻al黑| 国产成人免费无遮挡视频| 国产av精品麻豆| 午夜福利影视在线免费观看| 国产成人精品久久二区二区91| 咕卡用的链子| 天堂中文最新版在线下载| 国产日韩欧美视频二区| 99久久国产精品久久久| 老司机靠b影院| 亚洲免费av在线视频| 成年动漫av网址| 脱女人内裤的视频| 免费观看人在逋| 久热爱精品视频在线9| av福利片在线| 国产精品久久久av美女十八| 国产成人啪精品午夜网站| 男女床上黄色一级片免费看| 国产精品一区二区免费欧美| 午夜福利免费观看在线| 天天躁日日躁夜夜躁夜夜| 性高湖久久久久久久久免费观看| 国产亚洲精品一区二区www | 国产区一区二久久| 欧美黑人精品巨大| 国产成人一区二区三区免费视频网站| 国产日韩欧美在线精品| 亚洲情色 制服丝袜| 欧美激情高清一区二区三区| 一边摸一边做爽爽视频免费| 18禁黄网站禁片午夜丰满| 嫩草影视91久久| 国产极品粉嫩免费观看在线| 久久狼人影院| 午夜免费鲁丝| 国内毛片毛片毛片毛片毛片| 侵犯人妻中文字幕一二三四区| bbb黄色大片| 久久精品aⅴ一区二区三区四区| 黄色毛片三级朝国网站| 男人舔女人的私密视频| 777米奇影视久久| 一二三四社区在线视频社区8| av不卡在线播放| 如日韩欧美国产精品一区二区三区| 91成人精品电影| 亚洲成人免费av在线播放| 妹子高潮喷水视频| 欧美精品啪啪一区二区三区| 久久99热这里只频精品6学生| 黄频高清免费视频| 国产亚洲一区二区精品| 又大又爽又粗| 久久人人97超碰香蕉20202| 国产精品欧美亚洲77777| 国产精品av久久久久免费| 美女视频免费永久观看网站| 国产激情久久老熟女| 久久亚洲精品不卡| 亚洲专区字幕在线| 操出白浆在线播放| 久久天躁狠狠躁夜夜2o2o| 精品国产一区二区三区久久久樱花| 久久av网站| 精品一区二区三卡| 韩国精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 99riav亚洲国产免费| 久久精品熟女亚洲av麻豆精品| 咕卡用的链子| 亚洲人成电影免费在线| 国产欧美日韩一区二区三| 色婷婷久久久亚洲欧美| 精品一区二区三区四区五区乱码| 深夜精品福利| 一区二区三区国产精品乱码| 亚洲精品自拍成人| 丁香欧美五月| 亚洲av片天天在线观看| 女警被强在线播放| 亚洲男人天堂网一区| 国产精品99久久99久久久不卡| 久久精品亚洲av国产电影网| 久久久精品免费免费高清| 99国产精品一区二区三区| 操出白浆在线播放| 人妻一区二区av| 在线 av 中文字幕| 欧美大码av| 精品熟女少妇八av免费久了| 露出奶头的视频| 国产精品久久久av美女十八| 亚洲欧美精品综合一区二区三区| 最新美女视频免费是黄的| 别揉我奶头~嗯~啊~动态视频| 国产av精品麻豆| 日韩 欧美 亚洲 中文字幕| 大码成人一级视频| 亚洲欧美色中文字幕在线| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 欧美 日韩 精品 国产| 91老司机精品| tocl精华| 免费少妇av软件| 久久国产亚洲av麻豆专区| 桃红色精品国产亚洲av| 欧美变态另类bdsm刘玥| 不卡av一区二区三区| 99riav亚洲国产免费| 在线永久观看黄色视频| 欧美国产精品va在线观看不卡| 热99re8久久精品国产| 亚洲第一av免费看| 国产在视频线精品| 国产成人免费观看mmmm| 国产精品一区二区在线观看99| 五月天丁香电影| 日本wwww免费看| 俄罗斯特黄特色一大片| 涩涩av久久男人的天堂| 丰满少妇做爰视频| e午夜精品久久久久久久| 老司机影院毛片| 精品亚洲成a人片在线观看| av电影中文网址| 国产97色在线日韩免费| 成人免费观看视频高清| 欧美激情久久久久久爽电影 | 欧美av亚洲av综合av国产av| 一区二区三区乱码不卡18| 日韩大片免费观看网站| 免费在线观看完整版高清| 亚洲精品中文字幕在线视频| 久久青草综合色| 日本一区二区免费在线视频| 黄色视频不卡| 极品少妇高潮喷水抽搐| 国产亚洲精品一区二区www | 巨乳人妻的诱惑在线观看| 不卡av一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| kizo精华| 亚洲中文字幕日韩| 亚洲精品久久午夜乱码| 欧美日韩福利视频一区二区| 亚洲人成77777在线视频| 国产精品久久久久久精品古装| 看免费av毛片| 日韩一卡2卡3卡4卡2021年| 国产成人一区二区三区免费视频网站| 日韩有码中文字幕| 人人妻人人澡人人看| 亚洲av日韩在线播放| 男女下面插进去视频免费观看| 亚洲男人天堂网一区| 狠狠精品人妻久久久久久综合| 亚洲中文日韩欧美视频| 狂野欧美激情性xxxx| 国产一区二区激情短视频| av天堂在线播放| 国产精品一区二区精品视频观看| 成人影院久久| 黄色片一级片一级黄色片| 久久人妻熟女aⅴ| 久久人人97超碰香蕉20202| 国产日韩欧美亚洲二区| 久久久久精品人妻al黑| 大型av网站在线播放| av福利片在线| 成人18禁在线播放| 国产aⅴ精品一区二区三区波| 久久ye,这里只有精品| 精品第一国产精品| 午夜久久久在线观看| 亚洲精品av麻豆狂野| 亚洲国产欧美在线一区| 国产精品久久久久久人妻精品电影 | 精品熟女少妇八av免费久了| 亚洲精品粉嫩美女一区| 脱女人内裤的视频| 操美女的视频在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久久久久精品吃奶| 肉色欧美久久久久久久蜜桃| 91大片在线观看| 国产精品二区激情视频| 免费av中文字幕在线| 日本vs欧美在线观看视频| 欧美日韩视频精品一区| 一区二区三区精品91| 一边摸一边抽搐一进一小说 | 亚洲国产成人一精品久久久| 亚洲色图av天堂| 亚洲一区中文字幕在线| 免费在线观看视频国产中文字幕亚洲| 91精品国产国语对白视频| 亚洲 欧美一区二区三区| 天堂8中文在线网| 制服人妻中文乱码| 欧美精品av麻豆av| 国产高清videossex| 国产精品 欧美亚洲| 久久久久久久久久久久大奶| 久久精品国产亚洲av高清一级| 久久中文字幕人妻熟女| 露出奶头的视频| 亚洲人成77777在线视频| 成人免费观看视频高清| 美女高潮到喷水免费观看| 欧美激情极品国产一区二区三区| 精品少妇黑人巨大在线播放| 国产91精品成人一区二区三区 | 黄色片一级片一级黄色片| 日本一区二区免费在线视频| 国产精品美女特级片免费视频播放器 | 亚洲精品粉嫩美女一区| 亚洲精品国产精品久久久不卡| 啦啦啦 在线观看视频| 五月开心婷婷网| 一个人免费在线观看的高清视频| 亚洲av电影在线进入| 亚洲av电影在线进入| 黑丝袜美女国产一区| 俄罗斯特黄特色一大片| 91av网站免费观看| 最近最新中文字幕大全免费视频| 国产免费福利视频在线观看| 国产成人精品无人区| 久久久久久久久久久久大奶| 中亚洲国语对白在线视频| 99riav亚洲国产免费| 欧美乱妇无乱码| 日韩有码中文字幕| 日韩一卡2卡3卡4卡2021年| 久久久久久亚洲精品国产蜜桃av| av又黄又爽大尺度在线免费看| 亚洲综合色网址| www.熟女人妻精品国产| 国产精品一区二区免费欧美| 看免费av毛片| 亚洲av成人一区二区三| 国产精品久久久久成人av| www日本在线高清视频| a级毛片黄视频| 黄网站色视频无遮挡免费观看| 蜜桃国产av成人99| 日韩欧美国产一区二区入口| tube8黄色片| 两个人免费观看高清视频| 久久av网站| 亚洲av日韩精品久久久久久密| 俄罗斯特黄特色一大片| 亚洲黑人精品在线| 中文欧美无线码| 亚洲人成伊人成综合网2020| 自拍欧美九色日韩亚洲蝌蚪91| av一本久久久久| 免费在线观看日本一区| 精品国产国语对白av| 曰老女人黄片| 国产视频一区二区在线看| 一级a爱视频在线免费观看| 男男h啪啪无遮挡| 91九色精品人成在线观看| 亚洲全国av大片| 久久久精品免费免费高清| 午夜福利在线免费观看网站| 久久免费观看电影| 亚洲人成电影免费在线| 一级毛片女人18水好多| 婷婷成人精品国产| 亚洲精品美女久久av网站| 久久久久久久精品吃奶| 成人黄色视频免费在线看| 国产精品亚洲av一区麻豆| 久久影院123| 咕卡用的链子| 亚洲精品成人av观看孕妇| 纯流量卡能插随身wifi吗| 国产xxxxx性猛交| 国产欧美日韩精品亚洲av| 中文字幕人妻丝袜制服| 十八禁人妻一区二区| 精品久久久久久电影网| √禁漫天堂资源中文www| 日日爽夜夜爽网站| 国产精品久久久久久精品电影小说| 亚洲 国产 在线| 亚洲七黄色美女视频| 老司机福利观看| 777久久人妻少妇嫩草av网站| 久热这里只有精品99| 亚洲av国产av综合av卡| 国产精品一区二区在线不卡| 亚洲三区欧美一区| 99国产精品免费福利视频| 久久这里只有精品19| 制服人妻中文乱码| 99国产精品一区二区三区| 欧美亚洲日本最大视频资源| 国产日韩欧美在线精品| 国产一区二区三区综合在线观看| 性色av乱码一区二区三区2| 成人黄色视频免费在线看| 大片电影免费在线观看免费| 窝窝影院91人妻| 欧美精品亚洲一区二区| 亚洲全国av大片| 国产一卡二卡三卡精品| 俄罗斯特黄特色一大片| 欧美精品一区二区免费开放| 精品午夜福利视频在线观看一区 | 91九色精品人成在线观看| 动漫黄色视频在线观看| 国产在线观看jvid| 国产一区二区激情短视频| 色精品久久人妻99蜜桃| 亚洲国产av影院在线观看| 国产亚洲午夜精品一区二区久久| 成人精品一区二区免费| 一区二区三区激情视频| 亚洲国产看品久久| 国产精品一区二区在线观看99| 这个男人来自地球电影免费观看| 亚洲精品乱久久久久久| 国产免费视频播放在线视频| 九色亚洲精品在线播放| 日韩成人在线观看一区二区三区| 国产伦理片在线播放av一区| 久久精品国产亚洲av高清一级| 18禁黄网站禁片午夜丰满| 国产aⅴ精品一区二区三区波| 午夜91福利影院| 精品一区二区三区视频在线观看免费 | 夫妻午夜视频| 欧美性长视频在线观看| av线在线观看网站| 捣出白浆h1v1| 这个男人来自地球电影免费观看| 欧美成人免费av一区二区三区 | 黄频高清免费视频| 国产麻豆69| 日日爽夜夜爽网站| 丰满饥渴人妻一区二区三| 久久精品熟女亚洲av麻豆精品| 国内毛片毛片毛片毛片毛片| 视频区图区小说| 我的亚洲天堂| 欧美成人午夜精品| 久久国产精品男人的天堂亚洲| av有码第一页| 欧美日韩一级在线毛片| 国产成人av激情在线播放| 麻豆乱淫一区二区| 久久久精品免费免费高清| 亚洲欧美日韩另类电影网站| 国产在线观看jvid| 欧美激情高清一区二区三区| 亚洲人成伊人成综合网2020| 真人做人爱边吃奶动态| 久久久久久久大尺度免费视频| 欧美黑人欧美精品刺激| 丝袜人妻中文字幕| 正在播放国产对白刺激| 水蜜桃什么品种好| avwww免费| 69av精品久久久久久 | 黑人巨大精品欧美一区二区mp4| 欧美黑人精品巨大| 一级片免费观看大全| www.精华液| 精品福利观看| 久久久精品免费免费高清| 国产精品久久久久久人妻精品电影 | 午夜久久久在线观看| 下体分泌物呈黄色| 欧美黑人精品巨大| av超薄肉色丝袜交足视频| 老司机福利观看| 亚洲午夜精品一区,二区,三区| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 成人三级做爰电影| 手机成人av网站| 亚洲精品久久成人aⅴ小说| 亚洲黑人精品在线| 男女高潮啪啪啪动态图| 国产日韩欧美视频二区| 午夜福利,免费看| 欧美人与性动交α欧美精品济南到| 精品少妇久久久久久888优播| 亚洲精品国产区一区二| 亚洲精品国产精品久久久不卡| 国产成人av教育| 女人爽到高潮嗷嗷叫在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲精品久久久久5区| 婷婷丁香在线五月| 欧美日韩精品网址| 亚洲第一青青草原| 国产欧美日韩一区二区精品| 啦啦啦视频在线资源免费观看| 韩国精品一区二区三区| 香蕉丝袜av| 法律面前人人平等表现在哪些方面| 麻豆成人av在线观看| 国产在线视频一区二区| 久久久精品国产亚洲av高清涩受| 欧美日韩精品网址| 肉色欧美久久久久久久蜜桃| 亚洲成人手机| 日本精品一区二区三区蜜桃| 欧美另类亚洲清纯唯美| 人人妻,人人澡人人爽秒播| 国产色视频综合| 亚洲一区中文字幕在线| a级片在线免费高清观看视频| 日本一区二区免费在线视频| 亚洲专区国产一区二区| 岛国毛片在线播放| 婷婷丁香在线五月| 精品一区二区三区视频在线观看免费 | 亚洲精品中文字幕在线视频| 色94色欧美一区二区| 亚洲综合色网址| 久久久精品区二区三区| 国产男女超爽视频在线观看| 国产精品免费大片| 精品国产超薄肉色丝袜足j| 老司机影院毛片| svipshipincom国产片| 免费在线观看影片大全网站| 欧美激情高清一区二区三区| 伦理电影免费视频| 久久人人97超碰香蕉20202| 久久久久久久精品吃奶| 国产视频一区二区在线看| 国产精品 国内视频| 激情在线观看视频在线高清 | 国产一区二区三区在线臀色熟女 | 人人妻人人爽人人添夜夜欢视频| 免费久久久久久久精品成人欧美视频| 午夜福利视频在线观看免费| 国产av又大| 免费在线观看视频国产中文字幕亚洲| 蜜桃在线观看..| 两人在一起打扑克的视频| 又紧又爽又黄一区二区| 国产男女超爽视频在线观看| 久久久久精品国产欧美久久久| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 精品国产国语对白av| 亚洲成人国产一区在线观看| av片东京热男人的天堂| 黄网站色视频无遮挡免费观看| 99国产综合亚洲精品| 两个人看的免费小视频| 国产一区有黄有色的免费视频| 成年人黄色毛片网站| 老司机亚洲免费影院| 国产在线免费精品| svipshipincom国产片| 国产在线一区二区三区精| 这个男人来自地球电影免费观看| 在线观看舔阴道视频| 午夜福利,免费看| 亚洲综合色网址| 欧美+亚洲+日韩+国产| 最新的欧美精品一区二区| 欧美日韩国产mv在线观看视频| 欧美黄色片欧美黄色片| 中文字幕色久视频| 亚洲 国产 在线| 桃红色精品国产亚洲av| 不卡av一区二区三区| 美女高潮到喷水免费观看| av一本久久久久| 国产一区有黄有色的免费视频| 中文字幕人妻熟女乱码| 女同久久另类99精品国产91| 俄罗斯特黄特色一大片| 99精品在免费线老司机午夜| 久久毛片免费看一区二区三区| 99久久精品国产亚洲精品| 美女视频免费永久观看网站| 男女边摸边吃奶| 不卡av一区二区三区| 久久精品人人爽人人爽视色| 两个人看的免费小视频| 国产又色又爽无遮挡免费看| 日本撒尿小便嘘嘘汇集6| 女人被躁到高潮嗷嗷叫费观| 这个男人来自地球电影免费观看| 欧美 亚洲 国产 日韩一| 久久久国产精品麻豆| 怎么达到女性高潮| 黄片小视频在线播放| 免费久久久久久久精品成人欧美视频| videos熟女内射| 无人区码免费观看不卡 | 51午夜福利影视在线观看| 真人做人爱边吃奶动态| 国产精品免费一区二区三区在线 | 免费在线观看完整版高清| 午夜成年电影在线免费观看| 男女午夜视频在线观看| 视频区图区小说| 每晚都被弄得嗷嗷叫到高潮| 久久性视频一级片| 国产视频一区二区在线看| 啪啪无遮挡十八禁网站|