北京農(nóng)業(yè)信息技術研究中心

國家農(nóng)業(yè)信息化工程技術研究中心

南京農(nóng)業(yè)大學大豆研究所/國家大豆改良中心

大豆是世界上最重要的糧食作物之一，最常用來做豆制品、榨取豆油、釀造醬油和提取蛋白質。而近年來，中國的大豆進口量逐年增加，從側面反映了中國的大豆目前仍具有較大的高產(chǎn)潛力。利用高光譜技術可以進行作物生長監(jiān)測和產(chǎn)量估測，為實現(xiàn)簡單、快速、高效和無損地估產(chǎn)及監(jiān)測作物生長提供了有效的途徑。

為探討無人機載高光譜空間尺度對大豆產(chǎn)量預測精度的影響，以山東嘉祥圣豐大豆為研究對象，設計以多旋翼無人機為平臺搭載 Cubert UHD185（其詳細技術參數(shù)見表1） 成像高光譜傳感器的無人機遙感農(nóng)情監(jiān)測系統(tǒng)，光譜垂直向下測量，并獲取了大豆多個生育期的無人機高光譜數(shù)據(jù)。

表1 UHD185技術參數(shù)

利用 Cubert-Pilot 軟件對每張高光譜圖像和同時采集的相應全色圖像進行融合，得到融合后的高光譜影像；其次使用Agisoft PhotoScan軟件，借助全色圖像的點云數(shù)據(jù)完成高光譜圖像的拼接。以小區(qū)幾何中心為中心，由拼接好的高光譜影像提取冠層光譜反射率。如圖1所示為得到每個小區(qū) 21 個不同空間尺度梯度所對應的冠層光譜反射率，可以看出不同空間采樣面積所對應的冠層光譜反射率在可見光波段無明顯差異，但在近紅外區(qū)域的差異較為明顯提取冠層光譜反射率。

圖1 21個不同采樣面積所對應的冠層光譜反射率

在已有研究成果基礎上， 本文選取 4 種植被指數(shù)建立大豆產(chǎn)量與無人機成像高光譜遙感監(jiān)測模型，植被指數(shù)計算方法如表2所示。

表2 植被指數(shù)計算公式

基于 21 組空間尺度范圍下的大豆冠層高光譜數(shù)據(jù)分別與實測產(chǎn)量建立 PLSR 回歸方程， 對比不同空間尺度下的模型預測精度，進而確定最優(yōu)空間尺度。產(chǎn)量和 4 個植被指數(shù)的 PLSR 模型相關系數(shù) r 隨空間尺度面積的變化趨勢曲線見圖 2所示。相關系數(shù) r 介于 0.795～0.812 之間，大豆冠層光譜植被指數(shù)與大豆產(chǎn)量極顯著相關，在空間取樣面積小于 8 m2 時，建模精度隨著空間取樣面積逐漸增大而增大，且隨著取樣面積的增大，建模精度逐漸趨于平緩，當光譜取樣面積大于 10.13 m2后，建模精度隨取樣面積的增大而減小。光譜空間尺度為 4.25 m×2.125m和4.5 m×2.25 m時，即9.03和10.13 m2時，相關系數(shù)r達到最大值，最大相關系數(shù)約為 0.8117?？梢?，在基于無人機高光譜影像輔助大豆估產(chǎn)研究時，所取光譜的空間范圍長、寬與小區(qū)總長、寬之間的比例應該介于 4.25:5 和 4.5:5，這樣不僅能夠有效降低小區(qū)周圍土壤背景噪聲和邊緣陰影光譜的影響，又能夠避免大豆冠層信息的缺失。

圖2 產(chǎn)量和4個植被指數(shù)的PLSR模型相關系數(shù)和RMSE隨空間采樣窗口面積的變化曲線

以4.5×2.25m為空間取樣范圍提取的冠層光譜為最終光譜，建立產(chǎn)量和植被指數(shù)之間的PLSR模型，實際產(chǎn)量和預測產(chǎn)量的散點分布如圖3所示。所建立大豆估產(chǎn)模型為：

產(chǎn)量=68.8GNDVI+79.2NDVI+26.2RVI+293.7MSAVI2-227.87

使用無人機遙感平臺進行大豆估產(chǎn)，對實現(xiàn)大規(guī)模大豆產(chǎn)量高通量快速獲取具有重要意義，能夠幫助育種學家快速篩選大豆優(yōu)良品種。本研究基于無人機 UHD185 光譜成像儀獲取的大豆不同時期的冠層光譜對最優(yōu)空間尺度進行了分析， 結果表明，基于不同空間尺度計算的波譜數(shù)據(jù)在近紅外波段差異較大，且品系間各植被指數(shù)在不同空間尺度處的顯著性差異水平不同?；谠几吖庾V圖像的空間分辨率計算各空間取樣面積處的光譜信息，進而對空間尺度進行優(yōu)化選擇，在未來的研究中，可以對最優(yōu)空間取樣面積下的光譜進行重采樣處理，分析不同分辨率下光譜對大豆估產(chǎn)的影響。同時也可考慮基于不同試驗條件下的不同作物，探索無人機高光譜影像的空間尺度規(guī)律。另外，高光譜具有較高的光譜分辨率，可針對波譜信息進行深入分析對光譜空間尺度進行優(yōu)化選擇。

注：由于Cubert廠家2016年產(chǎn)品升級， UHD185已更名為S185。