孫世豪，符添瑋

(自然資源部海南基礎(chǔ)地理信息中心，海南 海口 570203)

0 引言

山火作為林地主要災(zāi)害之一，近些年來備受關(guān)注。大興安嶺山火，大涼山山火等嚴(yán)重的山火對(duì)人民生命、財(cái)產(chǎn)安全以及林地資源造成了嚴(yán)重的威脅。山火的防控十分重要，火災(zāi)后對(duì)于山火的影響情況以及后期恢復(fù)情況的監(jiān)測(cè)也十分重要，監(jiān)測(cè)以及各項(xiàng)措施有利于災(zāi)后重建工作以及未來的防控工作。利用遙感技術(shù)范圍覆蓋廣、獲取影像及時(shí)以及技術(shù)較為成熟等特點(diǎn)，可為以上工作提供有效技術(shù)手段[1]。其中，植被指數(shù)作為反應(yīng)植被覆蓋的主流指標(biāo)，用以反演地表林地覆蓋情況。現(xiàn)有的植被指數(shù)有幾十種，因地表覆蓋復(fù)雜，不同監(jiān)測(cè)區(qū)內(nèi)植被與水體種類、分布不同，要根據(jù)不同植被指數(shù)的特點(diǎn)和實(shí)驗(yàn)區(qū)域的實(shí)際情況，選擇相應(yīng)植被指數(shù)[2]。

1 研究區(qū)域與數(shù)據(jù)來源

研究區(qū)域位于海邊的山林，于2015年12月25日發(fā)生山火，植被被嚴(yán)重毀壞，同時(shí)也造成了當(dāng)?shù)厝嗣褙?cái)產(chǎn)、自然資源的嚴(yán)重?fù)p失。隨后幾年，植被逐年恢復(fù)，最終基本恢復(fù)到災(zāi)前水平。這一過程正適用于植被指數(shù)的山火后林地受災(zāi)以及恢復(fù)情況的研究。

實(shí)驗(yàn)數(shù)據(jù)來源于Landsat 8 OLI/TRIS C1 Level-2，其參數(shù)[3]如表1所示。數(shù)據(jù)范圍包含了實(shí)驗(yàn)區(qū)全境，時(shí)相為2015年12月10日、2016年1月11日、2016年9月7日、2017年8月9日、2018年7月27日，其中一個(gè)時(shí)相為山火發(fā)生之前，四個(gè)時(shí)相為山火發(fā)生之后，云量全部小于30%。

表1 Landsat 8 OLI/TRIS C1 Level-2傳感器參數(shù)Tab.1 Landsat 8 OLI/TRIS C1 Level-2 sensor parameters

2 研究方法

本文以Landsat 8 OLI/TRIS C1 Level-2數(shù)據(jù)為基礎(chǔ)，分析山火過后森林資源覆蓋情況以及恢復(fù)情況，主要流程為：

1)數(shù)據(jù)收集與處理；

2)植被指數(shù)的選擇；

3)閾值劃分；

4)計(jì)算受災(zāi)面積；

5)恢復(fù)情況分析。

2.1 數(shù)據(jù)收集與處理

利用Landsat 8 OLI/TRIS C1 Level-2的數(shù)據(jù)，選取相應(yīng)期數(shù)并盡量保證研究區(qū)域無云。將下載后的數(shù)據(jù)進(jìn)行裁剪，減小數(shù)據(jù)量，提升處理速度。根據(jù)影像特點(diǎn)，使用低通濾波和拉普拉斯濾波進(jìn)行影像增強(qiáng)，并進(jìn)行大氣校正。

2.2 植被指數(shù)選擇

根據(jù)本研究區(qū)域特征，選取歸一化植被指數(shù)(NDVI)和比值植被指數(shù)(RVI)進(jìn)行受災(zāi)情況分析和恢復(fù)情況分析。

1)歸一化植被指數(shù)計(jì)算公式：

NDVI=(NIR-R)/(NIR+R)

(1)

式中：NIR為近紅外波段的反射值；R為紅光波段的反射值。 NDVI的取值范圍是，-1<=NDVI<=1，當(dāng)NDVI<0時(shí)，表示地面覆蓋為云、水、雪等，對(duì)可見光高反射；當(dāng)NDVI=0時(shí)表示有巖石或裸土等，NIR和R近似相等；當(dāng)NDVI>0時(shí)，表示有植被覆蓋，且隨覆蓋度增大而增大。NDVI能反映出植物冠層的背景影響，如土壤、潮濕地面、雪、枯葉、粗糙度等，且與植被覆蓋有關(guān)。NDVI可作為標(biāo)準(zhǔn)化方式對(duì)健康植被進(jìn)行量化，并且準(zhǔn)確度較高。對(duì)于經(jīng)過山火的植被，其NDVI值與健康植被的NDVI值相比差距會(huì)較為明顯，可以準(zhǔn)確、高效確定受災(zāi)植被范圍，計(jì)算出受災(zāi)的林地面積。但是，NDVI用非線性拉伸的方式增強(qiáng)了NIR和R的反射率的對(duì)比度，NDVI對(duì)高植被區(qū)具有較低的靈敏度[4-7]。

2)比值植被指數(shù)計(jì)算公式：

RVI=NIR/R

(2)

式中：NIR為近紅外波段的反射值；R為紅光波段的反射值。綠色健康植被覆蓋地區(qū)的RVI遠(yuǎn)大于1，而無植被覆蓋的地面的RVI在1附近。植被的RVI通常大于2；其中，植被覆蓋度是RVI的影響因素之一，當(dāng)植被覆蓋度較高時(shí)，RVI對(duì)植被十分敏感；當(dāng)植被覆蓋度<50%時(shí)，這種敏感性顯著降低；由于研究區(qū)域處于大洋路沿線，氣候條件良好，植被覆蓋度較高。對(duì)于災(zāi)后恢復(fù)植被變化不是很明顯的情況，利用RVI可以彌補(bǔ)NDVI對(duì)于高植被覆蓋敏感度低的不足。同時(shí)，RVI也受大氣條件影響，大氣效應(yīng)大大降低對(duì)植被監(jiān)測(cè)的靈敏度，所以在計(jì)算前需要進(jìn)行大氣校正，或用反射率計(jì)算RVI[6-7]。

2.3 閾值劃分

閾值范圍劃分首先要對(duì)研究區(qū)域的植被指數(shù)進(jìn)行計(jì)算，再根據(jù)不同區(qū)域不同地物種類以及大量實(shí)地?cái)?shù)據(jù)進(jìn)行劃分。在閾值劃分的過程中，需要通過大量實(shí)地?cái)?shù)據(jù)對(duì)其分類的結(jié)果進(jìn)行驗(yàn)證，將目視解譯的結(jié)果與使用植被指數(shù)劃分后的結(jié)果進(jìn)行對(duì)比，如分類結(jié)果不理想，則需要重新劃分閾值或選擇指數(shù)[8]。當(dāng)分類結(jié)果達(dá)到一定的準(zhǔn)確性時(shí)，即可用于后續(xù)研究。

2.4 計(jì)算受災(zāi)面積

利用火災(zāi)后第一時(shí)相的數(shù)據(jù)對(duì)實(shí)驗(yàn)區(qū)域NDVI值計(jì)算結(jié)束后，將計(jì)算結(jié)果進(jìn)行重分類，根據(jù)劃分好的閾值確定山火范圍，通過像素個(gè)數(shù)與像素面積計(jì)算受災(zāi)面積[10]。計(jì)算公式為：受災(zāi)面積=像素個(gè)數(shù)*單個(gè)像素面積。

2.5 恢復(fù)情況分析

為了保證結(jié)果的準(zhǔn)確性，同時(shí)使用NDVI和RVI進(jìn)行林地恢復(fù)情況分析。分別計(jì)算研究區(qū)域火災(zāi)前一時(shí)相和后多時(shí)相數(shù)據(jù)的NDVI與RVI，分析其變化特征，從而反應(yīng)受災(zāi)過后林地恢復(fù)的情況[9]。

3 結(jié)果與分析

3.1 受災(zāi)面積

根據(jù)數(shù)據(jù)源，單個(gè)像素面積為900 m2，像素個(gè)數(shù)為12 701，代入公式可得受災(zāi)面積為11.43 km2。受災(zāi)前后影像及受災(zāi)范圍如圖1所示，其中白色部分為云層，黃色部分為受災(zāi)面積。

圖1 受災(zāi)范圍展示Fig.1 Images before and after disaster & display of affected area

由于除山火影響外其他影響因素(例如人工砍伐)也會(huì)導(dǎo)致林地覆蓋率的變化，因此此方法算出的面積可能存在一定的誤差，需要用目視解譯的方式進(jìn)行確認(rèn)。

3.2 林地恢復(fù)情況分析

山火發(fā)生前，研究區(qū)域的NDVI值較高(0.704-0.947)，大部分都處于0.8以上，這表明此區(qū)域的林地覆蓋率很高；研究區(qū)域的RVI幾乎都大于2(2-37)，這表明此區(qū)域的林地覆蓋率很高。

研究區(qū)域NDVI與RVI展示(山火前，2015年12月)，如圖2所示。

圖2 研究區(qū)域NDVI與RVI展示(山火前，2015年12月)Fig.2 Display of NDVI and RVI in study area before wildfire in December 2015

山火發(fā)生后，受災(zāi)區(qū)域NDVI值明顯降低，部分區(qū)域接近于0(0.028-0.393)，這表明林地破壞嚴(yán)重，近乎裸地。山火發(fā)生后四個(gè)時(shí)相(2016年1月、2016年9月、2017年8月、2018年7月)受災(zāi)區(qū)域的NDVI值范圍是0.028-0.393、0.357-0.585、0.519-0.996、0.666-0.988，NDVI值隨著年份的增加而增加。在山火發(fā)生后的一年之內(nèi)(2016年)，近一半的林地恢復(fù)；接下來的兩年內(nèi)(2017年、2018年)，受災(zāi)林地基本完全恢復(fù)，2018年的NDVI值也同2015年山火發(fā)生前的RVI值接近；山火發(fā)生后，受災(zāi)區(qū)域RVI值明顯降低，大部分區(qū)域處于1-2之間，表明覆蓋率很低，林地破壞嚴(yán)重。山火發(fā)生后四個(gè)時(shí)相(2016年1月、2016年9月、2017年8月、2018年7月)受災(zāi)區(qū)域的RVI值范圍是1-2、1-9、2-24、2-37，RVI值隨著年份的增加而增加。在山火發(fā)生后的三年內(nèi)(2016年、2017年、2018年)，受災(zāi)林地基本完全恢復(fù)，2018年的NDVI值也同2015年山火發(fā)生前的RVI值接近(見圖3)，同NDVI分析結(jié)論一致。

圖3 研究區(qū)域NDVI與RVI展示(山火后，2016年1月、2016年9月、2017年8月、2018年7月)Fig.3 Display of NDVI and RVI in study area after wildfire in January 2016，September 2016，August 2017，and July 2018

在整個(gè)的研究區(qū)域中，NDVI計(jì)算值受云層影響，在圖中可以看到云層部分，但是受災(zāi)區(qū)域無云層覆蓋；雖然NDVI反映各時(shí)相林地覆蓋的變化，但是由于其自身對(duì)于高植被區(qū)靈敏度低的不足，結(jié)果仍需要通過計(jì)算RVI進(jìn)一步驗(yàn)證。

研究區(qū)域依據(jù)RVI值分為三個(gè)等級(jí)：(1，2]，(2，9]和(9，37]，分別對(duì)應(yīng)著低林地覆蓋、中等林地覆蓋和高林地覆蓋。山火發(fā)生后4個(gè)時(shí)相RVI值變化如圖4所示。

圖4 山火發(fā)生后4個(gè)時(shí)相(2016年1月、2016年9月、2017年8月、2018年7月)RVI值變化Fig.4 Changes in RVI values in four phases after wildfire in January 2016，September 2016，August 2017，and July 2018

根據(jù)RVI值變化分析可以得出：低林地覆蓋區(qū)域面積開始有緩慢下降，在經(jīng)過迅速下降后最后趨于平緩，維持在較低值；中等林地覆蓋區(qū)域面積增長(zhǎng)緩慢，維持在中等值；高林地覆蓋區(qū)域面積增長(zhǎng)迅速，最終達(dá)到較高值。這表明火災(zāi)過后，裸地面積或草地面積下降，最終保持在較低值，被林地覆蓋；林地面積逐年增加，最終保持較高值，覆蓋整個(gè)區(qū)域。其反應(yīng)的恢復(fù)情況符合規(guī)律：山火發(fā)生后，草地等首先恢復(fù)，林地逐漸恢復(fù)，最終覆蓋草地。

4 結(jié)束語

本文通過NDVI和RVI指數(shù)分析了山火后林地受災(zāi)與恢復(fù)情況，計(jì)算出了火災(zāi)范圍以及分析了林地的恢復(fù)規(guī)律。得出以下結(jié)論：

1)山火受災(zāi)面積為11.43 km2，受災(zāi)區(qū)域林地基本全部被燒毀，林地資源損失嚴(yán)重。

2)林地恢復(fù)速度較快，災(zāi)后一年恢復(fù)一半，隨后裸地面積下降迅速，第三年基本恢復(fù)原有水平。

3)實(shí)驗(yàn)方法有待提高，需要引入更多適應(yīng)研究區(qū)域特點(diǎn)的指數(shù)，更全面地進(jìn)行分析；在計(jì)算火災(zāi)面積時(shí)，考慮可以提高精度的方法，排除其他影響因素；在分析林地恢復(fù)情況時(shí)，應(yīng)納入降水、溫度等氣候因素。