許童羽 楊佳欣 白駒馳 金忠煜 郭忠輝 于豐華

(1.沈陽農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院, 沈陽 110866; 2.遼寧省農(nóng)業(yè)信息化工程技術(shù)研究中心, 沈陽 110866)

0 引言

氮素是植物生長必需的營養(yǎng)物質(zhì)之一,為植物光合作用提供重要支持,對植物生長、產(chǎn)量和品質(zhì)有著重大影響[1]。及時、準(zhǔn)確地判斷水稻氮素營養(yǎng)狀況,進(jìn)而做出合理的施肥決策,對實現(xiàn)少氮肥高產(chǎn)量的目的尤為重要[2]。目前,無人機(jī)光譜遙感技術(shù)以其覆蓋面積大、速度快、無損壞等優(yōu)點被廣泛應(yīng)用于水稻、玉米、大豆、小麥等主要農(nóng)作物的生物量檢測中[3-5]。通常,水稻氮素含量可以從葉片和冠層兩種尺度進(jìn)行估算。冠層光譜參數(shù)測量作為作物氮素含量定量檢測的關(guān)鍵步驟,對監(jiān)測作物生長發(fā)育至關(guān)重要。光譜遙感技術(shù)通過獲取冠層的光譜反射特性間接反映作物葉片的營養(yǎng)成分,為作物營養(yǎng)診斷和合理施肥提供參考[6-7]。大量研究表明,水稻冠層光譜反射率與葉片含氮量呈顯著正相關(guān),且歸一化差值植被指數(shù)(NDVI)、綠波段歸一化植被指數(shù)(GNVI)、比值植被指數(shù)(RVI)等植被指數(shù)都是估測水稻葉片氮素含量的有效指標(biāo)[8-9]。傅友強(qiáng)等[10]提出美香占2號(MXZ2H)的葉片氮素含量與NDVI、比值植被指數(shù)(RVI)、差值植被指數(shù)(DVI)等呈正相關(guān);吉豐優(yōu)1002(JFY1002) 的葉片氮素含量與 RVI、EVI、改進(jìn)非線性植被指數(shù)(MNLI)、 DVI、NDVI呈正相關(guān)?？梢?利用植被指數(shù)進(jìn)行水稻冠層氮素含量的反演具有一定可行性。

上述研究中,利用無人機(jī)平臺進(jìn)行遙感監(jiān)測時,大多都基于垂直方向進(jìn)行反射強(qiáng)度測量,因此只能利用來自冠層頂部的光譜信息進(jìn)行氮素含量的反演。而對分蘗期的水稻,垂直方向上的測量數(shù)據(jù)中包含的水稻信息較少,且由于水稻冠層結(jié)構(gòu)復(fù)雜,受水體、葉傾角、太陽高度角、傳感器觀測角等[11]多種因素的影響,獲取的高光譜影像數(shù)據(jù)中包含鏡面反射分量和漫反射分量。但鏡面反射信息與植物葉片內(nèi)部的物質(zhì)組成無關(guān)[12]。因此,進(jìn)行氮素含量反演時,有必要考慮鏡面反射對反演精度的影響。LI等[13]基于偏振多光譜成像系統(tǒng),提出了一種新的鏡面反射去除植被指數(shù)(NSRVI)來消除植物的鏡面反射干擾,提高植物健康狀況檢測的準(zhǔn)確性。該方法通過構(gòu)建新的植被指數(shù)來消除鏡面反射對結(jié)果的影響。目前,在植物應(yīng)用研究方面,通常使用偏振測量技術(shù)對鏡面反射分量和漫反射分量進(jìn)行區(qū)分[14]。MICHIO等[15]采用不同觀測天頂角對小麥冠層進(jìn)行了偏振測量,發(fā)現(xiàn)采用偏振度表征鏡面反射分量時,能夠更好地描述小麥抽穗期開始時間。呂云峰[16]對抽穗前不同生長時期的玉米冠層進(jìn)行了垂直方向上的高光譜偏振觀測,發(fā)現(xiàn)偏振光在總反射光中的比例高達(dá)10%。國內(nèi)外學(xué)者將偏振遙感應(yīng)用于植物研究中,并取得了一定的成果,但大部分研究是利用偏振測量進(jìn)行角度效應(yīng)分析,缺乏偏振測量在水稻營養(yǎng)元素反演方面的應(yīng)用。

綜上,本文在進(jìn)行水稻冠層氮素含量反演時,利用無人機(jī)偏振遙感技術(shù),將鏡面反射信息從高光譜中剝離,結(jié)合多角度觀測,分析不同觀測天頂角下,去除鏡面反射信息前后的光譜信息與水稻冠層氮素含量的相關(guān)性,進(jìn)而得到無人機(jī)偏振測量的最優(yōu)觀測角度?；谧顑?yōu)觀測天頂角下偏振光譜數(shù)據(jù)(即去除鏡面反射信息的高光譜數(shù)據(jù)),構(gòu)建適用于偏振光譜數(shù)據(jù)的植被指數(shù)——偏振光譜植被指數(shù)(Polarisation spectrum vegetation index,PSVI),以期構(gòu)建一種高效、精準(zhǔn)的植被指數(shù),為水稻冠層氮素含量快速檢測提供監(jiān)測方法。

1 材料與方法

1.1 研究區(qū)概況

本研究試驗區(qū)位于遼寧省海城市耿莊鎮(zhèn)沈陽農(nóng)業(yè)大學(xué)精準(zhǔn)農(nóng)業(yè)航空科研基地(40°58′45.39″N,122°43′47.006 4″E),選用水稻品種為“沈農(nóng)9816”。在水稻試驗區(qū)(圖1)共設(shè)置5個氮肥梯度,即N0(0 kg/hm2)、N1(75 kg/hm2)、N2(150 kg/hm2)、N3(225 kg/hm2)和N4(300 kg/hm2),氮肥基追比為5∶3∶2。各個試驗小區(qū)之間采用田埂進(jìn)行劃分。每個梯度設(shè)置3個重復(fù)小區(qū),共設(shè)3×5=15個小區(qū)。每個小區(qū)面積為5 m×8 m=40 m2,除氮肥梯度外,每小區(qū)田間管理一致,磷鉀肥采用當(dāng)?shù)貥?biāo)準(zhǔn)進(jìn)行施用,磷肥標(biāo)準(zhǔn)施用量為144 kg/hm2,鉀肥標(biāo)準(zhǔn)施用量為192 kg/hm2,基追比為1∶1,其余田間管理同常規(guī)高產(chǎn)管理一致。田間采樣由分蘗期至抽穗期,采樣間隔7 d,每次采樣在各試驗小區(qū)選取具有代表性的3穴水稻進(jìn)行水稻氮素含量與地上干物質(zhì)量的獲取,最后結(jié)果取平均值作為該小區(qū)的水稻氮素含量與地上干物質(zhì)量。

圖1 研究區(qū)概況Fig.1 Overview of study area

1.2 數(shù)據(jù)獲取

1.2.1水稻冠層偏振光譜反射率獲取

采用耐司鍍膜偏振鏡,卡口尺寸為37 mm,厚度為3 mm。無人機(jī)是深圳大疆創(chuàng)新公司制造的 M600 PRO型六旋翼無人機(jī),搭載四川雙利合譜公司的GaiaSky-mini型高光譜成像系統(tǒng),將偏振片置于成像系統(tǒng)鏡頭前,偏振光譜波段為400～1 000 nm,分辨率為3 nm,單幅影像記錄時間為15 s,幀速率為162 f/s,設(shè)定無人機(jī)飛行高度為100 m,選取9個觀測角度(-60°、-45°、-30°、-15°、0°、15°、30°、45°和60°)進(jìn)行高光譜反射測量和偏振反射測量(其中“-”對應(yīng)后向散射方向,“0°”對應(yīng)垂直方向,“+”對應(yīng)前向散射方向),如圖2所示,2次測量在0.5 h內(nèi)完成。正午條件下,光照效果較好且風(fēng)速小,既可以確保光照的穩(wěn)定,又可以實現(xiàn)無人機(jī)平穩(wěn)飛行,故將無人機(jī)光譜數(shù)據(jù)采集時間確定在10:00—14:00。用工具軟件ENVI 5.3對獲取的偏振光譜影像進(jìn)行提取,通過計算每小區(qū)標(biāo)簽附近的平均光譜反射率來表征水稻冠層的偏振光譜反射率信息。

圖2 水稻冠層偏振光譜多角度觀測示意圖Fig.2 Multi-angle observation diagram of rice canopy polarization spectrum

1.2.2水稻冠層氮素含量測定

在每個試驗小區(qū)的取樣點對水稻進(jìn)行整穴的破壞性取樣,帶回實驗室后,將每穴水稻的所有新鮮葉片剪下置于干燥箱中,在120℃下殺青60 min后, 在80℃下干燥,達(dá)到恒定質(zhì)量后,稱量并粉碎。采用凱氏定氮法利用研磨好的粉末檢測葉片的氮素含量(mg/g)。最終共得到40組水稻冠層氮素含量的有效數(shù)據(jù)。之后,使用Kennard-Stone算法[17]按照訓(xùn)練集與驗證集比例7∶3將樣本劃分。氮素含量統(tǒng)計結(jié)果如表1所示。

表1 訓(xùn)練集與驗證集水稻葉片氮素含量統(tǒng)計結(jié)果Tab.1 Nitrogen content statistics of japonica rice leaves in training and validation sets

1.3 最優(yōu)觀測天頂角選擇

植被指數(shù)是對地表植被狀況簡單、有效和經(jīng)驗的度量,其在全球和區(qū)域土地覆蓋、植被分類和環(huán)境變化等方面都有廣泛應(yīng)用。此外,植被指數(shù)還可以轉(zhuǎn)換成葉冠生物的物理學(xué)參數(shù)。本研究結(jié)合現(xiàn)有文獻(xiàn),運用6種與水稻冠層氮素含量估算相關(guān)的植被指數(shù),即歸一化差值植被指數(shù)[18](Normalized difference vegetation index,NDVI)、比值植被指數(shù)[19](RVI)、差值植被指數(shù)[20](Difference vegetation index,DSI)、寬動態(tài)范圍植被指數(shù)[21](Wide dynamic range vegetation index,WDRVI)、再歸一化植被指數(shù)[22](Renormalized vegetation index,RDVI)以及轉(zhuǎn)換型植被指數(shù)[23](Conversion vegetation index,TVI)?；诓煌^測天頂角下的偏振光譜數(shù)據(jù),通過Matlab平臺對各植被指數(shù)針對的波段范圍進(jìn)行遍歷,找到各植被指數(shù)的最優(yōu)波長組合,在此基礎(chǔ)上構(gòu)建植被指數(shù),計算去除偏振信息前后,各觀測角度下,植被指數(shù)與對應(yīng)水稻冠層氮素含量的相關(guān)系數(shù)R。根據(jù)各天頂角下相關(guān)系數(shù)R的結(jié)果,對比得出研究所需要的最優(yōu)觀測角度。

1.4 水稻冠層偏振光譜反射率特征波長選擇

由于采集到的偏振光譜波段具有較高維度,存在大量無用和干擾波段,而構(gòu)建植被指數(shù)通常是對幾個特征波長進(jìn)行一定的組合而形成的一種數(shù)學(xué)表達(dá)式。因此,需要基于最優(yōu)觀測天頂角下的偏振光譜反射率數(shù)據(jù),提取與水稻氮素相關(guān)的特征波長,從而通過構(gòu)建植被指數(shù)反演水稻氮素含量。

本研究采用連續(xù)投影法[24](Successive projections algorithm,SPA)對最優(yōu)觀測角度下的偏振光譜反射率信息進(jìn)行特征波長的提取。連續(xù)投影算法是一種使矢量空間共線性最小化的前向波段選擇算法,它是將維度較高的光譜信息進(jìn)行初步壓縮,把篩選出的波段按照其對氮素含量的貢獻(xiàn)值做進(jìn)一步選擇,從而剔除不敏感波長,降低復(fù)雜度[25]。連續(xù)投影算法的優(yōu)勢在于可以提取全波段的幾個特征波長,有效降低變量間的共線性,建立冗余信息最少的波段組合,實現(xiàn)數(shù)據(jù)降維,大大減少建模所用波段的數(shù)量。本研究采用Matlab平臺實現(xiàn)連續(xù)投影算法提取特征波長。

1.5 水稻冠層氮素植被指數(shù)構(gòu)建

目前,大多數(shù)關(guān)于水稻氮素含量植被指數(shù)的構(gòu)建并未考慮鏡面反射信息的影響。而冠層光譜探測過程中,受觀測角度、植被冠層結(jié)構(gòu)等因素的影響,傳感器接收到的能量不僅包含進(jìn)入植物內(nèi)部多次散射后的能量(包含植物內(nèi)部的生物化學(xué)信息),還包含未進(jìn)入植物內(nèi)部器官,在植物葉片表面發(fā)生鏡面反射的能量(只與葉片表面結(jié)構(gòu)的光學(xué)特性相關(guān))[12]。

本研究利用偏振遙感手段去除光譜信息中的鏡面反射信息,采用1.4節(jié)中從最優(yōu)觀測天頂角下去除偏振信息的光譜數(shù)據(jù)中提取得到的特征波長來構(gòu)建植被指數(shù),構(gòu)建植被指數(shù)的具體步驟如下:

(1)通過連續(xù)投影法降維得到最優(yōu)觀測天頂角下水稻冠層偏振光譜的特征波長反射率R1、R2、…、Rn,n為光譜波段。

(3)由于氮素含量較大,均大于1 mg/g,而Qx較小,均小于1,為進(jìn)一步增強(qiáng)二者之間的相關(guān)性,故在Qx的基礎(chǔ)上,均加1得Bx。

(4)為提高回歸模型擬合優(yōu)度,在上述步驟的基礎(chǔ)上, 采用“平方變換”方法[26],最終得到最優(yōu)觀測天頂角下基于偏振光譜的植被指數(shù)(PSVI)

采用線性回歸法構(gòu)建水稻冠層氮素含量的反演模型,以均方根誤差(RMSE)和模型決定系數(shù)(R2)作為氮素反演模型的評價標(biāo)準(zhǔn)。

2 結(jié)果與分析

2.1 最優(yōu)觀測天頂角選擇結(jié)果

基于去除偏振信息前后的各觀測天頂角測量得到的光譜數(shù)據(jù),計算6種植被指數(shù)與水稻冠層氮素含量的相關(guān)性,為了便于對比去除偏振信息前后,各觀測角度下植被指數(shù)與氮素含量相關(guān)性的提升效果,本研究采用相關(guān)系數(shù)提升比率進(jìn)行直觀展示,如圖3所示。顯然,加入偏振片會對不同觀測天頂角下的植被指數(shù)與氮素含量的相關(guān)性造成不同的影響。觀測天頂角為-60°、0°、30°與60°時,各植被指數(shù)與水稻冠層氮素含量的相關(guān)性并沒有得到提升;觀測天頂角為-45°～-15°、15°與45°時,各植被指數(shù)與水稻冠層氮素含量的相關(guān)性有明顯提升,且在-15°時,提升效果最好。

圖3 去除偏振信息后不同觀測天頂角下各植被指數(shù)與氮素含量相關(guān)性提升結(jié)果Fig.3 Results of correlation enhancement between each vegetation index and nitrogen content at different observed zenith angles after removal of polarization information

為進(jìn)一步確定最優(yōu)觀測天頂角,基于各觀測天頂角測量得到的偏振光譜數(shù)據(jù)獲得的6種植被指數(shù)與水稻冠層氮素含量相關(guān)性如表2所示(表中ρ為光譜反射率,下角數(shù)字表示特征波長)。由表2可知,去除偏振信息后,不同觀測天頂角下的植被指數(shù)與水稻冠層氮素含量的相關(guān)性明顯存在差異。觀測天頂角的變化會對偏振光譜反射率所構(gòu)建的植被指數(shù)與水稻冠層氮素含量的相關(guān)性產(chǎn)生影響。在水稻冠層多角度偏振光譜實測數(shù)據(jù)中,當(dāng)垂直觀測時,植被指數(shù)NDVI、RVI、DSI、WDRVI、RDVI、TVI與水稻氮素含量的相關(guān)性最低,且相關(guān)系數(shù)最大為0.677;當(dāng)觀測天頂角變大時,即從垂直觀測到前向觀測15°和后向觀測15°時,各植被指數(shù)與氮素含量的相關(guān)性都有所提升,后向觀測15°時的相關(guān)性達(dá)到最高,且相關(guān)系數(shù)最小為0.802;當(dāng)觀測天頂角繼續(xù)增大時,即從前后向15°到前后向30°時,各植被指數(shù)與氮素含量的相關(guān)性又明顯下降,但仍高于垂直觀測時的相關(guān)性;當(dāng)觀測天頂角持續(xù)增大到前后60°時,各植被指數(shù)與氮素的相關(guān)性波動不大,但較前后30°時有所提升。綜上,-15°為最優(yōu)觀測天頂角。因此,本研究選擇觀測天頂角為-15°時的偏振光譜數(shù)據(jù)進(jìn)行后續(xù)植被指數(shù)的構(gòu)建。

表2 基于偏振光譜數(shù)據(jù)不同觀測天頂角下植被指數(shù)與水稻冠層氮素含量的相關(guān)系數(shù)Tab.2 Correlation coefficients between vegetation index and nitrogen content of rice canopy at different observed zenith angles based on polarization spectral data

2.2 偏振光譜特征波長選擇結(jié)果

水稻冠層偏振光譜反射率對水稻氮素的定量反演十分重要,偏振光譜的本質(zhì)是去除鏡面反射信息的高光譜,而連續(xù)波段的高光譜反射率信息中存在大量的冗余信息。因此,本研究應(yīng)用連續(xù)投影算法對全波段的水稻冠層偏振光譜進(jìn)行特征波長篩選,并利用校正集進(jìn)行內(nèi)部交叉驗證,當(dāng)內(nèi)部交叉驗證時,交叉驗證均方根誤差(Root mean square error of cross validation,RMSECV)呈現(xiàn)出持續(xù)下降趨勢且RMSECV達(dá)到最低,此時的RMSECV即為最優(yōu)波長集合包含的波長數(shù)量,也就是該時期的最優(yōu)子集波長數(shù)量。根據(jù)驗證結(jié)果的RMSECV篩選出6個水稻冠層氮素含量高光譜特征波長(圖4):500、566、663、691、736、763 nm。

圖4 水稻冠層偏振光譜特征波段SPA篩選結(jié)果Fig.4 Successive projections algorithm (SPA) screening results of hyperspectral characteristic bands of rice canopy

2.3 偏振光譜植被指數(shù)構(gòu)建結(jié)果

本研究將連續(xù)投影法篩選得到的6個特征波長的反射率兩兩組合,分別進(jìn)行PSVI指數(shù)運算,得到PSVI指數(shù)與水稻冠層氮素含量相關(guān)系數(shù)。圖5中最高時的組合為500、566 nm,相關(guān)系數(shù)絕對值為0.850。

圖5 PSVI指數(shù)與水稻冠層氮素含量的相關(guān)性Fig.5 Correlation between PSVI index and nitrogen content of rice canopy

2.4 水稻氮素含量反演結(jié)果

將本研究所構(gòu)建的PSVI新型植被指數(shù)作為模型的輸入,利用非線性關(guān)系建立水稻冠層氮素含量的反演模型,模型訓(xùn)練集R2為0.783 8, RMSE為 0.428 mg/g;驗證集RMSE為0.662 mg/g(圖6)。

圖6 水稻氮素含量反演結(jié)果Fig.6 Inversion results of rice nitrogen content

2.5 反演結(jié)果分析與評價

當(dāng)前,已有眾多研究人員將植被指數(shù)作為模型輸入,對水稻冠層氮素含量進(jìn)行反演,從IndexData-Base數(shù)據(jù)庫中選取6種常用于氮素反演的植被指數(shù)[17],用來與本研究所構(gòu)建的新型植被指數(shù)PSVI進(jìn)行對比。將6種植被指數(shù)同樣采用線性回歸的方式構(gòu)建基于最優(yōu)觀測角度下偏振光譜的水稻冠層氮素含量反演模型,如圖7所示,結(jié)果顯示,這6種指數(shù)所構(gòu)建的反演模型的決定系數(shù)R2均小于 0.783 8,且反演精度最高的R2為0.755。綜上,本研究所構(gòu)建的新型植被指數(shù)PSVI線性反演效果最好。

3 討論

目前,已有較多研究人員利用無人機(jī)遙感獲取的水稻冠層光譜反射率進(jìn)行氮素含量的反演建模,但大多數(shù)研究僅基于垂直觀測下的光譜數(shù)據(jù),通過數(shù)據(jù)驅(qū)動的方式,構(gòu)建反演模型以提升水稻冠層氮素含量的反演精度。為了進(jìn)一步提升水稻冠層氮素含量的反演精度,本研究結(jié)合多角度偏振測量技術(shù)與數(shù)據(jù)驅(qū)動方式,對無人機(jī)偏振遙感的水稻冠層氮素含量反演方法開展了研究,基于水稻冠層多角度偏振光譜數(shù)據(jù),構(gòu)建了新型植被指數(shù)PSVI定量估算水稻冠層葉片氮素含量。由于不同觀測天頂角下的偏振光譜存在一定的差異,本研究結(jié)合6種常用于植物氮素含量反演的植被指數(shù),通過分析不同觀測天頂角下各植被指數(shù)與水稻冠層氮素含量的相關(guān)性,得到最優(yōu)觀測天頂角為后向觀測15°,這一結(jié)果與HE等[27]研究結(jié)果一致。這是由于后向觀測主要來自具有較高反射率的水稻冠層上方陽光照射的葉片,而前向觀測包含更多來自遮陰葉片的信息,較少來自陽光照射冠層水稻葉片的信息,因此后向觀測15°在去除偏振信息的同時能夠更多地收集到來自水稻冠層的光譜信息,從而具有更高的相關(guān)性。此外,偏振光譜受植被冠層的幾何結(jié)構(gòu)、平均葉傾角等影響較大。不同觀測天頂角下去除鏡面反射信息后,相關(guān)性的提升效果存在一定的差異,后續(xù)還需進(jìn)一步研究。

在得到最優(yōu)觀測天頂角下的偏振光譜數(shù)據(jù)后,為盡可能減少該偏振光譜信息中的冗余信息,減小變量間的共線性,本研究采用連續(xù)投影法得到與水稻冠層氮素含量相關(guān)性最好的6個特征波長,分別為500、566、663、691、736、763 nm。由于該算法可以在確保不同波段之間共線性最小的同時,保留原始數(shù)據(jù)中的大部分特征,因此篩選到的波長分別處于與葉片氮素含量較高的綠光和紅光波段,與前期的研究結(jié)果相吻合[28-29]。本研究在構(gòu)建植被指數(shù)的過程中,基于表2的相關(guān)性結(jié)果可知,比值植被指數(shù)(RVI)的整體相關(guān)性較好,因此本研究將比值植被指數(shù)的比值形式作為指數(shù)構(gòu)建的基礎(chǔ)形式,為進(jìn)一步提升指數(shù)與氮素含量的相關(guān)性,利用數(shù)學(xué)變換的方法,構(gòu)建了新型植被指數(shù)PSVI。該指數(shù)基于后向觀測15°下的水稻冠層偏振光譜數(shù)據(jù)進(jìn)行計算,將計算結(jié)果作為模型輸入,采用線性回歸方式進(jìn)行反演,得到的PSVI植被指數(shù)具有較好的反演精度。 然而在對PSVI植被指數(shù)的反演效果進(jìn)行驗證時,雖然該指數(shù)的反演效果優(yōu)于其它6個植被指數(shù),但該指數(shù)用于其它數(shù)據(jù)集時可能出現(xiàn)反演精度下降的情況,因此存在普適性不足的問題,這也是數(shù)據(jù)驅(qū)動模型的弊端,缺乏農(nóng)學(xué)機(jī)理上關(guān)于氮素含量與波段之間關(guān)系的解釋。未來的科研中需要對上述內(nèi)容進(jìn)行更深入探討。但由于植被指數(shù)計算簡單,檢測裝置易于開發(fā)與集成,基于植被指數(shù)的水稻冠層氮素反演方法仍具有一定的研究價值[30]。本文采用線性回歸方法構(gòu)建了水稻葉片氮含量的反演模型。與其它神經(jīng)網(wǎng)絡(luò)算法相比,線性回歸建模方法思想簡單,易于實現(xiàn),建模迅速,對于小數(shù)據(jù)量、簡單的關(guān)系很有效,能夠?qū)崿F(xiàn)對水稻冠層氮素含量的反演。

4 結(jié)論

(1)基于多角度偏振光譜數(shù)據(jù),分析得到水稻冠層氮素含量反演的最優(yōu)觀測天頂角為后向觀測15°。

(2)通過連續(xù)投影法(SPA)篩選得到的特征波長為500、566、663、691、736、763 nm。將6個特征波長兩兩組合,計算得到的PSVI指數(shù)與水稻氮素含量相關(guān)性最高時的波長組合為500、566 nm,相關(guān)系數(shù)R為0.850。

(3)以PSVI為輸入,運用線性回歸法構(gòu)建的水稻冠層氮素含量反演模型,模型訓(xùn)練集R2為 0.783 8,RMSE為 0.428 mg/g;驗證集RMSE為0.662 mg/g,模型的擬合度優(yōu)于歸一化植被指數(shù)、比值植被指數(shù)等傳統(tǒng)植被指數(shù)所構(gòu)建的水稻氮素含量反演模型。