林 鵬

（甘肅省地質(zhì)礦產(chǎn)勘查開發(fā)局第一地質(zhì)礦產(chǎn)勘查院，甘肅 天水 740120）

1 概述

林火，對(duì)自然和經(jīng)濟(jì)社會(huì)都有多方面的影響。對(duì)于自然環(huán)境方面，過(guò)火后的森林生態(tài)環(huán)境受到影響[1]，動(dòng)植物不同程度的死亡，棲息地遭到破壞[2]，大量氣體排到空氣中。一方面有破壞作用，另一方面更新了森林生態(tài)系統(tǒng)，積累了新的養(yǎng)分，孕育了新的生命。對(duì)于社會(huì)經(jīng)濟(jì)方面，救火投入大量人力物力，造成一定經(jīng)濟(jì)損失，并且污染環(huán)境?；馃E地。

近些年來(lái)遙感技術(shù)因其成本低且可大范圍監(jiān)測(cè)過(guò)火區(qū)的火情被廣泛用于林火監(jiān)測(cè)[3]。陳本清、徐涵秋等人利用TM 影像選取三種不同的提取方法進(jìn)行精度比較并分析了影響精度的原因[4]。祖笑峰、覃先林等人利用高分一號(hào)衛(wèi)星及其計(jì)算出的光譜指數(shù)構(gòu)建了森林火燒跡地識(shí)別決策樹模型[5]。楊偉等人基于MODIS 數(shù)據(jù)建立了一套火燒跡地提取算法并對(duì)精度進(jìn)行了驗(yàn)證[6]。2020 年3 月30 日，西昌瀘山發(fā)生森林火災(zāi)，造成19 名滅火人員犧牲，火勢(shì)危及城區(qū)。本文以西昌瀘山為研究區(qū)，分別選擇決策樹分類法和NBR 提取火燒跡地，對(duì)比兩種提取方法的精度。

2 研究區(qū)概況與數(shù)據(jù)來(lái)源

2.1 研究區(qū)概況

瀘山位于四川省西昌市西南，海拔2317m，緯度27°46′51″-21°53′21″，經(jīng)度102°11′55″-102°17′55″，面積29km2，山脈走勢(shì)為北西—南東向，東南端與螺髻山接壤，西北端為安寧河谷，北東坡緊連邛海，與邛海一起構(gòu)成川西南一大景區(qū)。地形為北坡緩南坡陡。屬于亞熱帶高原季風(fēng)性氣候，平均降水量1020mm，平均氣溫17.5℃。研究區(qū)內(nèi)林業(yè)及資源豐富。如圖1 所示。

圖1 研究區(qū)地理位置示意圖（災(zāi)后遙感影像）

2.2 數(shù)據(jù)來(lái)源與預(yù)處理

Sentinel-2 衛(wèi)星作為歐洲SPOT 衛(wèi)星和美國(guó)陸地衛(wèi)星（Landsat-7）的延續(xù)衛(wèi)星，其擔(dān)任著生物物理變化制圖、監(jiān)測(cè)海岸和內(nèi)陸地區(qū)以及風(fēng)險(xiǎn)和災(zāi)害制圖等職責(zé)[7]。衛(wèi)星包括Sentinel-2A 和Sentinel-2B 兩顆衛(wèi)星，單顆衛(wèi)星的重返周期為10d，兩顆組合重返周期為5d。其攜帶的多光譜成像儀的成像幅寬為290km，光譜分辨率為15～180nm，空間分辨率分為10m、20m、60m[8]，其波段情況見表1。

表1 Sentine-2 衛(wèi)星波段情況

在歐空局下載兩景數(shù)據(jù)，災(zāi)前選擇S2B_MSIL2A_20200330T034529_N0214_R104_T48R TR_20200330T080145.SAFE,云量0.008709，災(zāi)后選擇S2B_MSIL2A_20200409T034529_N0214_R104_T4 8RTR_20200409T080223.SAFE，云量為9.9593，手動(dòng)選擇瀘山區(qū)域，數(shù)據(jù)基本無(wú)云覆蓋。在SNAP 中進(jìn)行重采樣，將20m、60m 的波段重采樣為10m 分辨率，每個(gè)地類選取20 個(gè)樣本，共100 個(gè)。70 個(gè)用于訓(xùn)練樣本，30 個(gè)用于驗(yàn)證樣本。

3 研究方法

3.1 典型地物光譜特征曲線分析

Sentinel-2 數(shù)據(jù)在紅邊范圍含有三個(gè)波段，可用于監(jiān)測(cè)植被健康狀況。從圖2 分析健康植被與火燒跡地，火燒跡地在紅（4）波段略高于健康植被，在紅邊波段（5.6.7）較健康植被大幅下降，在短波紅外（11.12）波段又呈現(xiàn)增長(zhǎng)趨勢(shì)，并且高于健康植被，這是因?yàn)榈乇砣紵龑?dǎo)致植被水分減少，短波紅外的反射率增強(qiáng)幅度要高于可見光波段的反射率。

圖2 典型地物光譜特征曲線圖

3.2 遙感指數(shù)提取

3.2.1 NBR

NBR 是衡量火災(zāi)嚴(yán)重程度的指標(biāo)之一，歸一化燃燒率(NBR)［9］計(jì)算公式如下所示:

式中:NIR、SWIR 分別為近紅外和短波紅外的反射率,即Sentinel-2 的第8、12 波段,中心波長(zhǎng)分別為0.865μm 和2.190μm.采用火災(zāi)前后兩幅影響的NBR 指數(shù)進(jìn)行差值處理能得到DNBR 指數(shù)。DNBR計(jì)算公式如下：

3.2.2 NDVI

歸一化植被指數(shù)（NDVI）是反映地表植被動(dòng)態(tài)及其空間狀況的重要指標(biāo)之一，它與不同地理區(qū)域和生態(tài)系統(tǒng)的生物量或年度地上凈初級(jí)生產(chǎn)力（ANPP）之間呈現(xiàn)正相關(guān)關(guān)系[10]，NDVI 可由近紅外波段的反射值與紅光波段的反射值之差比上兩者之和獲得，計(jì)算公式如下：

式中:NIR、R 分別為近紅外波段和紅光波段的反射率,即Sentinel-2 的第8、4 波段,中心波長(zhǎng)分別為0.865μm 和0.665 μm。

3.3 決策樹分類

提取火燒跡地傳統(tǒng)方法是利用NBR 指數(shù)，提取出來(lái)之后部分水體植被及裸地被錯(cuò)分到火燒跡地中，選取mndvi 指數(shù)，剔除水體，選取ndvi 指數(shù)，剔除植被，觀察裸地及火燒跡地的波段特征，選取波段8 進(jìn)行分類。將大于0.1 的歸為裸地，最終分出火燒跡地。決策樹分類圖如圖3 所示。

圖3 火燒跡地決策樹示意圖

4 結(jié)果與分析

4.1 指數(shù)提取

根據(jù)上述研究方法對(duì)遙感影像進(jìn)行指數(shù)提取，NBR 指數(shù)提取如圖4 所示，白色部分代表火燒跡地，但是影像右側(cè)（深色輪廓部分）被錯(cuò)分到火燒跡地中。NDVI 指數(shù)如圖5 所示，黑色采用部分代表水體。

圖4 DNBR 指數(shù)提取結(jié)果圖（白色部分為火燒跡地）

圖5 NDVI 提取結(jié)果圖（黑色部分為水體）

用OTSU 方法進(jìn)行閾值分割，反復(fù)實(shí)驗(yàn)后選擇閾值0.4，得到如圖所示的火燒跡地提取。但是仍有部分植被、裸地、水體被錯(cuò)分到了火燒跡地中。如圖6 所示。

圖6 DNBR 閾值分割后的提取結(jié)果

4.2 決策樹分類結(jié)果

綜合考慮地物錯(cuò)分采用圖7 的決策樹分類模型，得到如圖7 所示的分類結(jié)果。可以看出錯(cuò)分結(jié)果有所改善。

圖7 決策樹分類結(jié)果圖（紅色部分為火燒跡地）

4.3 精度檢驗(yàn)驗(yàn)證

采用混淆矩陣的方法對(duì)提取火燒跡地的結(jié)果進(jìn)行了定量估算。對(duì)于決策樹分類法，通過(guò)建立混淆矩陣，得出如表2 所示的結(jié)果，總體精度達(dá)到84.52%，Kappa 系數(shù)為0.87，通過(guò)驗(yàn)證。見表2。

表2 精度檢驗(yàn)

5 結(jié)論

本文通過(guò)對(duì)2020 年3 月30 日西昌瀘山火災(zāi)前后Sentinel-2 影像的分析，利用DNBR 指數(shù)和決策樹分類法對(duì)過(guò)火區(qū)域進(jìn)行火燒跡地面積提取。（1）通過(guò)目視驗(yàn)證，DNBR 指數(shù)提取之后再使用閾值分割可以提高提取精度，一些錯(cuò)分的地類可以被去除。（2）基于波段特征和植被指數(shù)的決策樹分類法通過(guò)了驗(yàn)證。精度達(dá)到84.52%，Kappa 系數(shù)為0.87。研究表明，結(jié)果表明，利用DNBR 指數(shù)和決策樹分類法對(duì)過(guò)火區(qū)域進(jìn)行火燒跡地面積提取是可行的。