基于多維高光譜植被指數(shù)的冬小麥葉面積指數(shù)估算

吾木提·艾山江 尼加提·卡斯木 陳 晨 買買提·沙吾提

(1.伊犁師范大學(xué)資源與生態(tài)研究所, 伊寧 835000; 2.伊犁師范大學(xué)生物與地理科學(xué)學(xué)院, 伊寧 835000；3.新疆大學(xué)地理科學(xué)學(xué)院， 烏魯木齊 830046)

0 引言

冬小麥?zhǔn)俏覈?guó)重要的糧食作物之一，其種植面積和產(chǎn)量?jī)H次于水稻[1]。LAI是與植被利用光能和冠層結(jié)構(gòu)密切相關(guān)，能夠很好地反映作物生長(zhǎng)趨勢(shì)和估算產(chǎn)量的重要作物參數(shù)之一[2]。因此，通過(guò)遙感準(zhǔn)確、快速地估算區(qū)域作物的LAI不僅有利于更好的作物生長(zhǎng)連續(xù)長(zhǎng)期監(jiān)測(cè)，且有利于在作物綜合管理和精確農(nóng)業(yè)中的應(yīng)用[3]。

隨著高光譜技術(shù)的發(fā)展，其數(shù)百或數(shù)千個(gè)窄波段為開發(fā)靈敏的植被指數(shù)提供了更多可能。諸多研究學(xué)者已經(jīng)完善了傳統(tǒng)光譜植被指數(shù)，并利用高光譜數(shù)據(jù)的光譜分辨率和連續(xù)性開發(fā)了新的數(shù)據(jù)分析技術(shù)[4]。目前，冬小麥葉面積指數(shù)的遙感估算大致從單波段反射率和波段組合植被指數(shù)著手，結(jié)合不同回歸算法來(lái)實(shí)現(xiàn)。對(duì)于單波段反射率而言，諸多學(xué)者利用不同的預(yù)處理數(shù)學(xué)方法。孟禹弛等[5]利用地面實(shí)測(cè)冬小麥生理參數(shù)和冠層光譜數(shù)據(jù)，通過(guò)一階微分預(yù)處理并計(jì)算植被指數(shù)，建立了冥函數(shù)模型和指數(shù)模型；梁棟等[6]通過(guò)對(duì)冬小麥高光譜曲線進(jìn)行連續(xù)小波變換，并與LAI建立支持向量機(jī)回歸模型，實(shí)現(xiàn)了冬小麥不同生育期LAI估算；付元元等[7]在高光譜遙感反演LAI時(shí)，利用主成分回歸(Principle component regression, PCR)或偏最小二乘回歸(Partial least squares regression, PLSR)結(jié)合的高光譜維數(shù)約簡(jiǎn)方法，降低基于全譜段建立的模型不穩(wěn)定，同時(shí)取得了較好的估算結(jié)果。而對(duì)于典型的植被指數(shù)而言，陳雪洋等[8]探討了歸一化植被指數(shù)(Normalized difference vegetation index, NDVI)、比值植被指數(shù)(Ratio vegetation index, RVI)、土壤調(diào)節(jié)植被指數(shù)(Soil-adjusted vegetation index, SAVI)、增強(qiáng)型植被指數(shù)(Enhanced vegetation index, EVI)和紅邊參數(shù)在建立高光譜模型反演冬小麥LAI的精度；除此之外，夏天等[9]、陳曉凱等[10]、田明璐等[11]、孫華林等[12]、陶惠林等[13]、吾木提·艾山江等[14]利用RVI、NDVI和差值植被指數(shù)(Difference vegetation index, DVI)、EVI、再歸一化植被指數(shù)(Resource difference vegetation index, RDVI)等典型的光譜指數(shù)，采用任意兩波段組合算法提取最優(yōu)窄波段組合植被指數(shù)，通過(guò)多元逐步回歸和偏最小二乘等線性回歸算法構(gòu)建了冬小麥LAI估算模型，取得了較好的效果。以上研究顯示以任意兩波段組合構(gòu)建植被指數(shù)能夠?qū)ふ易魑锢砘瘏?shù)的響應(yīng)位置，有利于提高模型的估算精度[15-16]。而結(jié)合高光譜數(shù)據(jù)的光譜分辨率優(yōu)勢(shì)，對(duì)典型的三波段植被指數(shù)的波段組合進(jìn)行優(yōu)化，并與單波段反射率、兩波段植被指數(shù)之間進(jìn)行系統(tǒng)的比較，尋求符合當(dāng)?shù)丨h(huán)境的LAI敏感波段組合植被指數(shù)和提高估算精度方面值得探討。

本文以冬小麥在拔節(jié)期的LAI和同步冠層高光譜為數(shù)據(jù)源，采用波段組合優(yōu)化方法，構(gòu)建基于不同類型數(shù)據(jù)下(包括原始光譜、一階導(dǎo)數(shù)、二階導(dǎo)數(shù))所有可能組合的兩波段植被指數(shù)和三波段植被指數(shù)，探索不同維度(單波段反射率、兩波段植被指數(shù)、三波段植被指數(shù))數(shù)據(jù)對(duì)LAI的敏感性變化，利用人工神經(jīng)網(wǎng)絡(luò)(ANN)、K近鄰(KNN)和支持向量回歸(SVR)構(gòu)建LAI估算模型，并進(jìn)行對(duì)比和分析篩選最優(yōu)模型，以其為精準(zhǔn)農(nóng)業(yè)估算作物參數(shù)提供理論基礎(chǔ)。

1 估算方法

1.1 研究區(qū)概況

研究區(qū)地理位置為89°13′～91°22′E，42°25′～45°29′N，位于新疆維吾爾自治區(qū)奇臺(tái)縣，年平均氣溫5.5℃，年平均無(wú)霜期153 d，年平均降水量269.4 mm，屬于沙漠干旱氣候[17]。該地區(qū)主要作物以冬小麥為主，播種時(shí)間為2017年9月25日。根據(jù)冬小麥的關(guān)鍵生長(zhǎng)期和當(dāng)?shù)靥鞖鈼l件，從2018年4—6月，在起身期、拔節(jié)期、開花期等階段采集樣品。地面觀測(cè)在東西方向每隔20 m和南北方向30 m，在130 m×420 m的地塊進(jìn)行。采用五點(diǎn)采樣法進(jìn)行采樣，共有78個(gè)采樣點(diǎn)，并用GPS記錄取樣點(diǎn)的坐標(biāo)(圖1)。

1.2 數(shù)據(jù)采集與處理

在起身期和拔節(jié)期，冬小麥冠層的高光譜數(shù)據(jù)由ASD FieldSpec3型光譜儀進(jìn)行測(cè)量，光譜測(cè)量于北京時(shí)間10:00—14:00無(wú)云照射中進(jìn)行。為每個(gè)樣本收集了10條光譜曲線，測(cè)量間隔為0.1 s，求平均值作為該樣本測(cè)量數(shù)據(jù)。高光譜數(shù)據(jù)預(yù)處理部分，因受邊緣帶350～399 nm的高光噪聲，將該區(qū)域進(jìn)行刪除。其余完整連續(xù)的光譜曲線400～1 300 nm作為數(shù)據(jù)源，利用Savitzky-Golay方法進(jìn)行平滑，再利用一階微分和二階微分等數(shù)學(xué)方法完成高光譜數(shù)據(jù)變換。

基于數(shù)字圖像處理技術(shù)的攝影方法可用于估算LAI[18]。本研究采樣區(qū)域面積為Ss，莖總數(shù)設(shè)為n。5根莖從樣本區(qū)取出，作為掃描的小麥樣本。將一張A4白紙放在平面寫字板上，面積為Sp(21 cm×29.7 cm)，用膠水將葉片粘在紙上，以確保葉片之間沒有重疊。利用ECOSYS FS-1125 MFP型多功能數(shù)字復(fù)合掃描儀掃描所有樣品。分辨率設(shè)置為600像素×600像素，掃描的圖像保存為TIFF格式。使用ENVI打開上述TIFF文件，使用“決策樹”的分類方法將綠葉和白紙分開，然后分別計(jì)算樣品占用的像素?cái)?shù)PZ和綠葉占用的像素?cái)?shù)PY。葉面積指數(shù)計(jì)算公式為

1.3 任意波段組合植被指數(shù)

探索冬小麥LAI數(shù)據(jù)與高光譜數(shù)據(jù)之間重要關(guān)系的最有效方法之一是對(duì)窄帶反射率因子進(jìn)行計(jì)算。該研究涉及到任意波段組合植被指數(shù)，分別為二維植被指數(shù)(Two-dimensional vegetation index, 2DVI)和三維植被指數(shù)(Three-dimensional vegetation index, 3DVI)。2DVI為兩個(gè)波段進(jìn)行任意組合構(gòu)建兩波段植被指數(shù)，如：比值植被指數(shù)(RVI)、歸一化植被指數(shù)(NDVI)和差值植被指數(shù)(DVI)；3DVI為3個(gè)波段進(jìn)行任意組合構(gòu)建三波段植被指數(shù)(表1)。通過(guò)研究團(tuán)隊(duì)在Java平臺(tái)上開發(fā)的軟件來(lái)實(shí)現(xiàn)兩波段組合計(jì)算過(guò)程，利用Matlab軟件完成三維植被指數(shù)的構(gòu)建與計(jì)算。

表1 波段組合植被指數(shù)Tab.1 Band combination vegetation index

1.4 模型建立與驗(yàn)證

采用的機(jī)器學(xué)習(xí)算法包括：人工神經(jīng)網(wǎng)絡(luò)(ANN)[22]、K近鄰(KNN)[23]和支持向量回歸(SVR)[24]。

以決定系數(shù)(Coefficient of determination,R2)和均方根誤差(Root mean square error, RMSE)來(lái)驗(yàn)證模型的穩(wěn)定性，相對(duì)分析誤差(Relative percent deviation, RPD)為計(jì)算過(guò)程標(biāo)準(zhǔn)偏差(Standard deviation, SD)與估計(jì)標(biāo)準(zhǔn)誤差(Standard error of prediction, SEP)之比，可驗(yàn)證模型的預(yù)測(cè)能力。RPD大于2.0，說(shuō)明該模型具有較好的預(yù)測(cè)能力；RPD在1.4～2.0之間，代表具有一般預(yù)測(cè)能力；RPD小于1.4，說(shuō)明其預(yù)測(cè)能力較差。

2 結(jié)果與分析

2.1 LAI數(shù)據(jù)的統(tǒng)計(jì)分析

對(duì)研究區(qū)冬小麥冠層LAI數(shù)據(jù)進(jìn)行了統(tǒng)計(jì)，結(jié)果見表2。本次采樣共78個(gè)，將數(shù)據(jù)劃分為兩個(gè)部分，一個(gè)子集用于模型建模(n=47)，另一個(gè)子集用于模型驗(yàn)證(n=31)。對(duì)應(yīng)于建模集和驗(yàn)證集的LAI平均值分別為6.38和6.46，標(biāo)準(zhǔn)偏差(SD)分別為1.28和0.96。所有采樣點(diǎn)的LAI平均值為6.42，變異系數(shù)為16.51%。

2.2 冠層高光譜數(shù)據(jù)與LAI相關(guān)性分析

通過(guò)整數(shù)階微分對(duì)高光譜進(jìn)行預(yù)處理，分析基于不同預(yù)處理方法的單波段以及兩波段植被指數(shù)對(duì)LAI相關(guān)性的影響。

表2 冬小麥LAI數(shù)據(jù)統(tǒng)計(jì)分析Tab.2 Statistical analysis of winter wheat LAI data

2.2.1單波段反射率與LAI關(guān)系

基于原數(shù)據(jù)(Original data, OD)、一階和二階微分的預(yù)處理，單波段反射率與LAI之間相關(guān)關(guān)系存在明顯的差異(圖2)。原數(shù)據(jù)光譜反射率與LAI的相關(guān)性在可見光(450～710 nm)、近紅外(750～910 nm、1 020～1 100 nm)部分0.01水平上呈顯著相關(guān)，相關(guān)系數(shù)在-0.4～0.3之間變化。整數(shù)階微分處理后，與LAI的相關(guān)系數(shù)有所提高，提高到0.6以上；一階微分和二階微分相比，一階微分預(yù)處理不僅提高了相關(guān)性，而且增加了通過(guò)0.01水平顯著性波段范圍。

圖2 基于不同預(yù)處理的高光譜數(shù)據(jù)與LAI相關(guān)性結(jié)果Fig.2 Correlation results between hyperspectral data and LAI based on different preprocessing

2.2.2兩波段植被指數(shù)與LAI關(guān)系

基于原數(shù)據(jù)、一階和二階微分的預(yù)處理，進(jìn)行兩波段植被指數(shù)(RVI、NDVI、DVI)與LAI之間相關(guān)性熱圖分析(圖3～5)。對(duì)原數(shù)據(jù)而言，所有兩波段組合的3種植被指數(shù)篩選發(fā)現(xiàn)(表3)，與LAI的相關(guān)系數(shù)絕對(duì)值最大的對(duì)應(yīng)波段組合為OD-RVI(728 nm,994 nm)、OD-NDVI(724 nm,987 nm)、OD-DVI(785 nm,860 nm)；對(duì)于一階導(dǎo)數(shù)變換而言，LAI敏感的波段組合植被指數(shù)為OD-RVI(627 nm，774 nm)、OD-NDVI(769 nm，974 nm)、OD-DVI(729 nm，811 nm)；對(duì)于二階導(dǎo)數(shù)變換而言，LAI敏感的波段組合植被指數(shù)為OD-RVI(735 nm，776 nm)、OD-NDVI(539 nm，762 nm)、OD-DVI(686 nm，764 nm)。從波段組合范圍以及相關(guān)程度分析發(fā)現(xiàn)，整數(shù)階微分預(yù)處理對(duì)3種植被指數(shù)與LAI之間的相關(guān)性有所提升，為估算模型的構(gòu)建提供了可選擇的變量。

圖3 基于原始光譜的兩波段植被指數(shù)與LAI的相關(guān)性Fig.3 Correlation between two-band vegetation index and LAI based on original spectral data

圖4 基于一階微分光譜的兩波段植被指數(shù)與LAI的相關(guān)性Fig.4 Correlation between two-band vegetation index and LAI based on the first derivative spectral data

圖5 基于二階微分光譜的兩波段植被指數(shù)與LAI的相關(guān)性Fig.5 Correlation between two-band vegetation index and LAI based on the second derivative spectral data

表3 最佳兩波段組合植被指數(shù)Tab.3 The best two-band combination vegetation index

圖6 基于原始光譜的三波段植被指數(shù)與LAI的相關(guān)性Fig.6 Correlation between three-band vegetation index and LAI based on original spectral data

2.2.3三波段植被指數(shù)與LAI關(guān)系

基于原數(shù)據(jù)、一階和二階微分的預(yù)處理，進(jìn)行三波段植被指數(shù)與LAI之間相關(guān)性分析(圖6～8)?；谠几吖庾V數(shù)據(jù)的三波段植被指數(shù)計(jì)算結(jié)果中，僅選取3種植被指數(shù)相關(guān)性結(jié)果，見圖6；分析發(fā)現(xiàn)3DVI-1、3DVI-3、3DVI-7分別與LAI呈顯著相關(guān)關(guān)系，相關(guān)系數(shù)絕對(duì)值介于0.4～0.6之間。其余4種指數(shù)(3DVI-2、3DVI-4、3DVI-5、3DVI-6)均與LAI的相關(guān)程度偏低，未參與建模過(guò)程。

圖7 基于一階微分光譜的三波段植被指數(shù)與LAI的相關(guān)性Fig.7 Correlation between three-band vegetation index and LAI based on the first derivative spectral data

圖8 基于二階微分光譜的三波段植被指數(shù)與LAI的相關(guān)性Fig.8 Correlation between three-band vegetation index and LAI based on the second derivative spectral data

基于一階微分的三波段植被指數(shù)計(jì)算結(jié)果中，與LAI呈顯著相關(guān)的植被指數(shù)有5種，分別是3DVI-1、3DVI-2、3DVI-4、3DVI-5、3DVI-7，相關(guān)系數(shù)絕對(duì)值介于0.6～0.9之間，相關(guān)性三維結(jié)果見圖7；其中3DVI-4與LAI之間相關(guān)系數(shù)達(dá)到最大?；诙A微分的三波段植被指數(shù)與LAI之間的相關(guān)性三維熱圖見圖8，7種三波段植被指數(shù)均呈顯著相關(guān)，相關(guān)系數(shù)絕對(duì)值介于0.6～0.7之間，見表4。

在不同光譜預(yù)處理下，對(duì)三波段植被指數(shù)波段組合進(jìn)行優(yōu)化，并進(jìn)行最佳波段組合篩選；若出現(xiàn)相同大小的最值，都被采集為最佳波段組合，并適用于估算模型的最佳參數(shù)。

2.3 模型建立

諸多研究表明，基于多種波段組合植被指數(shù)構(gòu)建的農(nóng)作物參數(shù)估算模型精度顯著高于敏感單波段反射率建立的模型，不同波段之間的運(yùn)算結(jié)果能夠有效加強(qiáng)農(nóng)作物參數(shù)的關(guān)聯(lián)程度，從而提高估算模型的精度[25-26]。本文對(duì)高光譜窄波段反射率及其不同形式組合運(yùn)算構(gòu)建的2DVI、3DVI植被指數(shù)為數(shù)據(jù)源，并以最佳變量構(gòu)建模型如下：以原數(shù)據(jù)、一階和二階微分?jǐn)?shù)據(jù)的KNN、ANN和SVR模型；對(duì)各變換光譜數(shù)據(jù)分別進(jìn)行最佳2DVI、3DVI植被指數(shù)的KNN、ANN和SVR模型；因參與模型變量不同，共構(gòu)建了27種估算模型，模型估算精度見圖9。

估算模型所需數(shù)據(jù)總體分為單波段反射率、2DVI、3DVI等3種類型(簡(jiǎn)稱：一維、二維、三維)，在3種類型光譜數(shù)據(jù)下，采用K近鄰(KNN)、人工神經(jīng)網(wǎng)絡(luò)(ANN)和支持向量回歸(SVR)建模對(duì)冬小麥LAI進(jìn)行估算；分析發(fā)現(xiàn)，對(duì)光譜數(shù)據(jù)進(jìn)行一階二階微分預(yù)處理，顯著提高了各維度數(shù)據(jù)構(gòu)建的估算模型精度。尤其是一階微分光譜處理在KNN模型各維度數(shù)據(jù)構(gòu)建模型中表現(xiàn)突出，決定系數(shù)R2分別為0.57、0.83、0.89；對(duì)于ANN和SVR模型而言，在一維、二維數(shù)據(jù)構(gòu)建的模型中表現(xiàn)突出；在三維數(shù)據(jù)構(gòu)建模型中，基于二階微分預(yù)處理的光譜數(shù)據(jù)結(jié)果比較偏高。從總體的估算模型對(duì)比分析，KNN構(gòu)建的估算模型精度在各維度數(shù)據(jù)中，決定系數(shù)均偏高，而均方根誤差均偏低。其中以一階微分處理運(yùn)算的最佳3DVI植被指數(shù)為自變量，建立冬小麥葉面積指數(shù)KNN估算模型表現(xiàn)最優(yōu)，其決定系數(shù)最高(R2=0.89)，均方根誤差最低(RMSE為0.31)。

依據(jù)模型穩(wěn)定性和預(yù)測(cè)能力評(píng)價(jià)指標(biāo)RPD對(duì)27種模型的估算能力進(jìn)行比較分析，結(jié)果見表5。本研究不具備估測(cè)樣本能力的模型(RPD小于1.5)主要是基于原數(shù)據(jù)單波段反射率所構(gòu)建的模型，其余部分模型的RPD均大于1.5，具備較好的估算能力。從總體分析發(fā)現(xiàn)，KNN模型的估算能力明顯高于ANN和SVR模型的估算能力；尤其是以最佳波段組合優(yōu)化植被指數(shù)為模型參數(shù)時(shí)，KNN模型RPD大于2.0的模型數(shù)量占比較大，共5種模型；而其余RPD大于2.0的模型數(shù)量各占3個(gè)。從光譜數(shù)據(jù)變換角度分析發(fā)現(xiàn)，在一階微分預(yù)處理后分別進(jìn)行最佳2DVI和3DVI植被指數(shù)的計(jì)算中，冬小麥LAI與最佳波段組合植被指數(shù)之間的相關(guān)性明顯提升，從而基于一階微分處理所構(gòu)建模型的估算能力顯著提高。

綜上所述，在不同的數(shù)據(jù)類型下，建立的機(jī)器學(xué)習(xí)估算模型存在一定的差異。冬小麥LAI估算模型(KNN-3DVI-FD)的實(shí)測(cè)值與預(yù)測(cè)值的空間分布見圖10，可以看出，以最佳3DVI植被指數(shù)變量，采用KNN算法構(gòu)建冬小麥LAI估算模型，其實(shí)測(cè)值和預(yù)測(cè)值空間上擬合效果更接近于1∶1，說(shuō)明在本研究中是最具有估算能力的模型。

3 討論

掌握農(nóng)作物的生長(zhǎng)狀況主要取決于作物的葉面積指數(shù)，LAI是能夠描述農(nóng)作物冠層結(jié)構(gòu)的重要參數(shù)之一[27-28]。在冬小麥整個(gè)生育期中，拔節(jié)期氮素積累過(guò)程較顯著，直接反映冬小麥健康生長(zhǎng)趨勢(shì)，對(duì)冬小麥拔節(jié)期LAI的精確、快速遙感監(jiān)測(cè)能夠?yàn)楹笃谔镩g尺度施肥管理提供理論依據(jù)[29-30]。

基于高光譜數(shù)據(jù)以及最優(yōu)植被指數(shù)的農(nóng)作物典型參數(shù)反演得到了諸多學(xué)者的研究，亞森江·喀哈爾等[31]、陳曉凱等[10]分別對(duì)春小麥的葉綠素含量和冬小麥的葉面積指數(shù)進(jìn)行高光譜監(jiān)測(cè)，探討了任意兩波段組合植被指數(shù)與作物參數(shù)之間的相關(guān)關(guān)系，表明基于最優(yōu)兩波段組合植被指數(shù)構(gòu)建的模型精度顯著提高。本研究基于光譜變換和最優(yōu)波段組合植被指數(shù)進(jìn)行冬小麥拔節(jié)期LAI估算研究，結(jié)果表明原數(shù)據(jù)、一階二階微分變換與LAI在紅邊波段750～910 nm、760～810 nm內(nèi)呈顯著相關(guān)(P<0.01)。通過(guò)分析發(fā)現(xiàn)，基于3種光譜變換所構(gòu)建的任意兩波段、三波段植被指數(shù)與LAI相關(guān)性顯著提高，顯著相關(guān)波段組合中包含紅邊波段(表4)，這與劉偉東等[32]、蘇偉等[33]研究成果相符，都反映了紅邊波段對(duì)作物L(fēng)AI的重要性，表明任意波段組合的植被指數(shù)應(yīng)用到農(nóng)作物參數(shù)反演是可行的。

表4 最佳三波段組合植被指數(shù)Tab.4 The best two-band combination vegetation index

在諸多研究中，經(jīng)常假設(shè)自變量和因變量呈線性關(guān)聯(lián)而構(gòu)建回歸模型，在實(shí)際中該條件較難滿足。陳曉凱等[10]研究表明最優(yōu)窄波段植被指數(shù)與冬小麥LAI的關(guān)系表現(xiàn)為非線性關(guān)系。因此，本文為擬合自變量和因變量之間的非線性關(guān)系，引入了KNN、ANN和SVR等3種機(jī)器學(xué)習(xí)算法；通過(guò)分析3種算法發(fā)現(xiàn)，K近鄰算法對(duì)數(shù)據(jù)的非線性擬合結(jié)果較好，尤其是基于一階微分光譜處理的9個(gè)3DVI指數(shù)(表4)與LAI之間構(gòu)建KNN模型相對(duì)于其它兩種算法明顯提高，更適合解決非線性擬合問(wèn)題。

4 結(jié)束語(yǔ)

以冬小麥冠層高光譜為數(shù)據(jù)源，在一階、二階微分預(yù)處理的基礎(chǔ)上，計(jì)算任意波段組合的2DVI指數(shù)和3DVI指數(shù)，利用人工神經(jīng)網(wǎng)絡(luò)(ANN)、K近鄰(KNN)和支持向量回歸(SVR)算法來(lái)構(gòu)建LAI估算模型，并進(jìn)行對(duì)比分析。結(jié)果表明：基于不同預(yù)處理光譜的任意波段組合植被指數(shù)在紅邊位置與LAI相關(guān)性顯著提高，尤其是基于一階微分預(yù)處理光譜的FD-3DVI-4(714 nm, 400 nm, 1 001 nm)指數(shù)與LAI相關(guān)系數(shù)達(dá)到0.93(P<0.01)，更能說(shuō)明充分利用海量光譜信息很有必要，能夠避免光譜信息的丟失和失效。對(duì)于3種模型構(gòu)建方法而言，以最優(yōu)波段組合3DVI為自變量，發(fā)現(xiàn)K近鄰算法數(shù)據(jù)擬合好，構(gòu)建的模型精度較好，其決定系數(shù)最高(R2=0.89)，均方根誤差最低(RMSE為0.31)，相對(duì)分析誤差為2.41，可為作物長(zhǎng)勢(shì)評(píng)估和估產(chǎn)提供科學(xué)依據(jù)。

圖9 冬小麥LAI高光譜估算模型精度對(duì)比Fig.9 Comparison of LAI hyperspectral estimation models for winter wheat

表5 冬小麥LAI的27種估算模型相對(duì)分析誤差Tab.5 Relative analysis errors of 27 estimation models of winter wheat LAI

圖10 實(shí)測(cè)值與預(yù)測(cè)值的散點(diǎn)圖Fig.10 Scatter plot of measured and predicted values