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

    融合無人機和衛(wèi)星影像的溫帶疏林草原木本和草本植物覆蓋度遙感估算

    2023-10-19 02:48:04李曉雅田昕段濤曹曉明楊凱捷盧琦王鋒
    遙感學報 2023年9期
    關鍵詞:特征模型

    李曉雅,田昕,段濤,曹曉明,楊凱捷,盧琦,王鋒

    1.中國林業(yè)科學研究院荒漠化研究所,北京 100091;2.中國林業(yè)科學研究院資源信息研究所,北京 100091;3.中國科學院微電子研究所,北京 100029

    1 引言

    溫帶疏林草原是中國半干旱區(qū)的重要生態(tài)系統(tǒng),是在中國北方森林、草原生態(tài)系統(tǒng)過渡區(qū)域發(fā)育于沙地的地帶性頂級群落(于順利,2011;史宇飛 等,2011)。榆樹疏林曾廣泛分布于松遼平原和內(nèi)蒙古高原,由于生境的喪失,目前主要分布于內(nèi)蒙古高原的渾善達克沙地、科爾沁沙地、呼倫貝爾沙地和毛烏素沙地等地,其中以渾善達克沙地和科爾沁沙地分布面積最廣,疏林草原景觀最為典型(于順利和陳宏偉,2007)。準確掌握疏林草原的分布及動態(tài),不僅有利于加深對中國溫帶疏林草原生態(tài)系統(tǒng)的理解,進一步也可為溫帶疏林草原的恢復和重建提供科學依據(jù)。

    植被覆蓋度是定量化表征地表植物群落的重要指標,同時反映著區(qū)域生態(tài)環(huán)境的變化(DeFries等,2007;Yang 等,2012)。定量獲取疏林草原植被覆蓋度、準確反演植被動態(tài)變化規(guī)律,是實現(xiàn)疏林草原動態(tài)監(jiān)測與評估的前提。獲取植被覆蓋度的方法主要包括地面測量和遙感解譯兩種途徑。野外觀測獲得的植被覆蓋度信息是實現(xiàn)覆蓋度遙感反演的基礎,為覆蓋度遙感反演提供訓練數(shù)據(jù)及驗證信息。目前使用遙感手段進行覆蓋度反演的方主要包括植被指數(shù)、經(jīng)驗模型、混合像元分解模型和機器學習方法,植被指數(shù)法通過建立植被指數(shù)與覆蓋度之間的關系反演植被覆蓋度(Harris等,2014;Wang 等,2020);經(jīng)驗模型即建立波段信息與覆蓋度之間的經(jīng)驗模型估算植被覆蓋度(岳健 等,2020),經(jīng)驗模型對數(shù)據(jù)集的正態(tài)分布特性和波段信息的獨立性要求較高,僅適用于特定的區(qū)域及特定值被,使其在大尺度植被覆蓋度估算中具有一定的局限性(賈坤 等,2013);混合像元分解方法同樣被廣泛應用于半干旱地區(qū)植被覆蓋度的估測中(Cui 等,2019),該方法利用像元分解模型分解像元中各端元的亮度信息貢獻,獲得端元植被豐度以估算值被覆蓋度(Adams 等,1986)?;旌舷裨纸夥椒ㄇ捌趯庾V端元的選取嚴重影響到像元分解的精度,植被覆蓋度估測能力不穩(wěn)定(李曉松 等,2010),且疏林草原具有較復雜的地表覆蓋情況,端元組合較為復雜。

    機器學習方法是通過構(gòu)建樣本數(shù)據(jù)集,結(jié)合地形、氣候、植被指數(shù)、光譜信息等預測變量,運用數(shù)據(jù)挖掘方法找出隱含在預測變量和分類樣本之間的相關關系,并自動建立估算模型,屬于遙感影像解譯中全監(jiān)督學習類方法(周培誠 等,2021)。由于機器方法更適合處理復雜海量的數(shù)據(jù)集(Brandt等,2016)并具有較高的計算效率與估算精度,在干旱半干旱區(qū)植被覆蓋類型圖的繪制以及植被覆蓋度反演等方面得到了廣泛的應用(Zhang等,2019;陳黔 等,2019;Hunter等,2020)。

    疏林草原具有喬木、灌木和草本植被混合生長的特點,在植物斑塊化分布、高異質(zhì)性景觀的疏林草原區(qū)域,利用中低分辨率衛(wèi)星遙感影像識別疏林草原,精度難以保證(Latifovic 和Olthof,2004)。近年來,覆蓋全球的高分辨率衛(wèi)星不斷涌現(xiàn),出現(xiàn)了覆蓋全球范圍的中高分辨率衛(wèi)星如哨兵(Sentinel)和國產(chǎn)高分(GF),為區(qū)域尺度的疏林草原植被覆蓋度估算提供了有利條件。

    無人機飛行時間、靈活、可搭載多種類型傳感器,具有獲取植被結(jié)構(gòu)特征、識別矮小植株的高空間分辨率優(yōu)勢(Kattenborn 等,2019),被廣泛應用到植被遙感分類和監(jiān)測中(韓東 等,2018;Wang等,2019;張順 等,2019;蔡棟 等,2019)。無人機近地遙感觀測技術(shù)的逐漸成熟,也帶動了遙感信息由過去地面監(jiān)測與衛(wèi)星遙感的“二維數(shù)據(jù)”到星—空—地“三維數(shù)據(jù)”的轉(zhuǎn)變,成為地面觀測與衛(wèi)星遙感之間重要的紐帶,為生態(tài)遙感的進一步發(fā)展提供了新的機遇(蔡宗磊 等,2019;郭慶華 等,2020;廖小罕 等,2019)。不過由于受到無人機續(xù)航時間的限制,通過無人機實現(xiàn)區(qū)域尺度的植被監(jiān)測還面臨較大的挑戰(zhàn)。

    在利用不同空間尺度遙感數(shù)據(jù)時,較大的空間尺度可以反應區(qū)域上的規(guī)律,較精細的空間尺度可以揭示復雜的地表特征(萬昌君 等,2019)。最近出現(xiàn)的這些新數(shù)據(jù)源為我們帶來了啟示,是否可同時利用近地遙感和高分辨率衛(wèi)星遙感兩種高空間分辨率數(shù)據(jù)源,基于機器學習算法融合兩種數(shù)據(jù)源在植被遙感監(jiān)測中各自的優(yōu)勢,實現(xiàn)區(qū)域尺度溫帶疏林草原不同類型植物的高效精準識別和動態(tài)監(jiān)測?

    因此,本研究探索了如何利用高分辨率衛(wèi)星(GF-6 和Sentinel-2)結(jié)合近地遙感獲取的超高分辨率影像,基于機器學習算法估算疏林草原上木本植物和草本植被的覆蓋度。具體目標為:(1)基于無人機正射影像,利用分類和回歸樹算法對疏林草原木本植物和草本植物進行自動精準分類;(2)利用隨機森林算法,實現(xiàn)從無人機正射影像到高分衛(wèi)星影像的疏林草原木本植物和草本植被覆蓋度的估算;(3)比較GF-6 和Sentinel-2 兩種高分辨率衛(wèi)星影像對疏林草原木本植物和草本植被覆蓋度估算的差異。

    2 研究區(qū)與數(shù)據(jù)

    2.1 研究區(qū)概況

    研究區(qū)位于中國內(nèi)蒙古自治區(qū)錫林郭勒盟正藍旗渾善達克沙地榆樹疏林草原長期定位觀測大樣地,地理位置為42°57′34″N—42°58′15″N,115°57′E—115°58′E,平均海拔高度約為1300 m,面積為1 km×1 km(圖1)。該區(qū)域?qū)儆诟珊?、半干旱大陸性氣候,年平均氣?.1 ℃,降水集中于7—9 月,年平均降水量為313 mm。疏林草原生態(tài)系統(tǒng)由喬木層、灌木層和草本層3 層垂直結(jié)構(gòu)組成(于順利 等,2007)。將研究區(qū)內(nèi)的唯一喬木種榆樹(Ulmus pumila)和伴生灌木如耬斗葉繡線菊(Spiraea aquilegifolia)、小紅柳(Salix microstachyavar.bordensis)、黃柳(Salix gordejevii)、柴樺(Betula fruticosa)等14種灌木種定義為木本植物。草本層則主要由多年生草本褐沙蒿(Artemisia intramongolica)和一年生草本叉分蓼(Polygonum divaricatum),羊草(Aneurolepidum chinense)等近190 種草本植物構(gòu)成(韓東 等,2018)。

    圖1 中國內(nèi)蒙古自治區(qū)渾善達克沙地榆樹疏林草原長期定位監(jiān)測大樣地Fig.1 Long-term monitoring plot in Elm(Ulmus pumila)sparse forest grassland ecosystem in the Otindag sandy land(ESFOGE-Plot)Inner Mongolia,China

    2.2 遙感數(shù)據(jù)獲取及預處理

    2.2.1 無人機影像數(shù)據(jù)

    為了建立精準的疏林草原植被覆蓋度訓練數(shù)據(jù)集,選用配置1 cm 精度導航定位和具有2000 萬像素CMOS傳感器的小型多旋翼高精度航測無人機(精靈4 RTK,大疆,中國深圳),采集榆樹疏林草原長期定位觀測大樣地的無人機正射影像。本研究根據(jù)研究區(qū)的邊界設計無人機飛行方案進行航空影像采集,影像拍攝時間為2020 年7 月20 日,拍攝時長為2 h,共采集到1334 張原始航拍影像,無人機飛行參數(shù)設置見表1。利用實景三維建模軟件Agisoft Metashape(V1.6.5,Agisoft,俄羅斯)對原始航拍影像進行校正和拼接處理,最終獲取樣地的數(shù)字正射影像圖DOM(空間分辨率為2.41 cm/pixel,圖2(a))。

    表1 榆樹疏林草原長期定位觀測大樣地無人機航拍飛行參數(shù)Table 1 Flight parameters of Unmanned Aerial Vehicle(UAV)aerial photogrammetry in the long-term monitoring plot in Elm Sparse Forest Grassland

    圖2 榆樹疏林草原長期定位觀測大樣地無人機正射影像和高分衛(wèi)星影像(框圖部分為分類精度驗證區(qū)域)Fig.2 High resolution UAV and satellite images of long-term monitoring plot in Elm sparse forest grassland ecosystem(The area covered by a black box is the sample for validating the classification accuracy)

    2.2.2 衛(wèi)星影像數(shù)據(jù)

    選用GF-6 和Sentinel-2 作為衛(wèi)星影像數(shù)據(jù)??紤]衛(wèi)星過境時間及云量,分別選取2020 年7 月22 日GF-6 影像和2020 年7 月1 日Sentinel-2 影像。GF-6影像預處理過程包括輻射定標、大氣校正等。GF-6 衛(wèi)星輻射定標采用中國資源衛(wèi)星應用中心的系數(shù),并采用光譜超立方體的快速視線大氣分析(FLAASH)方法進行大氣校正。由歐洲航天局網(wǎng)站(https://scihub.copernicus.eu/[2021-01-08])下載的Sentinel-2 數(shù)據(jù)是經(jīng)過輻射校正后的反射率產(chǎn)品,利用ESA 官方發(fā)布的大氣校正模型Sen2cor(V2.8,ESA)進行大氣校正。由于Sentinel-2 影像可見光波段與紅邊波段和紅外波段具有不同的空間分辨率,對10 m 分辨率的可見光波段進行重采樣,重采樣后的Sentinel-2影像空間分辨率為20 m。分別將無人機近地面遙感平臺獲取的數(shù)字正射影像重采樣到16 m 和20 m,以重采樣后的無人機正射影像為基準,基于一階多項式分別對GF-6 影像和Sentinel-2 影像進行精確的幾何配準,校正過后的誤差不超過一個像元,保證了后期處理中對無人機正射影像和衛(wèi)星影像疊加分析的準確性;最后,利用研究區(qū)邊界分別對影像進行裁剪,獲得兩景預處理后的研究區(qū)影像(圖2(b),2(c))。輻射校正、大氣校正與圖像裁剪等預處理在ENVI(V5.3)中完成。

    3 研究方法

    3.1 疏林草原木本和草本植物覆蓋度估算方法

    本研究發(fā)展了一套基于高分辨率衛(wèi)星影像和機器學習算法估算疏林草原木本植物和草本植物的新方法。首先,(1)基于高分辨率的無人機正射影像,利用分類與回歸樹算法識別木本植物和草本植物,分別計算對應GF-6 和Sentinel-2 衛(wèi)星像元尺度的木本植物和草本植被覆蓋度,作為疏林草原木本植物和草本植物覆蓋度估算模型的訓練樣本;(2)提取衛(wèi)星影像的波段信息、植被指數(shù)以及紋理特征等參數(shù),作為疏林草原木本植物和草本植物覆蓋度估算模型的特征因子;(3)利用(1)和(2)的訓練樣本數(shù)據(jù)集,基于隨機森林算法,針對不同衛(wèi)星數(shù)據(jù)源訓練模型建立木本植物和草本植物覆蓋度估算模型;(4)利用已構(gòu)建的疏林草原木本植物和草本植物覆蓋度估算模型,基于高分辨率衛(wèi)星影像計算疏林草原木本植物和草本植物覆蓋度(圖3)。

    圖3 榆樹疏林草原木本植物和草本植物覆蓋度估算方法流程圖Fig.3 Flowchart of woody and herbaceous Fractional Vegetation Cover(FVC)estimation in Elm sparse forest grassland

    3.1.1 分類訓練樣本的獲取

    (1)基于無人機影像估算木本和草本植物覆蓋度。機器學習是一種通過數(shù)據(jù)挖掘方法找出隱含在預測變量和分類樣本間的相關關系,并自動建立預測模型的方法。分類與回歸樹(CART)算法是由Breiman 等(1984)提出的一種非參數(shù)的有效監(jiān)督學習方法,即從一系列有特征和標簽的數(shù)據(jù)中總結(jié)出決策規(guī)則,并用樹狀圖結(jié)構(gòu)來呈現(xiàn)這些規(guī)則以解決分類和回歸問題的算法,能夠有效利用空間輔助信息,深度挖掘特征信息。利用分類與回歸樹算法對無人機正射影像進行地物分類,基于色彩空間拓展函數(shù),將正射影像各像元的RGB 信息進行擴展,將色彩空間拓展到包含RGB(紅,綠,藍)、HSV(色調(diào),飽和度,明度)、Lab(亮度,坐標軸從紅色到綠色,坐標軸從藍色到紅色)和XYZ(CIE1931標準色度系統(tǒng))(Reinhard等,2001)共12 種特征作為CART 分類算法的特征因子(Wang等,2019)。

    利用無人機正射影像RGB 信息建立色彩空間拓展函數(shù)作為分類與回歸樹模型的特征因子;基于地面調(diào)查數(shù)據(jù)和目視解譯方法,使用圖像編輯軟件(如ImageJ,Photoshop 等)分別選取2020 年正射影像中的木本、草本和裸地,建立研究區(qū)地表類型的樣本訓練集作為決策樹分類模型的分類樣本,生成決策樹分類模型。將地物類型劃分為木本植物,草本植物和沙地3 類,得到樣地木本、草本植物與裸沙地的地表類型圖。通過計算不同地物類型的像素數(shù)占植被分類圖中總像素數(shù)比例,獲取木本植物和草本植物的植被覆蓋度。

    基于研究區(qū)無人機影像的地表類型分布圖(圖4(b)),利用GF-6衛(wèi)星影像和Sentinel-2 衛(wèi)星影像分別對木本植物和草本植物進行掩膜,統(tǒng)計落在衛(wèi)星像元窗口內(nèi)的木本植物和草本植物像元總數(shù),并計算其占衛(wèi)星窗口像元總數(shù)的百分比(以下簡稱覆蓋度分類樣本),得到衛(wèi)星像元尺度的木本和草本植物覆蓋度作為疏林草原覆蓋度模型的分類樣本數(shù)據(jù)。每組GF-6影像為3903個樣本數(shù)據(jù),Sentinel-2影像為2497個樣本數(shù)據(jù)。

    圖4 基于無人機正射影像的衛(wèi)星像元尺度木本和草本植物覆蓋度計算方法Fig.4 Methods of woody and herbaceous FVC at pixel-scale of satellite images calculated from DOM obtained by UAV

    (2)分類結(jié)果驗證。為驗證分類與回歸樹模型對無人機正射影像的分類精度,在無人機正射影像中隨機選取4000 個像素點,將目視判別結(jié)果與模型分類結(jié)果進行對比,基于像素點建立混淆矩陣,并計算Kappa系數(shù)和總體誤差作為指標評價決策樹分類精度。

    3.1.2 模型特征因子的選取

    植被指數(shù)對植被的生物物理特征較為敏感。對于地表空間復雜的影像,提取影像的紋理特征可以獲取各種有用的信息,提高分類精度(劉龍飛等,2003)。本實驗選擇的特征包括3 類,分別是波段反射率、植被指數(shù)以及紋理特征(Haralick,1979)。首先,選擇衛(wèi)星像元波段反射率作為模型特征因子(表2);植被指數(shù)包括歸一化植被指數(shù)NDVI(Normalized Difference Vegetation Index),差值植被指數(shù)DVI(Difference Vegetation Index),歸一化綠度植被指數(shù)GNDVI(Green Normalized Difference Vegetation Index),綠波比指數(shù)GRI(Green Ratio Vegetation Index),紅邊指數(shù)(MTCI)和增強型植被指數(shù)EVI(Enhanced Vegetation Index);采用基于統(tǒng)計分析的灰度共生矩陣GLCM(Gray-Level-Co-occurrence Matrix)實現(xiàn)GF-6 和Sentinel-2衛(wèi)星各波段紋理特征的計算。植被指數(shù)和紋理特征的提取在ENVI(V5.3)中完成。GF-6衛(wèi)星模型特征因子選取8 個波段反射率、6 個植被指數(shù)和48個紋理特征,共計62個特征因子,Sentinel-2衛(wèi)星選取10 個波段反射率、6 個植被指數(shù)和60 個紋理特征,共計76個特征因子。

    表2 榆樹疏林草原木本和草本植物覆蓋度估算模型特征因子與描述Table 2 Features and description of woody and herbaceous FVC estimation model in Elm sparse forest grassland

    3.1.3 疏林草原木本和草本植物覆蓋度估算模型建立

    將影像波段反射率、6 種紋理特征和6 種植被指數(shù)作為特征因子,利用隨機森林算法建立疏林草原木本植物和草本植物覆蓋度估算模型。隨機森林是一個以分類與回歸樹算法為基礎,將決策樹作為弱模型的袋裝法集成算法,即通過構(gòu)建多個相互獨立的訓練樣本子集訓練弱模型,然后對各預測結(jié)果取平均值作為最終預測結(jié)果。隨機森林在利用袋裝法對獨立樣本子集進行模型訓練時,對每個特征因子出現(xiàn)的次數(shù)、深度以及模型的準確率進行統(tǒng)計,并賦予每個特征一定的權(quán)重值,該權(quán)重值即為各特征因子的重要性(Breiman,2001)。

    與復雜模型相比,優(yōu)先選擇特征因子較少且預測結(jié)果相對準確的模型(James 等,2013)。由于并非全部特征變量都與覆蓋度相關,利用所有特征因子建立隨機森林模型后,對特征因子進行優(yōu)選,以提高隨機森林的預測精度并簡化節(jié)點。特征因子優(yōu)選首先根據(jù)重要性得分對各特征因子進行排序,然后按照重要性排名依次疊加特征因子并訓練模型,利用10 倍交叉檢驗(Crossvalidation)(Breiman,2001)分別計算一定特征因子數(shù)量下隨機森林模型的相關系數(shù)范圍,并進行隨機森林穩(wěn)健性估計,當相關系數(shù)趨于穩(wěn)定后,排除對模型貢獻較小的變量,建立最優(yōu)性能下的隨機森林模型。

    確定模型最佳特征因子后,通過10 倍交叉檢驗查看模型預測精度的變化趨勢,調(diào)整隨機森林模型中樹的數(shù)量和樹的深度等參數(shù)以獲取模型參數(shù)的最優(yōu)值,進一步優(yōu)化模型、提高模型預測精度。應用建立的模型對研究區(qū)木本植物和草本植物的覆蓋度進行估算。以上圖像處理、模型訓練、模型構(gòu)建和圖像分類均利用Python(V3.7)語言完成。

    3.2 疏林草原木本和草本植物覆蓋度估算模型驗證

    為防止由空間自相關引起的空間過度擬合,建模之前在研究區(qū)內(nèi)隨機且均勻選取25%的獨立樣本數(shù)據(jù)作為精度驗證數(shù)據(jù),以評估模型對覆蓋度的估算精度。基于GF-6影像的數(shù)據(jù)集包含978組獨立驗證數(shù)據(jù),基于Sentinel-2 影像的數(shù)據(jù)集包含626組獨立驗證數(shù)據(jù)。比較無人機正射影像分類得到的覆蓋度與模型估算得到的覆蓋度值以評估模型精度,選擇的精度驗證指標包括模型樣本的決定系數(shù)R2(式(1))、均方根誤差RMSE(式(2))和估算精度EA(式(3))。

    式中,yi為基于正射影像分類得到的覆蓋度,為模型估算的第i個覆蓋度,N為樣本數(shù)。

    4 結(jié)果與分析

    4.1 疏林草原木本和草本植物覆蓋度估算結(jié)果

    4.1.1 無人機正射影像分類

    通過比較研究區(qū)內(nèi)部樣方的正射影像和分類圖,以目視判別方法評估分類的準確性。由圖5可見,木本植物、草本植物和裸沙地分類效果良好,草地和裸沙地中稀疏分布的灌木均被完整的提取,木本植物與草本植物邊界清晰,木本植物內(nèi)部完整,準確排除樹冠間隙和樹冠陰影。目視判別與機器分類之間的誤差主要是由于草本植物與裸沙地的混合像元分類產(chǎn)生。表3為基于驗證像素點建立的木本植物,草本植物和裸沙地的分類結(jié)果混淆矩陣,決策樹分類總體精度為0.89,Kappa 系數(shù)為0.83,分類精度較高。生產(chǎn)者精度和用戶精度的準確性充分證明了基于分類與回歸樹算法的影像分類能夠準確劃分無人機正射影像像元類型。

    表3 榆樹疏林草原長期定位觀測大樣地無人機正射影像像元分類精度評價Table 3 Accuracy evaluation in pixel-based classification of DOM in long-term monitoring plot in Elm sparse forest

    4.1.2 模型特征因子選擇

    不同特征因子數(shù)量下的隨機森林模型性能驗證結(jié)果如圖所示(圖6)。基于GF-6影像的木本植物覆蓋度估算模型在特征因子數(shù)量為14 時,模型趨于穩(wěn)健且預測精度較高,可作為平衡預測精度和簡化模型的最優(yōu)特征因子個數(shù),建立最優(yōu)性能的覆蓋度估算模型?;贕F-6 影像的草本植物覆蓋度模型最優(yōu)特征因子數(shù)量為17?;赟entinel-2影像的木本和草本植物覆蓋度模型最優(yōu)特征因子數(shù)量分別為11和17。

    圖6 不同特征因子數(shù)量下的隨機森林模型性能驗證結(jié)果Fig.6 Random forest model performance verification results under different number of feature factors

    表4 和表5 分別統(tǒng)計了基于GF-6 和Sentinel-2影像的疏林草原木本植物和草本植物覆蓋度估算模型重要性排名前5 的特征因子。基于GF-6 影像的疏林草原木本植物包含DVI、EVI、NDVI、MTCI和GRVI,均為植被指數(shù)(表4);基于Sentinel-2影像的疏林草原木本植物覆蓋度估算模型特征因子包括了NDVI、GNDVI、GRVI這3個植被指數(shù)、以及紅波段和植被紅邊1 波段?;贕F-6 的草本植物覆蓋度估算模型特征因子包含黃色波段、紅邊2波段、植被指數(shù)MTCI、紅邊2 波段紋理信息的均值和紅邊2波段紋理信息的方差;基于Sentinel-2的估算模型特征因子包含綠色波段、植被紅邊1、2波段和短波紅外1波段,紋理信息為短波紅外的均值。

    表4 基于GF-6影像的榆樹疏林草原木本和草本植物覆蓋度估算模型特征因子Table 4 Top five preferred features ranking of woody and herbaceous FVC estimation model in Elm sparse forest grassland based on GF-6 image

    表5 基于Sentinel-2影像的榆樹疏林草原木本和草本植物覆蓋度估算模型特征因子Table 5 Top five preferred features ranking of woody and herbaceous FVC estimation model in Elm sparse forest grassland based on Sentinel-2 image

    DVI 是基于GF-6 影像的木本植物覆蓋度估算模型中最重要的特征因子,而NDVI則是Sentinel-2影像的木本植物覆蓋度估算模型中最重要的特征因子。對兩個模型特征因子綜合分析可知,植被指數(shù)均對木本植物覆蓋度的估算起重要作用,基于Sentinel-2 影像的草本植物估算模型最佳預測因子為中心波長在705 nm 的植被紅邊1 波段,重要性達到0.60(表5),這與基于GF-6 的草本植物覆蓋度模型相似,除黃色波段外,模型特征因子的光譜特征和紋理特征同樣集中在紅邊波段(表4)。由此可知,GF-6 影像和Sentinel-2 影像在730—750 nm 處的反射率能有效反映草本植物植被覆蓋度。說明紅邊波段反射率對草本植物覆蓋度變化較為敏感。

    4.1.3 模型精度驗證

    將無人機影像分類得到的木本、草本植物覆蓋度樣本值與模型估算值進行對比(圖7)。基于GF-6 影像和Sentinel-2 影像估算的木本植物覆蓋度R2均為0.72,均方根誤差為6.76%和6.53%,估算精度分別為46.31%和54.30%?;贕F-6 影像比基于Sentinel-2 影像的木本植物覆蓋度估算值均方根誤差相比高0.23%,估算精低7.99%?;贕F-6影像和Sentinel-2影像估算的草本植物覆蓋度R2分別為0.66 和0.81,均方根誤差分別為10.69%和8.20%,估算精度分別為77.88%和83.17%。基于GF-6 影像比基于Sentinel-2 影像估算的草本植物覆蓋度R2低0.15,均方根誤差高2.49%,估算精度低5.29%。

    圖7 基于衛(wèi)星影像與基于無人機正射影像估算的榆樹疏林草原長期定位觀測大樣地木本和草本植物覆蓋度結(jié)果比較Fig.7 Comparison of woody and herbaceous FVC in Elm sparse forest grassland estimated by satellite images and DOM of UAV

    兩種衛(wèi)星影像對疏林草原木本植物覆蓋度估算精度總體差別較小,對草本植物覆蓋度估算精度Sentinel-2 影像略高于GF-6 影像。不論基于哪種衛(wèi)星影像,木本植物在0—20%區(qū)間的覆蓋度估算精度較高,而在20%—100%區(qū)間估算精度較低,原因可能由于木本植物覆蓋度在20%—100%區(qū)間的訓練樣本數(shù)量較少。

    4.2 研究區(qū)覆蓋度分析

    分別基于GF-6 影像和Sentinel-2 影像對研究區(qū)木本、草本植物覆蓋度進行制圖,并與基于無人機正射影像得到的木本、草本植物覆蓋度圖對比,如圖8所示。總體來看,基于高分辨率衛(wèi)星影像獲取的木本、草本植物覆蓋度,與基于無人機正射影像獲取的木本、草本植物覆蓋度空間分布一致?;贕F-6 衛(wèi)星像元尺度對正射影像分類結(jié)果進行統(tǒng)計,約3/4 的木本植物像元落在0—20%覆蓋度內(nèi)、草本植物像元落在20%—80%覆蓋度范圍內(nèi);基于Sentinel-2 衛(wèi)星像元尺度對正射影像分類結(jié)果進行統(tǒng)計,約4/5 的木本植物像元落在0—20%覆蓋度內(nèi)、草本植物像元落在20%—80%覆蓋度范圍內(nèi)。基于GF-6 影像估算得到的覆蓋度3/4 的木本植物集中在0—20%,草本植物在25%—75%;基于Sentinel-2影像得到的覆蓋度4/5的木本植物集中在0—20%,草本植物在25%—75%。

    圖8 基于高分衛(wèi)星影像和無人機正射影像的榆樹疏林草原長期定位觀測大樣地木本和草本植物覆蓋度空間分布圖Fig.8 Maps of woody and herbaceous FVC in Elm sparse forest grassland from satellite image(GF-6 and Sentinel-2)and classified DOM

    研究區(qū)木本植物覆蓋度在50%以上僅占研究區(qū)總像元數(shù)的2%—4%,且集中在研究區(qū)西部地勢較低的區(qū)域。研究區(qū)草本植物覆蓋度分布較為均衡,覆蓋度高于70%的區(qū)域集中分布在研究區(qū)的東部,覆蓋度低于30%的區(qū)域集中分布在地勢較高的沙丘區(qū)域。

    5 結(jié)論

    本研究探索了一種利用近地遙感和衛(wèi)星遙感影像,基于機器學習算法,估算疏林草原木本植物和草本植物覆蓋度的新方法。研究不僅實現(xiàn)了基于無人機正射影像不同地物類別的精準分類,獲取了清晰的不同植被類型生活型信息。同時利用重采樣升尺度的方法分別提取衛(wèi)星影像混合像元內(nèi)的疏林草原木本、草本植物覆蓋度,構(gòu)建了疏林草原木本、草本植物覆蓋度估算模型,實現(xiàn)了從無人機到衛(wèi)星的疏林草原木本植物和草本植被的覆蓋度精準估算。研究結(jié)果表明:(1)利用無人機遙感監(jiān)測平臺能夠?qū)崿F(xiàn)疏林草原景觀尺度的精準植被分類與覆蓋度計算;(2)構(gòu)建了基于GF-6 衛(wèi)星影像和Sentinel-2 衛(wèi)星影像的木本和草本植物覆蓋度估算模型,實現(xiàn)了基于高分辨率衛(wèi)星影像的疏林草原木本和草本植物覆蓋度估算;(3)基于Sentinel-2 衛(wèi)星影像的疏林草原木本和草本植物覆蓋度估算精度稍高于GF-6 衛(wèi)星,基于兩個衛(wèi)星影像的草本植物覆蓋度的估算精度都要顯著高于木本植物。

    近年來,覆蓋全球的高分辨率衛(wèi)星不斷涌現(xiàn),通過融合不同分辨率的多源衛(wèi)星遙感數(shù)據(jù),基于定量遙感尺度轉(zhuǎn)換模型,監(jiān)測熱帶地區(qū)稀樹草原上木本植被狀況和動態(tài)是稀樹草原研究的熱點(Brandt等,2016;Tian等,2017)。González-Roglich和Swenson(2016)利用中巴地球衛(wèi)星(CBERS)影像與陸地衛(wèi)星(Landsat)影像,基于隨機森林算法估計阿根廷中部稀樹草原樹木覆蓋度。Herrmann等(2013)在西非塞內(nèi)加爾使用軌道觀測衛(wèi)星(Orbview)高分辨率影像與中分辨率成像光譜儀(MODIS)的NDVI數(shù)據(jù),基于決策樹算法建立模型并反演木本覆蓋度。Higginbottom 等(2018)融合Landsat 影像與日本對地觀測衛(wèi)星(ALOS)高分辨率合成孔徑雷達數(shù)據(jù),基于隨機森林算法繪制出南非林波波省疏林草原地區(qū)木本植物覆蓋度圖。以上研究均應用機器學習方法結(jié)合多分辨率遙感影像實現(xiàn)了熱帶稀樹草原木本植物覆蓋度估算。中國溫帶疏林草原景觀與熱帶稀樹草原特征相似,在連續(xù)的草本層里鑲嵌有不連續(xù)的喬木和灌木植物,是中國半干旱區(qū)重要的生態(tài)系統(tǒng),也是世界上唯一分布于溫帶的疏林草原生態(tài)系統(tǒng)(于順利 等,2007)。本研究探索了將無人機近地遙感數(shù)據(jù)與高分辨率衛(wèi)星遙感數(shù)據(jù)相結(jié)合的新方法,以分布于中國內(nèi)蒙古高原上的溫帶疏林草原為對象,實現(xiàn)了疏林草原上木本植物和草本植物的同步監(jiān)測。新方法為實現(xiàn)區(qū)域尺度溫帶疏林草原不同生活型植物的高效精準監(jiān)測提供了方法基礎。

    本研究結(jié)果與同樣利用隨機森林方法,結(jié)合航空遙感數(shù)據(jù)與衛(wèi)星數(shù)據(jù),估算熱帶稀樹草原木本植物覆蓋度的驗證結(jié)果相比(Ludwig等,2019),具有更高的決定系數(shù)和更低的均方根誤差。Ludwig估測得到的木本植物覆蓋度估算模型決定系數(shù)R2為0.44,均方根誤差為11.6%;基于GF-6 影像和Sentinel-2 影像估算的木本植物覆蓋度估算模型決定系數(shù)R2均為0.72,均方根誤差為6.76%和6.53%。此外,基于隨機森林建立的植被覆蓋度模型預測值均高估了較低的植被覆蓋度值,而低估了高植被覆蓋度值,與隨機森林算法在預測結(jié)果時,對每棵回歸樹得到的預測變量求平均值作為預測結(jié)果的屬性相關(Prasad 等,2006)。隨機森林該屬性導致在覆蓋度預測中對于較低和較高的木本植物和草本植被覆蓋度估算誤差較大。

    Higginbottom 等(2018)融合了Landsat衛(wèi)星數(shù)據(jù)和PALSAR 合成孔徑雷達數(shù)據(jù),并通過組合旱季、雨季多季節(jié)指標訓練模型。得到的木本植物覆蓋度驗證結(jié)果中,決定系數(shù)最高達0.77。該結(jié)論對模型的改進具有一定的啟發(fā)性,可以結(jié)合輔助數(shù)據(jù)并加入生長季的不同時期影像作為預測指標,進一步提高疏林草原木本植物和草本植物覆蓋度估算精度。

    基于機器學習算法,融合無人機正射影像和高分辨率衛(wèi)星影像等多源數(shù)據(jù),充分發(fā)揮了近地遙感高空間分辨率和衛(wèi)星遙感區(qū)域大尺度覆蓋的特點。無人機正射影像為疏林草原木本、草本植被覆蓋度估算模型提供了高質(zhì)量的訓練數(shù)據(jù),衛(wèi)星影像為探測區(qū)域尺度植被狀況提供了基礎數(shù)據(jù)源,機器學習的方法為兩種類型數(shù)據(jù)的融合搭建了橋梁。利用這一方法,未來可以實現(xiàn)中國溫帶疏林草原分布范圍的精確制圖,進一步利用多時相的衛(wèi)星遙感數(shù)據(jù)也可監(jiān)測疏林草原上木本植物和草本植物的長期變化特征。

    猜你喜歡
    特征模型
    一半模型
    抓住特征巧觀察
    重要模型『一線三等角』
    新型冠狀病毒及其流行病學特征認識
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    如何表達“特征”
    不忠誠的四個特征
    當代陜西(2019年10期)2019-06-03 10:12:04
    抓住特征巧觀察
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    美女高潮到喷水免费观看| 亚洲熟妇熟女久久| 一级片免费观看大全| 日韩一卡2卡3卡4卡2021年| 精品国产乱子伦一区二区三区| 一级黄色大片毛片| 国产真人三级小视频在线观看| 国产av精品麻豆| 亚洲色图 男人天堂 中文字幕| 在线播放国产精品三级| 777久久人妻少妇嫩草av网站| 精品久久久久久,| 淫秽高清视频在线观看| 午夜精品在线福利| 日日爽夜夜爽网站| 国产精品亚洲一级av第二区| 不卡一级毛片| 1024香蕉在线观看| av有码第一页| 国产成人影院久久av| 亚洲国产中文字幕在线视频| 国产精品免费视频内射| 1024香蕉在线观看| 精品欧美国产一区二区三| 级片在线观看| 激情视频va一区二区三区| 两个人免费观看高清视频| 国产高清有码在线观看视频 | 国产又爽黄色视频| 国产野战对白在线观看| aaaaa片日本免费| 18美女黄网站色大片免费观看| 国产精品久久久久久人妻精品电影| 亚洲av美国av| 美女高潮到喷水免费观看| 国产精品一区二区三区四区久久 | 国产精品亚洲美女久久久| 欧美不卡视频在线免费观看 | 亚洲成a人片在线一区二区| 国产精品秋霞免费鲁丝片| 首页视频小说图片口味搜索| 国产麻豆成人av免费视频| 非洲黑人性xxxx精品又粗又长| 国产区一区二久久| 午夜久久久久精精品| 精品久久久久久久人妻蜜臀av | 国产精品,欧美在线| 叶爱在线成人免费视频播放| 两性夫妻黄色片| 欧美一级毛片孕妇| 50天的宝宝边吃奶边哭怎么回事| 曰老女人黄片| 村上凉子中文字幕在线| 十八禁网站免费在线| 亚洲五月天丁香| 18禁黄网站禁片午夜丰满| 国产精品 欧美亚洲| 国内精品久久久久久久电影| 极品人妻少妇av视频| 久久香蕉国产精品| 母亲3免费完整高清在线观看| 欧美激情高清一区二区三区| 变态另类成人亚洲欧美熟女 | 国产精品久久久av美女十八| 久久精品国产清高在天天线| 国产伦人伦偷精品视频| 色精品久久人妻99蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利成人在线免费观看| 日日摸夜夜添夜夜添小说| 老熟妇乱子伦视频在线观看| 日本免费一区二区三区高清不卡 | 国内久久婷婷六月综合欲色啪| 丁香六月欧美| 亚洲欧美精品综合一区二区三区| 婷婷精品国产亚洲av在线| 亚洲,欧美精品.| 香蕉丝袜av| 一区在线观看完整版| 欧美一级毛片孕妇| 精品国产国语对白av| 亚洲成人久久性| 婷婷六月久久综合丁香| 亚洲国产看品久久| 日韩欧美三级三区| 久久久久久国产a免费观看| 国产精品日韩av在线免费观看 | 热re99久久国产66热| 国产成+人综合+亚洲专区| 黄色成人免费大全| 久久久久久人人人人人| 美女大奶头视频| 国产精品二区激情视频| 国产av精品麻豆| 国内毛片毛片毛片毛片毛片| 久久精品成人免费网站| 国产熟女午夜一区二区三区| 欧美中文日本在线观看视频| 黑人欧美特级aaaaaa片| 中文字幕精品免费在线观看视频| 老司机靠b影院| 日韩三级视频一区二区三区| 国产麻豆69| 久久国产亚洲av麻豆专区| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久久精品吃奶| 欧美成人性av电影在线观看| 久久久久久久午夜电影| 男女下面插进去视频免费观看| 99riav亚洲国产免费| 女警被强在线播放| 亚洲欧美精品综合久久99| 多毛熟女@视频| 国产成人精品久久二区二区免费| 欧美乱色亚洲激情| 久热爱精品视频在线9| 老司机靠b影院| 波多野结衣av一区二区av| 亚洲精品美女久久久久99蜜臀| av福利片在线| 岛国在线观看网站| 亚洲精品美女久久久久99蜜臀| 国产欧美日韩综合在线一区二区| 91九色精品人成在线观看| 精品国产一区二区三区四区第35| 亚洲色图 男人天堂 中文字幕| 国产成人精品久久二区二区免费| 国产av一区二区精品久久| 成年人黄色毛片网站| 精品久久久久久久人妻蜜臀av | 国产欧美日韩一区二区精品| 在线国产一区二区在线| 在线永久观看黄色视频| 老司机午夜十八禁免费视频| 国产成+人综合+亚洲专区| 亚洲精品久久成人aⅴ小说| 一本久久中文字幕| 香蕉久久夜色| av天堂久久9| 国产一区二区激情短视频| 国产高清有码在线观看视频 | 老司机午夜十八禁免费视频| 亚洲欧美精品综合久久99| 欧美乱色亚洲激情| 久久久久亚洲av毛片大全| 一区二区日韩欧美中文字幕| 亚洲成av片中文字幕在线观看| 性色av乱码一区二区三区2| 成人手机av| 色av中文字幕| 久久国产精品人妻蜜桃| 又大又爽又粗| 97超级碰碰碰精品色视频在线观看| 色尼玛亚洲综合影院| 国产99久久九九免费精品| 热re99久久国产66热| 少妇 在线观看| 成人av一区二区三区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 岛国视频午夜一区免费看| 91av网站免费观看| 99国产精品一区二区蜜桃av| 无人区码免费观看不卡| xxx96com| 女人爽到高潮嗷嗷叫在线视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲电影在线观看av| 日韩免费av在线播放| 久久中文看片网| 成人免费观看视频高清| 一级片免费观看大全| 亚洲成a人片在线一区二区| 国产日韩一区二区三区精品不卡| 一级毛片精品| 成年版毛片免费区| 神马国产精品三级电影在线观看 | 淫妇啪啪啪对白视频| 欧美不卡视频在线免费观看 | av天堂久久9| 国产精品亚洲美女久久久| 色综合亚洲欧美另类图片| 18禁黄网站禁片午夜丰满| 欧美日韩瑟瑟在线播放| 亚洲熟妇熟女久久| svipshipincom国产片| 国产精品综合久久久久久久免费 | 一区二区三区国产精品乱码| 18禁国产床啪视频网站| 成年女人毛片免费观看观看9| 在线观看免费视频日本深夜| 一夜夜www| 91老司机精品| 亚洲片人在线观看| 日本撒尿小便嘘嘘汇集6| 成人国产综合亚洲| 俄罗斯特黄特色一大片| 一级毛片女人18水好多| avwww免费| 免费在线观看黄色视频的| 成人国产综合亚洲| av在线天堂中文字幕| 中文字幕最新亚洲高清| 97人妻精品一区二区三区麻豆 | 久久这里只有精品19| 9热在线视频观看99| 国产一区二区三区综合在线观看| 国产成人啪精品午夜网站| 99久久久亚洲精品蜜臀av| 最近最新中文字幕大全免费视频| 国产aⅴ精品一区二区三区波| 日韩欧美国产一区二区入口| 香蕉国产在线看| 免费在线观看视频国产中文字幕亚洲| 十八禁网站免费在线| 亚洲欧美精品综合久久99| 老汉色∧v一级毛片| 日本a在线网址| 50天的宝宝边吃奶边哭怎么回事| av视频免费观看在线观看| 国产精品秋霞免费鲁丝片| 精品不卡国产一区二区三区| 精品国产美女av久久久久小说| 久久久久久久午夜电影| 悠悠久久av| 日韩欧美国产一区二区入口| 国产精品,欧美在线| 男女之事视频高清在线观看| 国产精品久久久久久人妻精品电影| tocl精华| 国产精品久久电影中文字幕| 99久久99久久久精品蜜桃| 亚洲成人久久性| 精品人妻在线不人妻| 少妇裸体淫交视频免费看高清 | 欧美绝顶高潮抽搐喷水| 成人av一区二区三区在线看| 免费高清视频大片| 波多野结衣一区麻豆| 怎么达到女性高潮| 少妇粗大呻吟视频| 国产一区二区三区综合在线观看| 真人一进一出gif抽搐免费| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| 国产精品亚洲一级av第二区| 亚洲av美国av| 电影成人av| 日日爽夜夜爽网站| 国产亚洲av高清不卡| 国产一卡二卡三卡精品| 国产色视频综合| 国产av一区在线观看免费| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕一级| 国产激情久久老熟女| 日日爽夜夜爽网站| 国产精品久久久av美女十八| videosex国产| 久久久精品国产亚洲av高清涩受| 一卡2卡三卡四卡精品乱码亚洲| 美女高潮到喷水免费观看| 亚洲成av人片免费观看| 首页视频小说图片口味搜索| 久久久久久久久久久久大奶| 久久九九热精品免费| 国产精品电影一区二区三区| 欧美日韩黄片免| 亚洲精品久久成人aⅴ小说| 欧美激情高清一区二区三区| 成人国语在线视频| 精品高清国产在线一区| 老司机靠b影院| 少妇被粗大的猛进出69影院| 精品免费久久久久久久清纯| 性少妇av在线| 欧美一级a爱片免费观看看 | 成人手机av| 99国产精品一区二区三区| 国产一区二区在线av高清观看| 伦理电影免费视频| 一a级毛片在线观看| 精品国产一区二区久久| 久久久久久人人人人人| 久久香蕉激情| 国产一区二区三区视频了| 国产成人精品久久二区二区91| 久久影院123| 亚洲av五月六月丁香网| 午夜福利18| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频| 99国产精品一区二区三区| 国产亚洲精品第一综合不卡| 欧美午夜高清在线| 久久久久久久精品吃奶| 亚洲国产精品久久男人天堂| 亚洲人成网站在线播放欧美日韩| 美女高潮喷水抽搐中文字幕| 在线观看免费日韩欧美大片| 麻豆一二三区av精品| 国产又爽黄色视频| 淫秽高清视频在线观看| 桃色一区二区三区在线观看| 欧美乱码精品一区二区三区| 国产野战对白在线观看| 国产成人精品久久二区二区91| 最新美女视频免费是黄的| 亚洲精品国产区一区二| 欧美黄色片欧美黄色片| 欧美人与性动交α欧美精品济南到| 亚洲在线自拍视频| 国产午夜福利久久久久久| 国产亚洲欧美在线一区二区| 人人妻人人澡人人看| 日日爽夜夜爽网站| 韩国精品一区二区三区| 国产精品九九99| 精品一区二区三区四区五区乱码| 在线观看舔阴道视频| 国产真人三级小视频在线观看| 午夜福利免费观看在线| 黄色片一级片一级黄色片| 亚洲国产欧美网| 久久久久久人人人人人| 嫩草影院精品99| 夜夜夜夜夜久久久久| 午夜福利影视在线免费观看| 亚洲欧美日韩无卡精品| 亚洲欧美一区二区三区黑人| 日本三级黄在线观看| 国产精品久久久av美女十八| 婷婷丁香在线五月| 国产精品香港三级国产av潘金莲| 成人三级做爰电影| 一区福利在线观看| 多毛熟女@视频| 欧美日本亚洲视频在线播放| 亚洲成国产人片在线观看| 国产男靠女视频免费网站| 中文字幕人妻熟女乱码| 一进一出抽搐gif免费好疼| 男人的好看免费观看在线视频 | 精品卡一卡二卡四卡免费| 在线天堂中文资源库| 国产成人av激情在线播放| 日韩欧美三级三区| 久久午夜综合久久蜜桃| 岛国视频午夜一区免费看| 最近最新中文字幕大全电影3 | 黑丝袜美女国产一区| 黄色毛片三级朝国网站| 美女高潮喷水抽搐中文字幕| 色综合欧美亚洲国产小说| 精品电影一区二区在线| 超碰成人久久| av天堂在线播放| 久久亚洲精品不卡| 99riav亚洲国产免费| 最近最新中文字幕大全电影3 | 亚洲欧洲精品一区二区精品久久久| 午夜视频精品福利| 俄罗斯特黄特色一大片| 成人18禁在线播放| 国产精品亚洲一级av第二区| 女性被躁到高潮视频| 亚洲国产精品999在线| 长腿黑丝高跟| 亚洲成a人片在线一区二区| 欧美日韩精品网址| 午夜免费观看网址| 精品无人区乱码1区二区| netflix在线观看网站| 亚洲国产精品sss在线观看| 亚洲国产精品成人综合色| 欧美一级a爱片免费观看看 | 1024视频免费在线观看| 最新美女视频免费是黄的| 妹子高潮喷水视频| 视频区欧美日本亚洲| 国产亚洲欧美精品永久| 日韩免费av在线播放| 国产成人精品无人区| xxx96com| 精品电影一区二区在线| 成人免费观看视频高清| 国产精品永久免费网站| 日本五十路高清| 青草久久国产| 亚洲va日本ⅴa欧美va伊人久久| 色综合站精品国产| 日韩欧美国产在线观看| 男女午夜视频在线观看| a在线观看视频网站| 一区二区三区国产精品乱码| 欧美日本亚洲视频在线播放| 午夜福利视频1000在线观看 | 制服丝袜大香蕉在线| 一级毛片女人18水好多| 1024香蕉在线观看| 欧美国产日韩亚洲一区| 国产主播在线观看一区二区| 性少妇av在线| 在线十欧美十亚洲十日本专区| 亚洲激情在线av| 国内精品久久久久久久电影| 国产免费男女视频| 亚洲片人在线观看| 在线观看午夜福利视频| 国产精品野战在线观看| 啦啦啦免费观看视频1| 亚洲七黄色美女视频| 九色国产91popny在线| 嫩草影视91久久| 亚洲国产欧美一区二区综合| 欧美久久黑人一区二区| 国产一级毛片七仙女欲春2 | 黑丝袜美女国产一区| 国产精品永久免费网站| 1024香蕉在线观看| 日日摸夜夜添夜夜添小说| 亚洲人成77777在线视频| 国产精品亚洲美女久久久| 亚洲av熟女| 欧美精品啪啪一区二区三区| 欧美成人一区二区免费高清观看 | 桃色一区二区三区在线观看| 男女下面进入的视频免费午夜 | 午夜日韩欧美国产| 免费看美女性在线毛片视频| 欧美日韩精品网址| 久久午夜亚洲精品久久| av天堂在线播放| 精品欧美一区二区三区在线| 午夜免费观看网址| 国产精品野战在线观看| 一夜夜www| 一二三四在线观看免费中文在| 久久性视频一级片| 久久欧美精品欧美久久欧美| 精品午夜福利视频在线观看一区| 18禁国产床啪视频网站| 波多野结衣巨乳人妻| 国产亚洲精品久久久久5区| 亚洲一区中文字幕在线| 嫁个100分男人电影在线观看| 丝袜在线中文字幕| 夜夜夜夜夜久久久久| 一a级毛片在线观看| 国产人伦9x9x在线观看| 亚洲精品国产精品久久久不卡| 高清在线国产一区| 亚洲中文日韩欧美视频| av在线播放免费不卡| 欧美色欧美亚洲另类二区 | 51午夜福利影视在线观看| 国产99久久九九免费精品| 午夜福利欧美成人| 可以在线观看的亚洲视频| 亚洲精品中文字幕在线视频| 在线观看午夜福利视频| 91九色精品人成在线观看| www.999成人在线观看| 国产成人啪精品午夜网站| 日韩欧美国产一区二区入口| 久久精品成人免费网站| 超碰成人久久| 中文字幕另类日韩欧美亚洲嫩草| 黑人操中国人逼视频| 亚洲人成77777在线视频| 波多野结衣巨乳人妻| 亚洲精品粉嫩美女一区| 精品国产乱子伦一区二区三区| 亚洲,欧美精品.| 法律面前人人平等表现在哪些方面| ponron亚洲| 日韩高清综合在线| 老司机在亚洲福利影院| 一夜夜www| 咕卡用的链子| 久久久久久大精品| 免费在线观看日本一区| 国产欧美日韩一区二区三| √禁漫天堂资源中文www| 99re在线观看精品视频| 十分钟在线观看高清视频www| 男女下面进入的视频免费午夜 | 亚洲国产欧美网| 日韩av在线大香蕉| 国产高清视频在线播放一区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美黑人精品巨大| 亚洲国产精品久久男人天堂| 日韩欧美国产一区二区入口| 91成人精品电影| 欧美日韩一级在线毛片| 久久婷婷成人综合色麻豆| 香蕉国产在线看| 国产一级毛片七仙女欲春2 | 日日干狠狠操夜夜爽| 欧美乱色亚洲激情| 日本免费a在线| 又黄又粗又硬又大视频| 免费女性裸体啪啪无遮挡网站| 午夜福利影视在线免费观看| 国内毛片毛片毛片毛片毛片| 看黄色毛片网站| 日本a在线网址| 国产午夜福利久久久久久| 亚洲免费av在线视频| 1024香蕉在线观看| 咕卡用的链子| 嫁个100分男人电影在线观看| 国产精品影院久久| 级片在线观看| 国产成人精品在线电影| 国产精品98久久久久久宅男小说| 国产成人av教育| 一区二区日韩欧美中文字幕| 真人做人爱边吃奶动态| a级毛片在线看网站| 日本免费一区二区三区高清不卡 | 国产免费av片在线观看野外av| 看片在线看免费视频| 国产精品久久久久久人妻精品电影| 99久久久亚洲精品蜜臀av| 99久久综合精品五月天人人| 亚洲精品美女久久av网站| 欧美日韩瑟瑟在线播放| 中文字幕人妻熟女乱码| 亚洲精品久久成人aⅴ小说| 一级毛片精品| av超薄肉色丝袜交足视频| 18禁观看日本| 精品久久久久久成人av| 亚洲中文av在线| av免费在线观看网站| 国产成人一区二区三区免费视频网站| 午夜福利一区二区在线看| 中文字幕高清在线视频| 免费看十八禁软件| 男女午夜视频在线观看| 国产在线精品亚洲第一网站| 18禁美女被吸乳视频| 老司机靠b影院| 成人亚洲精品av一区二区| 欧美激情极品国产一区二区三区| 精品福利观看| 老司机午夜十八禁免费视频| 亚洲第一青青草原| 中文字幕精品免费在线观看视频| tocl精华| 一区福利在线观看| 午夜日韩欧美国产| 桃红色精品国产亚洲av| 国产在线观看jvid| 丁香六月欧美| 1024视频免费在线观看| 国产亚洲精品第一综合不卡| 我的亚洲天堂| av视频在线观看入口| 久久久久久人人人人人| 国内久久婷婷六月综合欲色啪| 亚洲狠狠婷婷综合久久图片| 日本黄色视频三级网站网址| 亚洲 国产 在线| 日本三级黄在线观看| 日韩视频一区二区在线观看| 日韩国内少妇激情av| 18美女黄网站色大片免费观看| 一边摸一边抽搐一进一出视频| 亚洲成人久久性| 一区在线观看完整版| 91字幕亚洲| 99久久综合精品五月天人人| 久久中文看片网| 看免费av毛片| 亚洲精品粉嫩美女一区| 久久国产亚洲av麻豆专区| 一个人免费在线观看的高清视频| 在线国产一区二区在线| 国产一区二区三区在线臀色熟女| 国产亚洲精品综合一区在线观看 | 中文字幕久久专区| 淫秽高清视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | svipshipincom国产片| 国产91精品成人一区二区三区| 中文字幕久久专区| 亚洲av电影在线进入| 国产精品免费视频内射| 女人被狂操c到高潮| 国产黄a三级三级三级人| 在线观看66精品国产| 国产成人欧美在线观看| 美女高潮到喷水免费观看| 国产亚洲精品久久久久久毛片| 久久久久国产一级毛片高清牌| 色综合站精品国产| 亚洲电影在线观看av| 熟女少妇亚洲综合色aaa.| 日韩视频一区二区在线观看| 亚洲精品中文字幕在线视频| 欧洲精品卡2卡3卡4卡5卡区| 国内毛片毛片毛片毛片毛片| 一进一出抽搐gif免费好疼| 国产av一区在线观看免费| 免费女性裸体啪啪无遮挡网站| 午夜福利高清视频| 最新美女视频免费是黄的| 两个人视频免费观看高清| 97碰自拍视频| 夜夜躁狠狠躁天天躁| 色精品久久人妻99蜜桃| 国产成人系列免费观看| 久久久久久久久中文|