殷曉飛，王春光，王海超，張文霞，宗哲英，張海軍

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 機電工程學(xué)院，呼和浩特 010018；2.呼和浩特職業(yè)學(xué)院 機電工程學(xué)院，呼和浩特 010051)

0 引言

草地生物量是單位面積內(nèi)牧草積累的物質(zhì)質(zhì)量，是草地生態(tài)系統(tǒng)的重要指標(biāo)之一。草地生物量常采用模型法和經(jīng)驗法進(jìn)行估測，但費時費力，不利于大面積開展且具有一定的破壞性[1]。高光譜遙感具有視野宏觀、動態(tài)的特點，是進(jìn)行大尺度草地動態(tài)監(jiān)測最為有效的工具之一，能夠?qū)崿F(xiàn)草地類型識別、草地退化監(jiān)測、植物生物化學(xué)參數(shù)檢測及植被覆蓋度、含水率、葉面積指數(shù)的估算。高光譜遙感已被廣泛應(yīng)用于草地生物量的估算，M. Mirik等對美國懷俄明州黃石國家公園的地面植物總量進(jìn)行了估算[2]；F. Fava等證明了高光譜數(shù)據(jù)在冠層水平上估算干草草生物量和植物多樣性指標(biāo)的潛力[3]；MA. Cho比較了偏最小二乘回歸、歸一化植被指數(shù)、紅邊位置對草地的生物量估算精度，結(jié)果表明偏最小二乘回歸的精度最高[4]；安海波等對內(nèi)蒙古烏蘭察布市四子王旗和呼和浩特市和林格爾縣天然和人工草地冠層光譜數(shù)據(jù)和生物量進(jìn)行分析，發(fā)現(xiàn)估算草地生物量最佳估算模型為優(yōu)化后的歸一化光譜指數(shù)[5]；胥慧等利用高光譜遙感對秋季干枯牧草生物量進(jìn)行研究，發(fā)現(xiàn)高光譜特征參數(shù)和草地生物量具有較高的相關(guān)性[6]；段敏杰等對西藏北部高寒草原的光譜反射率和植被生物量進(jìn)行測量分析，得出不同放牧強度下地上生物量和光譜植被指數(shù)之間存在顯著相關(guān)性[7]；孫成明等運用回歸分析方法分析歸一化植被指數(shù)與生物量的相關(guān)關(guān)系，構(gòu)建了南方草地地上生物量的光譜估算模型[8]；吳建付等利用高光譜遙感進(jìn)行了草地分類和生物量估測等研究[9]；劉占宇等對內(nèi)蒙古自治區(qū)錫林郭勒盟的天然草地進(jìn)行了高光譜遙感地面觀測[10]。

關(guān)于草地生物量的高光譜遙感估測模型因研究時間、區(qū)域、草地類型和光譜指數(shù)選擇不同而有所差別。從目前研究現(xiàn)狀來看，針對鄂爾多斯荒漠草原建立草地生物量高光譜遙感估測模型還未見相關(guān)報道。因此，利用近地面高光譜在不同退化梯度下估算鄂爾多斯荒漠草原草地生物量，能夠為開展草地生態(tài)效益評價、退牧還草效益檢測和航空遙感精細(xì)解譯提供技術(shù)支持，促進(jìn)高光譜遙感技術(shù)在草地定量檢測中的應(yīng)用，為合理制定草地載畜量提供科學(xué)依據(jù)。

1 試驗部分

1.1 研究區(qū)概況

研究區(qū)位于內(nèi)蒙古自治區(qū)鄂爾多斯市西部鄂托克旗境內(nèi)，東經(jīng)106°41′-108°54′，北緯38°18′-40°11′，平均海拔1 800m，年降水量為250mm，主要集中于7～9月，年均氣溫6.4℃；境內(nèi)地貌以波狀高原為主，西北高，東南低；天然草地以溫性荒漠草原為主，草地資源總面積為197.5萬hm2，可利用面積為174.7萬hm2，占草場總面積的88.43%。

1.2 試驗設(shè)計

根據(jù)李博提出的草地退化分級理論及相關(guān)專家多年實踐經(jīng)驗[11]，結(jié)合氣象資料，從鄂托克旗選出4個樣點，放牧強度基本一致，土壤類型均為棕鈣土。依草地干旱及利用程度由東向西劃分為對照(CK)、輕度退化(LD)、中度退化(MD)和重度退化(HD)[12]。 試驗區(qū)內(nèi)主要植被包括狹葉錦雞兒(Caraganastenophylla)、短花針茅(Stipabreviflora)、冷篙(Artemisiafrigida)、無芒隱子草(Cleistogenessongorica)、豬毛菜(Salsolacollina)和砂藍(lán)刺(Echinopsgmelini)等。

1.3 冠層光譜測量

采用美國ASD(Analytical Spectral Devices)公司便攜式高光譜儀器(350～1830nm)。其中，350～1 000nm光譜采樣間隔為1.4nm，光譜分辨率為3nm；1 000～1 830nm光譜采樣間隔為2nm，光譜分辨率為10nm。光譜測定在2014-2016年每年7月和8月選擇晴朗無云無風(fēng)天氣，于10：00-14:00進(jìn)行測定。采集數(shù)據(jù)前，先使用軟件程序進(jìn)行白板校正與去除暗電流操作，光譜采集探頭視場角為25°；采集冠層數(shù)據(jù)時探頭垂直于植被冠層上方，高度根據(jù)采樣植被冠層大小計算確定，每組樣品每間隔1s連續(xù)采集10次，將10次光譜平均后作為本組光譜數(shù)據(jù)。草地冠層反射光譜在草地類型、植被蓋度、植物水分、土壤類型和大氣狀況等多種因素作用下形成綜合的反射光譜。本文選取試驗樣地不同退化梯度下典型植物狹葉錦雞兒、短花針茅，退化指示植物無芒隱子草和冷蒿的草地類型進(jìn)行光譜測量和數(shù)據(jù)采集。

1.4 生物量測定

光譜數(shù)據(jù)采集后，在光譜測量區(qū)域正下方，將50cm×50cm樣方內(nèi)的植物鮮樣用剪刀齊地面剪下，按種類分別裝入自封袋，使用計量單位為1‰g天平立即其稱鮮重，并貼標(biāo)簽。

1.5 高光譜特征參數(shù)提取

高光譜特征參數(shù)是從原始光譜、一階微分光譜提取的基于高光譜位置變量、面積變量和植被指數(shù)變量等3種類型共19種特征參數(shù)，如表1所示。

2 結(jié)果與分析

2.1 草地生物量與高光譜變量相關(guān)分析

2.1.1 草地生物量與原始光譜相關(guān)分析

草地生物量與原始光譜相關(guān)分析如圖1所示。波長小于759nm時，原始光譜反射率與草地生物量呈負(fù)相關(guān)，在673nm相關(guān)系數(shù)達(dá)到最大，相關(guān)系數(shù)為-0.759，形成一個波谷；波長在759～943nm之間時，光譜反射率與草地生物量呈正相關(guān)，但相關(guān)性較差；波長在944～1 349nm和1 418～1 798nm之間時，原始光譜反射率與草地生物量呈負(fù)相關(guān)，波長在1 460nm時達(dá)到最大，相關(guān)系數(shù)為-0.709；波長在350～739nm和1 144～1 798nm之間時光譜反射率和草地生物量之間的相關(guān)系數(shù)均達(dá)到極顯著性檢驗水平。

表1 高光譜特征參數(shù)的定義和描述

續(xù)表1

————表示P≤0.05顯著相關(guān)(r=-0.079)水平線；——表示P≤0.01極顯著相關(guān)(r=-0.103)水平線。

圖1 草地生物量與原始光譜相關(guān)分析

Fig.1 Correlation Analysis of Grassland Biomass and Primitive Spectrum

2.1.2 草地生物量與一階微分光譜變量相關(guān)分析

反射率采用光譜微分技術(shù)處理提取一階微分，能夠有效地消除大氣效應(yīng)、陰影和土壤等對植被環(huán)境的影響，更好地反映植物本質(zhì)特征，增強光譜曲線在坡度上的細(xì)微變化。草地生物量與一階微分光譜變量相關(guān)分析如圖2所示。

————表示P≤0.05顯著相關(guān)(r=±0.080)水平線；——表示P≤0.01極顯著相關(guān)(r=±0.104)水平線。

圖2 草地生物量與一階微分光譜相關(guān)分析

Fig.2 The Correlation Analysis of Grassland Biomass and first derivative spectra

草地生物量與一階微分光譜之間在508～531nm和679～760nm呈正相關(guān)，達(dá)到極顯著性檢驗水平；在710nm，相關(guān)系數(shù)達(dá)0.812，且為全波段范圍內(nèi)相關(guān)系數(shù)最大值，位于紅邊范圍；草地生物量與一階微分光譜之間在421～505nm和535～673nm呈負(fù)相關(guān)，達(dá)到極顯著性檢驗水平；在568nm，相關(guān)系數(shù)達(dá)-0.789；草地生物量與一階微分光譜之間在其余光譜階段相關(guān)性較差。

2.1.3 草地生物量與高光譜特征參數(shù)相關(guān)分析

草地生物量與高光譜特征參數(shù)相關(guān)系數(shù)如表2所示。由表2可知：草地生物量與高光譜特征參數(shù)位置變量Dy、DR、Rg和Rr之間的相關(guān)系數(shù)較大，與Db、λy和λg之間的相關(guān)系數(shù)相對較小，上述變量均達(dá)到極顯著性檢驗水平。草地生物量與高光譜特征參數(shù)面積變量之間的相關(guān)系數(shù)中，與SDb相關(guān)系數(shù)較小，與SDy和SDR相關(guān)系數(shù)較大，且SDy呈負(fù)相關(guān)，SDR呈正相關(guān)，均達(dá)到極顯著性檢驗水平。草地生物量與高光譜特征參數(shù)植被指數(shù)變量之間的相關(guān)系數(shù)中，除SDR/SDy相關(guān)系數(shù)較小，其余相關(guān)系數(shù)都在0.7左右且都呈正相關(guān)，均達(dá)到極顯著性檢驗水平。在所有的高光譜特征參數(shù)中，SDR/SDb與草地生物量的相關(guān)系數(shù)最大，達(dá)到0.792。草地生物量與λb、λR、λr、SDR/SDy的相關(guān)關(guān)系未達(dá)到顯著檢驗水平，表明上述高光譜特征參數(shù)對草地生物量的變化并不敏感。

表2 草地生物量與高光譜特征參數(shù)相關(guān)系數(shù)

Table 2 Correlation coefficient between grassland biomass and hyperspectral characteristic parameters

特征參數(shù)相關(guān)系數(shù)特征參數(shù)相關(guān)系數(shù)特征參數(shù)相關(guān)系數(shù)Db0.418**λg-0.354**(Rg-Rr)/(Rg+Rr)0.720**λb-0.041Rr-0.752**SDR/SDb0.792**Dy-0.730**λr0.066SDR/SDy0.075λy0.509**SDb0.145**(SDR-SDb)/(SDR+SDb)0.699**DR0.608**SDy-0.780**(SDR-SDy)/(SDR+SDy)0.760**λR-0.130SDR0.785**Rg-0.701**Rg/Rr0.726**

*在0.05水平上顯著相關(guān)，**在0.01水平上顯著相關(guān)。

2.2 草地生物量高光譜遙感估測模型

為了對草地生物量進(jìn)行估算和精度檢驗，將在荒漠草原不同退化梯度下所獲取的600個樣本的試驗數(shù)據(jù)進(jìn)行隨機抽取，分為兩組。第1組為建模樣本(n=480)，用于進(jìn)行回歸分析；第2組為驗證樣本(n=120)，進(jìn)行模型檢驗。為了進(jìn)一步精確估算荒漠草原草地生物量，依據(jù)圖1、圖2和表2中相關(guān)系數(shù)分析，選取與草地生物量相關(guān)系數(shù)較大(大于0.7)且達(dá)到極顯著水平(P<0.01)的原始光譜敏感波段、一階微分光譜敏感波段和11種高光譜特征參數(shù)分別與草地生物量構(gòu)建單變量線性與非線性的回歸模型，如表3所示。

表3 草地生物量與原始光譜、一階微分光譜變量、高光譜特征參數(shù)的線性和非線性回歸方程決定系數(shù)

續(xù)表3

*在0.05水平上顯著相關(guān),**在0.01水平上顯著相關(guān)。

2.2.1 單變量線性與非線性回歸模型

選擇原始光譜、一階微分光譜變量與草地生物量相關(guān)系數(shù)最大波段R673、R′710，并從表2種選擇相關(guān)系數(shù)較大(大于0.7)的Dy、Rg、Rr、SDy、SDR、Rg/Rr、(Rg-Rr)/(Rg+Rr)、SDR/SDb和(SDR-SDy)/(SDR+SDy)等9種高光譜特征參數(shù)，進(jìn)行單變量線性與非線性回歸，建立草地生物量估測模型。

由表3可知：原始光譜R673、一階微分光譜R′710、9種高光譜特征參數(shù)與草地生物量的決定系數(shù)均通過0.01極顯著性檢驗水平。R673、Dy、Rg、Rr、SDy、Rg/Rr和(Rg-Rr)/(Rg+Rr)最佳估測模型為二次多項式，決定系數(shù)均大于0.5；R′710、SDR和SDR/SDb最佳估測模型為冪函數(shù)，決定系數(shù)均大于0.7；(SDR-SDy)/(SDR+SDy)最佳估測模型為指數(shù)函數(shù)，其回歸模型決定系數(shù)為0.729，且為單變量線性與非線性回歸分析中高光譜特征參數(shù)與草地生物量最佳估測模型。

2.2.2 逐步回歸模型

逐步回歸分析方法是從眾多自變量中有效選擇重要變量的一種統(tǒng)計方法[13]。確定與生物量相關(guān)的高光譜特征參數(shù)時，草地生物量和光譜變量值作為輸入變量，輸出結(jié)果是一系列包含不同光譜變量的多元線性方程及對應(yīng)的決定系數(shù)(R2)和回歸方程檢驗值(F檢驗)。逐步回歸方程表達(dá)式為

y=p0+p1x1+p2x2+…+pixi

其中，y為因變量；xi為第i個自變量；p0為常數(shù)項；pi為y對應(yīng)的偏回歸系數(shù)。

將Dy、Rg、Rr、SDy、SDR、Rg/Rr、(Rg-Rr)/(Rg+Rr)、SDR/SDb、(SDR-SDy)/(SDR+SDy)等9個高光譜特征參數(shù)作為自變量，草地生物量作為因變量，進(jìn)行逐步回歸分析。步進(jìn)條件為：當(dāng)F的概率≤0.05時，變量輸入；當(dāng)F的概率≥0.10時，變量去除。分析結(jié)果決定系數(shù)最大時在選定的9個高光譜特征參數(shù)中只有4個變量入選，用于草地生物量估測模型的逐步回歸方程如下：y=371.484+1713.954x1-2166.5x2+1207.404x3+5760.705x4。其中，x1為(Rg-Rr)/(Rg+Rr)、x2為Rg、x3為SDR、x4為SDy，其決定系數(shù)為0.709?？梢钥闯觯翰莸厣锪康奶卣髯兞扛鶕?jù)重要性由大到小依次是(Rg-Rr)/(Rg+Rr)、Rg、SDR和SDy。

2.2.3 BP神經(jīng)網(wǎng)絡(luò)回歸模型

采用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行回歸模型建立，網(wǎng)絡(luò)共有輸入層、中間層和輸出層3層。將相關(guān)系數(shù)較高的Dy、Rg、Rr、SDy、SDR、Rg/Rr、(Rg-Rr)/(Rg+Rr)、SDR/SDb、(SDR-SDy)/(SDR+SDy)等9個高光譜特征參數(shù)作為神經(jīng)網(wǎng)絡(luò)的輸入神經(jīng)元，將草地生物量作為輸出神經(jīng)元，經(jīng)過多次訓(xùn)練比較分析得到中間層神經(jīng)元的個數(shù)為5。結(jié)果表明：估算值與實測值的決定系數(shù)R2為0.766，RMSE為40.19g/m2，相對誤差為29.88%。模擬精度較逐步回歸分析法有所提高。

2.2.4 草地生物量的高光譜遙感估測模型精度驗證

上述表中挑選決定系數(shù)在0.7以上的回歸方程，作為估測模型進(jìn)行驗證。模型精度通過估算值與實測值之間的決定系數(shù)R2、均方根誤差RMSE(Root Mean Square Error)和相對誤差RE(Relative Error)進(jìn)行驗證。最優(yōu)的估測模型，應(yīng)該估算值與實測值具有較高的決定系數(shù)(R2)，較小的均方根誤差RMSE和相對誤差RE。

草地生物量估測模型的精度驗證如表5所示。從表5中可以看出：基于冪函數(shù)和指數(shù)函數(shù)單變量回歸模型決定系數(shù)小于基于多變量逐步回歸和BP神經(jīng)網(wǎng)絡(luò)。在預(yù)測的4種模型中，BP神經(jīng)網(wǎng)絡(luò)的多變量回歸預(yù)測模型最優(yōu)，此時草地生物量實測值與估算值決定系數(shù)最大，為0.766；均方根誤差RMSE和相對誤差RE值均最小，分別是40.19g/m2和29.88%。

表5 草地生物量估測模型的精度驗證

為更直觀反映各模型的預(yù)測效果，繪制草地生物量實測值與估算值1:1關(guān)系圖，如圖3所示。從圖3中可以看出：BP神經(jīng)網(wǎng)絡(luò)估算值與實測值分布最接近1:1黑色實線。這說明，該模型預(yù)測效果最好。

圖3 基于不同模型的草地生物量的實測值與估算值1:1分布圖

3 結(jié)論與討論

鄂爾多斯荒漠草原草地生物量與草地冠層原始光譜、一階微分光譜和高光譜特征參數(shù)都顯著相關(guān)。其中，在673nm波段，原始光譜與草地生物量相關(guān)系數(shù)為-0.759；在710nm波段，一階微分光譜與草地生物量相關(guān)系數(shù)為0.812，Dy、Rg、Rr、SDy、SDR、Rg/Rr、(Rg-Rr)/(Rg+Rr)、SDR/SDb和(SDR-SDy)/(SDR+SDy)等9個高光譜特征參數(shù)與草地生物量相關(guān)系數(shù)均大于0.7，其中SDR/SDb相關(guān)系數(shù)最大，相關(guān)系數(shù)為0.792。

通過線性、非線性、逐步回歸和BP神經(jīng)網(wǎng)絡(luò)等方法進(jìn)行回歸分析，選取最優(yōu)估測模型進(jìn)行驗證分析。其中，BP神經(jīng)網(wǎng)絡(luò)回歸模型決定系數(shù)最高，且均方根誤差和相對誤差最小，建立回歸模型最好，決定系數(shù)R2為0.766，均方根誤差RMSE為40.19g/m2,相對誤差RE為29.88%。其中，逐步回歸建立回歸模型決定系數(shù)大于冪函數(shù)和指數(shù)函數(shù)，均方根誤差小于冪函數(shù)和指數(shù)函數(shù)，但相對誤差較大。

采用原始光譜、一階微分光譜及高光譜特征參數(shù)等單變量構(gòu)建的線性與非線性回歸模型，在估算草地生物量時效果相對較差，而逐步回歸和BP神經(jīng)網(wǎng)絡(luò)充分利用多種高光譜特征參數(shù)構(gòu)建回歸模型，彌補單變量信息相對較少的缺點，估算精度相對較好。其中，BP神經(jīng)網(wǎng)絡(luò)在荒漠草原草地生物量估算結(jié)果為最優(yōu)。

利用草地冠層高光譜特征參數(shù)能夠精確估測鄂爾多斯荒漠草原草地生物量，為鄂爾多斯荒漠草原的動態(tài)監(jiān)測、草產(chǎn)估測、草地合理利用提供了科學(xué)有效的管理手段。