摘 要：為實(shí)現(xiàn)枸杞氮素含量的快速無損監(jiān)測(cè)，以“寧杞7號(hào)”為研究對(duì)象，同步測(cè)定枸杞葉片光譜反射率與葉片氮素含量，選用快速傅里葉變換（FFT）對(duì)測(cè)定的枸杞葉片光譜進(jìn)行平滑濾波處理，并獲取原始光譜（OS）。采用一階微分、二階微分及連續(xù)統(tǒng)去除法對(duì)原始光譜進(jìn)行變換處理，獲取枸杞葉片一階微分光譜（FDS）、二階微分光譜（SDS）及連續(xù)統(tǒng)去除光譜（CRS），將原始光譜及3種變式光譜分別與葉片氮素含量進(jìn)行相關(guān)性分析，進(jìn)而篩選出敏感波長(zhǎng)，并構(gòu)建預(yù)測(cè)枸杞葉片氮素含量的隨機(jī)森林回歸模型（RFRM）和多元線性回歸模型（MLRM）。結(jié)果表明，基于3種變換光譜構(gòu)建的RFRM和MLRM的預(yù)測(cè)精度均優(yōu)于基于OS構(gòu)建的模型；其中：基于FDS構(gòu)建的模型預(yù)測(cè)效果最優(yōu)，其次為基于SDS和基于CRS構(gòu)建的模型，基于OS構(gòu)建的模型預(yù)測(cè)精度最差。同時(shí)表明，RFRM的擬合度均優(yōu)于MLRM，基于原始光譜OS、一階微分光譜FDS、二階微分光譜SDS和連續(xù)統(tǒng)去除光譜CRS構(gòu)建的RFRM，同MLRM相比，其模型擬合度分別提高0.258、0.259、0.275和0.291，RMSE分別降低0.044、0.054、0.059、0.076，MAE分別降低0.045、0.043、0.066、0.059?；贔FT+FDS組合方法下篩選的敏感波長(zhǎng)構(gòu)建的RFRM，其建模數(shù)據(jù)集的擬合度、RMSE和MAE分別為0.897、0.071和0.058，檢驗(yàn)數(shù)據(jù)集的決定系數(shù)、RMSE和MAE分別為0.689、0.129和0.102，模型具有良好的精度和穩(wěn)定性，可作為枸杞葉片氮素的高光譜估測(cè)方法。

關(guān)鍵詞：氮素； 枸杞； 隨機(jī)森林回歸模型； 高光譜反射率； 一階微分； 二階微分； 連續(xù)統(tǒng)去除法

中圖分類號(hào)：S181； S567.1+9" " 文獻(xiàn)標(biāo)識(shí)碼：A" " " 文章編號(hào)：1002-204X（2024）06-0048-07

doi：10.3969/j.issn.1002-204x.2024.06.011

Predicting Leaf Nitrogen Content in Wolfberry Trees

Using Hyperspectral Reflectance

Abstract To realize rapid and non-destructive monitoring of nitrogen content in wolfberry tree， \"Ningqi No.7\" wolfberry was selected as the research object to synchronously measure the spectral reflectance and nitrogen content of wolfberry leaves. Fast fourier transform （FFT） was performed to smooth and filter the measured spectra to obtain the original spectra （OS）. Three types of mathematical transformations including first-derivative （FD）， second-derivative （SD） and continuum removal （CR） transformations were performed on the original spectra， and the corresponding spectral datasets including first-derivative spectra （FDS）， second-derivative spectra （SDS） and continuum removal spectra （CRS） were obtained. The correlation analysis between the spectra including OS， FDS， SDS and CRS and nitrogen of wolfberry leaves were performed to select sensitive wavelengths based on the value of correlation coefficients. Random forest regression models（RFRM） and multiple linear regression models （MLRM） were constructed using selected sensitive wavelengths to predict the nitrogen content of wolfberry leaves. The study indicated that the prediction accuracy of RFRM and MLRM constructed using the three types of transformation spectra were better than those constructed using OS. Among them， the models constructed using FDS had the best prediction performance， followed by the models constructed using SDS and the models constructed using CRS， the models constructed using OS had the worst prediction accuracy. Meanwhile， it is shown that the prediction accuracy of RFRM were superior to MLRM. Compared to MLRM， the fitting degree of RFRM constructed using OS， FDS， SDS and CRS increased by 0.258， 0.259， 0.275 and 0.291， the root mean square error （RMSE） decreased by 0.044， 0.054， 0.059 and 0.076， and the mean absolute error （MAE） decreased by 0.045， 0.043， 0.066 and 0.059， respectively. The RFRM model constructed using the sensitive wavelengths selected from FDS had best accuracy and stability， with the determination coefficient， RMSE and MAE of calibration set of 0.897， 0.071 and 0.058， respectively， with the determination coefficient， RMSE and MAE of validation set of 0.689， 0.129 and 0.102， respectively. It can be used as a hyperspectral estimation method for leaf nitrogen content

Key words Nitrogen content; Wolfberry; Random forest regression model; Hyperspectral reflectance; First-derivative; Second-derivative; Continuum removal

氮素是作物生長(zhǎng)發(fā)育必需的營(yíng)養(yǎng)元素，也是農(nóng)田生態(tài)系統(tǒng)變化的主控因子[1]。大量研究表明，植被氮素含量與光合作用密切相關(guān)，氮素通過參與光合作用，影響生態(tài)系統(tǒng)中生物量累積及碳水循環(huán)過程。因此，作物氮素豐缺診斷是實(shí)現(xiàn)農(nóng)業(yè)生產(chǎn)可持續(xù)發(fā)展和生態(tài)環(huán)境良性循環(huán)的重要途徑[2]。傳統(tǒng)作物氮素診斷方法主要是植物組織分析法，該方法需要進(jìn)行破壞性取樣、工作量大、操作繁瑣且推廣應(yīng)用困難[3]。目前快速、準(zhǔn)確地監(jiān)測(cè)作物氮素含量成為現(xiàn)代農(nóng)業(yè)研究的熱點(diǎn)[1]。由于高光譜遙感技術(shù)具有光譜分辨率高、波段連續(xù)性強(qiáng)及光譜信息量大等特點(diǎn)，可以對(duì)地物進(jìn)行微弱光譜差異的定量分析，在無損監(jiān)測(cè)作物營(yíng)養(yǎng)方面具有獨(dú)特的優(yōu)勢(shì)[4-5]，已成為監(jiān)測(cè)作物生長(zhǎng)狀況的一種新手段，在農(nóng)業(yè)養(yǎng)分管理中顯示出極好的應(yīng)用前景。

光譜信息測(cè)定時(shí)經(jīng)常會(huì)受到雜散光、噪聲、設(shè)備性能不穩(wěn)定等因素的干擾，從而對(duì)最終分析的結(jié)果造成影響，因此需要在建模前對(duì)原始光譜進(jìn)行預(yù)處理[6]。許多學(xué)者對(duì)原始光譜進(jìn)行多種預(yù)處理，然后基于預(yù)處理后的變換光譜構(gòu)建模型，根據(jù)建模性能的優(yōu)劣反過來選擇預(yù)處理方法，進(jìn)而實(shí)現(xiàn)對(duì)作物營(yíng)養(yǎng)狀況的良好監(jiān)測(cè)。王麗鳳等[7]應(yīng)用一階導(dǎo)數(shù)、二階導(dǎo)數(shù)及標(biāo)準(zhǔn)正態(tài)變換對(duì)原始光譜進(jìn)行變換，通過構(gòu)建偏最小二乘回歸模型（PLRM）選出最佳光譜預(yù)處理方法，實(shí)現(xiàn)玉米拔節(jié)期葉片氮含量的檢測(cè)。鄭一力等[8]基于原始反射率、一階微分、對(duì)數(shù)一階微分及二階微分的特征波段，構(gòu)建4種竹葉片氮含量估測(cè)模型，最終基于對(duì)數(shù)一階微分光譜的BP神經(jīng)網(wǎng)絡(luò)模型具有較好估測(cè)效果（R2=0.838）。李粉玲等[9]基于原始光譜、連續(xù)統(tǒng)去除光譜及吸收特征參數(shù)，構(gòu)建了預(yù)測(cè)冬小麥葉片全氮含量的一元回歸模型，最終獲取精度最高的模型為基于550～770 nm 波段的吸收峰總面積建立的指數(shù)回歸模型。RUBIO-DELGADO J等[10]采用SG平滑、標(biāo)準(zhǔn)正態(tài)變換，一階微分和二階微分對(duì)原始高光譜反射率進(jìn)行變換處理，在對(duì)變換光譜與葉片氮含量進(jìn)行相關(guān)分析的基礎(chǔ)上，構(gòu)建了PLSR對(duì)西班牙西南部橄欖樹葉片氮素含量進(jìn)行了較好的預(yù)測(cè)。易翔等[11]對(duì)原始光譜進(jìn)行多元散射校正（MSC）、SG平滑算法、變量標(biāo)準(zhǔn)化校正（SNV）和一階導(dǎo)數(shù)（FD）4種變換，采用隨機(jī)蛙跳和連續(xù)投影算法篩選特征波長(zhǎng)，構(gòu)建了估算棉花葉片全氮含量的PLSR，篩選出FD-SPA-PLSR組合方法下構(gòu)建的模型為最佳模型。

枸杞產(chǎn)業(yè)作為寧夏的“六特”產(chǎn)業(yè)之一，其高質(zhì)量發(fā)展備受關(guān)注。目前鮮見采用高光譜反射率監(jiān)測(cè)枸杞葉片氮素等生化參數(shù)的相關(guān)研究報(bào)道。本研究對(duì)枸杞葉片原始光譜進(jìn)行一階微分、二階微分及連續(xù)統(tǒng)去除3種變換處理，在對(duì)原始光譜、變換光譜與氮素定量分析的基礎(chǔ)上，構(gòu)建了估測(cè)葉片氮素含量的隨機(jī)森林回歸模型（random forest regression model， RFRM），通過構(gòu)建多元線性回歸模型對(duì)RFRM進(jìn)行評(píng)估，并采用驗(yàn)證數(shù)據(jù)對(duì)模型穩(wěn)定性進(jìn)行評(píng)價(jià)，以期探討枸杞葉片光譜的最佳預(yù)處理方法并優(yōu)選出估測(cè)枸杞葉片氮素的良好模型，為寧夏枸杞葉片氮素檢測(cè)與評(píng)估提供科學(xué)依據(jù)，為枸杞精準(zhǔn)施肥提供參考。

1 材料和方法

1.1 研究區(qū)概況

研究區(qū)位于寧夏中寧縣恩和鎮(zhèn)枸杞種植園。該地屬溫帶大陸性氣候，光照充足，熱量豐富，晝夜溫差大，干旱少雨，年平均氣溫9.5 ℃，年平均日照時(shí)間2 979.9 h，年平均降水量209 mm。園區(qū)地質(zhì)條件優(yōu)良，地貌平坦，土層厚，土壤質(zhì)地均勻，富含各種礦物質(zhì)和微量元素，屬于寧夏枸杞國(guó)家地理標(biāo)志產(chǎn)品保護(hù)區(qū)范疇。

1.2 試驗(yàn)材料

以“寧杞7號(hào)”枸杞為試驗(yàn)材料，栽培株行距為3 m×1 m，采集樣本時(shí)枸杞樹齡為5年。在開花期至果熟期，隨機(jī)選取長(zhǎng)勢(shì)不同的枸杞樹作為樣株，在樣株不同方位采集結(jié)果枝條上的健康葉片，作為待測(cè)樣本。將采集的鮮葉片樣本迅速裝入自封袋，編號(hào)，放入移動(dòng)式保鮮箱，帶回實(shí)驗(yàn)室。將每個(gè)樣本分成2組，一組用于光譜反射率的測(cè)定，一組用于氮素的測(cè)定。

1.3 數(shù)據(jù)測(cè)定

1.3.1 光譜反射率測(cè)定

光譜反射率測(cè)定在密閉暗室內(nèi)進(jìn)行。測(cè)定設(shè)備為美國(guó)ASD Field Spec Pro便捷式地物光譜儀，光纖視場(chǎng)角為25°，光譜范圍為350～2 500 nm。測(cè)定時(shí)，先將鮮葉片表面灰塵擦拭干凈，然后將葉片呈近似圓形平鋪在黑色絨布上，將光纖垂直對(duì)準(zhǔn)測(cè)定范圍中心位置，光纖距離樣本的垂直高度根據(jù)公式H=L/[2×tan（A/2）]進(jìn)行調(diào)整。其中：H為光纖與樣本的垂直距離，L為測(cè)定范圍直徑，A為光纖視場(chǎng)角。安裝50 W的鹵素?zé)裟M太陽光，且光線以與水平面成45°的恒定角度射入。每個(gè)樣本測(cè)定5條光譜，每次測(cè)定前對(duì)光譜儀進(jìn)行標(biāo)準(zhǔn)白板校正。

1.3.2 葉片氮素測(cè)定

將樣本放入烘箱，調(diào)溫至105 ℃殺青30 min，然后調(diào)溫至70 ℃烘干至恒重，將葉片研磨成粉末，通過0.25 mm的篩網(wǎng)，最后采用凱氏定氮法測(cè)定葉片氮素含量。

1.4 試驗(yàn)方法

1.4.1 光譜預(yù)處理

光譜測(cè)定中，由于光源波動(dòng)、外界干擾及設(shè)備自身等影響，使得光譜受到不同程度的干擾，因此，必須對(duì)光譜進(jìn)行降噪處理，以提高后期光譜分析結(jié)果的準(zhǔn)確性。本研究首先將光譜進(jìn)行比對(duì)，去除異常值，然后對(duì)測(cè)定的光譜進(jìn)行快速傅里葉變換（FFT），得到穩(wěn)定的原始光譜（OS）反射率。然后，采用微分和連續(xù)統(tǒng)去除法對(duì)原始光譜進(jìn)行數(shù)學(xué)變換處理，得到一階微分光譜（FDS）、二階微分光譜（SDS）和連續(xù)統(tǒng)去除光譜（CRS）。將原始光譜及3種變換光譜作為后期定量分析和模型建立的基礎(chǔ)數(shù)據(jù)集。

1.4.2 模型構(gòu)建

隨機(jī)森林（Random forest regression， RF）是以CART決策樹為基礎(chǔ)的一種集成機(jī)器學(xué)習(xí)算法，是將樣本不斷放回并進(jìn)行多次取樣形成訓(xùn)練集，通過決策樹的組合對(duì)預(yù)測(cè)結(jié)果求平均值來進(jìn)行預(yù)測(cè)[12]。本試驗(yàn)將光譜反射率作為輸入，枸杞葉片氮素作為輸出，采用十折交叉驗(yàn)證對(duì)模型進(jìn)行訓(xùn)練，模型的決策樹數(shù)量為100，樹的最大深度為10，分割節(jié)點(diǎn)數(shù)量為2，葉子節(jié)點(diǎn)的最大數(shù)量為50。為了評(píng)價(jià)隨機(jī)森林算法在枸杞葉片氮素預(yù)測(cè)上的表現(xiàn)，采用相同數(shù)據(jù)集，將光譜反射率作為自變量，枸杞葉片氮素作為因變量，構(gòu)建了多元線性回歸模型（multivariable linear regression model， MLRM）。

1.4.3 模型檢驗(yàn)

采用決定系數(shù)（determination coefficient， R2）、均方根誤差（root mean square error， RMSE）和平均絕對(duì)誤差（mean absolute error， MAE）對(duì)模型進(jìn)行獨(dú)立性檢驗(yàn)。

合度越高；RMSE和MAE數(shù)值越小，模型預(yù)測(cè)值與實(shí)測(cè)值間的誤差越小，表明模型的預(yù)測(cè)效果越好。

2 結(jié)果與分析

2.1 氮素正態(tài)檢驗(yàn)分析

將氮素?cái)?shù)據(jù)按照2∶1分為訓(xùn)練集和驗(yàn)證集，對(duì)全體樣本集、訓(xùn)練集和驗(yàn)證集進(jìn)行正態(tài)性檢驗(yàn)（表1）。由表1可見，3組數(shù)據(jù)集的峰度絕對(duì)值均小于10且偏度絕對(duì)值小于3，表明3組數(shù)據(jù)基本可接受為正態(tài)分布。

2.2 光譜反射率與氮素相關(guān)性分析

采用Pearson相關(guān)性分析，將原始光譜及變換光譜與氮素分別進(jìn)行定量分析（圖1）。可見，原始光譜與氮素的相關(guān)系數(shù)在-0.474～0.083之間變化。原始光譜729～860 nm波段與氮素在0.01水平呈顯著負(fù)相關(guān)，達(dá)到0.01水平顯著相關(guān)的波段有132個(gè)；722～728 nm和861～906 nm波段與氮素在0.05水平呈現(xiàn)顯著相關(guān)，達(dá)到0.05水平顯著相關(guān)的波段有54個(gè)；通過0.01和0.05水平顯著性檢驗(yàn)的波長(zhǎng)共有186個(gè)，主要集中在紅邊及近紅外波段。同原始光譜比，一階微分光譜和二階微分光譜與氮素的相關(guān)系數(shù)曲線波動(dòng)很大，其中一階微分光譜與氮素的相關(guān)系數(shù)在-0.550～0.571之間變化，與氮素在0.01水平顯著相關(guān)的波段有309個(gè)，在0.05水平顯著相關(guān)的有87個(gè)波長(zhǎng)，共有396個(gè)波長(zhǎng)通過0.01和0.05水平顯著性檢驗(yàn)。二階微分光譜與氮素的相關(guān)系數(shù)在-0.588～0.553之間，達(dá)到0.01水平顯著相關(guān)的波段有121個(gè)，達(dá)到0.05水平顯著相關(guān)的光譜有66個(gè)，共計(jì)187個(gè)波長(zhǎng)。一階微分光譜和二階微分光譜與氮素通過相關(guān)性檢驗(yàn)的波長(zhǎng)在可見光波段與近紅外波段均有分布。連續(xù)統(tǒng)去除光譜達(dá)到0.01水平顯著相關(guān)的波長(zhǎng)有437～530 nm、555～745 nm、834～848 nm，共有378個(gè)波長(zhǎng)；與氮素在0.05水平相關(guān)的波長(zhǎng)有93個(gè)，合計(jì)471個(gè)波段，主要分布于近紅外波段，其相關(guān)系數(shù)在-0.577～0.510之間變化。

根據(jù)相關(guān)系數(shù)絕對(duì)值最大原則，從4類光譜數(shù)據(jù)集中選取10個(gè)波長(zhǎng)作為敏感波長(zhǎng)（表2），表2中敏感波長(zhǎng)按照相關(guān)系數(shù)絕對(duì)值從大到小的順序排列?？梢?，原始光譜敏感波長(zhǎng)為近紅外768～777 nm波段內(nèi)連續(xù)分布的10個(gè)波長(zhǎng)。連續(xù)統(tǒng)去除光譜敏感波長(zhǎng)集中于777～938 nm的近紅外波段，同原始光譜相比，其敏感波長(zhǎng)的分布相對(duì)分散，所包含的光譜信息豐富。一階微分光譜的敏感波長(zhǎng)553 nm，位于550 nm附近的綠峰波段；745～748 nm的4個(gè)敏感波長(zhǎng)位于紅邊波段；825 nm、826nm、916 nm、936 nm及959 nm位于近紅外波段；可見一階微分光譜的敏感波長(zhǎng)分布最分散，包含的光譜信息最豐富。二階微分光譜的敏感波長(zhǎng)包含702～706 nm紅邊波段內(nèi)的連續(xù)5個(gè)波長(zhǎng)和747～749 nm紅邊波段內(nèi)的3個(gè)連續(xù)波長(zhǎng)，以及近紅外波段的837 nm和847 nm 2個(gè)波長(zhǎng)。

2.3 模型構(gòu)建

基于表2中的敏感波長(zhǎng)，構(gòu)建了預(yù)測(cè)枸杞葉片氮素的RFRM（表3）?？梢?，基于光譜變換構(gòu)建的模型預(yù)測(cè)效果均優(yōu)于基于原始光譜構(gòu)建的模型，其中，基于FDS構(gòu)建的RFRM擬合度最高（R2=0.897），同基于OS、SDS及CRS建立的RFRM比，其擬合度分別提高了0.137、0.027和0.034；RMSE和MAE分別降低了0.048、0.006和0.004，0.035、0.004和0.002。其次為基于SDS和基于CRS的RFRM，其擬合度分別為0.870和0.863，同OS構(gòu)建的RFRM比，擬合度分別提高了0.11和0.103；RMSE和MAE分別減小了0.042和0.044，0.031和0.033?；贠S構(gòu)建的RFRM預(yù)測(cè)精度最差（R2=0.760）。

為了評(píng)價(jià)隨機(jī)森林算法在枸杞葉片氮素預(yù)測(cè)上的表現(xiàn)，以表2中的敏感波長(zhǎng)作為自變量，構(gòu)建了MLRM（表3），可見，MLRM的擬合度高低依次為R2（FDS+MLRM）、R2（SDS+MLRM）、R2（CRS+MLRM）、R2（OS+MLRM），與RFRM表現(xiàn)規(guī)律一致。其中：基于FDS構(gòu)建的MLRM預(yù)測(cè)效果最優(yōu)，同其余3個(gè)MLRM比，其擬合度分別提高了0.136、0.043和0.066；RMSE和MAE分別降低了0.038、0.011和0.026，0.037、0.01和0.019；同MLRM相比，基于OS、FDS、SDS和CRS構(gòu)建的4種RFRM的擬合度分別提高了0.258、0.259、0.275和0.291，RMSE分別減小了0.044、0.054、0.059、0.076，MAE分別減小了0.045、0.043、0.066、0.059。

2.4 模型檢驗(yàn)

利用驗(yàn)證集對(duì)擬合度超過0.7的4個(gè)RFRM進(jìn)行預(yù)測(cè)精度檢驗(yàn)（表4），可見，模型的決定系數(shù)由大到小依次為R2（FDS+RFRM）、R2（SDS+RFRM）、R2（CRS+RFRM）、R2（OS+RFRM），RMSE由小到大依次為RMSE（FDS+RFRM）、RMSE（SDS+RFRM）、RMSE（CRS+RFRM）、RMSE（OS+RFRM），4個(gè)模型的MAE排序與RMSE一致。

通過對(duì)模型檢驗(yàn)可見，基于FDS構(gòu)建的RFRM預(yù)測(cè)精度及穩(wěn)定性最優(yōu)，其檢驗(yàn)R2、RMSE和MAE分別為0.689、0.129和0.102。基于OS構(gòu)建的RF回歸模型預(yù)測(cè)能力最差，其檢驗(yàn)R2僅為0.196，說明該模型只能解釋樣本變化的19.6%。

3 結(jié)論與討論

本研究從葉片尺度開展枸杞氮素含量估測(cè)研究，為了能夠準(zhǔn)確反映枸杞植株氮素狀況，以枸杞當(dāng)年生的結(jié)果枝條上的鮮葉片為研究對(duì)象，同步測(cè)定葉片光譜反射率與全氮含量，在定量分析光譜與氮素的定量關(guān)系的基礎(chǔ)上，建立了枸杞葉片氮素預(yù)測(cè)回歸模型，為枸杞精準(zhǔn)施肥提供科學(xué)參考。

已有研究表明，光譜變換可以突出光譜特征并提取植被生物化學(xué)信息[14-16]，已成為過濾噪聲、提高信噪比和提升建模精度的關(guān)鍵手段。快速傅里葉變換（FFT）[17]可以將信號(hào)從時(shí)域變換到頻域，使噪聲信號(hào)在頻域中更容易檢測(cè)和消除，從而提高去噪的效率。本研究采用FFT對(duì)測(cè)定的光譜進(jìn)行濾波平滑處理，獲取穩(wěn)定的原始光譜（OS），然后采用數(shù)學(xué)變換獲取3種變換光譜（FDS、SDS和CRS），并構(gòu)建了估測(cè)枸杞葉片氮素的隨機(jī)森林回歸模型和多元線性回歸模型。本研究表明基于3種變換光譜構(gòu)建的RFRM和MLRM的預(yù)測(cè)效果均優(yōu)于基于原始光譜構(gòu)建的模型，可見對(duì)枸杞葉片光譜進(jìn)行數(shù)學(xué)變換能提升模型預(yù)測(cè)精度，這與劉昌華等[18]得出的結(jié)果相一致?；贔DS的模型預(yù)測(cè)效果最優(yōu)，基于SDS的模型次之，第三為基于CRS的模型。這是由于采用Pearson相關(guān)性分析篩選枸杞葉片氮素的光譜敏感波長(zhǎng)時(shí)，F(xiàn)DS具有較好的效果，其敏感波長(zhǎng)分布于綠峰波段、紅邊波段及近紅外波段，涵蓋光譜信息更豐富。SDS的敏感波長(zhǎng)包含紅邊及近紅外波段。同F(xiàn)DS和SDS相比，CRS的敏感波長(zhǎng)分布相對(duì)集中。2種建模方法中，RFRM的預(yù)測(cè)精度明顯優(yōu)于MLRM。這是由于多元線性回歸模型在解決非線性問題時(shí)擬合能力較弱，而RF算法具有較強(qiáng)的全局搜索能力，能一定程度上改善信息缺失問題[19]，同時(shí)，該算法通過引入了兩個(gè)隨機(jī)性，使其具有較強(qiáng)的抗噪聲能力，也不易產(chǎn)生過度擬合現(xiàn)象[20]，在解決非線性問題時(shí)具有相對(duì)較高的精度。

今后可依托無人機(jī)平臺(tái)，采集枸杞高光譜圖譜信息，從枸杞冠層尺度，將圖譜、紋理與枸杞冠層結(jié)構(gòu)信息相結(jié)合，嘗試綜合更多的數(shù)學(xué)變換、光譜分析及建模方法，實(shí)現(xiàn)光譜數(shù)據(jù)的充分利用，進(jìn)一步提高枸杞氮素監(jiān)測(cè)精度，助力寧夏枸杞產(chǎn)業(yè)高質(zhì)量發(fā)展。

參考文獻(xiàn)：

[1] 李永梅，張立根. 作物氮素高光譜遙感監(jiān)測(cè)研究[J]. 寧夏農(nóng)林科技，2020，61（4）：55-59.

[2] DIACONO M， RUBINO P， MONTEMURRO F. Precision nitrogen management of wheat. A review[J]. Agronomy for Sustainable Development， 2013，33：219-241.

[3] 陳楊，樊明壽，李斐，等. 氮素營(yíng)養(yǎng)診斷技術(shù)的發(fā)展及其在馬鈴薯生產(chǎn)中的應(yīng)用[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2009，25（3）：66-71.

[4] 陳云浩，王思佳，趙逸飛，等. 農(nóng)業(yè)高光譜遙感研究進(jìn)展及發(fā)展趨勢(shì)[J]. 地理與地理信息科學(xué)，2019，35（5）：1-8.

[5] 李永梅，張立根，張學(xué)儉. 水稻葉片高光譜響應(yīng)特征及氮素估算[J]. 江蘇農(nóng)業(yè)科學(xué)，2017，45（23）：210-213.

[6] 第五鵬瑤，卞?；郏踝朔?，等. 光譜預(yù)處理方法選擇研究[J]. 光譜學(xué)與光譜分析，2019，39（9）：2800-2806.

[7] 王麗鳳，張長(zhǎng)利，趙越，等. 高光譜成像技術(shù)的玉米葉片氮含量檢測(cè)模型[J]. 農(nóng)機(jī)化研究，2017，39（11）：140-147.

[8] 鄭一力，趙燕東，董瑋，等. 竹葉片氮含量高光譜估測(cè)方法對(duì)比研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2018，49（S1）：393-400.

[9] 李粉玲，常慶瑞. 基于連續(xù)統(tǒng)去除法的冬小麥葉片全氮含量估算[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2017，48（07）：174-179.

[10] RUBIO-DELGADO J， Pérez C J， Vega-Rodríguez M A. Predicting leaf nitrogen content in olive trees using hyperspectral data for precision agriculture[J]. Precision Agriculture， 2020，22（1）：1-21.

[11] 易翔，呂新，張立福，等. 基于RF和SPA的無人機(jī)高光譜估算棉花葉片全氮含量[J]. 作物雜志，2023（2）：245-252.

[12] 陳曉凱，李粉玲，王玉娜，等. 無人機(jī)高光譜遙感估算冬小麥葉面積指數(shù)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36（22）：40-49.

[13] 谷振宇，陳聰，鄭家佳，等. 基于時(shí)空?qǐng)D卷積循環(huán)神經(jīng)網(wǎng)絡(luò)的交通流預(yù)測(cè)[J]. 控制與決策，2022，37（3）：645-653.

[14] 田安紅，熊黑鋼，趙俊三，等. 分?jǐn)?shù)階微分對(duì)鹽漬土野外光譜預(yù)處理精度提升的機(jī)理分析[J]. 光譜學(xué)與光譜分析，2019，39（8）：2495-2500.

[15] MOROS J. Vibrational spectroscopy provides a green tool for multi-component analysis[J]. Trends in Analytical Chemistry， 2009，29（7）：578–591.

[16] LUO S， HE Y， LI Q， et al. Assessment of unified models for estimating potato leaf area index under water stress conditions across ground-based hyperspectral data[J]. Journal of Applied Remote Sensing， 2020，14（1）：014517.

[17] XU P， CAO Y， XUE L， et al. Study on sparse representation and measurement matrices of compressive sensing of plant hyperspectral data for retrieving plant physiological and biochemical parameters[J]. Biosystems Engineering，2019，182：38-53.

[18] 劉昌華，方征，陳志超，等. ASD Field Spec3野外便攜式高光譜儀診斷冬小麥氮營(yíng)養(yǎng)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34（19）：162-169.

[19] 張春蘭，楊貴軍，李賀麗，等. 基于隨機(jī)森林算法的冬小麥葉面積指數(shù)遙感反演研究[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2018，51（5）：855-867.

[20] 高小梅，李燕麗，熊勤學(xué)，等. 受漬冬小麥不同葉位葉片的SPAD高光譜估算[J]. 生態(tài)學(xué)雜志，2024，43（3）：733-740.