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

    基于數(shù)字表面模型的冬小麥生物量估算

    2025-04-07 00:00:00郭燕賀佳位盼盼曾凱史舟葉粟楊秀忠鄭國(guó)清王來(lái)剛
    關(guān)鍵詞:遷移能力

    DOI:10.3969/j.issn.2095-1191.2025.01.005

    摘要:【目的】構(gòu)建冬小麥主要生育時(shí)期生物量估算模型,分析不同水處理和不同年份情景下估算模型的遷移能力,為冬小麥生物量快速估算、表型研究及制定作物水肥決策提供技術(shù)支撐。【方法】通過(guò)設(shè)置不同水氮處理,采用大疆M600 Pro無(wú)人機(jī)搭載安洲科技K6多光譜成像儀獲取冬小麥關(guān)鍵生育期影像,提取影像數(shù)字表面模型,基于無(wú)人機(jī)影像提取株高,通過(guò)BP神經(jīng)網(wǎng)絡(luò)構(gòu)建并改進(jìn)冬小麥生物量估算模型?!窘Y(jié)果】水氮耦合自然狀態(tài)條件下冬小麥實(shí)測(cè)株高的變化較小,但在氮充足條件下灌溉可增加冬小麥實(shí)測(cè)株高。無(wú)人機(jī)提取株高與實(shí)測(cè)株高的線性決定系數(shù)(R2)為0.81,即無(wú)人機(jī)提取株高可解釋81%的株高變異?;跓o(wú)人機(jī)遙感影像提取株高構(gòu)建的冬小麥生物量估算模型,R2、均方根誤差(RMSE)、相對(duì)分析誤差(RPD)分別為0.58、4528.23 kg/ha和1.25,說(shuō)明該模型可對(duì)冬小麥生物量進(jìn)行快速估算,但模型穩(wěn)健性較差(RPDlt;1.4),估算值(16198.27 kg/ha)較實(shí)測(cè)值(16960.23 kg/ha)偏小,且估算值較分散。通過(guò)數(shù)據(jù)轉(zhuǎn)換,基于生物量/無(wú)人機(jī)提取株高比值構(gòu)建的冬小麥生物量估算模型R2、RMSE、RPD分別為0.88、2291.90 kg/ha和2.75,改進(jìn)后的模型穩(wěn)健性較強(qiáng)(RPDgt;2.0),估算值(17478.21 kg/ha)與實(shí)測(cè)值(17222.59 kg/ha)較接近,模型估算精度提高了51.72%。經(jīng)驗(yàn)證,改進(jìn)的冬小麥生物量估算模型在不同水處理和不同年份情景下具有較強(qiáng)的遷移能力,遷移估算模型的R2均在0.85以上,能實(shí)現(xiàn)對(duì)冬小麥生物量的精準(zhǔn)快速估算?!窘Y(jié)論】利用無(wú)人機(jī)影像提取株高信息,通過(guò)數(shù)據(jù)轉(zhuǎn)換,能有效提高冬小麥生物量估算模型的估算精度。改進(jìn)的冬小麥生物量估算模型在不同水處理和不同年份情景下均表現(xiàn)出較強(qiáng)的遷移能力,但在不同氮水平情景下的遷移能力存在差異,因此,模型遷移利用前應(yīng)對(duì)不同情景數(shù)據(jù)集進(jìn)行直方圖特征分析,并綜合考慮多種影響因素以提升模型的泛化能力和魯棒性。

    關(guān)鍵詞:冬小麥;生物量;株高;數(shù)字表面模型(DSM);遷移能力

    中圖分類(lèi)號(hào):S127文獻(xiàn)標(biāo)志碼:A文章編號(hào):2095-1191(2025)01-0053-11

    Above-ground biomass estimation of winter wheat based on digital surface model

    GUO Yan1,2,HE Jia1,2,WEI Pan-pan1,2,ZENG Kai1,SHI Zhou3,YE Su3,YANG Xiu-zhong1,2,ZHENG Guo-qing1,WANG Lai-gang1,2*

    (1Institute of Agricultural Information Technology,Henan Academy of Agricultural Sciences/Key Laboratory of Huang-Huai-Hai Smart Agricultural Technology,Ministry of Agriculture and Rural Affairs,Zhengzhou,Henan 450002,China;2Henan Engineering Research Center of Crop Planting Monitoring and Warning,Zhengzhou,Henan 450002,China;3College of Environment and Resources,Zhejiang University,Hangzhou,Zhejiang 310058,China)

    Abstract:【Objective】To construct a biomass estimation model for the key growth stages of winter wheat and analyze the transferability of the estimation model under different water treatments and in different years scenarios,which could provide technical support for the rapid estimation of winter wheat above-ground biomass,phenotypic research,and crop water and fertilizer decision-making.【Method】In this study,by setting different water and nitrogen treatments,theDJI M600 Pro unmanned aerial vehicle(UAV)equipped with the Anzhou Technology K6 multispectral imager was used to acquire images of winter wheat during the key growth stages.The digital surface model(DSM)of the images was ex-tracted,and the plant height was extracted based on the UAV images.The winter wheat above-ground biomass estimation model was constructed and improved through the BP neural network method.【Result】Under the natural condition of water-nitrogen coupling,the change in the measured plant height of winter wheat was relatively small,but irrigation un-der nitrogen-sufficient conditions could increase the measured plant height of winter wheat.The linear determination coef-ficient(R2)between the plant height extracted by the UAV and the measured plant height was 0.81,indicating that the plant height extracted by the UAV could explain 81%of the plant height variation.For the winter wheat above-ground biomass estimation model constructed based on the plant height extracted from UAV remote sensing images,R2,root-mean-square error(RMSE)and relative performance deviation(RPD)were 0.58,4528.23 kg/ha and 1.25 respectively.This showed that the model could rapidly estimate the winter wheat above-ground biomass,but the model had poor ro-bustness(RPDlt;1.4).The estimated value(16198.27 kg/ha)was smaller than the measured value(16960.23 kg/ha),and the estimated values were relatively scattered.Through data transformation,for the winter wheat biomass estimation model constructed based on the ratio(above-graund biomass/plant height extracted by UVA ration,R2,RMSE and RPD were 0.88,2291.90 kg/ha and 2.75 respectively.The improved model had strong robustness(RPDgt;2.0).The estimated value(17478.21 kg/ha)was close to the measured value(17222.59 kg/ha),and the model estimation accuracy has in-creased by 51.72%.It was verified that the improved winter wheat above-ground biomass estimation model had strong transferability under different water treatments and in different years.R2 of the transfer estimation model was above 0.85,achieving accurate and rapid estimation of winter wheat above-ground biomass.【Conclusion】Extracting plant height infor-mation using UAV images and improving the winter wheat above-ground biomass estimation model through data transfor-mation can effectively improve the estimation accuracy of wheat biomass estimation model.The improved winter wheat above-ground biomass estimation model shows strong transfer ability under different water treatments and in different years scenarios.However,there are differences in its transferability under different nitrogen-level scenarios.Therefore,before applying the model for transfer estimation,histogram feature analysis should be carried out on the datasets of different scenarios,and various influencing factors should be comprehensively considered to enhance the generalization ability and robustness of the model.

    Key words:winter wheat;above-ground biomass;plant height;digital surface model(DSM);transfer ability

    Foundation items:National Key Research and Development Program of China(2022YFD2001105);Henan Central Government Guiding Local Scienc and Technology Development Fund Project(Z20231811179);Remote Sensing Innova-tion Team of Henan Academy of Agricultural Sciences Agricultural(2024TD28)

    0引言

    【研究意義】生物量對(duì)光能利用、干物質(zhì)及產(chǎn)量的形成均有重要影響,其精確估算是陸地生態(tài)學(xué)研究重點(diǎn)之一。一方面,作物生物量在全球碳循環(huán)中發(fā)揮重要作用(Piao etal.,2009;Bar-On etal.,2018;黃宏勝等,2024),準(zhǔn)確估算作物生物量有助于掌握氣候變化與農(nóng)田生態(tài)系統(tǒng)相互間的反饋?zhàn)饔?;另一方面,生物量與作物最終產(chǎn)量的形成密切相關(guān)(張傳波等,2022;張寧丹等,2022)。因此,及時(shí)準(zhǔn)確地進(jìn)行作物生物量估算,對(duì)科學(xué)管理與合理利用農(nóng)田、保護(hù)和增強(qiáng)其碳匯功能具有重要意義,同時(shí)為解決糧食安全問(wèn)題提供重要數(shù)據(jù)參考?!厩叭搜芯窟M(jìn)展】近年來(lái),國(guó)內(nèi)外應(yīng)用遙感技術(shù)估算作物生物量的研究越來(lái)越多(Bendig et al.,2015;賀佳等,2017;Adar et al.,2022;王晗等,2023)。不同的遙感數(shù)據(jù)源估算和反演精度存在明顯差異,以衛(wèi)星遙感影像為數(shù)據(jù)源的生物量估算與反演存在空間時(shí)間分辨率低、易受大氣狀況影響等問(wèn)題,一定程度上制約了模型估算精度(鄭陽(yáng)等,2017;Imran et al.,2020;Guo et al.,2022c;張科謙等,2023);近地高光譜影像光譜分辨率高,但耗時(shí)費(fèi)力,且監(jiān)測(cè)范圍受限(賀佳等,2017)。相對(duì)于衛(wèi)星遙感和近地面高光譜遙感,低空無(wú)人機(jī)可依據(jù)研究目的搭載RGB相機(jī)、高光譜成像儀及多光譜相機(jī)等不同類(lèi)型的傳感器,具有操作簡(jiǎn)便、靈活性強(qiáng)等優(yōu)勢(shì),是田塊尺度精準(zhǔn)觀測(cè)的重要手段,在現(xiàn)代化農(nóng)業(yè)定量精細(xì)化監(jiān)測(cè)中發(fā)揮重要作用(陶惠林等,2019;Fu etal.,2021;Guo et al.,2022b)。目前,農(nóng)業(yè)上采用無(wú)人機(jī)遙感進(jìn)行生物量估算主要是利用RGB和多光譜數(shù)字正射影像(Digital ortho-photo map,DOM)提取光譜參數(shù)、計(jì)算植被指數(shù)等構(gòu)建模型。楊俊等(2019)利用無(wú)人機(jī)數(shù)碼影像估算小麥生物量,發(fā)現(xiàn)圖像顏色指數(shù)與紋理特征結(jié)合可提高小麥生物量估算精度,生物量與光譜指數(shù)的相關(guān)系數(shù)最高為0.91。但無(wú)人機(jī)獲取冠層信息時(shí)受作物種植區(qū)域、作物類(lèi)型及作物生育期等多種因素的交互作用,植被指數(shù)存在飽和現(xiàn)象,進(jìn)而影響生物量估算模型的精度和遷移能力。數(shù)字表面模型(Digitalsurface model,DSM)是超高分辨率無(wú)人機(jī)影像的另一關(guān)鍵信息,與作物高度直接關(guān)聯(lián)(Xie et al.,2021;樊意廣等,2022;Qiuetal.,2022)。Chang等(2017)、牛慶林等(2018)、Guo等(2022a)分別利用無(wú)人機(jī)獲取玉米、高粱等作物的數(shù)字遙感影像,基于數(shù)字地面模型(Digital terrain model,DTM)和DSM估算得到作物株高,結(jié)果發(fā)現(xiàn)估算株高與實(shí)測(cè)株高的均方根誤差(RMSE)在0.33~0.88 m。此外,隨著實(shí)時(shí)動(dòng)態(tài)(Real-time kinematic,RTK)和事后動(dòng)態(tài)(Post processed kinematic,PPK)差分技術(shù)的興起,使得無(wú)人機(jī)在快速、準(zhǔn)確獲取株高方面具有顯著優(yōu)勢(shì),越來(lái)越多學(xué)者利用無(wú)人機(jī)影像提取株高并進(jìn)行生物量估算(Niu etal.,2019;劉楊等,2021;Wang et al.,2022)。Niu等(2019)通過(guò)使用精靈Phantom 4直接從點(diǎn)云中提取玉米高度,發(fā)現(xiàn)提取高度與實(shí)測(cè)高度的決定系數(shù)(R2)達(dá)0.90,結(jié)合植被指數(shù)對(duì)玉米生物量進(jìn)行估算的R2為0.82;Wang等(2022)研究表明,采用PPK和RTK校正DSM獲得的作物表面模型(CSM)可有效提高株高估算準(zhǔn)確性;郭燕等(2023a,2023b)研究表明,DSM提取的株高信息結(jié)合植被指數(shù)對(duì)提升作物生物量估算精度具有重要意義?!颈狙芯壳腥朦c(diǎn)】盡管上述研究在各自方向上已取得一定進(jìn)展,但直接利用DSM進(jìn)行生物量估算及模型遷移能力的研究尚無(wú)文獻(xiàn)報(bào)道?!緮M解決的關(guān)鍵問(wèn)題】通過(guò)無(wú)人機(jī)搭載多光譜相機(jī),結(jié)合RTK和地面控制點(diǎn)獲取小麥關(guān)鍵生育期的影像數(shù)據(jù)和高程信息,提取主要生育期的小麥株高,利用BP神經(jīng)網(wǎng)絡(luò)法構(gòu)建小麥主要生育時(shí)期生物量估算模型,并分析不同水處理和不同年份情景下估算模型的遷移能力,為冬小麥生物量快速估算、表型研究及制定作物水肥決策提供技術(shù)支撐。

    1數(shù)據(jù)來(lái)源與研究方法

    1.1研究區(qū)概況及田間試驗(yàn)設(shè)計(jì)

    于2020—2021年和2021—2022年,在河南省周口市商水縣國(guó)營(yíng)農(nóng)場(chǎng)(33°33′N(xiāo),114°37′E,海拔49 m)開(kāi)展冬小麥種植試驗(yàn),小區(qū)面積約42m2。該區(qū)域地處暖溫帶南部,為季風(fēng)半濕潤(rùn)氣候,全年溫度適宜。供試土壤為砂姜黑土,播種前0~20 cm土層基礎(chǔ)養(yǎng)分情況:pH 7.3,有機(jī)質(zhì)25.0 g/kg,全氮1.69 g/kg,速效磷22.6 mg/kg,速效鉀139.6 mg/kg。采用隨機(jī)區(qū)組試驗(yàn)設(shè)計(jì),設(shè)5個(gè)氮水平和2個(gè)水處理。氮水平分別為N0(0 kg/ha)、N6(90 kg/ha)、N12(180 kg/ha)、N18(270 kg/ha)和N24(360 kg/ha),每個(gè)水平重復(fù)3次;水處理為自然降水(W0)和因需灌溉(W1),總計(jì)30個(gè)小區(qū)。氮肥50%作為底肥施入,50%于拔節(jié)期追施;磷肥用量120 kg/ha,鉀肥用量90 kg/ha,均為播種前一次性施入。

    1.2研究技術(shù)路線

    總體研究技術(shù)路線如圖1所示。首先,通過(guò)大疆M600 Pro無(wú)人機(jī)獲取冬小麥拔節(jié)期、開(kāi)花期和灌漿期的遙感影像,利用Agisoft PhotoScan Professio-nal 12.0對(duì)影像進(jìn)行拼接及輻射校正等預(yù)處理,在3D點(diǎn)云建立過(guò)程中創(chuàng)建DSM,以獲得研究區(qū)植被和土壤的DSM數(shù)據(jù)。植被DSM數(shù)據(jù)減去土壤DSM數(shù)據(jù)即獲得冬小麥無(wú)人機(jī)提取株高(Hdsm)。然后,基于無(wú)人機(jī)提取株高、實(shí)測(cè)株高與生物量的關(guān)系,確認(rèn)株高對(duì)生物量的貢獻(xiàn),通過(guò)BP神經(jīng)網(wǎng)絡(luò)采用無(wú)人機(jī)提取株高構(gòu)建冬小麥生物量估算模型。同時(shí),采用生物量與株高的比值改進(jìn)生物量估算模型,并在不同水處理和不同年份情景下進(jìn)行模型遷移能力分析與應(yīng)用,探究影響模型遷移能力的原因,并給出研究建議。

    1.3數(shù)據(jù)獲取與處理

    1.3.1無(wú)人機(jī)數(shù)據(jù)獲取與處理在冬小麥生長(zhǎng)周期的關(guān)鍵生育期(拔節(jié)期、開(kāi)花期和灌漿期),分別選取天氣晴朗的10:00—14:00,采用大疆M600 Pro無(wú)人機(jī)搭載安洲科技K6多光譜成像儀進(jìn)行影像獲取。飛行前地面鋪設(shè)標(biāo)準(zhǔn)反射率為5%、20%、40%和70%的4塊輻射定標(biāo)板,同時(shí)使用RTK接收機(jī)測(cè)量試驗(yàn)區(qū)4個(gè)頂點(diǎn)及中間點(diǎn)的地面控制點(diǎn)(Ground control point,GCP)坐標(biāo),用于飛行位置和海拔高度精準(zhǔn)校正。飛行高度50m,獲取影像空間分辨率為0.02 m。影像獲取后按格式轉(zhuǎn)換、影像篩選、影像拼接、正射校正及輻射定標(biāo)等5個(gè)步驟進(jìn)行處理。

    1.3.2冬小麥無(wú)人機(jī)株高提取冬小麥無(wú)人機(jī)株高提取步驟如圖2所示。首先將無(wú)人機(jī)數(shù)據(jù)導(dǎo)入Agisoft PhotoScan Professional 12.0中,檢查無(wú)人機(jī)影像圖片質(zhì)量,生成點(diǎn)云并對(duì)齊照片;然后建立格網(wǎng)并輸入研究區(qū)GCP進(jìn)行幾何校正,重建DSM后輸出CSM;最后將得到的CSM數(shù)據(jù)導(dǎo)入到ArcGIS 10.6中,采用柵格計(jì)算得到無(wú)人機(jī)提取株高。此外,于2020年10月5日采集裸地影像建立研究區(qū)的地形數(shù)字高程模型(Digital elevation model,DEM)。值得注意的是,在影像拼接建立CSM的過(guò)程中,應(yīng)保證控制點(diǎn)處的高程與實(shí)測(cè)高程最大偏差小于0.01 m。根據(jù)小區(qū)范圍繪制矩形感興趣區(qū)域(Region of interest,ROI),再通過(guò)對(duì)標(biāo)準(zhǔn)化處理后的CSM進(jìn)行分區(qū)統(tǒng)計(jì),即可得到各小區(qū)的冬小麥無(wú)人機(jī)提取株高。

    1.3.3株高和生物量測(cè)定地面數(shù)據(jù)獲取與無(wú)人機(jī)影像的獲取同步進(jìn)行。田間樣本采集時(shí)使用北斗智能手持終端(UG908)獲取采樣點(diǎn)地理位置信息,在每小區(qū)各選取3處長(zhǎng)勢(shì)均勻的10株冬小麥測(cè)量株高,然后將0.5 m雙行的冬小麥地上植株樣本裝入密封袋。送回實(shí)驗(yàn)室,對(duì)莖、葉進(jìn)行分離,105℃下殺青,80℃烘干至恒重,稱(chēng)重莖、葉干重之和即為生物量,然后換算為kg/ha,共計(jì)獲取180個(gè)樣本數(shù)據(jù)。

    1.4 BP神經(jīng)網(wǎng)絡(luò)

    BP神經(jīng)網(wǎng)絡(luò)是一種采用誤差逆向傳播進(jìn)行算法訓(xùn)練的多層前饋網(wǎng)絡(luò),是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型(Du etal.,2021;Lan etal.,2023)。BP神經(jīng)網(wǎng)絡(luò)常用激活函數(shù)包括identity、sigmoid及ReLU等20多種,當(dāng)函數(shù)identity處于活動(dòng)狀態(tài)時(shí),節(jié)點(diǎn)的輸入等于輸出,最適合于潛在行為呈線性(類(lèi)似線性回歸)的任務(wù)(李航,2022),因此本研究采用identity作為訓(xùn)練模型的激活函數(shù)。本研究共設(shè)3層網(wǎng)絡(luò)結(jié)構(gòu),通過(guò)準(zhǔn)牛頓方法族優(yōu)化器(lbfgs)提高運(yùn)行速度,同時(shí)為防止過(guò)擬合,而引入學(xué)習(xí)率及正則化等參數(shù)對(duì)模型進(jìn)行優(yōu)化(周志華,2022)。

    1.5模型評(píng)價(jià)指標(biāo)

    冬小麥生物量估算模型評(píng)價(jià)指標(biāo)3個(gè),分別是決定系數(shù)(R2)、RMSE和相對(duì)分析誤差(RPD)。R2與PRD越大,RMSE越小,模型估算和遷移能力越好。當(dāng)估算值與實(shí)測(cè)值完全一致時(shí),RMSE為0,此時(shí)可視為完美模型;RMSE越大,模型誤差越大,其精度較低;反之,RMSE越小,模型精確度越高。RPD≥2.0時(shí)表明模型穩(wěn)健性較好,1.4≤RPDlt;2.0時(shí)表明模型穩(wěn)健性中等,RPDlt;1.4時(shí)表明模型穩(wěn)健性較差(Changetal.,2001;郭燕等,2023a)。

    2冬小麥生物量模型構(gòu)建及遷移估算

    2.1冬小麥實(shí)測(cè)株高和生物量統(tǒng)計(jì)結(jié)果

    不同水氮處理下冬小麥的實(shí)測(cè)株高和生物量變化見(jiàn)圖3。W0和W1處理均表現(xiàn)為冬小麥實(shí)測(cè)株高和生物量隨著氮水平的升高而呈逐漸增加趨勢(shì),但增加幅度存在明顯差異。當(dāng)?shù)接蒒0升至N12,W0和W1處理冬小麥實(shí)測(cè)株高的增幅分別為23.22和21.33 cm;氮水平由N12升至N24,對(duì)應(yīng)的增幅分別為2.11和0.89 cm。在冬小麥生物量方面,當(dāng)?shù)接蒒0升至N18,W0和W1處理的生物量增幅分別為15675.24和13912.85 kg/ha;氮水平由N18升至N24,W0處理的生物量基本無(wú)明顯變化,而W1處理的生物量仍然呈快速增加趨勢(shì)。

    從不同氮水平冬小麥實(shí)測(cè)株高和生物量的變化趨勢(shì)可看出,隨著氮肥施用量的增加,冬小麥實(shí)測(cè)株高和生物量均呈明顯的上升趨勢(shì),但增加幅度不同,冬小麥實(shí)測(cè)株高從N12開(kāi)始趨于穩(wěn)定,生物量則呈持續(xù)上升趨勢(shì)。單因素方差分析(One-way ANOVA)結(jié)果顯示,冬小麥實(shí)測(cè)株高在N0、N6和N12水平下存在顯著差異(Plt;0.05,下同),生物量在N0、N6、N12和N18水平下存在顯著差異。綜上所述,在水氮耦合過(guò)程中,水分對(duì)冬小麥實(shí)測(cè)株高產(chǎn)生的脅迫效應(yīng)大于生物量,肥料需要依靠水分才能被作物吸收,即充足水分是保證群體株高和生物量協(xié)調(diào)的基礎(chǔ)。

    2.2冬小麥無(wú)人機(jī)提取株高與實(shí)測(cè)株高的關(guān)系

    冬小麥無(wú)人機(jī)提取株高與實(shí)測(cè)株高的關(guān)系見(jiàn)圖4。無(wú)人機(jī)提取株高與實(shí)測(cè)株高的相關(guān)性隨氮水平的升高而下降,尤其在氮水平由N6升至N12時(shí)快速下降,R2從0.59降至0.08,究其原因是相同灌溉條件下,冬小麥株高受品種基因型的影響,灌漿期的株高趨于穩(wěn)定,氮水平的影響減弱(王聲鋒等,2010)。此外,氮水平由N12升至N24時(shí),無(wú)人機(jī)提取株高分布一致性較強(qiáng),但在N0和N6水平下無(wú)人機(jī)提取株高的分布相對(duì)分散,說(shuō)明氮水平在一定范圍可有效增加冬小麥株高,但超過(guò)某一閾值后對(duì)株高的影響變?nèi)酢?/p>

    不同水處理下,無(wú)人機(jī)提取株高與實(shí)測(cè)株高的相關(guān)性表現(xiàn)為W0處理優(yōu)于W1處理,R2分別為0.81和0.79。結(jié)合圖3可知,水氮耦合自然狀態(tài)條件下冬小麥株高的變化相對(duì)較小,但在氮充足條件下灌溉可增加冬小麥株高。無(wú)人機(jī)提取株高與實(shí)測(cè)株高的線性R2為0.81,即無(wú)人機(jī)提取株高可解釋81%的株高變異,不僅為冬小麥株高測(cè)量提供了簡(jiǎn)便快捷途徑,還為冬小麥生物量精準(zhǔn)估算提供了支撐。

    2.3冬小麥無(wú)人機(jī)提取株高與生物量的關(guān)系

    無(wú)人機(jī)提取冬小麥株高與生物量的關(guān)系見(jiàn)圖5。無(wú)人機(jī)提取株高與生物量的相關(guān)性隨氮水平升高的變化規(guī)律不明顯,但氮的效應(yīng)影響逐漸顯現(xiàn)出來(lái)。不同水處理下,無(wú)人機(jī)提取株高與生物量的相關(guān)性表現(xiàn)為W0處理優(yōu)于W1處理,R2分別為0.69和0.50,與部分株高的變化特征具有一致性。水分對(duì)冬小麥生物量和株高的影響存在明顯差異,主要是冬小麥株高變化趨于穩(wěn)定時(shí),灌溉使得生物量的增加更傾向于葉片和莖稈生物量的增加,對(duì)株高的影響相對(duì)較小??傮w來(lái)看,無(wú)人機(jī)提取株高與生物量建立簡(jiǎn)單線性關(guān)系模型的R2為0.57,即無(wú)人機(jī)提取株高可解釋57%的生物量變異,說(shuō)明利用無(wú)人機(jī)遙感影像提取的株高可用于冬小麥生物量快速估算。

    2.4冬小麥生物量估算模型構(gòu)建與改進(jìn)

    為實(shí)現(xiàn)對(duì)冬小麥生物量進(jìn)行快速估算,首先采用BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練無(wú)人機(jī)提取株高與生物量的關(guān)系模型(R2為0.57),然后基于訓(xùn)練模型對(duì)冬小麥生物量進(jìn)行估算,生物量的實(shí)測(cè)值和估算值關(guān)系見(jiàn)圖6。模型的R2、RMSE、RPD分別為0.58、4528.23 kg/ha和1.25,說(shuō)明該模型可對(duì)冬小麥生物量進(jìn)行快速估算,但模型穩(wěn)健性較差(RPDlt;1.4),估算值(16198.27 kg/ha)較實(shí)測(cè)值(16960.23 kg/ha)小,且估算值較分散。綜上所述,雖然采用無(wú)人機(jī)提取株高可對(duì)冬小麥生物量進(jìn)行估算,但由于模型的魯棒性較弱,可能會(huì)造成估算結(jié)果不準(zhǔn)確。這在圖6中95%置信區(qū)間覆蓋的數(shù)據(jù)上也得到驗(yàn)證,置信區(qū)間數(shù)據(jù)分散程度較高,且部分?jǐn)?shù)據(jù)并未包含在置信區(qū)間內(nèi),估算結(jié)果不確定性大。

    通過(guò)數(shù)據(jù)轉(zhuǎn)換發(fā)現(xiàn),冬小麥生物量/無(wú)人機(jī)提取株高比值隨氮水平的升高呈近似平行的遞增趨勢(shì)(圖7)。為此,提出一種基于比值的生物量快速估算方法,以期提高模型的估算精度。采用與上述相同的方法構(gòu)建模型對(duì)冬小麥生物量進(jìn)行估算,訓(xùn)練模型R2為0.89,生物量實(shí)測(cè)值和估算值的關(guān)系見(jiàn)圖8。模型的R2、RMSE、RPD分別為0.88、2291.90 kg/ha和2.75,改進(jìn)后的模型穩(wěn)健性明顯增強(qiáng)(RPDgt;2.0),說(shuō)明該模型可對(duì)冬小麥生物量進(jìn)行精準(zhǔn)估算,估算值(17478.21 kg/ha)與實(shí)測(cè)值(17222.59 kg/ha)較接近,即基于比值構(gòu)建的模型魯棒性顯著提升。與直接采用無(wú)人機(jī)提取株高建立冬小麥生物量估算模型相比,基于比值構(gòu)建的模型估算精度提高了51.72%,且模型估算置信區(qū)間較集中,能實(shí)現(xiàn)冬小麥生物量的簡(jiǎn)潔高效估算。

    假設(shè)冬小麥生物量(y)與比值(x)的改進(jìn)模型表示為y=f(x),將數(shù)據(jù)代入公式后得到模型為y=kx-4004.8。隨著拔節(jié)期、開(kāi)花期和灌漿期的推進(jìn),生物量不斷增加,采用多次截取生物量數(shù)據(jù)集的方法,形成相對(duì)應(yīng)的不同生物量與無(wú)人機(jī)提取株高數(shù)據(jù)集,采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建冬小麥生物量估計(jì)模型,k與無(wú)人機(jī)提取株高間存在k≈Hmax的關(guān)系,Hmax表示每個(gè)數(shù)據(jù)集中無(wú)人機(jī)提取株高的最大值。因此,將改進(jìn)模型進(jìn)行轉(zhuǎn)換,得到冬小麥生物量估算模型:AGB=4004.80×Hdsm/(Hmax-Hdsm)。該模型可對(duì)冬小麥生物量進(jìn)行精準(zhǔn)估算。

    2.5不同情景下改進(jìn)冬小麥生物量估算模型的遷移能力

    基于W0和W1處理數(shù)據(jù)集分別構(gòu)建冬小麥生物量估算模型,并對(duì)W1和W0處理下的生物量進(jìn)行遷移估算。將W0處理數(shù)據(jù)集構(gòu)建的模型遷移到W1處理,模型的R2、RMSE、RPD分別為0.90、2407.44 kg/ha和2.72;而將W1處理數(shù)據(jù)集構(gòu)建的模型遷移到W0處理,模型的R2、RMSE、RPD分別為0.88、2496.97 kg/ha和2.54。同時(shí),對(duì)模型在不同年份下進(jìn)行遷移估算,設(shè)定2020—2021年冬小麥生長(zhǎng)季節(jié)為n,2021—2022年冬小麥生長(zhǎng)季節(jié)為n+1,基于n和n+1年份數(shù)據(jù)集分別構(gòu)建冬小麥生物量估算模型,對(duì)n+1年份和n年份情境下的生物量進(jìn)行遷移估算,將n年份數(shù)據(jù)集構(gòu)建的模型遷移到n+1年份,模型的R2、RMSE、RPD分別為0.88、2257.36 kg/ha和2.75;將n+1年份數(shù)據(jù)集構(gòu)建的生物量模型遷移到n年份,模型的R2、RMSE、RPD分別為0.86、2239.37 kg/ha和2.93。綜合R2、RMSE和RPD,改進(jìn)的冬小麥生物量估算模型在不同水處理和不同年份情景下均表現(xiàn)出較好的魯棒性和遷移能力,能實(shí)現(xiàn)對(duì)冬小麥生物量的精準(zhǔn)快速估算。

    不同情景下的模型遷移能力存在差異,不同年份間的遷移估算效果優(yōu)于不同水處理,究其原因可能是由于數(shù)據(jù)特征存在差異。對(duì)不同水處理和不同年份用于訓(xùn)練模型的數(shù)據(jù)集進(jìn)行統(tǒng)計(jì)分析,結(jié)果(圖9)顯示,W0和W1處理,n年份和n+1年份實(shí)測(cè)冬小麥生物量與生物量/無(wú)人機(jī)提取株高數(shù)據(jù)集的平均值、變程、偏度系數(shù)和峰度系數(shù)較相近,且擬合曲線特征具有相似性,是不同情景條件下模型具有較強(qiáng)遷移能力的保證。不同水處理和不同年份間的數(shù)據(jù)特征也存在差異,W0、W1處理的生物量/無(wú)人機(jī)提取株高平均值分別為222.31和206.91,二者相差15.40;n、n+1年份的生物量/無(wú)人機(jī)提取株高平均比值分別為214.09和215.14,二者相差1.05;即不同年份情景下的差值小于不同水處理。同樣,變程、偏度系數(shù)和峰度系數(shù)也存在相似趨勢(shì),使得不同年份情境下的模型遷移效果優(yōu)于不同水處理。

    3討論

    本研究在分析實(shí)測(cè)株高和無(wú)人機(jī)提取株高與冬小麥生物量之間關(guān)系的基礎(chǔ)上,基于無(wú)人機(jī)提取株高,采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建并改進(jìn)冬小麥生物量估算模型,模型估算效果良好,且在不同水處理和不同年份情景下表現(xiàn)出較強(qiáng)的遷移能力。與實(shí)測(cè)株高相比,無(wú)人機(jī)提取株高整體上偏低,主要是冬小麥冠層影像中存在一些土壤背景,且部分冬小麥三維結(jié)構(gòu)空間信息可能缺失,而造成無(wú)人機(jī)提取株高相對(duì)偏低(陶惠林等,2019;劉楊等,2021)。這種情況在一定程度上會(huì)影響生物量估算模型的精度,為此,本研究通過(guò)對(duì)數(shù)據(jù)進(jìn)行比值變換,有效提高了冬小麥生物量估算模型的精度,且改進(jìn)的模型在不同水處理和不同年份情景下能獲得較好的遷移估算效果,與Guo等(2022b)對(duì)不同水處理植株氮含量模型進(jìn)行遷移能力分析得到的結(jié)果具有一致性。遷移估算是將一個(gè)環(huán)境的知識(shí)應(yīng)用到另一環(huán)境中,理想的遷移估算模型要求原訓(xùn)練集不僅有較大的數(shù)據(jù)量,還對(duì)數(shù)據(jù)質(zhì)量要求較高,模型在遷移至目標(biāo)模型時(shí)將包含更多的特征信息,因此,進(jìn)行遷移估算時(shí)首先要對(duì)數(shù)據(jù)的分布特征進(jìn)行分析,是保證模型能遷移并獲得較好的估算效果(李晨光等,2022;周志華,2022)。本研究采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建模型時(shí),訓(xùn)練數(shù)據(jù)集綜合了研究區(qū)冬小麥生物量的多個(gè)關(guān)鍵生育時(shí)期信息,且訓(xùn)練數(shù)據(jù)和測(cè)試數(shù)據(jù)的統(tǒng)計(jì)分析指標(biāo)擬合度高,是保證冬小麥生物量估算模型精度和魯棒性均較優(yōu)的重要原因。在實(shí)際應(yīng)用中,根據(jù)數(shù)據(jù)的基本特征有針對(duì)性地選擇預(yù)訓(xùn)練模型進(jìn)行遷移,更有助于提高模型估算效果(Zhu et al.,2021)。

    本研究基于無(wú)人機(jī)提取株高構(gòu)建的冬小麥生物量估算模型,雖然在某些情況下低于綜合利用植被指數(shù)、紋理特征等數(shù)據(jù)進(jìn)行生物量估算的模型精度(陶惠林等,2019;Luo et al.,2022),但該模型僅需無(wú)人機(jī)遙感影像的提取株高信息,節(jié)省了指數(shù)計(jì)算、紋理特征提取等步驟,具有簡(jiǎn)便快捷的特點(diǎn)。Steduto和Albrizio(2005)、韓文霆等(2021)研究認(rèn)為,大田作物絕大部分生育期的水分等資源利用情況不變,基于無(wú)人機(jī)遙感影像提取單一指標(biāo)可對(duì)其生物量進(jìn)行估算。不同的建模方法對(duì)模型估算結(jié)果也存在一定影響(Xu etal.,2022)。本研究采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建冬小麥生物量估算模型,采用的設(shè)備參數(shù)為Intel Core i7-9700K CPU和64 GB RAM,模型訓(xùn)練用時(shí)為0.044 s,實(shí)測(cè)值和預(yù)測(cè)值間的R2為0.88。與Wang等(2022)的研究相比,本研究構(gòu)建的冬小麥生物量估算模型精度優(yōu)于直線模型(R2=0.46~0.93)和PLSR模型(R2=0.50~0.75),但低于RF模型(R2=0.95~0.97)。這主要是由于RF是一種集成機(jī)器學(xué)習(xí)算法,進(jìn)行模型構(gòu)建時(shí)多種指標(biāo)的參與不斷進(jìn)行挑選,促使構(gòu)建的生物量估算模型精度隨之得到有效提升,但在進(jìn)行多種指標(biāo)提取時(shí)顯然要花費(fèi)更多時(shí)間,且增加了過(guò)程的復(fù)雜性。

    本研究利用無(wú)人機(jī)影像提取株高信息,通過(guò)數(shù)據(jù)變換,改進(jìn)冬小麥生物量估算模型,模型估算精度(R2=0.88)較原模型(R2=0.58)提高了0.30,且在不同水處理和不同年份情景下表現(xiàn)出較強(qiáng)的遷移能力,但在不同品種或面對(duì)極端氣候、病蟲(chóng)害侵襲時(shí),模型估算精度可能會(huì)出現(xiàn)波動(dòng)。因此,后續(xù)研究還需從以下2個(gè)方面進(jìn)行完善:(1)冬小麥生物量估算模型遷移利用前應(yīng)對(duì)不同情景數(shù)據(jù)集進(jìn)行直方圖特征分析。數(shù)據(jù)集特征差異是導(dǎo)致生物量估算模型遷移能力不同的關(guān)鍵因子,數(shù)據(jù)集直方圖特征分析為模型遷移應(yīng)用提供了前置條件,也為模型精度的提升奠定了基礎(chǔ)。(2)綜合考慮多種因素以提升模型的泛化能力和魯棒性。作物生物量的構(gòu)成和估算受多種因素影響,除了水分、氮水平和年份等因子外,生物量估算模型可能還受品種及建模方法等因素的影響,建議多角度、多維度綜合考慮冬小麥生物量估算的影響因素,進(jìn)一步改善模型的泛化能力和魯棒性,以提升模型的適用性及實(shí)用性。

    4結(jié)論

    利用無(wú)人機(jī)影像提取株高信息,通過(guò)數(shù)據(jù)變換,能有效提高冬小麥生物量估算模型的估算精度。改進(jìn)后的冬小麥生物量估算模型在不同水處理和不同年份情景下均表現(xiàn)出較強(qiáng)的遷移能力,但在不同氮水平情景下的遷移能力存在差異,因此,模型遷移利用前應(yīng)對(duì)不同情景數(shù)據(jù)集進(jìn)行直方圖特征分析,并綜合考慮多種影響因素以提升模型的泛化能力和魯棒性。

    參考文獻(xiàn)(References):

    樊意廣,馮海寬,劉楊,邊明博,孟煬,楊貴軍.2022.基于冠層光譜特征和株高的馬鈴薯植株氮含量估算[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),53(6):202-208.[Fan Y G,F(xiàn)eng H K,Liu Y,Bian M B,Meng Y,Yang G J.2022.Estimation of potato plant nitrogen content based on canopy spectral characteristics and plant height[J].Transactions of the Chinese Society for Agricultural Machinery,53(6):202-208.]doi:10.6041/j.issn.1000-1298.2022.06.021.

    郭燕,井宇航,賀佳,張會(huì)芳,賈德偉,王來(lái)剛.2023a.基于無(wú)人機(jī)影像多源信息的冬小麥生物量與產(chǎn)量估算[J].河南農(nóng)業(yè)科學(xué),52(12):149-161.[Guo Y,Jing Y H,He J,Zhang H F,Jia D W,Wang L G.2023a.Winter wheat aboveground biomass and yield estimation based on multi-source information from UAV imagery[J].Journal of Henan Agricultural Sciences,52(12):149-161.]doi:10.15933/j.cnki.1004-3268.2023.12.017.

    郭燕,井宇航,王來(lái)剛,黃競(jìng)毅,賀佳,馮偉,鄭國(guó)清.2023b.基于無(wú)人機(jī)影像特征的冬小麥植株氮含量預(yù)測(cè)及模型遷移能力分析[J].中國(guó)農(nóng)業(yè)科學(xué),56(5):850-865.[Guo Y,Jing Y H,Wang L G,Huang J Y,He J,F(xiàn)eng W,Zheng G Q.2023b.UAV multispectral image-based nitrogen con-tent prediction and the transferability analysis of the mo-delsin winter wheat plant[J].Scientia Agricultura Sinica,56(5):850-865.]doi:10.3864/j.issn.0578-1752.2023.05.004.

    韓文霆,湯建棟,張立元,牛亞曉,王彤華.2021.基于無(wú)人機(jī)遙感的玉米水分利用效率與生物量監(jiān)測(cè)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),52(5):129-141.[Han W T,Tang J D,Zhang LY,Niu Y X,Wang T H.2021.Maize water use efficiency and bio-mass estimation based on unmanned aerial vehicle remote sensing[J].Transactions of the Chinese Society for Agri-cultural Machinery,52(5):129-141.]doi:10.6041/j.issn.1000-1298.2021.05.014.

    賀佳,劉冰峰,郭燕,王來(lái)剛,鄭國(guó)清,李軍.2017.冬小麥生物量高光譜遙感監(jiān)測(cè)模型研究[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),23(2):313-323.[He J,Liu B F,Guo Y,Wang L G,Zheng G Q,Li J.2017.Biomass estimation model of winter wheat(Triticum aestivum L.)using hyperspectral reflec-tances[J].Journal of Plant Nutrition and Fertilizer,23(2):313-323.]doi:10.11674/zwyf.16173.

    黃宏勝,張馨月,居輝,韓雪.2024.大氣CO2濃度升高背景下冬小麥冠層光譜特征和地上生物量估算[J].作物學(xué)報(bào),50(4):991-1003.[Huang H S,Zhang X Y,Ju H,Han X.2024.Spectral characteristics of winter wheat canopy and estimation of aboveground biomass under elevated atmo-spheric CO2 concentration[J].Acta Agronomica Sinica,50(4):991-1003.]doi:10.3724/SP.J.1006.2024.31041.

    李晨光,張波,趙騫,陳小平,王行甫.2022.基于遷移學(xué)習(xí)的文本共情預(yù)測(cè)[J].計(jì)算機(jī)應(yīng)用,42(11):3603-3609.[Li C G,Zhang B,Zhao Q,Chen X P,Wang X F.2022.Empa-thy prediction from texts based on transfer learning[J].Journal of Computer Applications,42(11):3603-3609.]doi:10.11772/j.issn.1001-9081.2021091632.

    李航.2022.機(jī)器學(xué)習(xí)方法[M].北京:清華大學(xué)出版社.[Li H.2022.Machine learning methods[M].Beijing:Tsing-hua University Press.]

    劉楊,馮海寬,黃玨,孫乾,楊福芹,楊貴軍.2021.基于無(wú)人機(jī)高光譜影像的馬鈴薯株高和地上生物量估算[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),52(2):188-198.[Liu Y,F(xiàn)eng H K,Huang J,Sun Q,Yang F Q,Yang G J.2021.Estimation of potato plant height and above-ground biomass based on UAV hyper-spectral images[J].Transactions of the Chinese Society for Agricultural Machinery,52(2):188-198.]doi:10.6041/j.issn.1000-1298.2021.02.017.

    牛慶林,馮海寬,楊貴軍,李長(zhǎng)春,楊浩,徐波,趙衍鑫.2018.基于無(wú)人機(jī)數(shù)碼影像的玉米育種材料株高和LAI監(jiān)測(cè)[J].農(nóng)業(yè)工程學(xué)報(bào),34(5):73-82.[Niu Q L,F(xiàn)eng H K,Yang G J,Li C C,Yang H,Xu B,Zhao Y X.2018.Moni-toring plant height and leaf area index of maize breeding material based on UAV digital images[J].Transactions ofthe Chinese Society of Agricultural Engineering,34(5):73-82.]doi:10.11975/j.issn.1002-6819.2018.05.010.

    陶惠林,徐良驥,馮海寬,楊貴軍,楊小冬,苗夢(mèng)珂,代陽(yáng).2019.基于無(wú)人機(jī)數(shù)碼影像的冬小麥株高和生物量估算[J].農(nóng)業(yè)工程學(xué)報(bào),35(19):107-116.[Tao H L,Xu L J,F(xiàn)eng H K,Yang G J,Yang X D,Miao M K,Dai Y.2019.Estimation of plant height and biomass of winter wheat based on UAV digital image[J].Transactions of the Chi-nese Society of Agricultural Engineering,35(19):107-116.]doi:10.11975/j.issn.1002-6819.2019.19.013.

    王晗,向友珍,李汪洋,史鴻棹,王辛,趙笑.2023.基于無(wú)人機(jī)多光譜遙感的冬油菜地上部生物量估算[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),54(8):218-229.[Wang H,Xiang Y Z,Li W Y,Shi H Z,Wang X,Zhao X.2023.Estimation of winter rapeseed above-ground biomass based on UAV multi-spectral remote sensing[J].Transactions of the Chinese Society for Agricultural Machinery,54(8):218-229.]doi:10.6041/j.issn.1000-1298.2023.08.021.

    王聲鋒,段愛(ài)旺,徐建新.2010.冬小麥株高和葉面積指數(shù)變化動(dòng)態(tài)分析及模擬模型[J].灌溉排水學(xué)報(bào),29(4):97-100.[Wang S F,Duan A W,Xu J X.2010.Dynamic changes and simulation model of plant height and leaf areaindex of winter wheat[J].Journal of Irrigation and Drai-nage,29(4):97-100.]doi:10.13522/j.cnki.ggps.2010.04.027.

    楊俊,丁峰,陳晨,劉濤,孫成明,丁大偉,霍中洋.2019.小麥生物量及產(chǎn)量與無(wú)人機(jī)圖像特征參數(shù)的相關(guān)性[J].農(nóng)業(yè)工程學(xué)報(bào),35(23):104-110.[Yang J,Ding F,Chen C,Liu T,Sun C M,Ding D W,Huo Z Y.2019.Correlation of wheat biomass and yield with UAV image characteristic parameters[J].Transactions of the Chinese Society of Agri-cultural Engineering,35(23):104-110.]doi:10.11975/j.issn.1002-6819.2019.23.013.

    張傳波,李衛(wèi)國(guó),張宏,李偉,馬廷淮,張琤琤,陳華.2022.遙感光譜指標(biāo)和神經(jīng)網(wǎng)絡(luò)結(jié)合的冬小麥地上部生物量估測(cè)[J].麥類(lèi)作物學(xué)報(bào),42(5):631-639.[Zhang C B,Li W G,Zhang H,Li W,Ma T H,Zhang C C,Chen H.2022.Estimation of winter wheat aboveground biomass based on remote sensing spectral index and neural network[J].Jour-nal of Triticeae Crops,42(5):631-639.]doi:10.7606/j.issn.1009-1041.2022.05.14.

    張科謙,程鋼,吳微,宋向陽(yáng),張子謙,姚順,吳帥.2023.融合主被動(dòng)遙感影像的冬小麥種植面積提取研究[J].河南農(nóng)業(yè)科學(xué),52(6):160-171.[Zhang K Q,Cheng G,Wu W,Song X Y,Zhang Z Q,Yao S,Wu S.2023.Extraction of winter wheat planting area based on fused active and pas-sive remote sensing images[J].Journal of Henan Agricul-tural Sciences,52(6):160-171.]doi:10.15933/j.cnki.1004-3268.2023.06.017.

    張寧丹,任建強(qiáng),吳尚蓉.2022.基于花后累積地上生物量比例的冬小麥動(dòng)態(tài)收獲指數(shù)估算[J].農(nóng)業(yè)工程學(xué)報(bào),38(7):189-199.[Zhang N D,Ren J Q,Wu S R.2022.Esti-mating the dynamic harvest index of winter wheat usingthe fraction of accumulated aboveground biomass afterflowering[J].Transactions of the Chinese Society of Agri-cultural Engineering,38(7):189-199.]doi:10.11975/j.issn.1002-6819.2022.07.021.

    鄭陽(yáng),吳炳方,張淼.2017.Sentinel-2數(shù)據(jù)的冬小麥地上干生物量估算及評(píng)價(jià)[J].遙感學(xué)報(bào),21(2):318-328.[Zheng Y,Wu B F,Zhang M.2017.Estimating the above ground biomass of winter wheat using the Sentinel-2 data[J].Jour-nal of Remote Sensing,21(2):318-328.]doi:10.11834/jrs.20176269.

    周志華.2016.機(jī)器學(xué)習(xí)[M].北京:清華大學(xué)出版社.[Zhou Z H.2016.Machine learning[M].Beijing:Tsinghua Uni-versity Press.]

    Adra S,Sternberg M,Paz-Kagan T,Henkin Z,Dovrat G,Zaady E,Argaman E.2022.Estimation of aboveground biomass production using an unmanned aerial vehicle(UAV)and VENμS satellite imagery in Mediterranean and semiarid rangelands[J].Remote Sensing Applications:Society and Environment,26:100753.doi:10.1016/j.rsase.2022.100 753.

    Bar-On Y M,Phillips R,Milo R.2018.The biomass distribu-tion on earth[J].Proceedings of the National Academy ofSciences of the United States of America,1159(25):6506-6511.doi:10.1073/pnas.1711842115.

    Bendig J,Yu K,Aasen H,Bokten A,Bennerta S,Broscheit J,Gnyp M L,Bareth G.2015.Combining UAV-based plant height from crop surface models,visible,and near infrared vegetation indices for biomass monitoring in barley[J].International Journal of Applied Earth Observation and Geoinformation,39:79-87.doi:10.1016/j.jag.2015.02.012.

    Chang A J,Jung J H,Maeda M M,Landivar J.2017.Crop height monitoring with digital imagery from Unmanned Aerial System(UAS)[J].Computers and Electronics in Agriculture,141:232-237.doi:10.1016/j.compag.2017.07.008.

    Chang C W,Laird D A,Mausbach M J,Hurburgh C R.2001.Near-infrared reflectance spectroscopy-principal compo-nents regression analyses of soil properties[J].Soil Science Society of America Journal,65(2):480-490.doi:10.2136/sssaj2001.652480x.

    Du B,Lund P D,Wang J,Kolhe M,Hu E.2021.Comparative study of modelling the thermal efficiency of a novel straight through evacuated tube collector with MLR,SVR,BP and RBF methods[J].Sustainable Energy Technologies and Assessments,44:101029.doi:10.1016/j.seta.2021.10 1029.

    Fu YY,Yang G J,Song X Y,Li Z H,Xu X G,F(xiàn)eng H K,Zhao C J.2021.Improved estimation of winter wheat aboveg-round biomass using multiscale textures extracted from UAV-based digital images and hyperspectral feature analy-sis[J].Remote Sensing,13(4):581.doi:10.3390/rs 1304 0581.

    Guo Y H,Xiao Y,Li M W,Hao F H,Zhang X,Sun H Y,de Beurs K,F(xiàn)u Y S,He Y H.2022a.Identifying crop pheno-logy using maize height constructed from multi-sources images[J].International Journal of Applied Earth Observa-tion and Geoinformation,115:103121.doi:10.1016/j.jag.2022.103121.

    Guo Y,He J,Huang J Y,Jing Y H,Xu S B,Wang L G,Li S M,Zheng G Q.2022b.Effects of the spatial resolution of UAV images on the prediction and transferability of nitro-gen content model for winter wheat[J].Drones,6(10):299.doi:10.3390/drones6100299.

    Guo Y,Jing Y H,He J,Rong Y S,Wang L G,F(xiàn)eng W,Cheng Y Z,Liu T,Zheng G Q.2022c.Estimation of the height and above-ground biomass of winter wheat based on DOM and DSM of UAV images[J].International Conference on Optics and Machine Vision,12173:198-205.doi:10.1117/12.2634386.

    Imran A B,Khan K,Ali N,Ahmad N,Ali A,Shan K.2020.Narrow band based and broadband derived vegetation indi-ces using Sentinel-2 imagery to estimate vegetation bio-mass[J].Global Journal of Environmental Science andManagement,6(1):97-108.doi:10.22034/GJESM.2020.01.08.

    Lan W Q,Yang X,Gong T S,Xie J.2023.Predicting the shelf life of Trachinotus ovatus during frozen storage using a back propagation(BP)neural network model[J].Aqua-culture and Fisheries,8(5):544-550.doi:10.1016/j.aaf.2021.12.016.

    Luo S J,Jiang X Q,He Y B,Li J P,Jiao W H,Zhang S L,Xu F,Han Z C,Sun J,Yang J P,Wang X Y,Ma X T,Lin Z R.2022.Multi-dimensional variables and feature parameter selection for aboveground biomass estimation of potato based on UAV multispectral imagery[J].Frontiers in Plant Science,13:948249.doi:10.3389/fpls.2022.948249.

    Niu Y X,Zhang L Y,Zhang H H,Han W T,Peng X S.2019.Estimating above-ground biomass of maize using features derived from UAV-based RGB imagery[J].Remote Sen-sing,11(11):1261.doi:10.3390/rs 11111261.

    Piao S L,F(xiàn)ang J Y,Ciais P,Petlin P,Huang Y,Sitch S,Wang T.2009.The carbon balance of terrestrial ecosystems in China[J].Nature,458:1009-1013.doi:10.1038/nature 07944.

    Qiu R C,Zhang M,He Y.2022.Field estimation of maize plant height at jointing stage using an RGB-D camera[J].The Crop Journal,10(5):1274-1283.doi:10.1016/j.cj.2022.07.010.

    Steduto P,Albrizio R.2005.Resource use efficiency of field-grown sunflower,sorghum,wheat and chickpea:II.Water use efficiency and comparison with radiation use efficiency[J].Agricultural and Forest Meteorology,130(3-4):269-281.doi:10.1016/J.agrformet.2005.04.003.

    Wang F L,Yang M,Ma L F,Zhang T,Qin W L,Li W,Zhang Y H,Sun Z C,Wang Z M,Li F,Yu K.2022.Estimation of above-ground biomass of winter wheat based on consumer-grade multi-spectral UAV[J].Remote Sensing,14(5):1251.doi:10.3390/rs 14051251.

    Xie T J,Li J J,Yang C H,Jiang Z,Chen Y H,Guo L,Zhang J.2021.Crop height estimation based on UAV images:Me-thods,errors,and strategies[J].Computers and Electronics in Agriculture,185:106155.doi:10.1016/j.compag.2021.106155.

    Xu W K,Zhao L G,Li J,Shang X P,Ding X P,Wang T W.2022.Detection and classification of tea buds based on deep learning[J].Computers and Electronics in Agricul-ture,192:106547.doi:10.1016/j.compag.2021.106547.

    Zhu W B,Braun B,Chiang L H,Romagnoil J A.2021.Investi-gation of transfer learning for image classification and impact on training sample size[J].Chemometrics and In-telligent Laboratory Systems,211:104269.doi:10.1016/j.chemolab.2021.104269.

    (責(zé)任編輯:蘭宗寶)

    猜你喜歡
    遷移能力
    試析在高中生數(shù)學(xué)學(xué)習(xí)中,如何養(yǎng)成遷移能力
    初中化學(xué)遷移能力的培養(yǎng)
    高中數(shù)學(xué)的學(xué)習(xí)方法分析
    論小學(xué)語(yǔ)文閱讀教學(xué)探索路向與啟示
    淺談數(shù)學(xué)教學(xué)中對(duì)學(xué)生能力的培養(yǎng)
    提高高中生物理知識(shí)遷移能力的方法
    基于生活的科學(xué)探究,源于基礎(chǔ)的能力提升
    高中化學(xué)學(xué)生知識(shí)遷移學(xué)習(xí)力研究
    大學(xué)生就業(yè)路徑分析
    人民論壇(2016年5期)2016-03-24 22:13:11
    談?wù)劵瘜W(xué)教學(xué)中學(xué)生遷移能力的培養(yǎng)
    国产男靠女视频免费网站| 欧美精品啪啪一区二区三区| 黑人欧美特级aaaaaa片| 黄片播放在线免费| 日日摸夜夜添夜夜添小说| 国产免费视频播放在线视频| 国产成人av激情在线播放| 看免费av毛片| 国产1区2区3区精品| 午夜日韩欧美国产| 我要看黄色一级片免费的| 国产精品98久久久久久宅男小说| 国产精品久久久久久精品电影小说| 又黄又粗又硬又大视频| 黄色成人免费大全| 久久久精品区二区三区| 日本黄色视频三级网站网址 | 大片免费播放器 马上看| netflix在线观看网站| 老司机深夜福利视频在线观看| 在线 av 中文字幕| 少妇猛男粗大的猛烈进出视频| 国产日韩一区二区三区精品不卡| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦在线免费观看视频4| 怎么达到女性高潮| 日日夜夜操网爽| 成年女人毛片免费观看观看9 | 国产又爽黄色视频| 欧美变态另类bdsm刘玥| 高清视频免费观看一区二区| kizo精华| 欧美激情极品国产一区二区三区| 国产福利在线免费观看视频| 青草久久国产| 国产精品亚洲av一区麻豆| 久久人人爽av亚洲精品天堂| 精品少妇内射三级| 国产老妇伦熟女老妇高清| 久久久久久久大尺度免费视频| www.自偷自拍.com| 国产日韩欧美在线精品| 蜜桃在线观看..| 女性生殖器流出的白浆| 日本a在线网址| 久久天躁狠狠躁夜夜2o2o| 欧美乱妇无乱码| 性少妇av在线| 欧美一级毛片孕妇| 精品第一国产精品| 王馨瑶露胸无遮挡在线观看| 久久精品亚洲熟妇少妇任你| 国产精品一区二区精品视频观看| 国产无遮挡羞羞视频在线观看| 亚洲伊人色综图| 国产主播在线观看一区二区| 久久中文看片网| 亚洲中文av在线| 久久精品人人爽人人爽视色| 男女免费视频国产| 在线观看免费日韩欧美大片| 亚洲伊人久久精品综合| 91麻豆av在线| 欧美久久黑人一区二区| 王馨瑶露胸无遮挡在线观看| 水蜜桃什么品种好| 12—13女人毛片做爰片一| 精品人妻1区二区| 精品人妻在线不人妻| 日本精品一区二区三区蜜桃| 成人亚洲精品一区在线观看| 天天操日日干夜夜撸| 香蕉国产在线看| 国产精品亚洲av一区麻豆| 久久午夜综合久久蜜桃| 美女福利国产在线| 久久久久网色| 精品一区二区三区视频在线观看免费 | 久久 成人 亚洲| 免费日韩欧美在线观看| 亚洲精品国产区一区二| 在线观看人妻少妇| 一二三四社区在线视频社区8| 欧美日韩中文字幕国产精品一区二区三区 | 涩涩av久久男人的天堂| 欧美激情 高清一区二区三区| 国产精品二区激情视频| 18禁美女被吸乳视频| 久久九九热精品免费| 18禁裸乳无遮挡动漫免费视频| 最近最新中文字幕大全电影3 | 国产有黄有色有爽视频| 18在线观看网站| aaaaa片日本免费| 熟女少妇亚洲综合色aaa.| 日本av手机在线免费观看| 国产成+人综合+亚洲专区| 女同久久另类99精品国产91| 国产成人欧美| 免费看a级黄色片| 考比视频在线观看| 久久久国产精品麻豆| 国产精品久久久久久人妻精品电影 | √禁漫天堂资源中文www| 中文字幕av电影在线播放| 一级a爱视频在线免费观看| 国产精品久久久久久精品古装| 波多野结衣av一区二区av| 欧美老熟妇乱子伦牲交| 五月开心婷婷网| 99国产精品一区二区蜜桃av | 国产区一区二久久| 搡老乐熟女国产| 精品高清国产在线一区| 男女高潮啪啪啪动态图| 久久 成人 亚洲| 多毛熟女@视频| 男男h啪啪无遮挡| 国产精品免费大片| 亚洲欧美日韩高清在线视频 | 午夜福利在线免费观看网站| 一夜夜www| 青草久久国产| 人妻久久中文字幕网| 日日摸夜夜添夜夜添小说| 97人妻天天添夜夜摸| 午夜精品国产一区二区电影| 巨乳人妻的诱惑在线观看| 亚洲熟妇熟女久久| 久久亚洲真实| 久久精品国产99精品国产亚洲性色 | 蜜桃国产av成人99| 欧美精品高潮呻吟av久久| 深夜精品福利| 亚洲七黄色美女视频| 69精品国产乱码久久久| 欧美日韩亚洲国产一区二区在线观看 | 99国产精品免费福利视频| 大型av网站在线播放| 欧美日韩一级在线毛片| 亚洲avbb在线观看| 热re99久久国产66热| 在线观看66精品国产| 免费看a级黄色片| 99国产精品免费福利视频| 欧美精品一区二区大全| 色婷婷av一区二区三区视频| 免费av中文字幕在线| 欧美日韩国产mv在线观看视频| 曰老女人黄片| 午夜免费成人在线视频| 欧美激情极品国产一区二区三区| 黄片大片在线免费观看| 亚洲第一欧美日韩一区二区三区 | 亚洲精品在线观看二区| 欧美激情极品国产一区二区三区| a级毛片在线看网站| 成人黄色视频免费在线看| 正在播放国产对白刺激| 咕卡用的链子| 亚洲国产欧美日韩在线播放| 国产91精品成人一区二区三区 | 在线 av 中文字幕| videos熟女内射| 亚洲av日韩精品久久久久久密| 青草久久国产| 天堂8中文在线网| 激情在线观看视频在线高清 | 80岁老熟妇乱子伦牲交| 大型av网站在线播放| 一本—道久久a久久精品蜜桃钙片| 91精品国产国语对白视频| 一边摸一边做爽爽视频免费| 亚洲国产看品久久| 老司机影院毛片| 国产av又大| 最近最新中文字幕大全电影3 | 视频区欧美日本亚洲| 亚洲色图 男人天堂 中文字幕| 一级黄色大片毛片| av线在线观看网站| 男女床上黄色一级片免费看| 十八禁人妻一区二区| 国产精品免费大片| 欧美国产精品va在线观看不卡| 波多野结衣av一区二区av| 色94色欧美一区二区| 久久香蕉激情| 久久青草综合色| 热re99久久精品国产66热6| 亚洲,欧美精品.| 久久天躁狠狠躁夜夜2o2o| 国产精品欧美亚洲77777| 777久久人妻少妇嫩草av网站| 在线观看一区二区三区激情| 日韩视频一区二区在线观看| 免费在线观看黄色视频的| 高潮久久久久久久久久久不卡| 国产一区二区三区视频了| 国产伦理片在线播放av一区| 午夜福利在线观看吧| 在线播放国产精品三级| 欧美性长视频在线观看| 亚洲美女黄片视频| 黑丝袜美女国产一区| 午夜老司机福利片| 国产成人免费无遮挡视频| 成人av一区二区三区在线看| 亚洲成av片中文字幕在线观看| 十八禁人妻一区二区| 男女高潮啪啪啪动态图| 国产精品久久久久久精品电影小说| 波多野结衣av一区二区av| 国产一区二区三区在线臀色熟女 | 亚洲少妇的诱惑av| 国产成人精品在线电影| 国产欧美日韩综合在线一区二区| 99精品在免费线老司机午夜| 悠悠久久av| 美女高潮喷水抽搐中文字幕| 亚洲人成电影观看| 中文字幕最新亚洲高清| 欧美成人免费av一区二区三区 | 桃花免费在线播放| 色播在线永久视频| 老熟妇仑乱视频hdxx| 岛国毛片在线播放| 精品熟女少妇八av免费久了| 国产精品一区二区免费欧美| 国产精品久久久久久精品电影小说| 91成人精品电影| 中文字幕人妻丝袜制服| 日本撒尿小便嘘嘘汇集6| 国产成人精品久久二区二区免费| 在线亚洲精品国产二区图片欧美| 19禁男女啪啪无遮挡网站| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 亚洲少妇的诱惑av| 日本黄色日本黄色录像| 91国产中文字幕| 两个人看的免费小视频| 久久久精品区二区三区| 人人妻人人爽人人添夜夜欢视频| 一区二区三区精品91| 国产aⅴ精品一区二区三区波| 欧美日韩av久久| avwww免费| 久久久久久久久久久久大奶| 欧美成人免费av一区二区三区 | 91麻豆av在线| 日韩免费高清中文字幕av| 久久精品国产a三级三级三级| 露出奶头的视频| 午夜老司机福利片| 王馨瑶露胸无遮挡在线观看| 新久久久久国产一级毛片| 一个人免费在线观看的高清视频| 18禁黄网站禁片午夜丰满| 精品少妇内射三级| 三上悠亚av全集在线观看| 日韩熟女老妇一区二区性免费视频| 在线观看免费午夜福利视频| 一级毛片精品| 久久久久久久大尺度免费视频| 国产精品九九99| 亚洲人成电影免费在线| 精品高清国产在线一区| 老熟女久久久| 久久人妻福利社区极品人妻图片| 色婷婷av一区二区三区视频| 9191精品国产免费久久| 丝瓜视频免费看黄片| 成人免费观看视频高清| 欧美日韩一级在线毛片| 我的亚洲天堂| 老司机靠b影院| 日本精品一区二区三区蜜桃| 国产色视频综合| 我的亚洲天堂| 日韩一区二区三区影片| 亚洲第一青青草原| 嫁个100分男人电影在线观看| 国产男女内射视频| 成年人午夜在线观看视频| 日本欧美视频一区| 久久人妻福利社区极品人妻图片| 一区二区av电影网| 日本av手机在线免费观看| 久久国产精品影院| 国产av一区二区精品久久| 欧美精品人与动牲交sv欧美| 人妻久久中文字幕网| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久精品国产亚洲av高清涩受| 国产免费视频播放在线视频| 午夜精品国产一区二区电影| 视频区欧美日本亚洲| 麻豆国产av国片精品| 国产精品电影一区二区三区 | 国产精品99久久99久久久不卡| 成年动漫av网址| 中亚洲国语对白在线视频| 国产精品免费大片| 国产精品98久久久久久宅男小说| av视频免费观看在线观看| 精品国产一区二区三区久久久樱花| 国产97色在线日韩免费| 香蕉久久夜色| 一个人免费看片子| 纵有疾风起免费观看全集完整版| 久久国产精品人妻蜜桃| 精品一区二区三区视频在线观看免费 | 亚洲午夜理论影院| 亚洲成人免费av在线播放| 狠狠狠狠99中文字幕| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| av视频免费观看在线观看| 最新的欧美精品一区二区| 国产高清videossex| 妹子高潮喷水视频| 欧美激情极品国产一区二区三区| 黄色片一级片一级黄色片| 一本综合久久免费| 亚洲美女黄片视频| 69av精品久久久久久 | videos熟女内射| 久久久久久久精品吃奶| 久久人妻av系列| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 亚洲国产欧美在线一区| 在线观看66精品国产| 亚洲精品一二三| 国产熟女午夜一区二区三区| 精品国产乱子伦一区二区三区| videos熟女内射| 久久天躁狠狠躁夜夜2o2o| 美女福利国产在线| 天天躁狠狠躁夜夜躁狠狠躁| av超薄肉色丝袜交足视频| 久久青草综合色| 国产一区有黄有色的免费视频| 少妇粗大呻吟视频| 狠狠婷婷综合久久久久久88av| 人妻久久中文字幕网| 变态另类成人亚洲欧美熟女 | 18在线观看网站| 9191精品国产免费久久| 国产主播在线观看一区二区| 一区在线观看完整版| 美国免费a级毛片| 国产精品免费一区二区三区在线 | 欧美精品一区二区大全| 一级毛片女人18水好多| 涩涩av久久男人的天堂| 国产免费现黄频在线看| 一本—道久久a久久精品蜜桃钙片| 久久国产精品影院| 久久精品人人爽人人爽视色| 一级片'在线观看视频| 中文字幕制服av| 在线看a的网站| 99精品在免费线老司机午夜| 国产又爽黄色视频| 香蕉国产在线看| 免费在线观看视频国产中文字幕亚洲| 日韩欧美三级三区| 成人国产av品久久久| 国产一卡二卡三卡精品| 国产成人一区二区三区免费视频网站| 9热在线视频观看99| 在线天堂中文资源库| 亚洲avbb在线观看| 每晚都被弄得嗷嗷叫到高潮| 欧美精品人与动牲交sv欧美| 精品国内亚洲2022精品成人 | 丰满饥渴人妻一区二区三| 免费女性裸体啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 美女视频免费永久观看网站| 国产av精品麻豆| 精品第一国产精品| 欧美在线黄色| 操出白浆在线播放| 少妇 在线观看| 少妇裸体淫交视频免费看高清 | 亚洲国产毛片av蜜桃av| 汤姆久久久久久久影院中文字幕| 欧美乱码精品一区二区三区| 国产免费av片在线观看野外av| 99久久国产精品久久久| 黄色视频在线播放观看不卡| videosex国产| 国产人伦9x9x在线观看| 嫩草影视91久久| 国产成人影院久久av| 一区二区三区乱码不卡18| 欧美黑人欧美精品刺激| 桃花免费在线播放| 妹子高潮喷水视频| 狠狠狠狠99中文字幕| 男男h啪啪无遮挡| 国产精品免费一区二区三区在线 | 男女边摸边吃奶| 女人爽到高潮嗷嗷叫在线视频| 人人妻人人爽人人添夜夜欢视频| 久久久久国产一级毛片高清牌| 悠悠久久av| 久久亚洲真实| 在线观看舔阴道视频| 高清视频免费观看一区二区| 中文字幕色久视频| 精品一区二区三区视频在线观看免费 | 免费观看av网站的网址| 欧美精品av麻豆av| 日韩欧美国产一区二区入口| 免费人妻精品一区二区三区视频| 欧美人与性动交α欧美软件| 捣出白浆h1v1| 午夜成年电影在线免费观看| 老熟女久久久| 成人手机av| 免费在线观看视频国产中文字幕亚洲| 日本黄色日本黄色录像| 高清在线国产一区| 在线观看人妻少妇| 午夜精品久久久久久毛片777| 黑人巨大精品欧美一区二区mp4| 51午夜福利影视在线观看| 啦啦啦在线免费观看视频4| 午夜福利欧美成人| 精品一区二区三区四区五区乱码| 另类精品久久| 日韩大码丰满熟妇| 黄色毛片三级朝国网站| 日本五十路高清| 国产男女超爽视频在线观看| 亚洲人成电影免费在线| 亚洲成国产人片在线观看| svipshipincom国产片| 老熟妇乱子伦视频在线观看| 亚洲av电影在线进入| 欧美大码av| 日韩精品免费视频一区二区三区| 老司机午夜十八禁免费视频| 在线观看www视频免费| 麻豆av在线久日| 五月天丁香电影| 久久人妻福利社区极品人妻图片| 精品国内亚洲2022精品成人 | 12—13女人毛片做爰片一| 丝袜人妻中文字幕| 亚洲午夜精品一区,二区,三区| 精品久久久久久久毛片微露脸| 日日爽夜夜爽网站| a级片在线免费高清观看视频| 大片电影免费在线观看免费| 人人妻,人人澡人人爽秒播| 午夜老司机福利片| 老司机靠b影院| 嫁个100分男人电影在线观看| 国产精品98久久久久久宅男小说| 精品免费久久久久久久清纯 | 亚洲精品一卡2卡三卡4卡5卡| 人人妻人人澡人人爽人人夜夜| 黄色毛片三级朝国网站| 美女高潮喷水抽搐中文字幕| 欧美国产精品一级二级三级| 伊人久久大香线蕉亚洲五| 亚洲天堂av无毛| 亚洲欧美日韩另类电影网站| 精品免费久久久久久久清纯 | 在线十欧美十亚洲十日本专区| 国产在线免费精品| 黄色怎么调成土黄色| 亚洲av美国av| 国产片内射在线| 国产免费av片在线观看野外av| 国产精品 欧美亚洲| 免费少妇av软件| 国产97色在线日韩免费| 丁香六月天网| 操美女的视频在线观看| 纯流量卡能插随身wifi吗| 国产色视频综合| 亚洲三区欧美一区| 十八禁人妻一区二区| 麻豆成人av在线观看| 一夜夜www| 91精品国产国语对白视频| 日本精品一区二区三区蜜桃| 高潮久久久久久久久久久不卡| www.999成人在线观看| svipshipincom国产片| av在线播放免费不卡| 大片电影免费在线观看免费| av欧美777| 熟女少妇亚洲综合色aaa.| 国产又色又爽无遮挡免费看| 美女扒开内裤让男人捅视频| 久久精品熟女亚洲av麻豆精品| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| 欧美另类亚洲清纯唯美| 啦啦啦免费观看视频1| 精品欧美一区二区三区在线| 久久国产精品男人的天堂亚洲| 首页视频小说图片口味搜索| 久久久国产欧美日韩av| 免费不卡黄色视频| 女人高潮潮喷娇喘18禁视频| 热99re8久久精品国产| 久久久久久久国产电影| 国产熟女午夜一区二区三区| 成人av一区二区三区在线看| 69精品国产乱码久久久| 日本黄色视频三级网站网址 | 99精品久久久久人妻精品| 亚洲国产欧美日韩在线播放| 久久香蕉激情| 久久久久久久国产电影| 无人区码免费观看不卡 | 久久这里只有精品19| 久久人人97超碰香蕉20202| 18禁黄网站禁片午夜丰满| 在线亚洲精品国产二区图片欧美| 18禁黄网站禁片午夜丰满| 老司机靠b影院| 国产亚洲精品久久久久5区| 欧美久久黑人一区二区| 国产又色又爽无遮挡免费看| 自拍欧美九色日韩亚洲蝌蚪91| av在线播放免费不卡| 欧美国产精品一级二级三级| 欧美老熟妇乱子伦牲交| 久久久久久久久免费视频了| 十八禁网站网址无遮挡| 中文字幕av电影在线播放| 午夜福利免费观看在线| 大片免费播放器 马上看| 肉色欧美久久久久久久蜜桃| 热99国产精品久久久久久7| 人成视频在线观看免费观看| 国产精品1区2区在线观看. | 国产一区二区三区综合在线观看| 免费在线观看日本一区| 欧美日本中文国产一区发布| 久久国产亚洲av麻豆专区| 一个人免费在线观看的高清视频| 大片免费播放器 马上看| 极品教师在线免费播放| 一本大道久久a久久精品| 51午夜福利影视在线观看| 曰老女人黄片| 欧美亚洲 丝袜 人妻 在线| 国产一区二区三区视频了| 99精品久久久久人妻精品| 久久人人97超碰香蕉20202| 99久久精品国产亚洲精品| 无遮挡黄片免费观看| 精品午夜福利视频在线观看一区 | 超色免费av| 国产高清视频在线播放一区| 黄色片一级片一级黄色片| 国产精品影院久久| 欧美成人免费av一区二区三区 | 免费在线观看影片大全网站| 视频区图区小说| 在线观看人妻少妇| 亚洲国产欧美在线一区| av欧美777| av天堂久久9| 黄色视频不卡| 亚洲av成人一区二区三| 久久人妻福利社区极品人妻图片| 亚洲一区中文字幕在线| 久久精品91无色码中文字幕| 亚洲中文字幕日韩| 欧美人与性动交α欧美精品济南到| 亚洲专区中文字幕在线| 成人国语在线视频| 两个人免费观看高清视频| 亚洲情色 制服丝袜| 亚洲中文av在线| 又黄又粗又硬又大视频| 最近最新免费中文字幕在线| 999精品在线视频| 成人国语在线视频| 在线永久观看黄色视频| 啦啦啦中文免费视频观看日本| 精品国产一区二区三区久久久樱花| 岛国在线观看网站| 国产黄频视频在线观看| 大型黄色视频在线免费观看| 国产区一区二久久| 热99久久久久精品小说推荐| 一级毛片女人18水好多| 国产欧美亚洲国产| 欧美人与性动交α欧美软件| 日韩人妻精品一区2区三区| 国产精品久久久久久精品古装| 精品久久久久久电影网| 久久国产精品男人的天堂亚洲| 成人影院久久| 亚洲精品粉嫩美女一区| 99九九在线精品视频| 母亲3免费完整高清在线观看| 动漫黄色视频在线观看| 极品人妻少妇av视频| 精品视频人人做人人爽| 久久久久久久久免费视频了| 亚洲av日韩精品久久久久久密| 热99国产精品久久久久久7|