背景：

SPAD（Soil plant analysis development）即葉片綠度，反映葉片中葉綠素的相對(duì)含量。葉綠素是作物光合作用最重要的色素，其濃度變化直接影響作物的健康狀況。傳統(tǒng)葉綠素體外測(cè)量方法不適合大面積的實(shí)時(shí)監(jiān)測(cè)。近年來(lái)，遙感技術(shù)已成為作物生長(zhǎng)監(jiān)測(cè)的重要工具。特別是基于無(wú)人機(jī)的低空遙感探測(cè)技術(shù)，具有操作方便、快速靈活、效率高、空間分辨率理想等優(yōu)點(diǎn)，在精確農(nóng)業(yè)中得到了廣泛應(yīng)用。搭載在無(wú)人機(jī)中的高光譜傳感器具有較高的光譜分辨率和較強(qiáng)的波段連續(xù)性，可以準(zhǔn)確獲取冬小麥的光譜信息。

前人研究中，可以通過(guò)植被指數(shù)提取小麥生長(zhǎng)信息，或者采用支持向量回歸、隨機(jī)森林、反向傳播神經(jīng)網(wǎng)絡(luò)等方法監(jiān)測(cè)農(nóng)學(xué)參數(shù)。但是現(xiàn)有算法基于全局建模思想，無(wú)法最大限度地提高不同土壤肥力條件下小麥光譜特征差異。此外，在同一生育期，不同施肥處理下的冬小麥光譜特征也存在顯著差異。因此，為了更好地探索不同土壤肥力條件下的冬小麥SPAD濃度，必須使用聚類神經(jīng)網(wǎng)絡(luò)等差異化模型或地理加權(quán)回歸模型等局部模型。差異化模型的目標(biāo)是根據(jù)光譜角距離（Spectral angle distance, SAD）對(duì)小麥光譜進(jìn)行聚類，以區(qū)分不同肥力下的小麥光譜，然后對(duì)每個(gè)光譜聚類建立回歸模型。

因此，本研究采用K-means作為聚類方法，采用兩種集成學(xué)習(xí)算法作為回歸方法：基于Boosting集成策略的極端梯度提升（XGBoost），以及基于Bagging集成策略的隨機(jī)森林（RF）。比較了非聚類回歸模型（RF、XGBoost）和聚類回歸模型（cluster-RF、cluster-XGBoost）的性能，評(píng)估了土壤有機(jī)質(zhì)（OM）和土壤全氮（TN）對(duì)小麥SPAD遙感估算的影響。

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

安徽省農(nóng)業(yè)科學(xué)院花可可副研究員團(tuán)隊(duì)利用搭載有GaiaSky-mini2-VN高光譜相機(jī)（江蘇雙利合譜公司）的DJI M600無(wú)人機(jī)平臺(tái)，獲取了不同施肥處理下的小麥冠層高光譜影像。高光譜傳感器的波段范圍為400 ~ 1000 nm，光譜分辨率為3.5 nm，包含176個(gè)波段。獲取對(duì)應(yīng)采樣點(diǎn)的SPAD值以及土壤OM和TN。

研究中選擇的聚類方法為K-means，回歸算法為RF和XGBoost。整個(gè)算法流程分為兩步，第一步：使用K-means將小麥光譜劃分為不同的簇。對(duì)176個(gè)波段的光譜反射率進(jìn)行綜合劃分，并將SAD作為算法的判別距離。簇的最大數(shù)量設(shè)置為30，最終選擇的簇?cái)?shù)為2。圖1和2顯示了6個(gè)案例的實(shí)驗(yàn)流程圖，6個(gè)案例的設(shè)置如表1所示。第二步：針對(duì)小麥光譜的每個(gè)聚類，分別建立每個(gè)光譜聚類的回歸模型，并基于RF和XGBoost算法估測(cè)SPAD值。因此，每個(gè)光譜聚類都有相應(yīng)的回歸模型。

圖1 冬小麥SPAD估測(cè)流程圖

圖2 案例5和6的SPAD估測(cè)流程圖。

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

結(jié)論

與RF相比，XGBoost模型在所有田塊中表現(xiàn)出更好的估測(cè)性能（圖3）。RF模型和XGBoost模型的R²值分別為0.763和0.781，XGBoost的RMSE和MAPE均低于RF模型（表2）。RMSE和MAPE評(píng)價(jià)指標(biāo)表明，模型估計(jì)誤差較小且不顯著。

圖3 非聚類RF（a）和非聚類XGBoost（b）模型的預(yù)測(cè)散點(diǎn)圖

表2 非聚類回歸模型的準(zhǔn)確性

從圖4可以看出，無(wú)人機(jī)高光譜技術(shù)可以實(shí)現(xiàn)田塊尺度的SPA監(jiān)測(cè)。XGBoost和RF估測(cè)模型的估測(cè)結(jié)果相似，SPAD空間分布基本相似。XGBoost預(yù)測(cè)的SPAD值比RF分布更廣，并且XGBoost預(yù)測(cè)的SPAD平均值更接近于實(shí)測(cè)值。但兩種模型確定的標(biāo)準(zhǔn)差與真實(shí)數(shù)據(jù)略有差異。因此，有必要引入局部建模，如聚類回歸模型，以進(jìn)一步提高估測(cè)性能。

圖4 基于RF（a）和XGBoost（b）的SPAD空間分布

圖5為聚類后的小麥光譜。將小麥光譜聚為兩類，兩類光譜存在顯著差異，因此聚類結(jié)果是合理的。當(dāng)波長(zhǎng)小于732 nm時(shí)，Cluster1的平均光譜反射率低于Cluster2；當(dāng)波長(zhǎng)大于732 nm時(shí)，Cluster1的平均光譜反射率高于Cluster2。

圖5 高光譜反射率光譜曲線聚類結(jié)果

圖6顯示了聚類回歸模型預(yù)測(cè)的SPAD值與實(shí)測(cè)值之間的散點(diǎn)圖。cluster-XGBoost模型總體上優(yōu)于cluster-RF。cluster-XGBoost的R²（0.925）高于cluster-RF（0.901），RMSE（1.444）比cluster-RF（1.813）低0.369，MAPE比cluster-RF低25%。cluster-XGBoost比cluster-RF預(yù)測(cè)的SPAD范圍更大，并且cluster-XGBoost具有更高的標(biāo)準(zhǔn)差。兩種聚類回歸模型預(yù)測(cè)SPAD的空間分布如圖7所示，兩者表現(xiàn)出相似的分布。

圖6 cluster-RF（a）和cluster-XGBoost（b）模型的預(yù)測(cè)散點(diǎn)圖

圖7 基于cluster-RF（a）和cluster-XGBoost（b）的SPAD空間分布

前人研究表明，長(zhǎng)期施肥處理后土壤養(yǎng)分的變化會(huì)影響小麥的生長(zhǎng)，包括小麥葉綠素含量。如圖8所示，光譜反射率與小麥SPAD高度相關(guān)，176個(gè)波段的相關(guān)系數(shù)絕對(duì)值均大于0.68，說(shuō)明利用光譜信息可以建立小麥SPAD估測(cè)模型。土壤有機(jī)質(zhì)和全氮與小麥SPAD呈顯著正相關(guān)，相關(guān)系數(shù)分別為0.44和0.45。因此，可將土壤有機(jī)質(zhì)和全氮引入小麥SPAD遙感估測(cè)。如表3所示，通過(guò)添加土壤OM和TN，RF和XGBoost模型的精度都得到了提高。對(duì)于RF，R²提高9.95%，RMSE降低4.25%，MAPE降低4.66%；對(duì)于XGBoost，R²提高14.23%，RMSE降低21.74%，MAPE降低13.43%。因此，結(jié)合土壤養(yǎng)分和光譜信息估測(cè)小麥SPAD比單獨(dú)使用光譜信息更準(zhǔn)確。此外，RF和XGBoost模型均能利用光譜和土壤信息準(zhǔn)確檢測(cè)SPAD，且XGBoost模型在施用有機(jī)肥的田中表現(xiàn)優(yōu)于RF模型。

圖8 光譜反射率與小麥SPAD的相關(guān)系數(shù)

表3 RF和XGBoost模型的精度（包括土壤養(yǎng)分）

作者信息

花可可，博士，安徽省農(nóng)業(yè)科學(xué)院土壤肥料研究所副研究員。

主要研究方向：土壤有機(jī)質(zhì)提升過(guò)程與調(diào)控機(jī)制、活性氮污染與阻控技術(shù)。

參考文獻(xiàn)：

Yang, X., Yang, R., Ye, Y., Yuan, Z.R., Wang, D.Z., & Hua, K.K. (2021). Winter wheat SPAD estimation from UAV hyperspectral data using cluster-regression methods. International Journal of Applied Earth Observations and Geoinformation, 105, 102618.