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

    融合光譜和紋理特征的玉米產(chǎn)量預(yù)測研究

    2024-11-02 00:00:00馬元花汪樂印張禎鑫鄭大圣葉玉瀾崔志峰杜冰笑寸玉潔李軍王瑞

    摘 要 本研究采集夏玉米拔節(jié)、抽雄、灌漿和成熟4個(gè)時(shí)期的無人機(jī)可見光和多光譜影像,提取和篩選植被指數(shù)與紋理特征參數(shù),構(gòu)建植被指數(shù)與紋理特征融合變量,采用反向傳播神經(jīng)網(wǎng)絡(luò)、隨機(jī)森林和支持向量機(jī)3種機(jī)器學(xué)習(xí)方法構(gòu)建玉米產(chǎn)量預(yù)測模型。結(jié)果表明:相較于單一類型參數(shù),融合植被指數(shù)和紋理特征進(jìn)行產(chǎn)量預(yù)測模型精度更高;3種機(jī)器學(xué)習(xí)方法中以隨機(jī)森林構(gòu)建的玉米產(chǎn)量預(yù)測模型效果最好,且最佳預(yù)測時(shí)期為灌漿期(籽粒水泡期);綜合評價(jià)建模和驗(yàn)證結(jié)果,基于多光譜影像植被指數(shù)與紋理特征融合變量和隨機(jī)森林方法構(gòu)建的模型玉米產(chǎn)量預(yù)測效果最佳。遙感信息的多特征融合與機(jī)器學(xué)習(xí)方法的搭配能夠挖掘和利用更多信息并提高玉米產(chǎn)量預(yù)測的精度和魯棒性。

    關(guān)鍵詞 玉米;無人機(jī)遙感;可見光;多光譜;植被指數(shù);紋理特征;產(chǎn)量

    玉米作為三大谷物之首,是當(dāng)今全球種植范圍廣泛、產(chǎn)量最大的谷類作物,同時(shí)也是重要的糧食、經(jīng)濟(jì)和飼料兼用作物。近年來,中國玉米的播種面積和產(chǎn)量總體呈穩(wěn)定增長的態(tài)勢,及時(shí)準(zhǔn)確地獲得玉米產(chǎn)量信息,對于保障國家糧食安全以及農(nóng)業(yè)結(jié)構(gòu)的優(yōu)化有著重要意義[1]。由于受地理環(huán)境、自然氣候、品種和生產(chǎn)技術(shù)差異等因素影響,作物產(chǎn)量高效、準(zhǔn)確估測是一項(xiàng)具有挑戰(zhàn)性的任務(wù)[2-3]。

    傳統(tǒng)作物產(chǎn)量的估算方法在很大程度上依賴于實(shí)地調(diào)查,成本高昂、費(fèi)力耗時(shí)[4]。遙感監(jiān)測技術(shù)具有實(shí)時(shí)、無損和高通量的特性,可用于植物的長勢[5-6]、脅迫狀態(tài)[7-8]和產(chǎn)量[9-10]等監(jiān)測,已成為現(xiàn)代作物產(chǎn)量調(diào)查的主要方法之一[11]。得益于無人機(jī)平臺和機(jī)載傳感器領(lǐng)域的快速發(fā)展,無人機(jī)遙感可獲得較高的空間、光譜和時(shí)間分辨率數(shù)據(jù)[12],在農(nóng)田精細(xì)尺度和動(dòng)態(tài)連續(xù)監(jiān)測方面較衛(wèi)星遙感具有顯著優(yōu)勢[13-16]。無人機(jī)作物遙感監(jiān)測的主要方法是通過搭載可見光、多光譜和高光譜等傳感器,獲取作物光譜信息,提取植被指數(shù)等信息用于反演作物長勢。植被指數(shù)是從植被冠層的光譜信息中提取的一種獨(dú)特的光學(xué)參數(shù)[17]。利用植被指數(shù)進(jìn)行產(chǎn)量估算已在糧食作物[18-19]、經(jīng)濟(jì)作物[20]、飼用作物[21]等方面取得諸多成果,但預(yù)測效率與精度仍有較大提升空間。除植被指數(shù)等光譜信息外,用于表達(dá)空間細(xì)節(jié)差異、突出植物冠層結(jié)構(gòu)和幾何特征的紋理特征在產(chǎn)量預(yù)測方面也有較好的應(yīng)用潛力。前人研究發(fā)現(xiàn)相較于僅用植被指數(shù)進(jìn)行棉花葉面積指數(shù)監(jiān)測,通過增加遙感影像的紋理特征信息,可提高監(jiān)測精度[22]。

    遙感估產(chǎn)模型大致分為兩類:線性模型和非線性模型。使用基于植被指數(shù)、紋理指數(shù)等估測量和產(chǎn)量觀測值的傳統(tǒng)線性回歸模型,精度往往不足,且易出現(xiàn)過度擬合現(xiàn)象[23-24]。機(jī)器學(xué)習(xí)算法的發(fā)展為農(nóng)藝信息與遙感數(shù)據(jù)相結(jié)合并開發(fā)靈活的產(chǎn)量預(yù)測模型提供了新的機(jī)會。機(jī)器學(xué)習(xí)可高效編譯各種數(shù)據(jù)集和提取有用信息,提高模型準(zhǔn)確性,減少過擬合現(xiàn)象,已成功應(yīng)用于解決大型非線性問題[25],神經(jīng)網(wǎng)絡(luò)、隨機(jī)森林和支持向量機(jī)等算法日益廣泛用于產(chǎn)量預(yù)測問題的研究[26-28]。

    玉米產(chǎn)量的形成是一個(gè)從播種到籽粒成熟的漫長過程,產(chǎn)量波動(dòng)較大,不同時(shí)期植被指數(shù)和特征因子與產(chǎn)量的相關(guān)程度不同,遙感產(chǎn)量估算方案需要基于玉米生長規(guī)律篩選最佳生育時(shí)期、估產(chǎn)參數(shù)和建模方法。因此,本研究獲取夏玉米的拔節(jié)、抽雄、灌漿和成熟4個(gè)關(guān)鍵生育時(shí)期的無人機(jī)可見光和多光譜影像數(shù)據(jù),篩選與產(chǎn)量相關(guān)的植被指數(shù)和紋理特征參數(shù),基于多源遙感數(shù)據(jù)融合分別建立反向傳播神經(jīng)網(wǎng)絡(luò)、隨機(jī)森林和支持向量機(jī)的夏玉米產(chǎn)量估測模型并進(jìn)行精度對比,篩選最佳估測產(chǎn)量的時(shí)期、建模特征參數(shù)和方法,為夏玉米無人機(jī)遙感產(chǎn)量預(yù)測提供優(yōu)化方案。

    1 材料與方法

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

    試驗(yàn)于2022年在陜西省楊凌區(qū)西北農(nóng)林科技大學(xué)曹新莊農(nóng)場進(jìn)行(34°29′N,108°06′E,海拔454.8 m),該地位于陜西關(guān)中平原,海拔454.8 m,屬暖溫帶半濕潤季風(fēng)氣候,夏玉米生長季平均降水為304.9 mm,平均氣溫為26.0 ℃,總?cè)照諘r(shí)間為1 213 h,供試土壤為[HT5,6SS]土[KG-*3][HT5,6SS]婁土,本研究基于2009年起實(shí)施的冬小麥-夏玉米一年兩熟制長期定位氮磷施肥試驗(yàn)實(shí)施。

    供試玉米品種分別為‘鄭單958’‘先玉335’‘京科968’和‘陜單609’,田間試驗(yàn)采取品種、施氮量和施磷量三因素裂區(qū)設(shè)計(jì),夏玉米于2022年6月15日播種,9月30日收獲。5個(gè)氮肥(N)梯度:0、75、150、225和300 kg/hm2,4個(gè)磷肥(P2O5)梯度:0、60、120和180 kg/hm2,即20個(gè)施肥處理。設(shè)置3次重復(fù)。肥料全部基施。種植密度為67 500株/hm2,株距24.7 cm,行距 60 cm。其他田間管理措施同當(dāng)?shù)卮筇铩?/p>

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

    1.2.1 無人機(jī)影像獲取與處理 本研究使用大疆經(jīng)緯M300 RTK無人機(jī)搭載MicaSense Altum-PT多光譜鏡頭采集夏玉米拔節(jié)期(V6,7月24日)、抽雄期(VT,8月10日)、灌漿期(籽粒水泡期R2,9月1日)和成熟期(R6,9月28日)4個(gè)時(shí)期的遙感數(shù)據(jù),該相機(jī)可同時(shí)采集6通道(全色、藍(lán)光、綠光、紅光、紅邊、紅外)多光譜影像。利用DJI Pilot APP規(guī)劃航線,制定飛行任務(wù)。每次拍攝在晴朗無云天氣10:00—12:00進(jìn)行,無人機(jī)飛行高度設(shè)定為20 m,航向重疊率為80%,旁向重疊率為75%,每個(gè)生育時(shí)期的飛行計(jì)劃和相機(jī)設(shè)置保持一致。飛行前拍攝校準(zhǔn)反射面板,用于無人機(jī)影像的輻射校正。

    使用Pix4D mapper軟件對航拍獲得的影像進(jìn)行輻射校正和圖像拼接。采用ArcGIS軟件對合成后的可見光(RGB)影像和多光譜(Multi-spectral,MS)影像進(jìn)行試驗(yàn)區(qū)的選擇、地圖代數(shù)運(yùn)算和分區(qū)統(tǒng)計(jì)處理。

    1.2.2 產(chǎn)量數(shù)據(jù)獲取與處理 在玉米成熟期,每個(gè)小區(qū)選取行長5 m長勢均勻的玉米各3行,統(tǒng)計(jì)穗數(shù),在行內(nèi)選取20個(gè)均勻果穗,取3次重復(fù),風(fēng)干后于室內(nèi)考種,測定穗粒數(shù)和百粒質(zhì)量,按14%含水量折算籽粒產(chǎn)量[29]。

    1.3 植被指數(shù)和紋理特征的選取

    植被指數(shù)(Vegetation index,VI)是兩個(gè)或多個(gè)波段地物反射率的組合運(yùn)算,以增強(qiáng)植被某一特性或者細(xì)節(jié),是作物生長分析的重要參數(shù)。本文基于前人研究選取可見光波段構(gòu)建的歸一化紅光(r)、歸一化綠光(g)、歸一化藍(lán)光(b)、綠藍(lán)比指數(shù)(GBRI)、紅綠比指數(shù)(RGRI)、超藍(lán)指數(shù)(ExB)、超綠指數(shù)(ExG)、超紅指數(shù)(ExR)、超綠超紅差分指數(shù)(ExGR)、植被顏色提取指數(shù)(CIVE)、植被指數(shù)(VEG)、改進(jìn)型綠紅植被指數(shù)(MGRVI)、歸一化綠紅差分指數(shù)(NGRDI)、可見光波段差異植被指數(shù)(VDVI)、綜合指數(shù)1(COM1)和綜合指數(shù)2(COM2)16種可見光植被指數(shù)[30-32];基于多光譜波段構(gòu)建的耐大氣植被指數(shù)(ARVI)、差值植被指數(shù)(DVI)、增強(qiáng)植被指數(shù)(EVI)、綠色歸一化差異植被指數(shù)(GNDVI)、葉面葉綠素指數(shù)(LCI)、歸一化差異紅邊植被指數(shù)(NDRE)、歸一化差異植被指數(shù)(NDVI)、優(yōu)化的土壤調(diào)節(jié)植被指數(shù)(OSAVI)、比值植被指數(shù)(RVI)和土壤調(diào)節(jié)植被指數(shù)(SAVI)10種多光譜植被指數(shù)[33-34]進(jìn)行估測建模,共計(jì)26種植被指數(shù)。[JP]

    紋理特征(Textural features,T)是一種全局特征,反映圖像中同質(zhì)現(xiàn)象的視覺特征,體現(xiàn)物體表面具有緩慢變換或周期性變化的表面組織結(jié)構(gòu)排列屬性,對作物產(chǎn)量估測有很高的應(yīng)用潛力。本研究利用ENVI軟件中的灰度共生矩陣(Gray-level co-occurrence matrix,GLCM)提取出每個(gè)植被指數(shù)的8種紋理特征[35]:對比度(Contrast,CON)、相關(guān)性(Correlation,COR)、相異性(Dissimilarity,DIS)、熵(Entropy,ENT)、協(xié)同性(Homogeneity,HOM)、均值(Mean,MEA)、二階矩(Second moment,SEM)和方差(Variance,VAR)。分析窗口尺寸為3×3,空間相關(guān)矩陣偏移X和Y默認(rèn)為1。

    使用皮爾遜相關(guān)系數(shù)法(Pearson Correlation Coefficient,PCC)分析產(chǎn)量與植被指數(shù)和紋理特征的相關(guān)系數(shù)(r),并以其絕對值(|r|)的大小為標(biāo)準(zhǔn)對構(gòu)建的光譜指數(shù)進(jìn)行變量篩選并排序[36]。

    1.4 模型構(gòu)建

    本研究將各個(gè)時(shí)期植被指數(shù)和紋理特征作為輸入?yún)?shù),分別與產(chǎn)量數(shù)據(jù)進(jìn)行相關(guān)性分析和顯著性檢驗(yàn),篩選出最優(yōu)參數(shù)后,選取反向傳播神經(jīng)網(wǎng)絡(luò)(Back Propagation Neural Network,BPNN)、隨機(jī)森林(Random Forest,RF)和支持向量機(jī)(Support Vector Machine,SVM)3種機(jī)器模型進(jìn)行產(chǎn)量估測。本試驗(yàn)共采集產(chǎn)量數(shù)據(jù)240個(gè),按3∶1比例劃分訓(xùn)練集和驗(yàn)證集,即隨機(jī)選取其中75%作為建模集,25%作為訓(xùn)練集樣本。取各小區(qū)植被指數(shù)和紋理特征平均值作為建模數(shù)據(jù)。使用Matlab軟件進(jìn)行數(shù)據(jù)分析。

    1.5 模型精度評價(jià)指標(biāo)

    選取決定系數(shù)(Coefficient of determination,R2)、均方根誤差(Root mean square error,RMSE)和平均絕對誤差(Mean absolute error,MAE)評價(jià)模型精度。R2用于評價(jià)估測值和實(shí)測值的擬合程度,其值越接近1則模型擬合效果越好;RMSE用于評價(jià)估測值和實(shí)測值的偏差程度,其值越小則模型擬合效果越好;MAE用于評價(jià)估測值和實(shí)測值的實(shí)際偏差,其值越小則模型擬合效果越好[5]。

    2 結(jié)果與分析

    2.1 可見光與多光譜植被指數(shù)及其紋理特征與玉米籽粒產(chǎn)量的相關(guān)性

    將各生育時(shí)期的植被指數(shù)和紋理特征分別與夏玉米產(chǎn)量進(jìn)行相關(guān)性分析。根據(jù)皮爾遜相關(guān)系數(shù)界值表,n=80時(shí),|r|>0.217表示P<0.05水平顯著,|r|>0.283表示P<0.01水平顯著,|r|>0.357表示P<0.001水平顯著。對植被指數(shù)篩選發(fā)現(xiàn)(表1、表2),玉米V6、VT、R2和R6的植被指數(shù)與產(chǎn)量相關(guān)系數(shù)絕對值|r|> 0.217(P<0.05)的可見光和多光譜植被指數(shù)分別為25、24、22和20個(gè),且多光譜植被指數(shù)與產(chǎn)量的相關(guān)性大于可見光植被指數(shù)。4個(gè)生育時(shí)期與產(chǎn)量相關(guān)性最高的植被指數(shù)分別為NDRE、NDRE、RVI和GNDVI。以R2時(shí)期植被指數(shù)與產(chǎn)量相關(guān)性最高。

    為便于表示,用“_”連接植被指數(shù)與紋理特征,如r_CON表示從r植被指數(shù)提取的對比度(Contrast,CON)。對可見光和多光譜植被指數(shù)的紋理特征篩選發(fā)現(xiàn)(圖1),可見光紋理特征的在4個(gè)生育時(shí)期分別以r_MEA、r_MEA、RGRI_MEA和RGRI_COR與產(chǎn)量的相關(guān)性最好,|r|介于0.657~0.763;多光譜紋理特征的|r|在4個(gè)生育時(shí)期分別以NDRE_MEA、ARVI_COR、RVI_MEA和NDRE_VAR相關(guān)性最好,|r|介于0.779~0.843,均表現(xiàn)出極顯著相關(guān)關(guān)系(P< 0.001)。通過統(tǒng)計(jì)篩選4個(gè)生育時(shí)期最優(yōu)可見光與多光譜在的紋理特征類型,以MAE特征頻率最高,占比為52.5%,其次是COR特征,占比為15%。

    綜上,基于多光譜波段構(gòu)建的植被指數(shù)和紋理特征與產(chǎn)量的相關(guān)性要優(yōu)于可見光,紋理特征與產(chǎn)量的相關(guān)性略優(yōu)于植被指數(shù)。4個(gè)生育時(shí)期中以R2時(shí)期提取的有關(guān)參量與產(chǎn)量的相關(guān)性表現(xiàn)最好。

    2.2 建模參數(shù)的優(yōu)選

    依據(jù)皮爾遜相關(guān)系數(shù)法對特征進(jìn)行排序,分別針對RGB和MS數(shù)據(jù)篩選各生育時(shí)期相關(guān)系數(shù)絕對值最高的5個(gè)植被指數(shù)(VI)和紋理特征(T)作為建模的特征參數(shù),其中VI+T為相對應(yīng)參數(shù)的疊加,形成RGB_VI、RGB_T、RGB_VI+T、MS_VI、MS_T、MS_VI+T 6種參數(shù)組合類型(表3)。根據(jù)表1、表2和圖1結(jié)果,篩選的特征參數(shù)與在V6、VT、R2、R6 4個(gè)生育時(shí)期|r|范圍表現(xiàn)為RGB_VI:0.627~ 0.659、0.629~0.708、0.616~0.757、0.398~ 0.563;MS_VI:0.722~0.782、0.714~0.821、 0.855~0.878、0.675~ 0.757,RGB_T:0.627~0.657、0.679~0.755、 0.662~0.763和0.503~0.658,MS_T:0.733~0.783、0.798~0.843、 0.862~0.880和0.722~0.779,均表現(xiàn)出極強(qiáng)相關(guān)性(P<0.001)。

    2.3 基于機(jī)器學(xué)習(xí)的玉米產(chǎn)量預(yù)測模型的構(gòu)建與驗(yàn)證

    分別使用BPNN、RF和SVM 3種機(jī)器學(xué)習(xí)方法進(jìn)行建模和驗(yàn)證(圖2、圖3和圖4),基于不同生育時(shí)期估測產(chǎn)量,均以籽粒灌漿期(R2)模型的精度最優(yōu);在參數(shù)類型選擇方面,均以多光譜植被指數(shù)與紋理特征疊加(MS_VI+T)類型參數(shù)建模與驗(yàn)證精度最高,多光譜紋理特征(MS_T)類型參數(shù)次之。結(jié)合建模和驗(yàn)證效果,3種機(jī)器學(xué)習(xí)建模方法中以RF的總體建模和驗(yàn)證精度 最優(yōu)。

    基于BPNN模型中,最優(yōu)模型為MS_VI+T(建模集R2=0.763,RMSE=466.890 kg/hm2,MAE=360.531 kg/hm2;驗(yàn)證集R2=0.864,RMSE=446.387 kg/hm2,MAE=353.725 kg/hm2)。當(dāng)僅利用單一類型參數(shù)建模時(shí),以MS_T模型結(jié)果最優(yōu)(模型R2=0.822,RMSE=482.861 kg/hm2,MAE=386.144 kg/hm2)。相比于MS_VI模型,MS_T模型驗(yàn)證效果R2提升0.002,RMSE降低99.786 kg/hm2,MAE降低78.448 kg/hm2。表明,基于單一紋理特征參數(shù)建模效果優(yōu)于單一植被指數(shù)模型。

    基于RF模型中,最優(yōu)模型為MS_VI+T(建模集R2=0.859,RMSE=400.288 kg/hm2,MAE=308.907 kg/hm2;驗(yàn)證集R2=0.900,RMSE=295.291 kg/hm2,MAE=249.082 kg/hm2)。其中,相比于RGB_VI+T模型,MS_VI+T模型驗(yàn)證效果R2提升0.224,RMSE降低344.799 kg/hm2,MAE降低268.906 kg/hm2。由此,基于多光譜參數(shù)類型構(gòu)建的模型優(yōu)于可見光參數(shù)類型的模型。

    基于SVM模型中,最優(yōu)模型為MS_VI+T(建模集R2=0.775,RMSE=506.824 kg/hm2,MAE=346.592 kg/hm2;驗(yàn)證集R2=0.904,RMSE=291.351 kg/hm2,MAE=217.139 kg/hm2)。相比于MS_VI和MS_T模型,MS_VI+T模型驗(yàn)證效果R2提升0.017和0.038,RMSE降低142.053和88.274 kg/hm2,MAE降低112.016和79.655 kg/hm2。因此,融合植被指數(shù)與紋理特征參數(shù)建模能夠提升模型的精度,其中對于植被指數(shù)提升空間高于紋理特征參數(shù) 模型。

    綜上,融合植被指數(shù)和紋理特征構(gòu)建的模型較單一類型變量模型的建模與驗(yàn)證精度均有所提升,基于紋理特征構(gòu)建的產(chǎn)量預(yù)測模型的精度總體優(yōu)于植被指數(shù)模型,基于隨機(jī)森林方法構(gòu)建模型的精度總體優(yōu)于反向神經(jīng)網(wǎng)絡(luò)和支持向量機(jī)方法。綜合建模和驗(yàn)證精度,本研究中最優(yōu)玉米產(chǎn)量預(yù)測模型為融合多光譜植被指數(shù)和紋理特征的基于隨機(jī)森林方法在灌漿期(R2)構(gòu)建的模型。

    2.4 玉米產(chǎn)量分布反演

    基于篩選的最優(yōu)玉米產(chǎn)量預(yù)測進(jìn)行產(chǎn)量反演,得到玉米產(chǎn)量示意圖(圖5)。不同的氮磷水平處理下夏玉米產(chǎn)量存在差異,產(chǎn)量隨施氮量的增加呈現(xiàn)先增加后降低的趨勢,即N3施肥處理下產(chǎn)量達(dá)到最高,且產(chǎn)量隨施磷量的增加而增加,在所有處理中產(chǎn)量介于3 709~7 322 kg/hm2,各小區(qū)不同品種玉米產(chǎn)量的實(shí)際產(chǎn)量和預(yù)測產(chǎn)量較為接近。

    3 討 論

    3.1 玉米產(chǎn)量預(yù)測生育時(shí)期選擇

    遙感信息會隨著作物的生長發(fā)生變化,本文基于玉米拔節(jié)期、抽雄期、灌漿期和成熟期對產(chǎn)量進(jìn)行遙感估測的精度存在較大差異。玉米苗期(出苗至拔節(jié))根系發(fā)展迅速,莖葉生長緩慢,前期地上部長勢及遙感信息與產(chǎn)量的相關(guān)性較小。穗期(拔節(jié)至開花)是營養(yǎng)生長向生殖生長并進(jìn)的時(shí)期,決定穗數(shù),莖節(jié)間迅速生長,葉片增長,根系繼續(xù)擴(kuò)展,干物質(zhì)迅速增加,是玉米生長發(fā)育最旺盛的階段。玉米粒期(開花至成熟)以生殖生長為中心,經(jīng)過開花受精而進(jìn)入了籽粒產(chǎn)量形成時(shí)期,此時(shí)主要功能葉片是植株的中上層葉片,是決定粒數(shù)和粒質(zhì)量的重要時(shí)期。玉米成熟期中下部葉片會變黃,基部葉片變干枯,會降低利用光譜信息預(yù)測的精度。本研究篩選出的玉米最佳預(yù)測產(chǎn)量時(shí)期是灌漿前期,同前人研究一致[37-38]。作物產(chǎn)量的形成是一個(gè)長期不斷積累的過程,基于單一的生育時(shí)期進(jìn)行玉米產(chǎn)量預(yù)測可能代表性不足,而綜合利用多生育時(shí)期的信息更能夠全面、科學(xué)合理地反映產(chǎn)量形成過程,因此,如何綜合利用多時(shí)相的作物遙感信息進(jìn)行產(chǎn)量預(yù)測,還需有待進(jìn)一步研究。

    3.2 玉米產(chǎn)量預(yù)測參數(shù)類型的選擇

    本研究利用無人機(jī)搭載的傳感器獲得田塊尺度的可見光和多光譜的影像,對比可見光和多光譜影像產(chǎn)量預(yù)測的精度,以多光譜傳感器產(chǎn)量預(yù)測的精度更高。作物生長狀況在可見光(被葉綠素吸收)、近紅外(受葉片內(nèi)部構(gòu)造影響)[39]和紅邊波段(植物色素狀態(tài))呈現(xiàn)不同特征,而多光譜數(shù)據(jù)相比于可見光包含的信息更多,能夠?yàn)楫a(chǎn)量預(yù)測精度的提高提供潛力。通過植被指數(shù)單一類型作為參數(shù)能夠?qū)崿F(xiàn)產(chǎn)量的預(yù)測[40],但玉米產(chǎn)量和生長過程是多因素共同作用的結(jié)果,利用多類型(植被指數(shù)+紋理特征)的數(shù)據(jù)融合進(jìn)行作物產(chǎn)量估算[41],能夠充分利用不同數(shù)據(jù)的優(yōu)點(diǎn),從而達(dá)到提高模型精度和泛化能力。本研究利用多源影像數(shù)據(jù)計(jì)算了植被指數(shù)和紋理特征兩種參數(shù)對作物產(chǎn)量進(jìn)行預(yù)測,其中將植被指數(shù)與紋理特征進(jìn)行參數(shù)融合后的模型精度明顯高于單一類型參數(shù),與前人研究結(jié)果一致[42-44]。多傳感器數(shù)據(jù)融合可以增強(qiáng)光譜信息,提高估產(chǎn)的精度,但精度提高并不明顯,且多傳感器的數(shù)據(jù)融合不僅成倍增加了遙感作業(yè)的工作量,也增加了設(shè)備投入成本,性價(jià)比降低[9]。因此相比于融合多傳感器數(shù)據(jù),可以考慮融合同一傳感器的不同特征類型的數(shù)據(jù)。而作物的圖像特征除紋理特征外,還包括能描述圖像或圖像區(qū)域所對應(yīng)的景物的表面性質(zhì)的顏色特征、針對物體的外邊界圖像的輪廓和區(qū)域關(guān)系的形狀區(qū)域,以及圖像中分割出來的多個(gè)目標(biāo)之間相互的空間關(guān)系等特征。因此如何利用有限的光譜信息中獲取更多維度的信息,融合更多包含作物信息的輔助特征,從而進(jìn)一步提高模型精度有待更深入的研究。

    3.3 玉米產(chǎn)量預(yù)測機(jī)器學(xué)習(xí)方法的選擇

    機(jī)器學(xué)習(xí)算法可以充分利用遙感信息,在作物監(jiān)測方面表現(xiàn)出較好應(yīng)用潛力。但采用算法不同、監(jiān)測作物種類不同可能導(dǎo)致預(yù)測效果存在較大差異[21,23,26,45]。為了探究不同的機(jī)器學(xué)習(xí)算法對產(chǎn)量的預(yù)測效果,本研究選擇了BPNN、RF和SVM 3種方法用于比較夏玉米產(chǎn)量預(yù)測的潛力?;贐PNN模型可以學(xué)習(xí)和逼近任意復(fù)雜的非線性函數(shù)關(guān)系,較好的并行分布能力,但調(diào)參數(shù)目過多,因此對結(jié)果的解釋欠佳[46];RF模型可以集成多個(gè)決策樹來減少過擬合,在處理各種類型的數(shù)據(jù)時(shí),通常具有較高的準(zhǔn)確度[18,47];基于RF模型和SVM模型都能夠?qū)?shù)據(jù)產(chǎn)生的噪聲有較好的魯棒性,從而減少模型過擬合的風(fēng)險(xiǎn),但SVM由于需要計(jì)算核函數(shù)和對參數(shù)的設(shè)置調(diào)節(jié),參數(shù)的選擇可能對模型的性能產(chǎn)生較大的影響[48]。在本研究中基于RF構(gòu)建的訓(xùn)練集和測試集模型性能最佳。此外,深度學(xué)習(xí)算法已逐漸作為估測產(chǎn)量的新方法,其通過捕捉數(shù)據(jù)中復(fù)雜的模式和關(guān)系,在自動(dòng)學(xué)習(xí)特征和生成可靠的作物產(chǎn)量預(yù)測方面具有強(qiáng)大的潛力[49-50]。然而深度學(xué)習(xí)對數(shù)據(jù)及算力的要求較高,需要投入大量的人力物力和時(shí)間來開發(fā)新的算法和模型,相關(guān)研究和應(yīng)用案例還 較少。

    本研究期限較短,數(shù)據(jù)集較小,在后期的研究與應(yīng)用中還需要更多的數(shù)據(jù)積累驗(yàn)證與模型優(yōu)化。面對大尺度產(chǎn)量預(yù)測的應(yīng)用,如何將低空無人機(jī)遙感的細(xì)致性和高空衛(wèi)星遙感的宏觀性進(jìn)行結(jié)合從而實(shí)現(xiàn)不同尺度的作物產(chǎn)量估測還需更多的研究。產(chǎn)量由多因素影響,針對不同品種和氣候條件,未來可以考慮進(jìn)行時(shí)序特征分析和動(dòng)態(tài)指數(shù)的構(gòu)建,融合結(jié)構(gòu)特征(株高、葉面積指數(shù)等)、生長動(dòng)態(tài)參數(shù)(生長速率等)以及土壤、氣象等多源數(shù)據(jù)繼續(xù)優(yōu)化模型,構(gòu)建更加精準(zhǔn)和泛化能力強(qiáng)的玉米生長監(jiān)測與產(chǎn)量預(yù)測模型。

    4 結(jié) 論

    本研究基于無人機(jī)可見光和多光譜影像分別計(jì)算植被指數(shù)和紋理特征,并進(jìn)行了特征融合,豐富了反演建模的維度,利用相關(guān)系數(shù)篩選夏玉米拔節(jié)期、抽雄期、灌漿期、成熟期最優(yōu)建模參數(shù),比較BPNN、RF和SVM 3種機(jī)器學(xué)習(xí)對夏玉米產(chǎn)量預(yù)測的差異。相比于單一維度參數(shù)建模,將植被指數(shù)與紋理特征融合建模預(yù)測效果更佳,且基于多光譜數(shù)據(jù)構(gòu)建的模型優(yōu)于可見光,基于多光譜植被指數(shù)與紋理特征(MS_VI+T)的組合作為最佳建模參數(shù)。篩選最佳估產(chǎn)時(shí)期為夏玉米灌漿期(R2),以RF方法建模精度最好。通過多源數(shù)據(jù)的對比分析以及融合多特征數(shù)據(jù)和采用機(jī)器學(xué)習(xí)方法能夠較好地實(shí)現(xiàn)夏玉米產(chǎn)量預(yù)測。

    參考文獻(xiàn) Reference:

    [1]WART V J,KERSEBAUM C K,PENG S,et al.Estimating crop yield potential at regional to national scale[J].Field Crops Research,2013,143:34-43.

    [2]MUELLER N D,GERBER J S,JOHNSTO M,et al.Closing yield gaps through nutrient and water management[J].Nature,2012,490 (7419):254-7.

    [3]DENG L,MAO Z,LI X,et al.UAV-based multispectral remote sensing for precision agriculture:a comparison between different cameras [J].ISPRS Journal of Photogrammetry and Remote Sensing,2018,146:124-136.

    [4]MESHESHA T D,ABEJE M.Developing crop yield forecasting models for four major Ethiopian agricultural commodities [J].Remote Sensing Applications:Society and Environment,2018,11:83-93.

    [5]徐云飛,程 琦,魏祥平,等.變異系數(shù)法結(jié)合優(yōu)化神經(jīng)網(wǎng)絡(luò)的無人機(jī)冬小麥長勢監(jiān)測[J].農(nóng)業(yè)工程學(xué)報(bào),2021, 37(20):71-80.

    XU Y F,CHENG Q,WEI X P,et al.Monitoring of winter wheat growth under UAV using variation coefficient method and optimized neural network[J].Transactions of the Chinese Society of Agricultural Engineering,2021, 37(20):71-80.

    [6]CHIRANJIT S,CHANDRA K S.Rice crop growth monitoring with sentinel-1 SAR data using machine learning models in google earth engine cloud[J/OL].Remote Sensing Applications:Society and Environment,2023,32:1-2.https://doi.org/10.1016/j.rsase.2022.101029.

    [7]齊文棟,何黎明,王安鵬,等.低溫冷害脅迫下的多熟水稻產(chǎn)量遙感監(jiān)測[J].遙感技術(shù)與應(yīng)用,2023,38(3):558-565.

    QI W D,HE L M,WANG A P,et al. Remote sensing monitoring of yield loss of multiple cropping paddy caused by low temperature [J].Remote Sensing Technology and Application,2023,38(3):558-565.

    [8]ZHOU J,KHOT R L,BOYDSTON A R,et al. Low altitude remote sensing technologies for crop stress monitoring:a case study on spatial and temporal monitoring of irrigated pinto bean[J].Precision Agriculture,2018,19(3)555-569.

    [9]宋成陽,耿洪偉,費(fèi)帥鵬,等.基于多源數(shù)據(jù)的小麥品種產(chǎn)量估測研究[J].光譜學(xué)與光譜分析,2023,43(7):2210-2219.

    SONG CH Y,GENG H W,F(xiàn)EI SH P,et al. Study on yield estimation of wheat varieties based on multi-source data [J].Spectroscopy and Spectral Analysis,2023,43(7):2210-2219.

    [10] QU G,SHUAI Y,SHAO C,et al.County scale corn yield estimation based on Multi-Source Data in Liaoning province[J].Agronomy,2023,13(5):1428.

    [11]CHENG T,YANG Z,INOUE Y,et al.Preface:recent advances in remote sensing for crop growth monitoring[J].Remote Sensing,2016,8(2):116.

    [12]BO D,SHENGHUI F,YAN G,et al.Remote estimation of grain yield based on UAV data in different rice cultivars under contrasting climatic zone[J/OL].Field Crops Research,2021,267:1-11.https://doi.org/10.1016/j.fcr.2021.108148.

    [13]孫 陽,吳 瓊,馬 鋼.基于無人機(jī)影像的農(nóng)作物生長反演方法及應(yīng)用[J].測繪與空間地理信息,2021,44(9):165-166,170.

    SUN Y,WU Q,MA G.Crop growth inversion method and application based on UAV image[J].Geomatics & Spatial Information Technology,2021,44(9):165-166,170.

    [14]LEW T T S,SAROJAM R,JANG izT08UK0ApBZzUk6P24gtIvg6G+dUZc5vEz8Vk9VM50=I C,et al. Species-independent analytical tools for next-generation agriculture.[J].Nature Plants,2020,6(12):1408-1417.

    [15]FUHONG C,WEN L,WUXIONG S,et al.A mobile device-based imaging spectrometer for environmental monitoring by attaching a lightweight small module to a commercial digital camera[J/OL].Scientific Reports,2017, 7(1):1-9.doi:10.1038/s41598-017-15848-x.

    [16]KEITA H,NAOTO O,MOTOKI H,et al.High-dynamic-range spectral imaging system for omnidirectional scene capture[J].Journal of Imaging,2018,4(4):53.

    [17]TEAL K R,TUBANA B,GIRMA K,et al.In-season prediction of corn grain yield potential using normalized difference vegetation index[J].Agronomy Journal,2006,98(6):1488-1494.

    [18]王晶晶,李長碩,卓 越,等.基于多時(shí)相無人機(jī)遙感生育時(shí)期優(yōu)選的冬小麥估產(chǎn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2022,53(9):197-206.

    WANG J J,LI CH SH,ZH Y,et al.Yield estimation of winter wheat based on optimization of growth stages by multi-temporal UAV remote sensing[J].Transactions of the Chinese Society for Agricultural Machinery,2022, 53(9):197-206..

    [19]穆 軍,許 可,鐘曼麗,等.基于CiteSpace的農(nóng)作物遙感估產(chǎn)知識圖譜分析[J].現(xiàn)代農(nóng)業(yè)裝備,2020,41(2):41-48.

    MU J,XU K,ZHONG M L,et al.Knowledge structure analysis of crop yield estimation using remote sensing based on citespace [J].Modern Agricultural Equipment,2020,41(2):41-48.

    [20]祁佳峰,劉 笑,杜文玲,等.基于無人機(jī)高清影像的棉花單產(chǎn)預(yù)測[J].棉花學(xué)報(bào),2022,34(4):286-298.

    QI J F,LIU X,DU W L,et al.Cotton yield estimation based on UAV high-resolution images [J].Cotton Science,2022,34(4):286-298.

    [21]嚴(yán)海軍,卓 越,李茂娜,等.基于機(jī)器學(xué)習(xí)和無人機(jī)多光譜遙感的苜蓿產(chǎn)量預(yù)測[J].農(nóng)業(yè)工程學(xué)報(bào),2022,38(11):64-71.

    YAN H J,ZHUO Y,LI M N,et al.Alfalfa yield prediction using machine learning and UAV multispectral remote sensing[J].Transactions of the Chinese Society of Agricultural Engineering,2022,38(11):64-71.

    [22]邵亞杰,湯秋香,崔建平,等.融合無人機(jī)光譜信息與紋理特征的棉花葉面積指數(shù)估測[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2023, 54(6):186-196.

    SHAO Y J,TANG Q X,CUI J P,et al. Cotton leaf area index estimation combining UAV spectral and textural features [J].Transactions of the Chinese Society for Agricultural Machinery,2023,54(6):186-196.

    [23]ENHUI C,BING Z,DAILIANG P,et al.Wheat yield estimation using remote sensing data based on machine learning approaches[J/OL].Frontiers in Plant Science,2022,13:1-16.https://doi.10.3389/fpls.2022.1090970.

    [24]JOHNSON D M,HSIEH W W,CANNON J A,et al.Crop yield forecasting on the Canadian Prairies by remotely sensed vegetation indices and machine learning methods[J].Agricultural and Forest Meteorology,2016,218:74-84.

    [25]WENGANG Z,XIN G,LI H,et al.Comprehensive review of machine learning in geotechnical reliability analysis:Algorithms,applications and further challenges[J/OL].Applied Soft Computing Journal,2023,136:1-18.https://doi.org/10.1016/j.asoc.2023.110066.

    [26]競 霞,張 杰,王嬌嬌,等.水稻產(chǎn)量遙感監(jiān)測機(jī)器學(xué)習(xí)算法對比[J].光譜學(xué)與光譜分析,2022,42(5):1620-1627.

    JING X,ZHANG J,WANG J J,et al. Comparison of machine learning algorithms for remote sensing monitoring of rice yields [J].Spectroscopy and Spectral Analysis,2022,42(5):1620-1627.

    [27]于海洋,陳圣波,楊北萍,等.基于遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的玉米遙感估產(chǎn)方法[J].世界地質(zhì),2020,39(1):208-214.

    YU H Y,CHEN SH B,YANG B P,et al.Method of remote sensing estimation of corn yield based on genetic algorithm optimized BP neural network [J].Global Geology,2019,39(1):208-214.

    [28]ALEXIS P,L W M C,YANG C.Evaluation of random forests (RF) for regional and local-scale wheat yield prediction in southeast Australia[J].Sensors,2022,22(3):717-717.

    [29]劉 苗.不同水平氮磷配施對夏玉米產(chǎn)量和養(yǎng)分利用的影響[D].陜西楊凌:西北農(nóng)林科技大學(xué),2022.

    LIU M.Effects of different nitrogen and phosphorus application rate on yield,nutrient uptake and utilization of summer maize [D].Yangling Shaanxi:Northwest A&F University,2022.

    [30]王嘉盼,武紅旗,王德俊,等.基于無人機(jī)可見光影像與生理指標(biāo)的小麥估產(chǎn)模型研究[J].麥類作物學(xué)報(bào),2021, 41(10):1307-1315.

    WANG J P,WU H Q,WANG D J,et al.Research on wheat yield estimation model based on UAV visible light image and physiological index [J].Journal of Triticeae Crops,2021,41(10):1307-1315.

    [31]范軍亮,王 涵,廖振棋,等.基于紋理-顏色特征與植被指數(shù)融合的冬小麥LAI估測[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2023, 54(7):347-359.

    FAN J L,WANG H,LIAO ZH Q,et al.Winter wheat LAI estimation based on texture-color features and vegetation indices [J].Transactions of the Chinese Society for Agricultural Machinery,2023,54(7):347-359.

    [32]JOYCY S J,PRABAVATHY K A.Survey on automatic segmentation of relevant textures in agricultural images[J].Digital Image Processing,2012,4(3):123-129.

    [33]張 通,金 秀,饒 元,等.基于無人機(jī)多光譜的大豆旗葉光合作用量子產(chǎn)量反演方法[J].農(nóng)業(yè)工程學(xué)報(bào),2022,38(13):150-157.

    ZHANG T,JIN X,RAO Y,et al.Inversing photosynthesis quantum yield of the soybean flag leaf using a UAV-carrying multispectral camera [J].Transactions of the Chinese Society of Agricultural Engineering,2012, 38(13):150-157.

    [34]敖 登,楊佳慧,丁維婷,等.54種植被指數(shù)研究進(jìn)展綜述[J].安徽農(nóng)業(yè)科學(xué),2023,51(1):13-21,28.

    AO D,YANG J H,DING W T,et al.Review of 54 vegetati on indices [J].Journal of Anhui Agricultural Sciences,2023,51(1):13-21,28.

    [35]YUAN L,CHENWEI N,ZHEN Z, et al. Evaluating how lodging affects maize yield estimation based on UAV observations[J/OL].Frontiers in Plant Science,2023,13:1-16.https://doi:10.3389/fpls.2022.979103.

    [36]杜菲菲,安慧君,李賀新.應(yīng)用光譜指數(shù)和機(jī)器學(xué)習(xí)反演紫丁香葉片的葉綠素面密度[J].東北林業(yè)大學(xué)學(xué)報(bào),2024,52(2):75-83.

    DU F F,AN H J,LI H X.Inversion of chlorophyll contents in syringe oblata leaves using spectral index and MPA-TSVR [J].Journal of Northeast Forestry University,2024,52(2):75-83.

    [37]葉 強(qiáng),楊鳳海,劉煥軍,等.引入地形特征的田塊尺度玉米遙感估產(chǎn)與空間格局分析[J].科學(xué)技術(shù)與工程,2021, 21(24):10215-10221.

    YE Q,YANG F H,LIU H J,et al.Maize yield estimation and spatial pattern analysis based on topographic features by remote sensing at the field scale [J].Science Technology and Engineering,2021,21(24):10215-10221.

    [38]馬冠南,丁春雨.基于植被指數(shù)和最優(yōu)物候期的玉米產(chǎn)量預(yù)測研究[J].農(nóng)業(yè)與技術(shù),2017,37(9):17-18,25.

    MA G N,DING CH Y.Prediction of maize yield based on vegetation index and optimal phenological period [J].Agriculture and Technology,2017,37(9):17-18,25.

    [39]張芳毓,譚永毅,聶 婧,等.無人機(jī)多光譜在農(nóng)業(yè)中的應(yīng)用[J].智慧農(nóng)業(yè)導(dǎo)刊,2022,2(24):11-13.

    ZHANG F Y,TAN Y Y,NIE J,et al. Application of UAV multi-spectrum in agriculture [J].Journal of Smart Agriculture,2022,2(24):11-13.

    [40]韓文霆,彭星碩,張立元,等.基于多時(shí)相無人機(jī)遙感植被指數(shù)的夏玉米產(chǎn)量估算[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,51(1):148-155.

    HAN W T,PENG X SH,ZHANG L Y,et al.Summer maize yield estimation based on vegetation index derived from multi-temporal UAV remote sensing [J].Transactions of the Chinese Society for Agricultural Machinery, 2019,51(1):148-155.

    [41]劉欣誼,仲曉春,陳 晨,等.利用無人機(jī)圖像顏色與紋理特征數(shù)據(jù)在小麥生育前期對產(chǎn)量進(jìn)行預(yù)測[J].麥類作物學(xué)報(bào),2020,40(8):1002-1007.

    LIU X Y,ZHONG X CH,CHEN CH,et al. Prediction of wheat yield using color and texture feature data of UAV image at early growth stage [J].Journal of Triticeae Crops,2020,40(8):1002-1007.

    [42]QUAN W,ANDIKA N P,YI G,et al.Combining both spectral and textural indices for alleviating saturation problem in forest LAI estimation using Sentinel-2 data[J].Geocarto International,2022,37(25):10511-10531.

    [43]賈 丹,陳鵬飛.低空無人機(jī)影像分辨率對冬小麥氮濃度反演的影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2020,51(7):164-169.

    JIA D,CHEN P F.Effect of low-altitude UAV image resolution on inversion of winter wheat nitrogen concentration[J].Transactions of the Chinese Society for Agricultural Machinery,2020,51(7):164-169.

    [44]JIAYI Z,XIAOLEI Q,YUETING W,et al.Combining texture,color,and vegetation indices from fixed-wing UAS imagery to estimate wheat growth parameters using multivariate regression methods[J/OL].Computers and Electronics in Agriculture,2021,185:1-14.https://doi.org/101016/j.compag.2021.106138.

    [45]崔孟然.基于機(jī)器學(xué)習(xí)的馬鈴薯產(chǎn)量高光譜估算[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué),2022.

    CUI M R.Machine learning based hyperspectral estimation of potato yield[D].Hohhot:Inner Mongolia Agricultural University,2022.

    [46]樂夏唯.基于機(jī)器學(xué)習(xí)的傾斜基板激光增材制造工藝優(yōu)化研究[D].江蘇徐州:中國礦業(yè)大學(xué),2023.

    LE X W.Optimisation of laser additive manufacturing process for inclined substrates based on machine learning [D].Xuzhou Jiangsu :China University of Mining and Technology,2023.

    [47]JEONG H J,RESOP P J,MUELLER D N,et al. Random forests for global and regional crop yield predictions.[J].PLoS One,2017,11(6):e0156571.

    [48]LIU H,SHI T,CHEN Y,et al.Improving spectral estimation of soil organic carbon content through semi-supervised regression[J].Remote Sensing,2017,9(1):29-29.

    [49]DILLI P,ALLARD W D,HENDRIK B,et al. Interpretability of deep learning models for crop yield forecasting[J/OL].Computers and Electronics in Agriculture,2023,206:1-14.https://doi.org/10.1016/j.compay.2023.107663.

    [50]S K M,TIMOTHY D,MBULISI S,et al. Fine-scale characterization of irrigated and rainfed croplands at national scale using multi-source data,random forest,and deep learning algorithms[J].ISPRS Journal of Photogrammetry and Remote Sensing,2023,204:117-130.

    Remote Sensing Estimation for Maize Yield Integrating Spectral and Texture Features

    MA Yuanhua, WANG Leyin, ZHANG Zhenxin, ZHENG Dasheng, YE Yulan,CUI Zhifeng, DU Bingxiao, CUN Yujie, LI Jun and WANG Rui

    (College of Agronomy,Northwest A&F University,Yangling Shaanxi 712100,China)

    Abstract This study investigates the integration of spectral and texture features for pridicting maize yield using remote sensing techniques.Visible light and multispectral images captured by drone were obtained at four key growth stages of summer maize:jointing, tasseling, filling, and maturity.Vegetation indexes and texture feature parameters were extracted and used to formulate spectral and texture feature fusion variables.Three machine learning methods(a back propagation neural network, a random forest, and a support vector machine)were employed to develop predictive models for predicting maize yield.The results showed that the accuracy of maize yield prediction models incorporating both spectral and texture feature fusion was superior to that of models using only spectral or texture features.Among the three machine learning methods, the random forest yielded the best maize prediction model.Among the four growth stages, the model for the filling stage (blister) had the highest accuracy.An overall evaluation of the results showed that the model obtained using the fused vegetation indexes and texture feature variables extracted from the multispectral images using the random forest method predicted maize yield with the highest accuracy.The fusion of multiple remote sensing features and the application of machine learning methods allow the extraction and utilization of more information, thereby improving the accuracy and robustness of yield prediction.

    Key words Maize; UAV remote sensing; RGB; Multi-spectral; Vegetation index; Texture features; Yield

    Received 2023-10-16 Returned 2024-01-21

    Foundation item Natural Science Basic Research Program of Shaanxi Province(No.2023-JC-QN-0192).

    First author MA Yuanhua,female,master student.Research area:crop monitoring by remote sensing.E-mail:1162392907@qq.com

    Corresponding author WANG Rui, male, associate professor, master supervisor.Research area:crop monitoring by remote sensing.E-mail:rico@nwafu.edu.cn

    (責(zé)任編輯:史亞歌 Responsible editor:SHI Yage)

    基金項(xiàng)目:陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃(2023-JC-QN-0192)。

    第一作者:馬元花,女,碩士研究生,研究方向?yàn)樽魑镞b感監(jiān)測。E-mail:1162392907@qq.com

    通信作者:王 瑞,男,副教授,碩士生導(dǎo)師,研究方向?yàn)樽魑镞b感監(jiān)測。E-mail:rico@nwafu.edu.cn

    美女大奶头视频| 亚洲国产欧美网| 久久草成人影院| 亚洲成a人片在线一区二区| 最新在线观看一区二区三区| 校园春色视频在线观看| 老司机午夜福利在线观看视频| 99久久成人亚洲精品观看| 少妇的丰满在线观看| svipshipincom国产片| 国产色爽女视频免费观看| 99热这里只有是精品50| 日韩国内少妇激情av| 婷婷亚洲欧美| 久久精品人妻少妇| 国内揄拍国产精品人妻在线| 久久久久久国产a免费观看| 深爱激情五月婷婷| 成年人黄色毛片网站| 男人舔女人下体高潮全视频| www日本黄色视频网| 国产精品免费一区二区三区在线| 欧美日韩乱码在线| 国产精华一区二区三区| 老鸭窝网址在线观看| 午夜精品在线福利| 成年免费大片在线观看| 99热精品在线国产| 成人高潮视频无遮挡免费网站| 久久精品91蜜桃| 国产男靠女视频免费网站| 黄色女人牲交| 听说在线观看完整版免费高清| 老鸭窝网址在线观看| 亚洲人成伊人成综合网2020| 两个人看的免费小视频| 激情在线观看视频在线高清| 国产精品亚洲美女久久久| 白带黄色成豆腐渣| 精品人妻1区二区| 他把我摸到了高潮在线观看| 天美传媒精品一区二区| 可以在线观看的亚洲视频| 欧美成人免费av一区二区三区| 天堂网av新在线| 亚洲无线观看免费| 国产亚洲精品久久久久久毛片| 久久欧美精品欧美久久欧美| 18禁国产床啪视频网站| 岛国在线免费视频观看| 很黄的视频免费| 久久久久精品国产欧美久久久| 午夜免费男女啪啪视频观看 | 亚洲成人免费电影在线观看| 久久精品国产自在天天线| 午夜亚洲福利在线播放| 日韩欧美在线二视频| 国产精品一区二区三区四区免费观看 | 国产乱人视频| 听说在线观看完整版免费高清| 亚洲欧美激情综合另类| 中文字幕av在线有码专区| 免费在线观看影片大全网站| 淫秽高清视频在线观看| 狂野欧美激情性xxxx| 无限看片的www在线观看| 午夜影院日韩av| av在线天堂中文字幕| aaaaa片日本免费| 中文资源天堂在线| 久久国产乱子伦精品免费另类| 欧美性感艳星| 国产69精品久久久久777片| 亚洲国产日韩欧美精品在线观看 | 丝袜美腿在线中文| 免费观看的影片在线观看| 丰满的人妻完整版| 久久久久国内视频| av黄色大香蕉| 国产精品98久久久久久宅男小说| 欧美bdsm另类| 午夜亚洲福利在线播放| 中文字幕精品亚洲无线码一区| 丰满人妻熟妇乱又伦精品不卡| netflix在线观看网站| 亚洲成人免费电影在线观看| 日韩 欧美 亚洲 中文字幕| 国产男靠女视频免费网站| 亚洲成人免费电影在线观看| 免费av观看视频| 国产欧美日韩精品一区二区| 亚洲男人的天堂狠狠| 久久国产精品人妻蜜桃| 欧美日韩瑟瑟在线播放| 欧美三级亚洲精品| 少妇高潮的动态图| 在线观看舔阴道视频| av天堂中文字幕网| 精品国产三级普通话版| 观看美女的网站| 一进一出抽搐动态| 最新美女视频免费是黄的| 搞女人的毛片| 无遮挡黄片免费观看| 欧美一区二区国产精品久久精品| 久久九九热精品免费| 欧美丝袜亚洲另类 | 九九久久精品国产亚洲av麻豆| 九九热线精品视视频播放| 在线观看午夜福利视频| 一级作爱视频免费观看| 久久国产精品人妻蜜桃| 一级毛片高清免费大全| 亚洲成人久久性| 午夜福利欧美成人| 亚洲av电影在线进入| 性色avwww在线观看| 国产精品久久久久久久久免 | 日本一二三区视频观看| 国产视频内射| 在线播放无遮挡| 12—13女人毛片做爰片一| 丰满的人妻完整版| 亚洲黑人精品在线| 精品99又大又爽又粗少妇毛片 | 日韩有码中文字幕| 最近最新中文字幕大全电影3| 亚洲最大成人中文| 一级a爱片免费观看的视频| 国产精品日韩av在线免费观看| 高清毛片免费观看视频网站| 99热这里只有精品一区| 亚洲国产精品合色在线| 在线观看免费视频日本深夜| 精品国产美女av久久久久小说| 丝袜美腿在线中文| 99热这里只有精品一区| 国产毛片a区久久久久| av天堂在线播放| 老汉色av国产亚洲站长工具| 国产欧美日韩一区二区精品| 亚洲在线自拍视频| 亚洲国产欧美网| 亚洲va日本ⅴa欧美va伊人久久| 国产爱豆传媒在线观看| 在线免费观看的www视频| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 国内精品一区二区在线观看| 中文字幕人成人乱码亚洲影| 国产精品女同一区二区软件 | 国产伦精品一区二区三区四那| 亚洲欧美日韩高清在线视频| 一个人看视频在线观看www免费 | 亚洲七黄色美女视频| 一进一出抽搐动态| 亚洲不卡免费看| 亚洲在线观看片| 伊人久久精品亚洲午夜| 精品不卡国产一区二区三区| 一进一出好大好爽视频| 亚洲精品久久国产高清桃花| 免费在线观看成人毛片| 免费看日本二区| 蜜桃久久精品国产亚洲av| 国产精品乱码一区二三区的特点| 人妻夜夜爽99麻豆av| 久久婷婷人人爽人人干人人爱| 日本a在线网址| 一级作爱视频免费观看| 51午夜福利影视在线观看| 国产成人福利小说| 国产男靠女视频免费网站| 亚洲av成人不卡在线观看播放网| 国产精品一区二区免费欧美| 久久精品亚洲精品国产色婷小说| 婷婷精品国产亚洲av| 在线观看av片永久免费下载| 国产黄片美女视频| 免费av不卡在线播放| 亚洲国产日韩欧美精品在线观看 | 9191精品国产免费久久| 国产毛片a区久久久久| 精品国产亚洲在线| 亚洲五月婷婷丁香| 亚洲成a人片在线一区二区| 欧美另类亚洲清纯唯美| 九色国产91popny在线| 欧美乱码精品一区二区三区| 激情在线观看视频在线高清| 国产爱豆传媒在线观看| 五月伊人婷婷丁香| 香蕉av资源在线| 亚洲性夜色夜夜综合| 日本 欧美在线| 久久久久九九精品影院| 午夜激情欧美在线| 内地一区二区视频在线| 可以在线观看毛片的网站| 在线十欧美十亚洲十日本专区| 国产亚洲精品久久久久久毛片| 在线观看av片永久免费下载| 少妇人妻一区二区三区视频| 精品国产美女av久久久久小说| 亚洲国产欧美人成| 国产淫片久久久久久久久 | 色综合欧美亚洲国产小说| 国产一区二区三区在线臀色熟女| 精品不卡国产一区二区三区| 国产亚洲精品av在线| 欧美最新免费一区二区三区 | 18禁裸乳无遮挡免费网站照片| 哪里可以看免费的av片| 欧美性猛交黑人性爽| www.999成人在线观看| a级毛片a级免费在线| 十八禁网站免费在线| 听说在线观看完整版免费高清| 午夜亚洲福利在线播放| 天天躁日日操中文字幕| 高清毛片免费观看视频网站| 久99久视频精品免费| 亚洲av不卡在线观看| 99热这里只有是精品50| 国内久久婷婷六月综合欲色啪| 制服人妻中文乱码| 中文字幕av在线有码专区| 国内精品一区二区在线观看| 一级黄色大片毛片| 国产乱人视频| eeuss影院久久| 国产私拍福利视频在线观看| 国产精品三级大全| 亚洲av第一区精品v没综合| 日日干狠狠操夜夜爽| 国产精华一区二区三区| 欧美极品一区二区三区四区| 国产伦在线观看视频一区| 午夜福利在线在线| 日韩 欧美 亚洲 中文字幕| 男女床上黄色一级片免费看| 色综合站精品国产| 欧美激情在线99| 美女 人体艺术 gogo| 99热这里只有精品一区| 亚洲欧美日韩高清在线视频| 国产亚洲精品综合一区在线观看| 婷婷亚洲欧美| 日本黄大片高清| 亚洲自拍偷在线| 99在线人妻在线中文字幕| 最近最新中文字幕大全免费视频| 欧美日韩福利视频一区二区| 欧美极品一区二区三区四区| 丁香六月欧美| 日韩精品中文字幕看吧| 怎么达到女性高潮| 精品久久久久久久人妻蜜臀av| 色老头精品视频在线观看| 少妇人妻一区二区三区视频| 亚洲熟妇熟女久久| 精品熟女少妇八av免费久了| av黄色大香蕉| 久久久久国产精品人妻aⅴ院| 午夜福利在线观看免费完整高清在 | 亚洲欧美日韩东京热| 成年女人看的毛片在线观看| www日本黄色视频网| 久久久久久久午夜电影| 中亚洲国语对白在线视频| aaaaa片日本免费| 国产又黄又爽又无遮挡在线| netflix在线观看网站| 999久久久精品免费观看国产| 亚洲精品一区av在线观看| 国产精品,欧美在线| 日韩欧美精品免费久久 | 国模一区二区三区四区视频| 中出人妻视频一区二区| 免费在线观看日本一区| av福利片在线观看| 欧美av亚洲av综合av国产av| 色综合站精品国产| 99在线视频只有这里精品首页| 亚洲天堂国产精品一区在线| 国产亚洲精品综合一区在线观看| 国产一级毛片七仙女欲春2| 日韩欧美三级三区| 91在线精品国自产拍蜜月 | 身体一侧抽搐| 97超视频在线观看视频| 天天躁日日操中文字幕| 婷婷丁香在线五月| 国产精品野战在线观看| 少妇熟女aⅴ在线视频| 国产av在哪里看| 九色成人免费人妻av| 国产亚洲av嫩草精品影院| 成年女人毛片免费观看观看9| 中文字幕久久专区| 婷婷精品国产亚洲av| 日本与韩国留学比较| 国产成人福利小说| x7x7x7水蜜桃| 国产激情偷乱视频一区二区| 99视频精品全部免费 在线| 国产真人三级小视频在线观看| 90打野战视频偷拍视频| 好男人在线观看高清免费视频| 国产成人av激情在线播放| 精品不卡国产一区二区三区| 国内揄拍国产精品人妻在线| 国产精品女同一区二区软件 | 老司机在亚洲福利影院| 日本与韩国留学比较| 长腿黑丝高跟| 国产美女午夜福利| 91九色精品人成在线观看| 91av网一区二区| 久久久久久久午夜电影| 欧美乱妇无乱码| svipshipincom国产片| 日本熟妇午夜| 日本五十路高清| 综合色av麻豆| 一级作爱视频免费观看| 国产精品一区二区免费欧美| 在线十欧美十亚洲十日本专区| 男人舔奶头视频| 成年女人永久免费观看视频| 国产精华一区二区三区| 国产高清videossex| 国产毛片a区久久久久| 在线视频色国产色| 天天添夜夜摸| 国产精品美女特级片免费视频播放器| 精品人妻一区二区三区麻豆 | 热99re8久久精品国产| 国产乱人伦免费视频| a级毛片a级免费在线| 级片在线观看| 女人高潮潮喷娇喘18禁视频| 午夜福利视频1000在线观看| 亚洲成av人片在线播放无| 亚洲真实伦在线观看| 亚洲国产精品sss在线观看| 哪里可以看免费的av片| 女人被狂操c到高潮| 不卡一级毛片| 18禁国产床啪视频网站| 国产精品嫩草影院av在线观看 | 亚洲真实伦在线观看| 99国产精品一区二区蜜桃av| 精品一区二区三区av网在线观看| 男女下面进入的视频免费午夜| 国产中年淑女户外野战色| 国产老妇女一区| 国产精品自产拍在线观看55亚洲| 成年免费大片在线观看| 色播亚洲综合网| 在线播放国产精品三级| 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| 久久亚洲真实| 欧美日韩亚洲国产一区二区在线观看| 久久精品亚洲精品国产色婷小说| 国产一区二区在线观看日韩 | 日本 欧美在线| 观看免费一级毛片| 久久精品人妻少妇| 制服丝袜大香蕉在线| 亚洲第一电影网av| 99久久久亚洲精品蜜臀av| 亚洲av第一区精品v没综合| 国产精品 欧美亚洲| 国语自产精品视频在线第100页| 香蕉av资源在线| 久久这里只有精品中国| 国产一区二区三区在线臀色熟女| 欧美丝袜亚洲另类 | 19禁男女啪啪无遮挡网站| 18禁黄网站禁片午夜丰满| 天堂网av新在线| 国产爱豆传媒在线观看| 狠狠狠狠99中文字幕| 成人18禁在线播放| 动漫黄色视频在线观看| 成人特级av手机在线观看| 欧美一级毛片孕妇| 变态另类成人亚洲欧美熟女| 亚洲精品色激情综合| 精品一区二区三区视频在线 | 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 亚洲色图av天堂| 一二三四社区在线视频社区8| 无人区码免费观看不卡| 国产精品国产高清国产av| 亚洲激情在线av| 18禁黄网站禁片午夜丰满| 日本熟妇午夜| 美女黄网站色视频| 免费在线观看亚洲国产| 久久精品亚洲精品国产色婷小说| www.色视频.com| 三级毛片av免费| 757午夜福利合集在线观看| 波多野结衣高清作品| 色视频www国产| 女同久久另类99精品国产91| 美女高潮喷水抽搐中文字幕| 日本撒尿小便嘘嘘汇集6| 少妇高潮的动态图| 五月玫瑰六月丁香| 少妇人妻精品综合一区二区 | 老汉色av国产亚洲站长工具| 少妇人妻精品综合一区二区 | 国产精品三级大全| 法律面前人人平等表现在哪些方面| 国产精品久久久久久久电影 | 国产激情欧美一区二区| 成人三级黄色视频| 高清毛片免费观看视频网站| 18+在线观看网站| 久久99热这里只有精品18| 亚洲人成电影免费在线| 岛国在线观看网站| 丰满乱子伦码专区| 亚洲人成网站在线播| 亚洲中文字幕一区二区三区有码在线看| 国产精品自产拍在线观看55亚洲| 国产色爽女视频免费观看| 村上凉子中文字幕在线| 俺也久久电影网| 日日摸夜夜添夜夜添小说| 美女免费视频网站| 欧美乱妇无乱码| 韩国av一区二区三区四区| 高清日韩中文字幕在线| 三级毛片av免费| 精品福利观看| 国产精品自产拍在线观看55亚洲| 国产色婷婷99| 日韩欧美一区二区三区在线观看| av片东京热男人的天堂| 身体一侧抽搐| 性色av乱码一区二区三区2| 国产高清三级在线| 久久久久亚洲av毛片大全| 熟女少妇亚洲综合色aaa.| 欧美性感艳星| 亚洲人与动物交配视频| 国产精品乱码一区二三区的特点| 婷婷精品国产亚洲av在线| 女人十人毛片免费观看3o分钟| 国产又黄又爽又无遮挡在线| 人人妻,人人澡人人爽秒播| 日本 欧美在线| 午夜视频国产福利| 99久久精品热视频| 亚洲欧美激情综合另类| 精品一区二区三区人妻视频| 免费搜索国产男女视频| 国产激情欧美一区二区| 热99在线观看视频| 别揉我奶头~嗯~啊~动态视频| 搡老熟女国产l中国老女人| 18禁国产床啪视频网站| 狂野欧美激情性xxxx| 国产高潮美女av| 欧美成人性av电影在线观看| 免费观看精品视频网站| 成人三级黄色视频| a在线观看视频网站| 亚洲午夜理论影院| 男女做爰动态图高潮gif福利片| 精品一区二区三区视频在线 | 欧美午夜高清在线| e午夜精品久久久久久久| 国产黄片美女视频| or卡值多少钱| 久久久久精品国产欧美久久久| 99国产极品粉嫩在线观看| 男插女下体视频免费在线播放| 国产精品乱码一区二三区的特点| 国产成人福利小说| 欧洲精品卡2卡3卡4卡5卡区| 亚洲 国产 在线| 青草久久国产| 欧美+日韩+精品| 99国产精品一区二区蜜桃av| svipshipincom国产片| 国产三级中文精品| 国产欧美日韩一区二区精品| 黄色成人免费大全| 我要搜黄色片| 免费观看精品视频网站| 亚洲美女视频黄频| 欧美3d第一页| 91字幕亚洲| 十八禁人妻一区二区| 99久久九九国产精品国产免费| 噜噜噜噜噜久久久久久91| 黑人欧美特级aaaaaa片| 男人舔女人下体高潮全视频| 少妇人妻精品综合一区二区 | 十八禁网站免费在线| 欧美黄色淫秽网站| 国产免费一级a男人的天堂| 99国产精品一区二区蜜桃av| 成人欧美大片| 欧美黄色片欧美黄色片| 欧美中文综合在线视频| 国产av在哪里看| 亚洲成人久久爱视频| 夜夜看夜夜爽夜夜摸| 男女之事视频高清在线观看| 看片在线看免费视频| 国产精品免费一区二区三区在线| 亚洲男人的天堂狠狠| 精品久久久久久,| 18禁黄网站禁片午夜丰满| 国产精品久久久久久久电影 | 91九色精品人成在线观看| 19禁男女啪啪无遮挡网站| 精品一区二区三区av网在线观看| 午夜福利欧美成人| 老司机午夜福利在线观看视频| 日韩欧美 国产精品| 久久九九热精品免费| 欧美日韩精品网址| 91在线精品国自产拍蜜月 | 哪里可以看免费的av片| 亚洲av成人精品一区久久| 在线观看一区二区三区| 亚洲av不卡在线观看| 美女免费视频网站| 99久久精品国产亚洲精品| a级一级毛片免费在线观看| 99热只有精品国产| 法律面前人人平等表现在哪些方面| 国产真人三级小视频在线观看| 国产亚洲欧美98| 女同久久另类99精品国产91| av在线天堂中文字幕| 国产精品野战在线观看| 亚洲熟妇中文字幕五十中出| 级片在线观看| 一a级毛片在线观看| 两人在一起打扑克的视频| 日韩人妻高清精品专区| 狠狠狠狠99中文字幕| 国产精品野战在线观看| 天堂√8在线中文| 欧美最黄视频在线播放免费| 波多野结衣高清作品| 高潮久久久久久久久久久不卡| 五月玫瑰六月丁香| 尤物成人国产欧美一区二区三区| 精品国产三级普通话版| 色综合欧美亚洲国产小说| 午夜福利18| 中文字幕精品亚洲无线码一区| 国产爱豆传媒在线观看| av女优亚洲男人天堂| 精品国产亚洲在线| 日韩欧美精品免费久久 | 国产午夜精品论理片| 嫩草影院精品99| 亚洲国产欧美网| 男女之事视频高清在线观看| 成年版毛片免费区| a在线观看视频网站| 国产午夜福利久久久久久| 国产欧美日韩一区二区精品| 国产亚洲精品久久久久久毛片| 99久久无色码亚洲精品果冻| 内射极品少妇av片p| 观看美女的网站| 变态另类成人亚洲欧美熟女| 国产淫片久久久久久久久 | 日韩欧美国产在线观看| 黄片小视频在线播放| 尤物成人国产欧美一区二区三区| 午夜免费观看网址| 国产在线精品亚洲第一网站| 日韩欧美三级三区| 精品久久久久久久久久免费视频| 亚洲熟妇中文字幕五十中出| 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 久久精品国产综合久久久| 国产精品自产拍在线观看55亚洲| 国产伦在线观看视频一区| 听说在线观看完整版免费高清| 内射极品少妇av片p| av视频在线观看入口| 黄色日韩在线| 美女cb高潮喷水在线观看| 一区二区三区激情视频| 最近在线观看免费完整版| 亚洲 国产 在线| 免费看日本二区| 岛国视频午夜一区免费看| 国产av麻豆久久久久久久| 久久久久久大精品| 老熟妇乱子伦视频在线观看| tocl精华| 日本与韩国留学比较| 成人18禁在线播放| 成年女人永久免费观看视频| 国产熟女xx| 最近最新免费中文字幕在线| 国产精品综合久久久久久久免费| 性欧美人与动物交配|