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

    基于無人機(jī)獲取草地植物地上生物量的信息損失分析及修正

    2025-08-03 00:00:00李嘉欣靳瑰麗劉文昊王生菊陳夢(mèng)甜李文雄杜玟霖李超胡秀雯
    草地學(xué)報(bào) 2025年7期
    關(guān)鍵詞:植被指數(shù)波段修正

    中圖分類號(hào):TP751;S812 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1007-0435(2025)07-2206-13

    引用格式:,等.基于無人機(jī)獲取草地植物地上生物量的信息損失分析及修正[J].草地學(xué)報(bào),2025,33 (7) :2206-2218 LI Jia-xin,JINGui-li,LIU Wen-hao,et al.Analyzing and Modifyingof Grassland Plant InformationLoss ofAboveg round Biomass based onUAV Data[J].Acta Agrestia Sinica,2025,33(7) :2206-2218

    Analyzingand ModifyingofGrassland PlantInformation Lossof AbovegroundBiomassbasedonUAVData

    LI Jia-xin,JIN Gui-li*,LIU Wen-hao,WANG Sheng-ju,CHEN Meng-tian,LI Wen-xiong, DU Wen-lin,LI Chao,HU Xiu-wen

    (Colege ofGrassand Science,Xinjiang Agricultural University/Xinjiang KeyLaboratoryof GrasslandResources andEcology, Urumqi,Xinjiang 830005,China)

    Abstract:In order to explore the information lossin UAV(Unmanned aerial vehicle) monitoring of grassland biomass,this study compared the ground hyperspectral and on-site measured data with the multispectral images of the desert grassland community of Seriphidium transiliense,and used the vegetation index combined with linear regresson method to construct the inversion model of aboveground biomass of the plant community at different phenological periods.The error and rule of UAV acquisition of aboveground biomass were discussed and the modified model was established.The result shows that:(1)In April,June and September,the best UAV multi-spectral models were NDVI(Normalized diference vegetation index),RVI(Ration vegetation index) and NDVI,with accuracy of 65.61% , 48.26% and 61.59% ,respectively. Among all the ground hyperspectral data, NDVI is the most optimal, and the accuracy was 71. 77% , 53.63% and 67.58% ,respectively.(2) The ground hyperspectral retrieval ability was better than the UAV multispectral retrieval ability,the loss rates ofmultispectral data from UAV at each phenological stagewere ΔNDVIgt;ΔRVIgt;I OVI(Differencevegetation index),and the loss rates of optimal model were 8.58% , 9.89% and 8.87% ,respectively.(3)The UAV multi-spectral accuracy of all models was improved by 0.02%~3.74% ,The modified NDVI mode in April had the highest accuracy of 66.36%(y=403.431x+17.5936) .In summary,by screening the best inversion models,analyzing the information loss of diferent platforms,and modifying the UAV multi-spectral of the UAV can give full play to the advantages of diferent platforms and improve the estimation accuracy.It also proves the feasibility of using high-resolution remote sensing data to modify low-spatial-resolution data.

    Keywords: Unmanned aerial vehicle;Ground hyperspectral;Aboveground;Inversion;Informationloss;Modification

    地上生物量(Abovegroundbiomass,AGB)是反映草地生長狀況的關(guān)鍵指標(biāo),準(zhǔn)確高效獲取地上生物量信息對(duì)草地資源管理至關(guān)重要[1-2]。目前多采用地面測量和遙感監(jiān)測估算地上生物量,地面測量獲取數(shù)據(jù)信息準(zhǔn)確有效,但費(fèi)時(shí)費(fèi)力,受時(shí)間、天氣、交通等因素影響難以實(shí)現(xiàn)區(qū)域性動(dòng)態(tài)監(jiān)測[3]。光學(xué)遙感為大面積高效監(jiān)測植物群落結(jié)構(gòu)特征提供了可能。無人機(jī)(Unmannedaerialvehicle,UAV能夠捕獲地物的光譜特征,常用于植物冠層群落特征提取[4-5],其搭載傳感器具有靈活、高效、便捷的特點(diǎn),克服了衛(wèi)星和航空遙感的時(shí)空分辨率限制,彌補(bǔ)了地面遙感采集效率較低,空間覆蓋能力有限的不足,與地面調(diào)查數(shù)據(jù)更好地銜接,在地面觀測、衛(wèi)星和航空遙感之間搭建起信息推繹的橋梁,顯著提高了生態(tài)系統(tǒng)動(dòng)態(tài)監(jiān)測的時(shí)間頻度和空間尺度[6-7]。

    目前,大量研究通過光譜植被指數(shù)構(gòu)建多元線性回歸(Multiplelinearregression,MLR)模型對(duì)植被的葉面積指數(shù)、葉綠素和地上生物量等屬性進(jìn)行精確研究[8-10],并且利用優(yōu)化植被指數(shù)[11]、結(jié)合光譜變換等形式進(jìn)一步提高模型的估測精度,在預(yù)測森林蓄積量12、農(nóng)作物[13和草地生物量[14等研究中得到應(yīng)用。然而,無人機(jī)在數(shù)據(jù)采集過程中易受到外部環(huán)境干擾和下墊面影響,即使使用最先進(jìn)的無人機(jī)傳感器系統(tǒng),仍然可能遺漏大量的植被信息,從而降低植物群落特征的估計(jì)精度[15],因此了解無人機(jī)獲取草地植物群落地上生物量數(shù)據(jù)信息損失情況亟待進(jìn)一步探討,通過與地面高光譜進(jìn)行對(duì)比能夠分析不同數(shù)據(jù)源的精度差異,判斷無人機(jī)多光譜在植物地上生物量估測中是否存在信息損失。

    本研究以集中分布在我國新疆,具有重要生態(tài)和生產(chǎn)功能且處于退化階段的伊犁絹蒿(Seriphid-iumtransiliense)荒漠草地為研究對(duì)象,基于多光譜無人機(jī)獲取群落冠層影像,以地面高光譜及現(xiàn)場測量數(shù)據(jù)為對(duì)照,通過最佳指數(shù)因子法篩選特征波段,采用不同植被指數(shù)結(jié)合線性回歸法構(gòu)建反演模型,明確無人機(jī)獲取草地群落地上生物量的反演模型及精度,探討其誤差與規(guī)律,通過誤差關(guān)系式建立修正模型,以期提高無人機(jī)對(duì)草地群落的估算精度。

    1 材料與方法

    1. 1 研究區(qū)域概況

    研究區(qū)位于新疆烏魯木齊市米東區(qū)柏楊河鄉(xiāng)典型伊犁絹蒿荒漠草地分布區(qū)(圖1),地處 87° 5259′′E~87°5513\"E , 44°0016′′N~44°0120′′N 之間,年降水量為 250mm 左右,屬于溫帶大陸性半干旱半荒漠氣候。該地區(qū)建群種為伊犁絹蒿,伴生種為角果藜(Ceratocarpusarenarius)和叉毛蓬(Petrosimonia sibirica)等[16]。

    1. 2 試驗(yàn)設(shè)計(jì)

    伊犁絹蒿荒漠草地建群種伊犁絹蒿春季3—4月返青,5一6月營養(yǎng)生長期,7一8月進(jìn)入休眠期,8一9月進(jìn)入花果期,因此參照伊犁絹蒿荒漠草地物候節(jié)律以及研究區(qū)各類植物生長周期[17-18],于2022年4月、6月、9月采集無人機(jī)多光譜、地面高光譜及實(shí)測數(shù)據(jù)。實(shí)驗(yàn)選擇晴朗無云,地面能見度高,風(fēng)力不要超過3級(jí)的天氣,當(dāng)天時(shí)間在12一16點(diǎn)之間采集光譜數(shù)據(jù)。遵循植叢分布均勻的原則,在伊犁絹蒿優(yōu)勢(shì)度的草地內(nèi)設(shè)置3條長 550m 寬 50m ,相隔 200m 的樣帶用于無人機(jī)多光譜采集,每條樣線內(nèi)布設(shè)55個(gè)大小為 0.5m×0.6m 的小樣方用于地面高光譜和實(shí)測數(shù)據(jù)采集,共165個(gè)(圖2)。

    1.3數(shù)據(jù)采集及預(yù)處理

    1.3.1無人機(jī)多光譜使用大疆精靈4多光譜版無人機(jī)對(duì)地采集光譜數(shù)據(jù),設(shè)置飛行高度為 15m 飛行速度為 1.0m?s-2 ,空間分辨率為 0.8cm ,拍照模式為等時(shí)間拍照(間隔為2.0s),相機(jī)沿航線方向拍攝樣線,航向重疊率和旁向重疊率設(shè)置為75%[16] 。通過大疆智圖軟件預(yù)處理操作得到5幅單波段圖像,使用ENVI5.3軟件LayerStacking功能融合為多光譜圖像,RadiomatricCalibration和QuickAtmosphericCorrection(QUAC)功能進(jìn)行輻射定標(biāo)和大氣校正,提取群落感興趣區(qū)(Regionofinterest,ROI保存,BandMath功能計(jì)算紅光波段和近紅外波段的植被指數(shù)光譜影像[19-20],剔除因陣風(fēng)引起變形的圖像異常數(shù)據(jù)。

    1.3.2地面高光譜使用SOC710VP高光譜成像光譜儀采集光譜數(shù)據(jù),將高光譜掃描儀安裝在三腳架上,成像光譜儀鏡頭距離地面 1.5m ,垂直正射于樣方植物冠層自然狀態(tài)的表面,保證測定樣方中無陰影進(jìn)入,前方無物體遮擋,對(duì)樣方內(nèi)的植物群落進(jìn)行光譜成像。每次測量前均使用標(biāo)準(zhǔn)白板校準(zhǔn)反射率。在SRAnal710軟件中提取光譜數(shù)據(jù),將圖片添加RGB假彩色,選定圖片的最大范圍,生成光譜灰度值,再將光譜灰度值(Digitalnumber,DN)值轉(zhuǎn)為光譜反射率,選取所有像素點(diǎn)的光譜反射率平均值作為其的光譜反射值[21]

    1.3.3植物群落實(shí)測地上生物量在采集地面高光譜數(shù)據(jù)后,采用齊地劉割法齊地采集樣方內(nèi)所有植物群落地上生物量,利用電子天平稱量鮮重,將0.5m×0.6m 小樣方內(nèi)群落地上生物量單位換算轉(zhuǎn)化為 g?m-2 。在調(diào)查時(shí)記錄每個(gè)樣地經(jīng)緯度(WGS84坐標(biāo)系)和海拔等信息,同時(shí)用相機(jī)拍攝實(shí)際植被照片,以便后期將無人機(jī)圖片中的地物與

    實(shí)際地物進(jìn)行對(duì)比觀察。

    1.3.4特征波段篩選光譜特征波段篩選的核心為從全波段光譜中確定最具代表性的特征光譜波段,其能夠提高建模效率與模型性能[21]。因此本研究選擇。CHAVEZ等22]提出的針對(duì)整體組合的最佳指數(shù)因子法(Optimumindexfactor,OIF)[16.23]篩選兩個(gè)遙感平臺(tái)的特征波段。使用ENVI軟件ComputeStatistics功能計(jì)算遙感影像各個(gè)波段的標(biāo)準(zhǔn)差,按照光譜區(qū)間范圍將光譜所有的波段分成5個(gè)集合,選擇各集合中標(biāo)準(zhǔn)差最大的波段兩兩組合進(jìn)行相關(guān)性分析,計(jì)算最佳指數(shù)因子。多光譜無人機(jī)只有5個(gè)波段,直接計(jì)算標(biāo)準(zhǔn)差,相關(guān)系數(shù)及OIF值。最佳指數(shù)因子計(jì)算公式如下:

    式中,OIF為最佳指數(shù)因子, Si 為第i個(gè)波段的標(biāo)準(zhǔn)差, Rij 表示ij兩波段的相關(guān)系數(shù)且要取絕對(duì)值,n為在所有波段中選取的波段數(shù)。

    1.4模型建立與精度驗(yàn)證

    在ENVI軟件中,提取ROI范圍內(nèi)無人機(jī)的植被指數(shù)光譜影像,光譜影像上所有像素值的平均值作為群落的植被指數(shù),并對(duì)提取地面高光譜數(shù)據(jù)的反射率進(jìn)行植被指數(shù)計(jì)算,以此獲取不同遙感平臺(tái)的歸一化植被指數(shù)(Normalizeddifferencevegetation,NDVI)、比值植被指數(shù)(Rationvegetationindex,RVI)、差值植被指數(shù)(Differencevegetationindex, DVI)[24] ,采用SPSS軟件對(duì)兩個(gè)遙感平臺(tái)的植被指數(shù)與實(shí)測地上生物量進(jìn)行相關(guān)性分析,并選擇線性逐步回歸法建立光譜數(shù)據(jù)與群落地上生物量反演模型。數(shù)據(jù)集按照7:3的比例隨機(jī)進(jìn)行模型建立與精度驗(yàn)證。選用決定系數(shù)(Determinationcoefficients, R2 )、均方根誤差(Rootmean squareerror,RMSE)、精度(Accuary)3個(gè)指標(biāo)驗(yàn)證模型精度。計(jì)算公式如下:

    式中: bi 和 bi 分別為實(shí)測和估測生物量, 分別為實(shí)測和估測生物量的平均值。

    1. 5 信息損失分析

    使用Origin2021軟件將無人機(jī)多光譜估測值與實(shí)測值通過線性回歸模型擬合,選擇 R2 進(jìn)行精度驗(yàn)證,通過與地面高光譜對(duì)比,判斷不同數(shù)據(jù)源之間的精度差異[25-26]。在Excel2016計(jì)算多光譜無人機(jī)損失率,損失率若為正值表示地面高光譜優(yōu)于無人機(jī)多光譜,若為負(fù)值表示無人機(jī)多光譜優(yōu)于地面高光譜。將實(shí)測值按照從小到大順序依次排列,并對(duì)實(shí)測值、地面高光譜估測值、無人機(jī)多光譜估測值在Origin2021軟件中進(jìn)行行統(tǒng)計(jì),計(jì)算均值和標(biāo)準(zhǔn)差,繪制誤差棒圖,評(píng)估3者的誤差關(guān)系,以分析多光譜無人機(jī)在植物冠層的信息損失。損失率計(jì)算公式如下:

    1.6 模型修正

    在Exce12016中計(jì)算無人機(jī)多光譜估測值與地面高光譜估測值之間的差值,將無人機(jī)植被指數(shù)與差值擬合獲取誤差關(guān)系式,利用無人機(jī)反演模型減去誤差關(guān)系式構(gòu)建修正模型,計(jì)算得出修正后無人機(jī)多光譜估測值并驗(yàn)證精度。

    2 結(jié)果與分析

    2.1草地植物地上生物量反演模型的篩選

    2.1.1特征波段篩選無人機(jī)多光譜數(shù)據(jù)4月、6月和9月特征波段均為Red,RedEdge和NIR,OIF值分別為0.69,1.32和1.22。地面高光譜數(shù)據(jù)選擇的植被指數(shù)需要紅光波段和近紅外波段,因此選擇特征波段4月為448.90,622.60和 761.56nm ,6月為499.69,622.6和 761.56nm ,9月為499.69,622.60和 764.35nm ,OIF值分別為0.71,0.64和0.16。利用以上特征波段計(jì)算植被指數(shù)構(gòu)建地上生物量反演模型(表1)。

    表1不同波段組合的最佳指數(shù)因子

    Table1 OIF fordifferentband combinations

    2.1.2植被指數(shù)與群落地上生物量相關(guān)性分析植被指數(shù)與群落地上生物量關(guān)系如表2所示,相關(guān)系數(shù)均達(dá)到顯著性水平 (Plt;0.01) ,無人機(jī)多光譜數(shù)據(jù)的地上生物量在4月和9月與NDVI的相關(guān)系數(shù)最高,分別為0.500和0.364,6月與RVI的相關(guān)系數(shù)最高,為0.385;地面高光譜數(shù)據(jù)地上生物量3個(gè)月份均與NDVI的相關(guān)系數(shù)最高,分別為0.527,0.505和0.529。

    表2植被指數(shù)與群落地上生物量相關(guān)性分析Table2Correlation analysis between vegetation index and aboveground biomass
    注:*在0.05水平(雙側(cè))上顯著相關(guān),**在0.01水平(雙側(cè))上顯著相關(guān) Note: * Significant correlationatO.O5 level, ** significantcorrelationat O.Ol level

    2.1.3模型建立及精度評(píng)價(jià)群落地上生物量反演模型如表3,無人機(jī)多光譜4月和9月基于NDVI建立的反演模型精度最高,分別為 65.61% 和61.59% ,4月最佳模型為 y=487.7x+5.2916,9 月最佳模型為 y=175.35x+15.777;6 月RVI精度最高為 48.26% ,最佳模型為 y=112.97x-106.47 。地面高光譜3個(gè)月份均與NDVI精度最高,分別為71. 77% , 53.63% 和 67.58% ,4月最佳模型為 y= 434.81x-2.3008,6 月最佳模型為 y=334.19x- 7.2206,9月最佳模型為 $y = 5 7 . 6 3 7 x + 5 . 0 7 6 7 。$ 對(duì)比植被群落地上生物量最佳模型可知,地面高光譜模型反演精度優(yōu)于無人機(jī)多光譜,通過模型的篩選和精度差異對(duì)比能夠分析無人機(jī)獲取地物的信息損失,為模型修正奠定基礎(chǔ)。

    表3基于不同遙感平臺(tái)的群落地上生物量反演模型及精度評(píng)價(jià)Table3Inversionmodelandaccuracyevaluationofcommunityaboveground biomassbasedondiferentremotesensingplatfor

    2.2草地植物地上生物量反演模型的信息損失分析

    2.2.1擬合效果分析將所有樣本的無人機(jī)多光譜估測值與實(shí)測值進(jìn)行擬合,并通過與地面高光譜對(duì)比得到實(shí)測值與估測值的1:1線型圖(圖3),以評(píng)估無人機(jī)多光譜與地面高光譜對(duì)地上生物量估測值的準(zhǔn)確性。整體上地面高光譜反演模型得到的植被地上生物量估測值與實(shí)測值的擬合線更接近1:1線,4月、6月和9月NDVI的 R2 分別為0.5213,0.325和0.4948,RVI的 R2 分別為0.5091,0.323和0.4875;無人機(jī)多光譜4月、6月和9月NDVI的 R2 分別為0.3978,0.174和0.2188,RVI的 R2 分別為0.3977,0.1767和0.2171,對(duì)比不同平臺(tái)可知不同物候期NDVI和RVI的地上生物量模型 R2 均高于無人機(jī)多光譜。說明地面高光譜對(duì)植被覆蓋度的估算能力優(yōu)于無人機(jī)多光譜,進(jìn)一步證明了不同遙感平臺(tái)對(duì)地上生物量估測的精度差異。

    2.2.2損失率分析對(duì)反演模型估算精度分析可知(表4),在4月、6月和9月各植被指數(shù)的損失率均表現(xiàn)為NDVIgt;RVIgt;DVI,說明地上生物量損失率在不同物候期和植被指數(shù)之間具有差異性。其中,對(duì)最佳反演模型的計(jì)算結(jié)果為:4月無人機(jī)多光譜NDVI模型相比地面高光譜估算精度低 8.58% ,9月低 8.87% ,6月無人機(jī)多光譜RVI模型相比地面高光譜NDVI模型估算精度低 9.89% 。說明無人機(jī)多光譜在地上生物量反演過程中存在信息損失,需要根據(jù)損失規(guī)律和差異利用地面高光譜對(duì)其分月份和指標(biāo)修正。

    2.2.3誤差分析不同遙感平臺(tái)估測值與實(shí)測值之間有誤差,且并不均勻,需要進(jìn)行誤差分析。由圖4可知不同物候期和植被指數(shù)均表現(xiàn)為在較低和較高生物量時(shí)誤差較大,在中生物量誤差相對(duì)較小,其中4月、6月和9月生物量范圍分別在 69.13~ 83.87,11.02\~16.77和 19.33~29.67g?m-2 之間誤差較小,呈現(xiàn)出兩頭大中間小的趨勢(shì),因此需要根據(jù)物候期、植被指數(shù)誤差的差異性對(duì)地上生物量進(jìn)行修正。

    表4基于無人機(jī)多光譜反演地上生物量的損失率

    Table4 Loss rate of aboveground biomass inversion based onUAV multispectral

    2.3草地植物地上生物量反演模型的修正

    表5為地面高光譜修正無人機(jī)多光譜的地上生物量誤差關(guān)系式,將無人機(jī)估算模型減去誤差關(guān)系式構(gòu)建修正模型(表6),并利用實(shí)測值、無人機(jī)多光譜修正前后估測值繪制對(duì)比分析圖(圖4)。

    表5無人機(jī)多光譜的地上生物量修正誤差關(guān)系式Table5Aboveground biomass modification error relation of UAV multispectrum

    由表6可知,所有地上生物量修正模型均得到不同程度的提高,提高了 0.02%~3.74%, 4月以RVI的精度提升最多,為 1.04%,6 月和9月DVI提升最多,分別為 2.94% 和 3.74% ,在月份間和植被指數(shù)間無明顯規(guī)律性。修正后模型4月和9月以NDVI的精度最高,分別為 66.36% 和61. 61% ,6月以RVI的精度最高為 50.11% ,與修正前對(duì)比,地上生物量各月反演最佳模型沒有改變,但精度均得到提高。

    圖4地上生物量誤差關(guān)系圖Fig.4Aboveground biomass error diagram
    表6無人機(jī)多光譜的地上生物量修正模型及精度評(píng)價(jià)Table 6Aboveground biomass modification model and accuracy evaluation of UAV multispectrum
    注:提升精度 修正后精度一修正前精度 Note:Improved accuracy accuracyaftermodification-accuracybefore modification

    由圖5所示,各月份無人機(jī)多光譜修正后估測生物量呈現(xiàn)出NDVI,RVI隨實(shí)測生物量增加而更接近實(shí)測值的趨勢(shì),并且4月DVI修正后效果更加明顯,而9月DVI隨地上生物量增加無明顯變化。綜上,無人機(jī)多光譜在4月的修正效果最佳,修正后NDVI精度最高,修正模型為 y=403.431x+17.5936°

    3討論

    3.1反演模型的篩選

    本研究通過最佳指數(shù)因子法篩選無人機(jī)多光譜特征波段為藍(lán)光、紅邊、近紅外波段,地面高光譜4月為藍(lán)光、紅光、近紅外波段,6月和9月為綠光、紅光、近紅外波段。對(duì)植被指數(shù)與群落地上生物量進(jìn)行相關(guān)性分析,3個(gè)月份相關(guān)系數(shù)均達(dá)到了顯著性水平,這與Wang等2研究結(jié)果類似。利用特征波段構(gòu)建反演模型,可以發(fā)現(xiàn)相關(guān)系數(shù)較高的植被指數(shù),構(gòu)建的反演模型精度也較高,這與針對(duì)冬小麥、棉花(Gossypiumhirsutum)等[28-29]農(nóng)作物的生物量估算研究結(jié)果一致,并且通過篩選相關(guān)系數(shù)較大的敏感植被指數(shù)組合作為輸入特征,結(jié)合多種機(jī)器學(xué)習(xí)算法構(gòu)建模型均能獲得較高的估算精度[27.30-31],有效說明了相關(guān)性與模型估算精度之間的密切聯(lián)系。本研究利用3個(gè)植被指數(shù)構(gòu)建反演模型,整體上NDVI和RVI比DVI更適用于地上生物量的估測,與Jin等[32和牛亞曉等33學(xué)者篩選的植被指數(shù)各有不同。綜上,需要針對(duì)遙感平臺(tái)、物候期和地物等差異性需要對(duì)植被指數(shù)進(jìn)行專一性篩選。

    3.2 信息損失分析

    通過擬合效果分析可知地面高光譜在各物候期NDVI和RVI的 R2 均高于無人機(jī)多光譜,進(jìn)一步證明其對(duì)地上生物量的估算能力更優(yōu)。引起無人機(jī)多光譜精度差異的原因是地面高光譜連續(xù)的窄波段能夠捕捉地物更細(xì)微光譜差異[34],有利于實(shí)現(xiàn)地上生物量的高精度反演,而多光譜遙感較少的光譜通道和不連續(xù)的寬波段對(duì)荒漠草地植被參數(shù)提取存在一定困難[35];并且本研究中地面高光譜采集高度為 1.5m ,無人機(jī)飛行高度為 15m ,高度的增加也是引起空間分辨率的差異的原因之一,有研究表明距離目標(biāo)地物更近,受到的環(huán)境干擾相對(duì)更小,反映地物的光譜特征更加真實(shí)[36],說明適宜的空間分辨率有助于提高估測精度[37]。冠層信息損失在不同物候期也具有差異性,春季4月伊犁絹蒿處于返青期,角果藜處于苗期;夏季6月群落植物處于生長旺盛期,但由于干旱葉片水分含量降低,逐漸枯黃,7一8月伊犁絹蒿有休眠現(xiàn)象,秋季9月伊犁絹蒿進(jìn)入花期,角果藜枯黃[16.38],為此,在不同物候期的地物監(jiān)測中,需重視因植物生長形態(tài)變化而造成的反演精度信息損失的差異。無人機(jī)多光譜最佳模型的損失率分別為 8.58%,9.89% 和 8.87% ,誤差分析結(jié)果表明隨地上生物量增加呈現(xiàn)出兩頭大中間小的趨勢(shì),證明了無人機(jī)多光譜在植被地上生物量反演過程中存在信息損失。趙菡等[39]對(duì)研究樣地進(jìn)行立地質(zhì)量分級(jí),認(rèn)為中間立地等級(jí)的誤差最小,立地質(zhì)量越接近平均水平,生物量估計(jì)的相對(duì)誤差越小,與本文研究結(jié)果一致。魏高磊等40認(rèn)為中高植被覆蓋度時(shí)空間尺度效應(yīng)對(duì)植被指數(shù)影響顯著,這歸因于植被內(nèi)部光譜差異和冠層陰影的影響更為突出,而土壤光譜差異較小。因此需要考慮遙感平臺(tái)、物候期和植被指數(shù)間誤差的差異性對(duì)地上生物量進(jìn)行修正。

    3.3 反演模型的修正

    一般而言,僅采用單一的低、中等分辨率的遙感數(shù)據(jù)在大尺度范圍上對(duì)植被地上生物量的反演能力較為有限,例如部分研究者利用GF-1衛(wèi)星遙感構(gòu)建的反演模型的 R2 在 0.14~0.40[41-42] ,而馬儀等[43利用無人機(jī)構(gòu)建升尺度修正后GF-1衛(wèi)星遙感數(shù)據(jù)土壤含水率反演模型,葉靜蕓等44同樣實(shí)現(xiàn)了高分辨率WorldView-3數(shù)據(jù)的估算結(jié)果對(duì)低分辨率Landsat-OLI數(shù)據(jù)的修正,尤其是部分研究者融合了地面與無人機(jī)遙感信息[45],或構(gòu)建“空-地\"光譜融合構(gòu)建植被監(jiān)測模型并驗(yàn)證,發(fā)現(xiàn)兩種尺度融合模型精度優(yōu)于單一尺度模型[46-47],均證明了利用高分辨率遙感數(shù)據(jù)修正較低空間分辨率遙感數(shù)據(jù)的必要性。本研究通過地面高光譜修正無人機(jī)多光譜后所有模型估測精度均提高了 0.02%~ 3.74% ,但從修正前后對(duì)比分析圖中可知,地上生物量較低時(shí)修正后估測值高于修正前,較高時(shí)呈現(xiàn)相反的結(jié)果,在中間時(shí)更接近實(shí)測值,并且對(duì)DVI的修正效果不佳。有研究表明在低、高生物量等級(jí)的估測值與實(shí)測值之間誤差較大,并且不同等級(jí)對(duì)植被指數(shù)的選擇不同[48];俞靜等49]篩選出最優(yōu)反演模型對(duì)研究區(qū)域內(nèi)植被生物量進(jìn)行分級(jí)評(píng)估,發(fā)現(xiàn)大部分區(qū)域的生物量處于較低等級(jí),說明生物量的空間分布并不均勻;對(duì)地上生物量進(jìn)行分級(jí)反演,特別是針對(duì)中、高等級(jí)的精度有顯著提高[50],并且隨等級(jí)的增加,整體估測精度會(huì)逐步提高[51];考慮到草本覆蓋度背景的異質(zhì)性,采用分層方案能進(jìn)一步提高光譜指數(shù)與喬灌木地上部分生物量模型的反演精度[52]。因此在今后研究中需要根據(jù)物候期與植被指數(shù)的差異性對(duì)地上生物量分級(jí)建模與修正,有利于進(jìn)一步提高無人機(jī)的反演精度。

    本研究通過高精度的地面高光譜對(duì)無人機(jī)進(jìn)行修正,在一定程度上提高了其反演精度,使得地面高光譜的作用得以發(fā)揮。從結(jié)果來看仍存在修正后模型精度提升不顯著的問題,研究需進(jìn)一步結(jié)合對(duì)數(shù)、二次、冪、指數(shù)函數(shù),尤其是多種機(jī)器學(xué)習(xí)法探索多個(gè)組分之間的非線性關(guān)系;在建模時(shí)加大樣本數(shù)量,分級(jí)建立反演模型以降低樣本異質(zhì)性,提高模型精度和穩(wěn)定性。今后需優(yōu)化實(shí)驗(yàn)流程,加強(qiáng)地面光譜數(shù)據(jù)修正無人機(jī)的相關(guān)研究,以進(jìn)一步應(yīng)用無人機(jī)修正航天遙感,提升草地大面積遙感動(dòng)態(tài)監(jiān)測的精度。

    4結(jié)論

    (1)通過最佳指數(shù)因子法篩選無人機(jī)多光譜特征波段3個(gè)月份均為Red,RedEdge和NIR,而地面高光譜特征波段在不同物候期表現(xiàn)不一致。各物候期植被指數(shù)與群落地上生物量相關(guān)性達(dá)到顯著性水平( (Plt;0.01) ,無人機(jī)多光譜4月和9月基于NDVI建立的模型精度最高,分別為 65.61% 和61. 59% ,6月RVI精度最高為 48.26% ;地面高光譜各月均與NDVI精度最高,分別為 71.77% 53.63% 和 67.58% 。針對(duì)性地篩選模型能夠減少無人機(jī)獲取地物的信息損失,對(duì)比最佳模型地面高光譜反演精度均優(yōu)于無人機(jī)多光譜,為模型修正奠定基礎(chǔ)。

    (2)通過無人機(jī)與地面高光譜反演效果對(duì)比分析,整體上地面高光譜估測值與實(shí)測值的擬合線更接近1:1線,優(yōu)于無人機(jī)多光譜反演能力。無人機(jī)多光譜生物量損失率在各物候期均為NDVIgt;RVIgt; DVI,最佳模型損失率分別為 8.58% , 9.89% 和8.87% 。可見,相比地面高光譜平臺(tái),無人機(jī)多光譜在植物冠層反演時(shí)在不同月份和模型間存在不同程度的信息損失,應(yīng)分物候期和植被指數(shù)對(duì)其修正。

    (3)通過建立誤差關(guān)系式修正無人機(jī)多光譜反演模型,所有模型估測精度均得到不同程度的提高,提高了 0.02%~3.74% 。與修正前對(duì)比,修正后的無人機(jī)多光譜地上生物量各月反演最佳模型沒有改變,但精度均得到提高;以4月NDVI模型最高,精度為 66.36% ,修正模型為 y=403.431x+ 17.5936??梢?,利用地面高光譜平臺(tái)修正,提高了無人機(jī)對(duì)地監(jiān)測精度。

    參考文獻(xiàn)

    [1]LIU Y,F(xiàn)ENG HK,YUEJB,et al.Estimation of aboveg-round biomass of potatoes based on characteristic variablesextracted fromUAVhyperspectral imagery[J].RemoteSens-ing,2022,14(20):5121

    [2]LIW,DOUZG,WANGY,etal.Estimationof above-groundbiomass of reed(Phragmitescommunis)based onin situ hyperspectral data inBeijingHanshiqiaoWetland,China[J].Wet-landsEcologyand Management,2019,27(1):87-102

    [3] WANG D L,XIN XP,SHAO QQ,et al. Modeling aboveg-round biomass in Hulunber grassland ecosystem by usingunmanned aerialvehicle discrete lidar[J].Sensors,2O17,17(1):180

    [4] VILLOSLADAPECINAM,BERGAMOTF,WARDRD,etal.Anovel UAV-basedapproach forbiomasspredictionandgrassland structure assessment in coastal meadows[J].Ecologi-calIndicators,2021(122):107227

    [5] 楊棟誤,李亞強(qiáng),刀劍,等,基于無人機(jī)多光譜與地面高光譜遙感的土壤主要養(yǎng)分含量估測[J].江蘇農(nóng)業(yè)科學(xué),2022,50(2):178-186

    [6] 胡健波,張?。疅o人機(jī)遙感在生態(tài)學(xué)中的應(yīng)用進(jìn)展[J].生態(tài)學(xué)報(bào),2018,38(1):20-30

    [7]WUSB,WANGJ,YANZB,et al.Monitoringtree-crownscaleautumn leaf phenologyina temperate forest with aninte-gration of PlanetScope and drone remote sensing observations[J].ISPRS Journal of Photogrammetry and Remote Sensing,2021(171) :36-48

    [8]POTGIETER A B,GEORGE-JAEGGLI B,CHAPMAN SC,et al. Multi-spectral imaging from an unmanned aerialvehicle enables the assessment of seasonal leaf area dynamics ofSorghum breeding lines[J]. Frontiers in Plant Science,2017(8):1532

    [9]王建華.基于無人機(jī)多光譜遙感的玉米葉綠素含量監(jiān)測研究[D].太谷:山西農(nóng)業(yè)大學(xué),2023:29-35

    [10]李天馳,馮海寬,朱貝貝,等.基于無人機(jī)高光譜和數(shù)碼影像數(shù)據(jù)的冬小麥生物量反演[J].現(xiàn)代農(nóng)業(yè)科技,2020(20):1-5

    [11]安海波,李斐,趙萌莉,等.基于優(yōu)化光譜指數(shù)的牧草生物量估算[J].光譜學(xué)與光譜分析,2015,35(11):3155

    [12]霍寶民,職慶利,俞秀玲,等.基于無人機(jī)機(jī)載激光雷達(dá)構(gòu)建森林蓄積量模型反演及精度分析[J].河南林業(yè)科技,2023,43(3):13-16

    [13]YUE JB,YANG G J,LI C C,et al. Estimation of winterwheat above-ground biomass using unmanned aerial vehicle-based snapshot hyperspectral sensor and crop height improvedmodels[J].Remote Sensing,2017,9(7):708

    [14]王秀梅.內(nèi)蒙古典型草原植被地上生物量遙感反演——基于多源數(shù)據(jù)的整合[D].呼和浩特:內(nèi)蒙古大學(xué),2022:70-77

    [15]HU TY,SUN XL,SUYJ,et al.Development and perfor-mance evaluation of a very low-cost UAV-lidar system for for-estry applications[J]. Remote Sensing,2021,13(1) :77

    [16]馬建.基于多源遙感數(shù)據(jù)的伊犁絹蒿荒漠草地植物的識(shí)別與反演[D].烏魯木齊:新疆農(nóng)業(yè)大學(xué),2022:11-12

    [17]鄭偉.伊犁絹蒿的生態(tài)學(xué)特性及利用和保護(hù)研究進(jìn)展[J].草原與草坪,2013,33(1):68-75,83

    [18]武紅旗,范燕敏,靳瑰麗,等.伊犁絹蒿荒漠草地植物光譜特征[J].草業(yè)科學(xué),2019,36(7):1765-1773

    [19]張恒瑞,段喜明,魏征,等.基于無人機(jī)多光譜遙感的華北地區(qū)夏玉米LAI監(jiān)測[J].山西農(nóng)業(yè)科學(xué),2021,49(5):608-614

    [20]曹中盛,李艷大,潘玉霞,等.基于無人機(jī)RGB圖像的金沙柚葉片葉綠素含量監(jiān)測研究[J].江西農(nóng)業(yè)大學(xué)學(xué)報(bào),2022,44(6):1407-1418

    [21]鮑浩,張艷.基于改進(jìn)哈里斯鷹優(yōu)化算法的光譜特征波段選擇模型研究[J].光譜學(xué)與光譜分析,2024,44(1):148-157

    [22] CHAVEZ P,BERLIN G,SOWERS L. Statistical method forselecting landsat mss ratios[J].Journal of AppliedPhoto-graphic Engineering,1982,8(1):23-30

    [23]魯向暉,龔榮新,張海娜,等.基于無人機(jī)多光譜遙感的矮林芳樟光合參數(shù)估測[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2023,54(10):179-187

    [24]BANNARI A,MORIN D,BONNF,et al. A review of vegeta-tion indices[J].Remote Sensing Reviews,1995,13(1/2):95-120

    [25]陶惠林,徐良驥,馮海寬,等.基于無人機(jī)數(shù)碼影像的冬小麥株高和生物量估算[J].農(nóng)業(yè)工程學(xué)報(bào),2019,35(19):107-116

    [26] ZHAO X X,SU Y J,HU T Y,et al. Analysis of UAV lidarinformation loss and its influence on the estimation accuracy ofstructural and functional traits in a meadow steppe[J].Ecologi-cal Indicators,2022(135):108515

    [27]WANG Y W,XUE W X. Sustainable development early warn-ingand financing risk management of resource-based industrialclusters using optimization algorithms[J].Journal of EnterpriseInformation Management,2022,35(4/5):1374-1391

    [28]陶惠林,馮海寬,徐良驥,等.基于無人機(jī)高光譜遙感數(shù)據(jù)的冬小麥生物量估算[J].江蘇農(nóng)業(yè)學(xué)報(bào),2020,36(5):1154-1162

    [29]周萌,韓曉旭,鄭恒彪,等.基于參數(shù)化和非參數(shù)化法的棉花生物量高光譜遙感估算[J].中國農(nóng)業(yè)科學(xué),2021,54(20):4299-4311

    [30]馬彥鵬,邊明博,樊意廣,等.基于無人機(jī)RGB影像的馬鈴薯植株鉀含量估算[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2023,54(7):196-203,233

    [31]DONGJW,XIAO XM,WAGLEP,et al.Comparison offour EVI-based models for estimating grossprimary productionof maize and soybean croplands and tallgrass prairie undersevere drought[J]. Remote Sensing of Environment,2015(162):154-168

    [32]JINXL,YANGGJ,XUXG,etal.Combined multi-temporaloptical and radar parameters for estimating LAI and biomass inwinterwheat usingHJ and RADARSAR-2 data[J].RemoteSensing,2015,7(10):13251-13272

    [33]牛亞曉,張立元,韓文霆,等,基于無人機(jī)遙感與植被指數(shù)的冬小麥覆蓋度提取方法[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2018,49(4):212-221

    [34]孫偉偉,劉圍圍,王煜淼,等.2010—2022年全球濕地高光譜遙感研究進(jìn)展與展望[J].遙感學(xué)報(bào),2023,27(6):1281-1299

    [35]苗春麗,伏帥,劉潔,等.基于UAV成像高光譜圖像的高寒草甸地上生物量一以海北試驗(yàn)區(qū)為例[J].草業(yè)科學(xué),2022,39(10):1992-2004

    [36]王圓,基于無人機(jī)高光譜遙感及人工智能的荒漠化草原地物分類研究[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué),2023:2-6

    [37]張婧,劉詠梅,徐健,等,影像分辨率對(duì)植被覆蓋度提取的影響[J].水土保持研究,2014,21(3):120-124

    [38]孫宗玖,安沙舟,許鵬.伊犁絹蒿構(gòu)件動(dòng)態(tài)變化研究[J].草地學(xué)報(bào),2007,15(5):454-459

    [39]趙菡,雷淵才,符利勇.江西省不同立地等級(jí)的馬尾松林生物量估計(jì)和不確定性度量[J」.林業(yè)科學(xué),2017,53(8):81-93

    [40]魏高磊,吳太夏,王樹東,等.基于無人機(jī)數(shù)據(jù)的植被指數(shù)空間尺度效應(yīng)研究[J].地理空間信息,2021,19(4):4-9,40

    [41]張智韜,陳欽達(dá),黃小魚,等,基于加權(quán)算法的空-天遙感升尺度土壤含鹽量監(jiān)測模型[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2022,53(9):226-238,251

    [42]姚一飛,王爽,張珺銳,等.基于GF-1衛(wèi)星遙感的河套灌區(qū)土壤含水率反演模型研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2022,53(9):239-251

    [43]馬儀,黃組桂,賈江棟,等.基于無人機(jī)-衛(wèi)星遙感升尺度的土壤水分監(jiān)測模型研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2023,54(6):307-318

    [44]葉靜蕓,吳波,賈曉紅,等.極干旱區(qū)稀疏荒漠植被地上生物量遙感估算[J].干旱區(qū)地理,2022,45(2):478-487

    [45]朱航.基于多源遙感信息融合的作物營養(yǎng)狀況監(jiān)測與噴灑控制系統(tǒng)的研究[D].長春:吉林大學(xué),2011:104-109

    [46]印彩霞,基于“空-地\"光譜融合的棉花氮營養(yǎng)監(jiān)測研究[D].石河子:石河子大學(xué),2023:89-95

    [47]KAUFMAN JR,EISMANN M T,CELENK M. Assessmentof spatial-spectral feature-level fusion for hyperspectral targetdetection[J].IEEE Journal of Selected Topics in AppliedEarth Observations and Remote Sensing,20l5,8(6):2534-2544

    [48]王惠寧,靳瑰麗,范燕敏,等.不同蓋度下伊犁絹蒿荒漠草地光譜特征及蓋度反演精度研究[J].中國草地學(xué)報(bào),2019,41(2):51-57

    [49]俞靜,張世文,芮婷婷,等.基于無人機(jī)遙感的多特征組礦區(qū)草本植物地上生物量反演[J].草業(yè)科學(xué),2024,41(1):35-48

    [50]李嘉欣,王惠寧,劉文昊,等.伊犁絹蒿光譜反射率及反演精度對(duì)覆蓋度變化的響應(yīng)[J].草業(yè)科學(xué),2023,40(11):2775-2786

    [51]王艷芳,譚露,郭紅麗,等.基于多源遙感數(shù)據(jù)的溧陽市林地植被覆蓋度時(shí)空差異研究[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2023,47(6):183-191

    [52]石永磊,王志慧,李世明,等.基于光學(xué)遙感的稀疏喬灌木地上部分生物量反演方法[J].林業(yè)科學(xué),2022,58(2):13-22

    (責(zé)任編輯劉婷婷)

    猜你喜歡
    植被指數(shù)波段修正
    不同光譜變換形式對(duì)土壤有機(jī)質(zhì)偏最小二乘估算模型精度的影響
    涂改帶有話說
    基于無人機(jī)多光譜遙感和機(jī)器學(xué)習(xí)的棉花SPAD值預(yù)測
    運(yùn)用機(jī)會(huì)成本法進(jìn)行土地評(píng)估的理論探討
    輻射制冷涂料研究進(jìn)展:提升太陽光反射率的機(jī)理與實(shí)現(xiàn)途徑
    陽泉市植被覆蓋遙感監(jiān)測與時(shí)序分析
    基于Sentinel-2多光譜數(shù)據(jù)的紅棗葉面積指數(shù)反演
    基于單層二硫化鉬和金納米柱的可見光波段寬帶吸收器
    基于無人機(jī)多光譜數(shù)據(jù)的葉綠素a濃度反演模型研究
    基于無人機(jī)多光譜和改進(jìn)BPNN的煙草病毒病檢測
    舔av片在线| 男人舔女人下体高潮全视频| 国产极品精品免费视频能看的| av在线老鸭窝| 成人二区视频| 嫩草影院新地址| 日本成人三级电影网站| 色综合站精品国产| 搞女人的毛片| 不卡视频在线观看欧美| 麻豆国产97在线/欧美| 国内精品一区二区在线观看| 欧美zozozo另类| 亚洲中文字幕日韩| 成人性生交大片免费视频hd| 日日干狠狠操夜夜爽| 精品乱码久久久久久99久播| 色在线成人网| 久久99热这里只有精品18| 看片在线看免费视频| 精品久久久久久久人妻蜜臀av| 免费在线观看日本一区| 日本三级黄在线观看| 1024手机看黄色片| 制服丝袜大香蕉在线| 尾随美女入室| 欧美黑人巨大hd| 亚洲精品日韩av片在线观看| 在线观看舔阴道视频| 免费无遮挡裸体视频| 日韩欧美 国产精品| 深夜a级毛片| 88av欧美| 天美传媒精品一区二区| 亚洲欧美清纯卡通| 少妇人妻精品综合一区二区 | 国产精品,欧美在线| 在线观看免费视频日本深夜| 12—13女人毛片做爰片一| 精华霜和精华液先用哪个| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产亚洲av天美| 亚洲在线自拍视频| 小蜜桃在线观看免费完整版高清| 久久久精品大字幕| 国产成人影院久久av| 嫩草影视91久久| 欧美潮喷喷水| 热99re8久久精品国产| 最后的刺客免费高清国语| 日韩 亚洲 欧美在线| 日韩欧美在线二视频| 午夜日韩欧美国产| 欧美日韩综合久久久久久 | 亚洲专区中文字幕在线| 男人舔奶头视频| 午夜福利高清视频| 尤物成人国产欧美一区二区三区| 国产精品日韩av在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 一区二区三区四区激情视频 | 精品久久久久久久久亚洲 | 国产伦精品一区二区三区四那| 国产乱人伦免费视频| 在线观看av片永久免费下载| 久久精品久久久久久噜噜老黄 | 成人欧美大片| 国产黄色小视频在线观看| 尾随美女入室| 一区二区三区免费毛片| www日本黄色视频网| 国产成人一区二区在线| 在线观看舔阴道视频| 欧美精品啪啪一区二区三区| 日本-黄色视频高清免费观看| av国产免费在线观看| a级毛片a级免费在线| 一区二区三区四区激情视频 | 亚洲欧美日韩高清在线视频| 在线观看美女被高潮喷水网站| 在线天堂最新版资源| 人妻久久中文字幕网| a级毛片免费高清观看在线播放| 人妻夜夜爽99麻豆av| 国产久久久一区二区三区| 99国产精品一区二区蜜桃av| 国产精品久久久久久久电影| 国产单亲对白刺激| 黄色视频,在线免费观看| 亚洲真实伦在线观看| 久久精品人妻少妇| 一区福利在线观看| 丰满乱子伦码专区| 亚洲最大成人手机在线| 校园春色视频在线观看| 精品一区二区三区av网在线观看| 午夜免费男女啪啪视频观看 | 如何舔出高潮| 久久久久免费精品人妻一区二区| 亚洲人与动物交配视频| 日韩欧美免费精品| 日本与韩国留学比较| 亚洲av免费在线观看| 欧美激情在线99| 国内精品久久久久久久电影| 十八禁国产超污无遮挡网站| 免费黄网站久久成人精品| 久99久视频精品免费| 男女做爰动态图高潮gif福利片| 99视频精品全部免费 在线| 午夜久久久久精精品| 亚洲精华国产精华精| 女生性感内裤真人,穿戴方法视频| 两个人的视频大全免费| 色综合婷婷激情| 午夜免费男女啪啪视频观看 | 麻豆成人av在线观看| 国产精品福利在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 99精品久久久久人妻精品| 麻豆国产av国片精品| 桃色一区二区三区在线观看| 最近最新免费中文字幕在线| 干丝袜人妻中文字幕| 午夜视频国产福利| 一区二区三区激情视频| 国产高清激情床上av| 观看美女的网站| 久久精品国产亚洲网站| 69人妻影院| 国产黄a三级三级三级人| 日韩欧美国产一区二区入口| 免费看光身美女| 九色国产91popny在线| 日韩 亚洲 欧美在线| 午夜福利成人在线免费观看| 一进一出抽搐动态| 亚洲欧美激情综合另类| 在线播放国产精品三级| 国产在视频线在精品| 午夜爱爱视频在线播放| 久久精品国产亚洲网站| 性欧美人与动物交配| 亚洲成a人片在线一区二区| 亚洲在线自拍视频| av女优亚洲男人天堂| 日韩精品青青久久久久久| 97超视频在线观看视频| 少妇猛男粗大的猛烈进出视频 | 日本三级黄在线观看| 免费看光身美女| 最好的美女福利视频网| 超碰av人人做人人爽久久| 久久精品国产99精品国产亚洲性色| a级一级毛片免费在线观看| 国产熟女欧美一区二区| av国产免费在线观看| 最近中文字幕高清免费大全6 | 黄色一级大片看看| 国内揄拍国产精品人妻在线| 久久久久久久精品吃奶| 国产精品久久久久久久久免| 人妻丰满熟妇av一区二区三区| 亚洲专区国产一区二区| www日本黄色视频网| 我的女老师完整版在线观看| 亚洲va在线va天堂va国产| 久久人妻av系列| 免费av不卡在线播放| 97人妻精品一区二区三区麻豆| 免费电影在线观看免费观看| 少妇裸体淫交视频免费看高清| xxxwww97欧美| 91在线观看av| 99国产极品粉嫩在线观看| 成年女人永久免费观看视频| 简卡轻食公司| 老司机午夜福利在线观看视频| 精品久久久久久久人妻蜜臀av| 国产熟女欧美一区二区| 亚洲av中文av极速乱 | 51国产日韩欧美| 麻豆av噜噜一区二区三区| 国产高清三级在线| 国产黄片美女视频| 中文字幕av成人在线电影| 香蕉av资源在线| 久久久精品欧美日韩精品| 狠狠狠狠99中文字幕| 国产人妻一区二区三区在| 国内久久婷婷六月综合欲色啪| 尾随美女入室| 我的老师免费观看完整版| 亚洲av美国av| 18禁裸乳无遮挡免费网站照片| 国产在线男女| 久久久精品欧美日韩精品| 成年免费大片在线观看| 成人特级黄色片久久久久久久| 亚洲无线在线观看| 九色国产91popny在线| av在线观看视频网站免费| 免费大片18禁| 日韩大尺度精品在线看网址| 亚洲美女视频黄频| 日本黄色片子视频| 国产精品一区二区性色av| 制服丝袜大香蕉在线| 国产色爽女视频免费观看| 欧美又色又爽又黄视频| 桃红色精品国产亚洲av| 欧美不卡视频在线免费观看| 中文字幕熟女人妻在线| 精品一区二区三区av网在线观看| 草草在线视频免费看| 午夜精品在线福利| 免费无遮挡裸体视频| 久久久久久国产a免费观看| 免费黄网站久久成人精品| 国产乱人视频| 久久午夜亚洲精品久久| 亚洲av五月六月丁香网| 听说在线观看完整版免费高清| 一本精品99久久精品77| 黄片wwwwww| 熟女人妻精品中文字幕| 一级黄色大片毛片| 观看美女的网站| 精品久久久噜噜| 精品久久久久久久久亚洲 | 村上凉子中文字幕在线| 日韩 亚洲 欧美在线| 国语自产精品视频在线第100页| 欧美性猛交黑人性爽| 亚洲精品456在线播放app | 一级av片app| 欧美日韩黄片免| 岛国在线免费视频观看| 搡老熟女国产l中国老女人| 免费av毛片视频| x7x7x7水蜜桃| 亚洲中文字幕一区二区三区有码在线看| 精品免费久久久久久久清纯| 国产亚洲精品综合一区在线观看| 国产精品免费一区二区三区在线| 蜜桃亚洲精品一区二区三区| 久久欧美精品欧美久久欧美| 最近在线观看免费完整版| 99久久九九国产精品国产免费| 床上黄色一级片| 亚洲男人的天堂狠狠| 久久精品人妻少妇| 久久久色成人| 又紧又爽又黄一区二区| 日本爱情动作片www.在线观看 | 在线观看av片永久免费下载| 亚洲av免费高清在线观看| 国产免费男女视频| 国产精品爽爽va在线观看网站| 国产 一区 欧美 日韩| 少妇人妻一区二区三区视频| 简卡轻食公司| 日本黄大片高清| 亚洲av日韩精品久久久久久密| 午夜激情福利司机影院| 欧美日韩黄片免| 男女边吃奶边做爰视频| 国产单亲对白刺激| 人人妻人人澡欧美一区二区| 国产麻豆成人av免费视频| av福利片在线观看| 亚洲黑人精品在线| 国产精品人妻久久久久久| 亚洲av美国av| 特大巨黑吊av在线直播| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩高清专用| 亚洲中文字幕日韩| 亚洲精品久久国产高清桃花| 伦理电影大哥的女人| 国产一区二区三区在线臀色熟女| 男女那种视频在线观看| 久久99热这里只有精品18| 国产单亲对白刺激| 欧美又色又爽又黄视频| 精品欧美国产一区二区三| 亚洲精华国产精华液的使用体验 | 色噜噜av男人的天堂激情| 麻豆精品久久久久久蜜桃| 亚洲成人久久性| 最近最新免费中文字幕在线| 亚洲色图av天堂| 偷拍熟女少妇极品色| 啦啦啦观看免费观看视频高清| 国产免费av片在线观看野外av| 美女高潮的动态| 夜夜看夜夜爽夜夜摸| 免费搜索国产男女视频| 亚洲,欧美,日韩| 色哟哟·www| 最近最新中文字幕大全电影3| 成人三级黄色视频| 日本黄色片子视频| 最好的美女福利视频网| 亚洲成a人片在线一区二区| 久久久久久九九精品二区国产| videossex国产| 欧美激情国产日韩精品一区| 亚洲狠狠婷婷综合久久图片| 日日摸夜夜添夜夜添小说| 丰满人妻一区二区三区视频av| 亚洲av电影不卡..在线观看| 成人毛片a级毛片在线播放| 噜噜噜噜噜久久久久久91| 日本在线视频免费播放| 在线a可以看的网站| 久久久色成人| 亚洲国产精品久久男人天堂| 一进一出抽搐gif免费好疼| 性插视频无遮挡在线免费观看| 少妇的逼水好多| 亚洲av.av天堂| 非洲黑人性xxxx精品又粗又长| 伊人久久精品亚洲午夜| 搡女人真爽免费视频火全软件 | 欧美bdsm另类| 欧美成人性av电影在线观看| 亚洲人成网站在线播| 成年女人看的毛片在线观看| 美女xxoo啪啪120秒动态图| 男人狂女人下面高潮的视频| 亚洲成人免费电影在线观看| 男人舔女人下体高潮全视频| 久久久久久大精品| 极品教师在线视频| 欧美最新免费一区二区三区| 日日撸夜夜添| 欧美色视频一区免费| 国产亚洲精品综合一区在线观看| 亚洲精品成人久久久久久| 国产精品女同一区二区软件 | 成人三级黄色视频| 日本爱情动作片www.在线观看 | 国产探花在线观看一区二区| 中文字幕av成人在线电影| 欧美成人性av电影在线观看| 99热这里只有是精品50| 人人妻,人人澡人人爽秒播| 日本 欧美在线| 看免费成人av毛片| 亚洲国产精品成人综合色| 成人鲁丝片一二三区免费| 两个人的视频大全免费| 欧美xxxx黑人xx丫x性爽| 白带黄色成豆腐渣| 18禁在线播放成人免费| 老司机深夜福利视频在线观看| 中文字幕熟女人妻在线| 亚洲,欧美,日韩| 免费看美女性在线毛片视频| 国产精品99久久久久久久久| 亚洲av熟女| 人妻制服诱惑在线中文字幕| 亚洲欧美日韩卡通动漫| 亚洲经典国产精华液单| 啪啪无遮挡十八禁网站| 国产精品98久久久久久宅男小说| 中国美女看黄片| 国产高清不卡午夜福利| 毛片女人毛片| av在线蜜桃| 真人一进一出gif抽搐免费| 三级国产精品欧美在线观看| 一进一出抽搐动态| 全区人妻精品视频| 欧美成人性av电影在线观看| 国产av麻豆久久久久久久| 日韩欧美精品v在线| 午夜福利视频1000在线观看| 亚洲中文日韩欧美视频| 人人妻人人澡欧美一区二区| 两个人视频免费观看高清| 国产大屁股一区二区在线视频| 大型黄色视频在线免费观看| 久久6这里有精品| 淫妇啪啪啪对白视频| 午夜影院日韩av| 中国美女看黄片| 亚洲中文日韩欧美视频| 我的女老师完整版在线观看| 国产人妻一区二区三区在| h日本视频在线播放| 我要看日韩黄色一级片| 级片在线观看| 亚洲成人精品中文字幕电影| bbb黄色大片| 中文字幕av成人在线电影| 中文字幕精品亚洲无线码一区| 中文字幕av成人在线电影| 一级黄片播放器| 亚洲乱码一区二区免费版| 日本黄色片子视频| 男女下面进入的视频免费午夜| 亚洲人成网站在线播| 久久99热这里只有精品18| 亚洲人成网站在线播| 无遮挡黄片免费观看| 午夜福利在线观看吧| 国产精品电影一区二区三区| 国产免费av片在线观看野外av| 性欧美人与动物交配| 欧美+亚洲+日韩+国产| 麻豆国产av国片精品| 欧美区成人在线视频| 国产精品爽爽va在线观看网站| 十八禁国产超污无遮挡网站| 亚洲精品成人久久久久久| 在线免费十八禁| 国产在线男女| 亚洲av美国av| 日本撒尿小便嘘嘘汇集6| 国产午夜福利久久久久久| 噜噜噜噜噜久久久久久91| 黄色丝袜av网址大全| 久久欧美精品欧美久久欧美| 国产精华一区二区三区| 欧美激情国产日韩精品一区| 成年女人永久免费观看视频| 欧美最黄视频在线播放免费| 免费不卡的大黄色大毛片视频在线观看 | 天堂影院成人在线观看| 黄色欧美视频在线观看| 国产精品久久久久久久电影| 国产精品女同一区二区软件 | 精品国内亚洲2022精品成人| 亚洲精品一区av在线观看| 麻豆国产97在线/欧美| 国产高清三级在线| 一a级毛片在线观看| 麻豆av噜噜一区二区三区| 嫩草影视91久久| 最近在线观看免费完整版| 亚洲专区国产一区二区| 我要搜黄色片| 99久久精品一区二区三区| 日日撸夜夜添| 国产精品伦人一区二区| 久久精品人妻少妇| 国产亚洲精品久久久com| 小说图片视频综合网站| 色综合站精品国产| 联通29元200g的流量卡| 亚洲自拍偷在线| 在线免费观看不下载黄p国产 | 久久久色成人| 欧美区成人在线视频| 亚洲欧美日韩卡通动漫| 日韩国内少妇激情av| 中文字幕熟女人妻在线| 精品国内亚洲2022精品成人| 久久中文看片网| 精品日产1卡2卡| 日韩一本色道免费dvd| 性插视频无遮挡在线免费观看| 成人永久免费在线观看视频| 国国产精品蜜臀av免费| 久久人妻av系列| 久久久成人免费电影| 99久国产av精品| 精品久久久久久,| 久久精品国产亚洲网站| 嫩草影视91久久| 毛片女人毛片| 亚洲中文字幕日韩| 国产精品久久视频播放| 国产黄片美女视频| 99国产精品一区二区蜜桃av| 欧美人与善性xxx| 日韩强制内射视频| 小蜜桃在线观看免费完整版高清| 成年女人永久免费观看视频| 啦啦啦观看免费观看视频高清| 久久久久久久久久黄片| 亚洲不卡免费看| 在线观看一区二区三区| 中出人妻视频一区二区| 精品人妻视频免费看| 最新在线观看一区二区三区| 国内揄拍国产精品人妻在线| 啦啦啦韩国在线观看视频| 日韩欧美免费精品| 在线天堂最新版资源| 三级毛片av免费| 五月伊人婷婷丁香| 精品一区二区三区av网在线观看| 免费大片18禁| 全区人妻精品视频| xxxwww97欧美| 高清日韩中文字幕在线| 中文字幕高清在线视频| 日本爱情动作片www.在线观看 | 成人特级av手机在线观看| 亚洲无线观看免费| 中亚洲国语对白在线视频| 国产免费一级a男人的天堂| 国内精品久久久久久久电影| 最新中文字幕久久久久| 又黄又爽又刺激的免费视频.| 嫩草影院精品99| 国产精品国产高清国产av| 日韩亚洲欧美综合| 亚洲精品日韩av片在线观看| 亚洲av成人精品一区久久| 欧美性猛交╳xxx乱大交人| 免费av不卡在线播放| 亚洲成人久久爱视频| 国产一区二区在线观看日韩| 免费观看精品视频网站| а√天堂www在线а√下载| 在线天堂最新版资源| 在线看三级毛片| 久久亚洲精品不卡| 天堂av国产一区二区熟女人妻| 久久久久久久精品吃奶| 成年女人看的毛片在线观看| 欧美色视频一区免费| 老熟妇仑乱视频hdxx| 18禁黄网站禁片午夜丰满| 99国产极品粉嫩在线观看| 成人二区视频| 一区福利在线观看| 中文字幕人妻熟人妻熟丝袜美| 91av网一区二区| 亚洲欧美精品综合久久99| 精华霜和精华液先用哪个| 免费黄网站久久成人精品| 久久婷婷人人爽人人干人人爱| 午夜久久久久精精品| 欧美黑人欧美精品刺激| 中文字幕久久专区| 长腿黑丝高跟| 波多野结衣巨乳人妻| 国产精品女同一区二区软件 | 极品教师在线免费播放| 特级一级黄色大片| 丝袜美腿在线中文| 变态另类成人亚洲欧美熟女| 午夜免费激情av| 亚洲一区二区三区色噜噜| 成人性生交大片免费视频hd| 又粗又爽又猛毛片免费看| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 成人鲁丝片一二三区免费| 成年人黄色毛片网站| 黄色日韩在线| а√天堂www在线а√下载| 久久精品国产亚洲av涩爱 | 免费观看精品视频网站| 我要看日韩黄色一级片| 精品人妻一区二区三区麻豆 | 美女cb高潮喷水在线观看| 亚洲一区高清亚洲精品| 久久精品久久久久久噜噜老黄 | 两人在一起打扑克的视频| 一个人看的www免费观看视频| 国产一区二区在线av高清观看| 精品人妻一区二区三区麻豆 | 少妇人妻一区二区三区视频| 九九爱精品视频在线观看| 香蕉av资源在线| 淫妇啪啪啪对白视频| 露出奶头的视频| 欧美日韩瑟瑟在线播放| 亚洲精品久久国产高清桃花| 亚洲中文字幕一区二区三区有码在线看| 欧美xxxx黑人xx丫x性爽| 国产视频一区二区在线看| 亚洲精品影视一区二区三区av| 国产白丝娇喘喷水9色精品| 午夜亚洲福利在线播放| 亚洲国产精品合色在线| 性色avwww在线观看| 精品一区二区三区人妻视频| 国内精品久久久久久久电影| 亚洲乱码一区二区免费版| 成熟少妇高潮喷水视频| 久久久久久伊人网av| 午夜精品一区二区三区免费看| 国产一区二区三区在线臀色熟女| 亚洲,欧美,日韩| 夜夜爽天天搞| 国产老妇女一区| 精品一区二区三区视频在线| 国产高清视频在线播放一区| 成人性生交大片免费视频hd| 男女啪啪激烈高潮av片| 最新在线观看一区二区三区| 精品久久久久久久久亚洲 | 欧美日韩精品成人综合77777| 欧美国产日韩亚洲一区| 十八禁网站免费在线| 日日夜夜操网爽| 久久久久久大精品| 国产精品国产高清国产av| 男女视频在线观看网站免费| 国产白丝娇喘喷水9色精品| 性色avwww在线观看| 男人舔奶头视频| 久久99热6这里只有精品| 久久国产乱子免费精品| 色尼玛亚洲综合影院| 久久精品影院6| 中国美白少妇内射xxxbb| 嫩草影院新地址|