宋春雨，甘 淑， ，袁希平 ，胡 琳，李 雁，閆馨方

(1.昆明理工大學(xué) 國(guó)土資源工程學(xué)院，云南 昆明 650093；2.云南省高校高原山區(qū)空間信息測(cè)繪技術(shù)應(yīng)用工程研究中心，云南 昆明 650093；3.滇西應(yīng)用技術(shù)大學(xué)，云南 大理 671009)

植被的光譜特性對(duì)植被遙感解譯和定量分析具有重要意義[1]。掌握植被資源現(xiàn)狀及其變化，對(duì)于促進(jìn)林業(yè)發(fā)展和人類(lèi)生存環(huán)境可持續(xù)發(fā)展具有重要意義，因此加強(qiáng)植被資源的監(jiān)測(cè)已經(jīng)成為當(dāng)前的熱點(diǎn)話題[2-4]。采用傳統(tǒng)手段進(jìn)行植被監(jiān)測(cè)會(huì)浪費(fèi)大量時(shí)間和人力財(cái)力，多光譜遙感可以識(shí)別植被的大類(lèi)，但是難以區(qū)分更為詳細(xì)的植被種類(lèi)[5]。高光譜可以通過(guò)窄而連續(xù)的波段描述植被的光譜特征，從而區(qū)分識(shí)別不同的植被類(lèi)型[6-8]。使用高光譜技術(shù)分析植被的光譜特征不僅可以識(shí)別區(qū)分不同的植被類(lèi)型，還可以完善植被光譜數(shù)據(jù)庫(kù)，為處理衛(wèi)星遙感影像提供重要信息，這對(duì)于遙感監(jiān)測(cè)植被及其反演都具有重要意義[9-10]。側(cè)柏(Platycladusorientalis)、樟木(Cinnamomumlongepaniculatum)、柳杉(CryptomeriafortuneiHooibrenk)都是滇中典型的喬木植被，都可用于園林綠化和醫(yī)藥用途；此外，這3種喬木還是普遍使用的建筑和家具木材。因此，研究側(cè)柏、樟木、柳杉的高光譜特征具有重要意義。

測(cè)量分析植被的光譜特性是林業(yè)遙感的基礎(chǔ)。目前國(guó)內(nèi)外已經(jīng)有許多學(xué)者對(duì)植被的光譜特征進(jìn)行了大量的研究，并且取得了許多成果[11-13]。如加力戈等[14]通過(guò)分析桂柳、梭梭、白刺3種荒漠植被生長(zhǎng)期的光譜特征發(fā)現(xiàn)，這3種植被可見(jiàn)光波段的光譜特征變化與不同生長(zhǎng)期相關(guān)。王波等[15]使用高光譜技術(shù)分析東祁的6種灌木植被光譜，通過(guò)篩選敏感波段計(jì)算歸一化植被指數(shù)(NDVI)值和相對(duì)活力指數(shù)(RVI值)，實(shí)現(xiàn)了對(duì)6種灌木植被的分類(lèi)識(shí)別。馬東輝等[16]分析南京6種植被的冠層和落葉光譜特征，發(fā)現(xiàn)6種植被之間存在明顯差異。Schmidt等[17]利用連續(xù)統(tǒng)去除手段凸顯27種鹽沼植被的光譜差異，篩選出6個(gè)特征波段來(lái)區(qū)分識(shí)別這些鹽沼植被。使用光譜儀進(jìn)行實(shí)驗(yàn)設(shè)計(jì)，現(xiàn)場(chǎng)觀察不同植被的有效光譜信息，分析后選擇植被特征波段，可以為結(jié)合高光譜遙感圖像進(jìn)行植被識(shí)別和分類(lèi)提供依據(jù)[18]。目前植被光譜測(cè)量所采用的儀器基本都是ASD光譜測(cè)量?jī)x。相對(duì)ASD地物光譜儀，SOC710VP成像光譜儀能在光譜成像的圖像上自主選擇所需研究范圍，所得結(jié)果更為準(zhǔn)確[19-20]。參考許多文獻(xiàn)發(fā)現(xiàn)，雖然目前對(duì)植被的光譜特性研究較多，但研究時(shí)間多集中在夏季和秋季，而研究時(shí)間選擇冬季的相對(duì)較少，對(duì)滇中植被光譜特性的研究也比較少[21-24]。本文以SOC710VP所測(cè)的高光譜數(shù)據(jù)為數(shù)據(jù)源，使用ENVI5.3和SOC710自帶的處理軟件，以側(cè)柏、樟木、柳杉為研究對(duì)象，對(duì)這3種典型的滇中喬木植被光譜曲線進(jìn)行分析，并使用光譜微分和連續(xù)統(tǒng)去除凸顯其光譜特征，分析3種滇中典型喬木植被的光譜特性與差異性，為植被分類(lèi)和遙感監(jiān)測(cè)提供參考。

1 材料與方法

1.1 實(shí)驗(yàn)區(qū)概況

植被光譜監(jiān)測(cè)實(shí)驗(yàn)選在昆明市西山森林公園，公園的植被多而集中，且受人為活動(dòng)影響較小。西山的海拔為1 890～2 360 m，東臨滇池湖畔，屬于滇中高原典型地段。地理坐標(biāo)為102°37′E～102°38′E，24°59′N(xiāo)～27°57′N(xiāo)，北高南低。昆明西山公園森林資源豐富，植物以喬木為主，也有少量灌木和草本植物，共計(jì)1 086種植物。

1.2 數(shù)據(jù)采集

在試驗(yàn)區(qū)分別選取側(cè)柏、樟木、柳杉3種典型喬木(常綠植被)作為光譜測(cè)試對(duì)象(表1)。2019年12月中旬，在天氣晴、地面能見(jiàn)度高、風(fēng)力小(風(fēng)力等于或小于3級(jí))的情況下，在10：00～15：00利用SOC710VP高光譜儀對(duì)試驗(yàn)區(qū)選擇的樣本進(jìn)行測(cè)量并記錄相關(guān)參數(shù)。選擇長(zhǎng)勢(shì)良好的植被作為實(shí)驗(yàn)樣本。測(cè)量時(shí)應(yīng)保障植被周?chē)哂辛己玫耐ㄒ晽l件，無(wú)陰影遮擋，避免被測(cè)植被和參照板晃動(dòng)，保證所測(cè)植被和參考板在相同的條件下進(jìn)行測(cè)量。在實(shí)際測(cè)量中，記錄每個(gè)樣本的點(diǎn)號(hào)、植被類(lèi)型、位置和其他相關(guān)參數(shù)，并拍攝樣本照片。

表1 三種喬木植被及其特征Table 1 Three kinds of tree vegetation and their characteristics

傳統(tǒng)的光譜掃描儀獲取點(diǎn)數(shù)據(jù)，而SOC710VP則獲取兼具圖像信息和光譜信息的面狀高光譜數(shù)據(jù)。SOC710VP高光譜成像儀可以提供從可見(jiàn)光到近紅外128個(gè)波段的光譜數(shù)據(jù)，從單葉片測(cè)量植物光譜反射率和光合作用變化情況，適用于植被的測(cè)量。此次所測(cè)數(shù)據(jù)光譜范圍為400～1 000 nm，光譜分辨率為1.3 nm。

1.3 數(shù)據(jù)處理

1.3.1 數(shù)據(jù)預(yù)處理

鑒于所有數(shù)據(jù)都是實(shí)時(shí)現(xiàn)場(chǎng)采樣，會(huì)受到許多環(huán)境背景因素的影響，因此，有必要識(shí)別和篩選原始光譜數(shù)據(jù)，剔除有缺陷和無(wú)效的值(比如飽和的數(shù)據(jù)，因?yàn)閿?shù)據(jù)飽和會(huì)損失信息量)。采用S-Golay濾波方法平滑原始光譜曲線，可以在一定程度上抑制光譜抖動(dòng)、噪聲，且基本保留了真實(shí)的光譜特征，保證分析研究時(shí)的數(shù)據(jù)質(zhì)量。使用SOC710自帶的相關(guān)軟件進(jìn)行數(shù)據(jù)標(biāo)定，最后得到植被葉片的光譜反射率。

不同高度的植被葉片生化參數(shù)空間分布不均勻，導(dǎo)致不同位置葉片對(duì)應(yīng)的光譜反射率發(fā)生變化。此外，獲取的高光譜數(shù)據(jù)為面數(shù)據(jù)，為保證實(shí)驗(yàn)數(shù)據(jù)具有代表性，將對(duì)植被葉片進(jìn)行均勻抽樣，計(jì)算平均反射率。利用ENVI軟件，在側(cè)柏、樟木、柳杉不同高度和不同位置(包括冠層、中層和底層)的葉片上分別選取20個(gè)樣本，每個(gè)樣本包括10個(gè)像元，將這些區(qū)域的光譜反射率求平均值，將其作為該植被的光譜反射率。處理的方法流程如下所示(圖1)。

圖1 技術(shù)路線圖Fig.1 Technology roadmap

1.3.2 光譜分析方法

微分法：光譜微分可以在一定程度上消除光譜數(shù)據(jù)之間的系統(tǒng)誤差，減少背景噪聲(如大氣輻射、散射和吸收)對(duì)目標(biāo)光譜的影響，還可以突出顯示光譜的斜率，區(qū)分重疊光譜并提取可識(shí)別地物的光譜吸收峰參數(shù)[25-28]。計(jì)算過(guò)程如下[29]。

(1)

式中：FRDλt為波段t和波段t+1之間的光譜一階微分，Rλt和Rλt+1為相對(duì)應(yīng)的波段t和t+1處的原始光譜反射，λt和λt+1為每個(gè)波段的波長(zhǎng)，Δλ為波長(zhǎng)λt+1到λt的間隔。

連續(xù)統(tǒng)去除法：增強(qiáng)光譜曲線中的吸收和反射特征可以通過(guò)連續(xù)統(tǒng)去除變換有效的突出顯示，并且還能將反射率歸一化，連續(xù)統(tǒng)去除也叫包絡(luò)線去除[30-33]。計(jì)算過(guò)程如下。

(2)

式中：CRi為i處的連續(xù)統(tǒng)去除值，Ri為i處的光譜反射率，RHi為i處直線的光譜反射率。

2 結(jié)果與分析

2.1 原始曲線反射光譜特征分析

從圖2可以看出，側(cè)柏、樟木、柳杉3種植被的原始光譜特征曲線在400～1 000 nm大致符合綠色植被的光譜特征曲線趨勢(shì)，曲線上波峰和波谷、紅邊現(xiàn)象出現(xiàn)的位置基本相同。反射率閾值范圍均在0.03～0.80。側(cè)柏的反射率閾值在0.03～0.38，樟木的反射率閾值在0.03～0.69，柳杉的反射率閾值在0.05～0.80。

圖2 三種喬木植被反射光譜曲線Fig.2 Reflection spectra of three tree species

可見(jiàn)光范圍內(nèi)，葉綠素對(duì)藍(lán)光吸收強(qiáng)，410 nm附近出現(xiàn)第1個(gè)波谷(藍(lán)光吸收帶)，波谷不明顯。此時(shí)柳杉反射率曲線呈下降趨勢(shì)，區(qū)別于側(cè)柏和樟木。550 nm附近因葉綠素對(duì)綠光的反射作用形成反射峰，側(cè)柏、樟木、柳杉的反射率分別為0.12、0.14、0.13，樟木的反射率相對(duì)其他2種植被略高，說(shuō)明樟木具有最強(qiáng)的葉綠素反射能力。670～680 nm的位置因葉綠素對(duì)紅光的強(qiáng)吸收作用，出現(xiàn)紅光吸收帶，也就是第2個(gè)波谷，此時(shí)樟木的反射率最高。680～760 nm反射率急劇上升，出現(xiàn)一個(gè)明顯的“陡坡”，即紅邊特征，此時(shí)除柳杉外，其余植被皆達(dá)到反射率的最高值。

近紅外波段，側(cè)柏、樟木、柳杉的光譜曲線在走勢(shì)上表現(xiàn)出明顯差異。700～760 nm側(cè)柏反射率上升的幅度較為平緩，明顯區(qū)別于柳杉和樟木。760～1 000 nm側(cè)柏、樟木、柳杉的反射率值存在明顯差異。柳杉的反射率持續(xù)上升，側(cè)柏和樟木的反射率持續(xù)下降。柳杉的反射率保持在60%以上，明顯高于其他樹(shù)種，樟木其次，側(cè)柏的反射率值最低，保持在40%以下。

2.2 紅邊波段光譜一階微分分析

“紅邊”現(xiàn)象是植被與其他物質(zhì)光譜特征最大的不同，它是植被在670～760 nm反射率增高最快的點(diǎn)，也是植被光譜曲線最明顯的特征。通過(guò)一階微分可以反映光譜響應(yīng)的變化速率，因此，將側(cè)柏、樟木、柳杉3種植被在紅邊波段的一階微分曲線作為主要的研究對(duì)象。如圖3所示，3種植被的一階微分值代表紅邊波段的陡峭程度，側(cè)柏、樟木、柳杉的一階微分值分別在700、720、730 nm達(dá)到了極大值點(diǎn)。3種植被在極大值點(diǎn)的增速最快，其中，柳杉的增速最快，樟木次之，側(cè)柏的增速最慢。側(cè)柏在700 nm和750 nm分別有1個(gè)峰值，在730 nm處有1個(gè)波谷，具有“雙峰一谷”的特性。樟木在720 nm和765 nm處分別有1個(gè)峰值，在755 nm和775 nm處分別有1個(gè)波谷，具有“兩峰兩谷”的特性。柳杉在730 nm處有1個(gè)峰值，此處的一階微分值高于其他2種植被，在768 nm處有1個(gè)波谷，此處的一階微分值介于側(cè)柏和樟木之間，具有“單峰谷”的特性。

圖3 紅邊反射波段一階導(dǎo)數(shù)曲線Fig.3 The first derivative curve of the red edge reflection band

2.3 連續(xù)統(tǒng)去除光譜特征分析

基于側(cè)柏、樟木、柳杉3種喬木植被的原始光譜曲線進(jìn)行連續(xù)統(tǒng)去除處理。結(jié)果如圖4所示。在400～1 000 nm內(nèi)側(cè)柏、樟木、柳杉3種喬木植被的連續(xù)統(tǒng)去除值都介于0～1。側(cè)柏、樟木、柳杉在近紅外波段均沒(méi)有明顯差異，差異主要表現(xiàn)在可見(jiàn)光波段，明顯凸顯了原始光譜曲線的吸收和反射特征。原始光譜反射率曲線的可見(jiàn)光波段，側(cè)柏、樟木、柳杉3種植被的反射率沒(méi)有明顯的差異，3個(gè)明顯波峰波谷特征的反射率之差小于0.05。連續(xù)統(tǒng)去除處理后，側(cè)柏在可見(jiàn)光波段具有4個(gè)波谷、3個(gè)波峰的特征，其曲線走勢(shì)明顯區(qū)別于樟木和柳杉，較好區(qū)分。樟木和柳杉曲線走勢(shì)基本相似，連續(xù)統(tǒng)去除值在440～760、950～1 000 nm存在差異，在這2個(gè)波段內(nèi)，樟木的連續(xù)統(tǒng)去除值均略高于柳杉。

圖4 連續(xù)統(tǒng)去除變換光譜曲線Fig.4 Continuum removal transform spectrum curve

3 結(jié)論

本文針對(duì)滇中地區(qū)典型喬木植被，基于實(shí)測(cè)高光譜數(shù)據(jù)，分析了側(cè)柏、樟木和柳杉的葉片光譜特征差異，并借助一階微分和連續(xù)統(tǒng)去除法突出不同植被的可區(qū)分光譜，得到了以下結(jié)論：

側(cè)柏、樟木、柳杉3種喬木植被的光譜曲線具有相似的起伏特征，并具有綠色植被在400～1 000 nm的光譜特征，比如可見(jiàn)光波段的“雙峰一谷”和“紅邊”特征。與樟木和柳杉相比，側(cè)柏在波峰和波谷處的反射率最低，在“綠峰”和“紅谷”處樟木的反射率最大，在“藍(lán)谷”處柳杉的反射率最大。反射率在近紅外波段一直保持在一個(gè)較高的值，近紅外波段反射率依次為柳杉>樟木>側(cè)柏。

從3種植被原始光譜曲線分析得出，450～700 nm 3種植被區(qū)分度不大，結(jié)合光譜曲線和光譜反射率可以在400～420 nm和700～1 000 nm將3種植被區(qū)分開(kāi)。

紅邊波段的一階微分變換凸顯了側(cè)柏、樟木、柳杉的不同波峰波谷的差異特征。出現(xiàn)紅邊位置的先后順序依次為側(cè)柏、樟木、柳杉；極大值的大小順序依次為樟木、側(cè)柏、柳杉。從3種植被的一階微分形態(tài)得出，側(cè)柏、樟木、柳杉分別具有“兩峰一谷”“兩峰兩谷”“一峰一谷”的特性。其中，側(cè)柏的波峰波谷更為明顯，識(shí)別度更大。

連續(xù)統(tǒng)去除凸顯了側(cè)柏、樟木、柳杉3種植被原始光譜在可見(jiàn)光波段的吸收和反射特征差異。側(cè)柏的曲線走勢(shì)明顯區(qū)別于樟木和柳杉，能較好地區(qū)別側(cè)柏與另外2種植被的差異。在440～760 nm和950～1 000 nm，樟木和柳杉的連續(xù)統(tǒng)去除值差異明顯。400～440 nm和760～950 nm 3種植被連續(xù)統(tǒng)去除曲線混淆大，區(qū)分效果不理想。原始光譜中3種植被在440～700 nm區(qū)分度不大，但在連續(xù)統(tǒng)去除光譜中區(qū)分效果明顯。