基于隨機森林算法的冬小麥葉片氮含量遙感估算研究

韓玉杰　丁寧　馮海寬　張春蘭

摘要 ? ?精確、快速估算冬小麥葉片氮含量，對冬小麥長勢監測及田間管理指導具有重要的研究意義。為精確反演冬小麥葉片氮含量（leaf nitrogen content ，LNC），本文利用遙感方法，依托不同氮處理水平對冬小麥的影響試驗，基于獲取的高光譜遙感數據和LNC地面實測數據，對比分析光譜指數與隨機森林算法（random forest ，RF）反演冬小麥葉片氮含量的精度和穩健性。結果表明，以敏感波段496 nm和604 nm為自變量，利用隨機森林算法構建的LNC回歸模型精度較光譜指數法有了大幅提高，模型的建模精度為R2=0.922，均方根誤差為0.290;驗證精度為R2=0.873，均方根誤差為0.397，并且相對分析誤差值為2.220，表明將敏感波段與隨機森林算法組合構建的反演模型能較好地反演冬小麥葉片氮含量。

關鍵詞 ? ?冬小麥;葉片氮含量;遙感;隨機森林算法

中圖分類號 ? ?S127;S512.1 ? ? ? ?文獻標識碼 ? ?A

文章編號 ? 1007-5739（2019）13-0001-03 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 開放科學（資源服務）標識碼（OSID）

Remote ?Sensing ?Estimation ?of ?Nitrogen ?Content ?in ?Winter ?Wheat ?Leaves ?Based ?on ?Random ?Forest ?Algorithm

HAN Yu-jie 1 ? ?DING Ning 1 ? ?FENG Hai-kuan 2 ? ?ZHANG Chun-lan 2

（1 College of Surveying and Geo-Informatics，Shandong Jianzhu University，Jinan Shandong 250101;

2 National Engineering Research Center for Information Technology in Agriculture）

Abstract ? ?Precisely and rapidly monitoring of nitrogen content in winter wheat leaves was of great significance for monitoring the growth of winter wheat and guiding field management.By remote sensing，the accuracy and robustness of estimating nitrogen content in winter wheat leaves by spectral index and random forest algorithm were compared and analyzed based on hyperspectral data and LNC field data.The results showed that the precision of the LNC regression model constructed by the RF was much higher than that of the spectral index with sensitive bands 496 nm and 604 nm as independent variables.The accuracy of the model was R2=0.922，RMSE=0.290，the validation accuracy was R2=0.873，RMSE=0.397，and the RPD value was 2.220.It showed that the model constructed by the RF has good precision of wheat LNC.

Key words ? ?winter wheat;leaf nitrogen content（LNC）;remote sensing;random forest algorithm

作物的光合作用、生長發育等離不開礦物質元素的攝入，作為礦物質元素之一的氮元素是作物生長需求量與施用量較多的礦物質元素[1]。

當氮營養元素施用不足時，會嚴重影響作物的產量及品質;相反，當作物施氮過多時，又會對環境造成一定的污染。因此，精確估算冬小麥葉片氮素含量，對冬小麥長勢監測、提高氮肥利用效率以及深入開展精準農業研究具有重要的意義[2]。

冬小麥作為我國重要的糧食作物之一，其氮素營養診斷對于小麥長勢監測及田間技術指導具有重要的研究意義[3]。近年來，新興起的高光譜遙感技術具有諸多優勢，如信息量大、光譜分辨率高、波段連續性強等，已受到作物氮素營養診斷科研人員的青睞[4]。

目前，科研人員對作物氮素診斷開展了大量的研究。如楊寶華等[5]通過篩選敏感波段構建冬小麥氮含量反演模型，能夠較好地估算冠層氮含量;王仁紅等[4]利用高光譜數據建立氮營養指數（NNI）的經驗模型，定性診斷氮營養情況，說明利用光譜反射率估算作物冠層氮含量是具有可行性的;還有一些科研人員利用機器學習算法（如人工神經網絡法、隨機森林算法）估算作物的冠層氮含量[6-7]。

目前，前人在監測作物氮含量方面開展了大量的工作，并取得了豐碩的成果，但是不同的反演方法各有特點，如估算作物氮含量的線性回歸模型比較容易構建，而且還具有可視化的優點，但也存在對大數據集處理需要降維、參數過多容易出現過飽和現象等問題;而隨機森林算法作為一種新型學習方法，具有抗噪聲能力強、不易出現過飽和等優點，在數據分析與挖掘方面具有重要的應用價值，同時也可以很好地解決線性回歸模型的缺陷問題。因此，本文運用統計分析方法，將敏感波段與隨機森林算法組合估算冬小麥葉片氮含量（LNC），進而構建估算冬小麥葉片氮含量（LNC）的最優模型。

1 ? ?材料與方法

1.1 ? ?試驗地概況

冬小麥試驗在北京市小湯山國家精準農業示范基地（東經116.44°，北緯40.18°）開展。試驗區平均海拔約為36 m，年平均氣溫約14.0 ℃。

1.2 ? ?試驗設計

試驗共設置3個因素處理，即小麥品種為2個，分別是京9843、中麥175;氮素水平設置4個處理，分別為施純N 0、91.065、182.130、273.195 kg/hm2;水分供給設置3個水平，即雨養、正常水、2倍正常水。其他田間管理措施與當地正常管理情況一致。

1.3 ? ?數據采集

試驗數據于冬小麥4個關鍵生育期（拔節期、挑旗期、開花期及灌漿期）分別同步獲取冬小麥冠層高光譜反射率數據和LNC數據，其試驗數據的70%樣本用于構建冬小麥LNC反演模型，剩余30%用于模型的驗證。其中，冬小麥地面冠層高光譜數據利用美國ASD Field SpecFR Pro 2500光譜儀獲取，其光譜范圍為350～2 500 nm，采樣間隔為1 nm。冬小麥葉片氮含量采用凱氏定氮法測定。

1.4 ? ?數據分析方法

1.4.1 ? ?組合光譜指數法。由于冠層反射率數據受外界監測環境的影響，有必要篩選敏感波段，然后組合得到一種光譜參數即為光譜指數[8]。由代數形式組合得到的光譜指數能夠消除一些噪聲的影響[9-10]，從而可提高光譜指數對作物生化指標的監測精度。本文組合得到的光譜指數為歸一化光譜指數（Normalized difference spectral index，NDSI）、比值光譜指數（Ratio spectral index，RSI），其公式組合形式如下[11]：

NDSI（x，y）=（1）

RSI（x，y）=（2）

式（1）和式（2）中，x和y分別為一定波段范圍內任意2個波段所對應的光譜反射率。

1.4.2 ? ?隨機森林算法。隨機森林算法（random forest，RF）是由Breiman L[12]提出，該算法是基于集成學習的思想將多棵樹集成的一種算法，其基本單元是決策樹，本質屬于機器學習的一大分支。隨機森林算法是一種靈活性強、實用性高的學習算法，能有效地運用在大數據集上，而且不需要降維，能夠獲得精確的結果。因此，在數據分析與挖掘方面具有極其重要的應用價值。

1.5 ? ?模型評價指標

1.5.1 ? ?決定系數。決定系數（R2）是用來衡量擬合效果的指標，R2在0～1之間，值越大表明模型擬合精度越高，其計算公式為：

1.5.2 ? ?均方根誤差。均方根誤差（RMSE）是衡量模擬值偏離實測值程度大小的指標，其值越小，表示模擬值與實測值吻合性越好，反映構建的反演模型精度較高。其計算公式為：

RMSE=（4）

式（4）中，變量解釋同公式（3）。

1.5.3 ? ?相對分析誤差。相對分析誤差（RDP）是衡量模型質量的指標，當RPD<1.4時，表明構建的反演模型效果差，無法對樣本數據進行預測;當1.4≤RPD<2.0時，表明構建的反演模型效果一般，可以對樣本數據進行粗略評估;當RPD≥2.0時，表明構建的反演模型預測能力較好。其計算公式為：

其中，標準差SD=，式中變量解釋同公式（3）。

2 ? ?結果與分析

2.1 ? ?篩選敏感波段

由于作物的光譜反射率在350～400 nm及1 000 nm以后受水汽影響及環境噪聲影響較大，本文選取數據質量較好的400～1 000 nm的光譜反射率數據進行以下分析。篩選敏感波段的程序在Matlab內運行，通過在400～1 000 nm波段范圍篩選敏感波段，組合得到的NDSI和RSI對LNC的敏感性如圖1所示。結果表明，496 nm與604 nm為敏感波段，其組合得到的光譜指數RSI（R604，R496）和NDSI（R604，R496））與葉片氮含量（LNC）的決定系數R2分別為-0.871、-0.865，表明由敏感波段組合的光譜指數與LNC高度負相關，這是由于歸一化和比值的數學計算形式能夠降低或消除一些噪聲的影響[13]。

2.2 ? ?利用光譜指數法構建冬小麥LNC模型

由篩選的敏感波段（496、604 nm）組合得到的NDSI（R604，R496）、RSI（R604，R496）為自變量，分別構建冬小麥LNC的反演模型。由圖2可知，以RSI（R604，R496）為自變量構建的冬小麥LNC反演模型的決定系數R2為0.759，高于以RSI（R604，R496）為自變量構建的冬小麥LNC反演模型的決定系數R2。說明以比值形式組合的光譜指數RSI（R604，R496）為自變量構建的冬小麥LNC反演模型的估算效果較好。

為了進一步檢驗利用光譜指數法構建的冬小麥LNC估算模型的精度與穩定性，采用相同生育期采集的獨立樣本集對其精度進行驗證。由圖3可知，利用光譜指數RSI（R604，R496）構建的冬小麥LNC估算模型的驗證精度R2為0.793，并且其估算模型的RPD值為1.86，表明以光譜指數RSI（R604，R496）為自變量構建的冬小麥LNC估算模型的預測精度高于以光譜指數NDSI（R604，R496）為自變量構建的冬小麥LNC估算模型的預測精度。

2.3 ? ?基于隨機森林法構建冬小麥LNC模型

為了構建較優的冬小麥LNC估算模型，本文基于與光譜指數法相同的建模樣本集，利用隨機森林算法（RF）構建冬小麥LNC的RF回歸模型，并采用與光譜指數法相同的驗證樣本集對模型精度進行驗證。為了便于這2種方法間的對比與分析，利用隨機森林法構建冬小麥LNC估算模型的輸入因子為2.1篩選的敏感波段496 nm和604 nm處的光譜反射率，因變量為冬小麥LNC，在MATLAB軟件中運行隨機森林算法（RF）回歸試驗代碼。利用隨機森林算法（RF）構建的冬小麥LNC估算模型如圖4（a）所示，其決定系數R2為0.922，均方根誤差為0.290，利用隨機森林算法構建的冬小麥LNC估算模型的預測精度較高。采用獨立的驗證樣本集對利用隨機森林算法構建的冬小麥LNC估算模型的精度及穩定性進行驗證。由圖4（b）可知，該模型的擬合精度較高，可以對冬小麥LNC進行精確算。

與2.2利用光譜指數法構建的冬小麥LNC估算模型相比，利用隨機森林法構建的冬小麥LNC反演模型的估算精度（R2=0.922）和預測精度（R2=0.873）均較高。利用光譜指數RSI（R604，R496）構建的LNC估算模型的RPD值為1.86，而利用RF構建的冬小麥LNC估算模型的RPD值為2.220，表明利用RF構建的冬小麥LNC估算模型反演效果較好，能用于冬小麥LNC的精確估算。

3 ? ?結論與討論

目前，已有較多研究者利用遙感技術開展冬小麥葉片氮含量估算研究工作，但由于線性回歸模型存在受外界生長環境影響較大的問題，尋找高效的冬小麥葉片氮含量反演方法仍需進一步探討。由于監測環境中影響因素較多，本文首先進行了敏感波段的篩選，通過篩選發現496 nm與604 nm波段為敏感波段，而且由敏感波段組合的光譜指數對LNC的敏感性較好，R2分別為-0.871、-0.865;然后分別利用線性回歸和隨機森林法構建冬小麥LNC估算模型，結果表明，隨機森林法構建的LNC估算模型的建模精度（R2=0.922，RMSE=0.290）和驗證精度（R2=0.873，RMSE=0.397）均較高，這是因為線性回歸模型存在對大數據集需要降維、擬合過程中參數過多容易出現過飽和等缺陷;新興起的隨機森林算法具有抗噪能力強、適應能力好、對大數據集處理精度高且不需要降維等多種優勢[14-15]，可以很好地解決線性回歸模型的缺陷。因此，通過篩選敏感波段，利用隨機森林算法構建冬小麥LNC估算模型可以作為氮素田間快速診斷的有效方法。

4 ? ?參考文獻

[1] 田永超，朱艷，姚霞，等.基于光譜信息的作物氮素營養無損監測技術[J].生態學雜志，2007（9）：1454-1463.

[2] 李永梅，張學儉，張立根.水稻氮素高光譜遙感估測模型研究[J].江蘇農業科學，2016，44（8）：435-439.

[3] 尚艷，常慶瑞，劉秀英，等.關中地區小麥冠層光譜與氮素的定量關系[J].西北農林科技大學學報（自然科學版），2016，44（5）：38-44.

[4] 王仁紅，宋曉宇，李振海，等.基于高光譜的冬小麥氮素營養指數估測[J].農業工程學報，2014，30（19）：191-198.

[5] 楊寶華，陳建林，陳林海，等.基于敏感波段的小麥冠層氮含量估測模型[J].農業工程學報，2015，31（22）：176-182.

[6] ZHANG C H，KOVACS J M，WACHOWIAK M P，et al.Relationship between hyperspectral measurements and mangrove leaf nitrogen conce-ntrations[J].Remote Sensing，2013，5（2）：891-908.

[7] YI Q X，HUANG J F，WANG F M，et al.Monitoring rice nitrogen status using hyperspectral reflectance and artificial neural network[J].Environ-mental Science & Technology，2007，41（19）：6770-6775.

[8] 童慶禧，張兵，鄭蘭芬.高光譜遙感[M].北京：高等教育出版社，2006：121.

[9] HUNT E R，DORAISWAMY P C，MCMURTREY J E，et al.A visible band index for remote sensing leaf chlorophyll content at the canopy scale[J].International Journal of Applied Earth Observation and Geoinf-ormation，2013：103-112.

[10] 姜海玲，楊杭，陳小平，等.利用光譜指數反演植被葉綠素含量的精度及穩定性研究[J].光譜學與光譜分析，2015，35（4）：975-981.

[11] INOUE Y，GU？魪RIF M，BARET F，et al.Simple and robust methods for remote sensing of canopy chlorophyll content：A comparative analysis of hyperspectral data for different types of vegetation[J].Plant，Cell and En-vironment，2016，39：2609-2623.

[12] BREIMAN L.Random forests[J].Machine Learning，2001，45（1）：5-32.

[13] RICHARDSON A J，WEIGAND C L.Distinguishing vegetation from soil background information[J].Photogrammetric Engineering and Remote Sensing，1977，43（12）：1541-1552.

[14] 張春蘭.基于綜合長勢參數指標的冬小麥長勢遙感監測[D].西安：西安科技大學，2018.

[15] 韓兆迎，朱西存，房賢一，等.基于SVM與RF的蘋果樹冠LAI高光譜估測[J].光譜學與光譜分析，2016，36（3）：800-805.