基于遙感數(shù)據(jù)的林火監(jiān)測與過火面積估算方法研究

1、    基于遙感數(shù)據(jù)的林火監(jiān)測與過火面積估算方法研究    郟江杰+王有志+穆振娟+張宇涵+戴憶如摘要： 野外森林大火每年造成巨大的經(jīng)濟(jì)損失和生態(tài)損失，利用遙感技術(shù)對森林火災(zāi)進(jìn)行監(jiān)測與災(zāi)后受災(zāi)面積的測算已成為重要的技術(shù)手段，本文使用landsat8高分辨率的遙感影像對昆明市西山區(qū)2014年4月16日林火發(fā)生前后的植被覆蓋進(jìn)行了對比，運(yùn)用了三種不同的分類方法對西山火災(zāi)進(jìn)行了研究，統(tǒng)計出各種方法算得的過火面積，通過對比分析三種方法的優(yōu)勢與存在的問題，并對可能存在的誤差進(jìn)行分析。abstract： every year， wild forest fires c

2、aused huge economic losses and ecological loss. using remote sensing technology to monitor forest fire and estimate disaster affected area has become an important method nowadays. in this paper， landsat 8 high spatial resolution remote sensing image were used to compare vegetation coverage before an

3、d after the forest fire on kunming xishan district in april 16， 2014. three different remote sensing classification methods are used to estimate the burned area of xishan fire， through the comparison， problems and advantages of the three methods are analyzed， and the possible reasons of errors are d

4、iscussed in this paper.關(guān)鍵詞： 過火面積；遙感；火災(zāi)key words： burned area；remote sensing；wild fire：s762 ：a ：1006-4311（2017）12-0219-040 引言森林是陸地生態(tài)系統(tǒng)的重要組成部分之一，對全球氣候和生態(tài)環(huán)境有重大的影響作用。森林可以極大地改善人類生存環(huán)境與維護(hù)自然生態(tài)平衡。然而，中國的森林資源十分匱乏，人均占有量不足的全球平均水平的1/4，且我國森林資源區(qū)域分布極度不均勻。除此之外，我國森林資源還經(jīng)常遭受氣候災(zāi)害、病蟲害的困擾。野外森林火災(zāi)是造成森林破壞的主要原因之一，它不僅會給社會造成巨大的經(jīng)

5、濟(jì)損失，而且會對土壤、野生動物和微生物造成巨大的破壞，造成生態(tài)環(huán)境的持續(xù)惡化。據(jù)統(tǒng)計，全球森林火災(zāi)造成的經(jīng)濟(jì)損失約占全球經(jīng)濟(jì)的0.1%，部分火災(zāi)頻發(fā)國家損失甚至超過了經(jīng)濟(jì)收入的0.5%。因此，森林火災(zāi)因其原因復(fù)雜、防治難度大等原因已成為森林資源和自然環(huán)境破壞的主要威脅。及時、準(zhǔn)確、有效地監(jiān)測野生森林火災(zāi)，及時準(zhǔn)確地評估其損失，對于維護(hù)生態(tài)平衡，保護(hù)森林資源具有重要的現(xiàn)實意義。隨著遙感技術(shù)的逐步成熟和快速發(fā)展，遙感技術(shù)因其具有高時空分辨率、覆蓋面廣、成本低等特點，已被廣泛應(yīng)用于野外火災(zāi)監(jiān)測和災(zāi)害評估中。不僅如此，它已成為我國大面積野外火災(zāi)監(jiān)測與過火面積估算的重要手段。本文對比了三種常用的遙感方法

6、：支持向量機(jī)監(jiān)督分類（svm）、pc3閾值，和ndvi閾值，利用landsat 8多光譜遙感數(shù)據(jù)提取昆明市西山區(qū)森林過火面積的信息（研究對象為2014年4月16日昆明市西山區(qū)森林火災(zāi)）。通過對遙感圖像的序列分析和火災(zāi)前后圖像的對比，本文分析了該區(qū)域火災(zāi)的程度并估算了森林過火面積。1 研究對象西山區(qū)位于昆明市西部，地理坐標(biāo)為東經(jīng)102°21°至102°45°，北緯24°'到25°'，東西寬36.8公里，南北寬為72公里。西山區(qū)靠近五華區(qū)和官渡區(qū)，毗鄰安寧市。西山森林火災(zāi)發(fā)生在2014年4月16日。本文使用兩幅landsa

7、t8多光譜遙感圖像數(shù)據(jù)，分別拍攝于西山森林火災(zāi)之前和之后（2013年4月20日和2014年4月23日），作為研究的基礎(chǔ)材料。這些如圖1所示。landsat8是美國nasa于2013年發(fā)射升空的陸地衛(wèi)星（landsat）計劃發(fā)射的第8顆衛(wèi)星，陸地衛(wèi)星主要是以地球陸地的環(huán)境與資源為探測目標(biāo)，為農(nóng)業(yè)、林業(yè)、地質(zhì)礦產(chǎn)、環(huán)境等國民經(jīng)濟(jì)各部門源源不斷地提供信息和資料。 landsat- 8上攜帶有兩個主要載荷：oli和tirs。oli陸地成像儀包括9個波段，空間分辨率為30m，其中包括一個15米的全色波段，成像寬幅為185×185km。2 研究方法2.1 支持向量機(jī)監(jiān)督分類（svm）由于land

8、sat 8多光譜圖像具有11條波段，因此尋找最佳的波段組合來表示森林火災(zāi)在數(shù)據(jù)處理中顯得尤為重要。經(jīng)過多次實驗，我們最終將波段5（nir，0.85-0.88m），波段4（red，0.64-0.67m）和波段3（green，053-0.59m）組合為假彩色影像地圖，反映森林火災(zāi)之前和森林火災(zāi)后的地表用地變化。然后通過目視解譯選擇陸地衛(wèi)星圖像的訓(xùn)練樣本和測試樣本，利用谷歌地球上更高分辨率的圖像進(jìn)行協(xié)助?；蛟诎l(fā)生前，圖像主要有四個土地覆蓋類別：水，植被，不透水和裸露的土地，森林火災(zāi)之后，我們添加了另外的類別過火區(qū)域。經(jīng)過支持向量機(jī)的監(jiān)督分類，得到了西山的過火區(qū)域，然后對分類結(jié)果進(jìn)行精度評價，得到分類

9、結(jié)果的可靠性。2.2 pc3閾值法陳等人2用pc3閾值的方法提取森林火燒跡地信息在火災(zāi)后的tm 5羅源縣形象。主成分分析（pca）方法又稱k-l變換，它將高相關(guān)多波段圖像中的有用信息轉(zhuǎn)化為一些無關(guān)的主成分圖像，從而大大減少了總數(shù)據(jù)量，增強(qiáng)了圖像信息量。 在對波段1，波段2，波段3，波段4，波段5，波段7進(jìn)行 k-l變換后，陳等人發(fā)現(xiàn)第三主成分（pc3）可以增強(qiáng)過火區(qū)的信息。pc3成分表達(dá)式的原始6個波段的非線性方程：值=0.104 *波段 1 + 0.321 *波段2 + 0.072 * 波段3 + 0.788 *波段4-0.004 *波段5-0.509 *波段7在我們的研究中，我們利用lan

10、dsat 8 oli圖像為數(shù)據(jù)源，與tm 5圖像波段的名字是不同的。所以我們研究的方程式是：值=0.104 *波段2 + 0.321 *波段3 + 0.072 *波段4 + 0.788 *波段5 - 0.004 *波段6 - 0.509 *波段7公式表明，波段5和波段7的系數(shù)較高，波段5的系數(shù)為正，而波段7系數(shù)為負(fù)。從原始圖像中可以看出，燃燒區(qū)域的反射率在波段5是低的，而在波段7是相對高的。因此，在pc3閾值在確定過火區(qū)域是有利的，因為它可以提高燃燒區(qū)的信息，盡量減少干擾其他特征信息。2.3 ndvi閾值法envi軟件可用于計算圖像的ndvi指數(shù)（歸一化植被指數(shù)）。歸一化植被指數(shù)的范圍從1到1

11、，負(fù)值表明地面上覆蓋著云、水、雪，0值代表的地區(qū)都是巖石或裸露的土壤，由于近紅外波段和紅光波段大致與裸土和巖石類似。當(dāng)ndvi值為正，該地區(qū)被植被所覆蓋，隨著植被覆蓋的增加，ndvi值也會增加。ndvi圖像如圖2所示。把原始圖像轉(zhuǎn)換成植被指數(shù)ndvi圖像之后，進(jìn)行歸一化處理，在林區(qū)過火區(qū)域植被指數(shù)明顯不同與植被覆蓋較好的地區(qū)。然后應(yīng)用決策樹方法將這兩幅土地圖像劃分為兩類土地利用（植被和非植被）。分類后，變化檢測工具是用來比較兩個圖像之間的差異，對過火面積的的統(tǒng)計也在這一過程中進(jìn)行計算。3 研究結(jié)果3.1 支持向量機(jī)分類監(jiān)督分類結(jié)果圖3顯示了西山森林大火前后圖像的分類結(jié)果。中間深色區(qū)域是被燒毀的

12、區(qū)域，在這個圖像中與其他類明顯不同。結(jié)果表明，2013年4月20日的分類結(jié)果總體準(zhǔn)確率為96.6252%，kappa系數(shù)為0.9480。2014年4月23日分類結(jié)果總體準(zhǔn)確率為97.4164%，kappa系數(shù)為0.9652。的混淆矩陣顯示，不透水表面分類結(jié)果的準(zhǔn)確性相對低于其他類別，植被和過火面積分類的結(jié)果的準(zhǔn)確性更好并接近100%。支運(yùn)用持向量機(jī)監(jiān)督分類計算過火面積，從直方圖的分類結(jié)果來看，過火區(qū)域面積為1089.18公頃。3.2 pc3閾值法分類結(jié)果在圖4左圖中，過火區(qū)域在pc3圖像中顯示為暗像元， pc3數(shù)據(jù)的值范圍從4464.94到22893.44。經(jīng)過不斷調(diào)整，pc3的閾值設(shè)為768

13、7，低于7687的像元表示過火區(qū)域，高于7678的像元則為其他區(qū)域，圖4右圖展示出分類結(jié)果。根據(jù)pc3閾值法結(jié)果的統(tǒng)計，燒傷區(qū)面積為1051.2公頃。3.3 ndvi閾值法分類結(jié)果根據(jù)不同的土地覆蓋類型顯示不同的ndvi值，實驗使用決策樹進(jìn)行分類，通過反復(fù)調(diào)整ndvi值達(dá)到最佳的分類結(jié)果，最后確定0.232是第一個圖像的最佳閾值。當(dāng)ndvi大于0.232，土地覆蓋類型分為植被。同樣，第二個圖像的閾值為0.248。圖5是火災(zāi)前后決策樹分類的結(jié)果。對兩個結(jié)果進(jìn)行變化檢測后發(fā)現(xiàn)，由于季節(jié)變化。一些周邊植被成為其他類型的土地利用和一些轉(zhuǎn)向植被。因此，我們的地圖大小調(diào)整為170 * 170像素地圖完全覆

14、蓋整個燒傷面積，用于排除干擾。以下圖像顯示了重新裁剪后的分類結(jié)果（見圖6）。經(jīng)過變化檢測后結(jié)果顯示（見表4），過火面積占整個圖像的14363像素，實際過火面積為1292.67公頃。4 結(jié)果分析pc3閾值方法性能總體更優(yōu)，用pc3閾值法得到區(qū)域的面積相當(dāng)接近于支持向量機(jī)監(jiān)督分類的結(jié)果。然而，仍然存在一些問題。首先，比較pc3閾值法的分類結(jié)果與原始圖像可以發(fā)現(xiàn)，過火面積的邊緣容易被誤分為未燃區(qū)。森林火災(zāi)在邊緣區(qū)域不如中部地區(qū)那么嚴(yán)重，所以過火區(qū)域的像元可能會與未著火的區(qū)域像元有一定程度的混合。其次，一些其他像元誤被歸類為過火面積。根據(jù)監(jiān)督分類結(jié)果，這些錯分像元一般零星分布位于水域和植被（林）區(qū)。這

15、些零星像素的植被被歸類為燒傷面積的原因可能是，是在這些地點的植被正遭受嚴(yán)重的植物病蟲害，這種情況可能會導(dǎo)致植被葉綠素含量變低，從而在近紅外波段（波段5）的反射率將低于健康植被。對于被劃分為過火面積的水體的像元，在近紅外波段（5波段）的反射率與過火區(qū)域相當(dāng)。因此，使用p c3組件為代表的原有的6波段的線性無關(guān)的方程后，對水的像素值和燒傷面積的差異可能很小，有些水域可誤分為過火區(qū)域。對ndvi閾值法的結(jié)果，利用ndvi得到的的過火區(qū)域比監(jiān)督分類和pc3閾值結(jié)果相對較高。一個可能的原因是一些土地的樹木分布稀疏或像素結(jié)合植被和非植被區(qū)，所以他們的ndvi不夠高，這些地區(qū)可能很容易被歸類為非植被區(qū)。另一

16、種假設(shè)的原因是，雖然面積小，但也有一些地區(qū)是由于季節(jié)變化，原來的植被覆蓋逐漸變成裸露的土壤，而不是受森林火災(zāi)的影響。這兩個原因可以導(dǎo)致植被到非植被的像元增加。5 結(jié)論利用支持向量機(jī)監(jiān)督分類方法、pc3閾法和ndvi閾值法從landsat 8多光譜遙感數(shù)據(jù)中提取信息并估算昆明西山區(qū)森林火災(zāi)的過火面積在，結(jié)果分別為1089.18公頃， 1051.2公頃和1292.67公頃。pc3閾值方法獲得的過火區(qū)域的面積相當(dāng)接近支持向量機(jī)監(jiān)督分類的結(jié)果，而用ndvi閾值法獲得的過火區(qū)域面積比監(jiān)督分類和pc3閾值結(jié)果相對較高。比較不同類型的森林覆蓋面積估算方法，有助于提高森林植被調(diào)查結(jié)果的準(zhǔn)確性，維護(hù)生態(tài)平衡和保

17、護(hù)森林資源。 參考文獻(xiàn)：1jiao lin lin， et al. using tm to validate burnt area of large forest fires based on modis： a study in heilongjiang province.2chen， ben-qing， and han-qiu xu. the application of remote sensing in extraction of forest burnt scar information j. journal of fuzhou university （natural sciences

18、 edtion） ，2001（2）： 006.3張兵，韋銳，宋昆侖&趙彥博.基于 tm8 遙感影像的森林火災(zāi)監(jiān)測與過火面積測算j.礦山測量，2015（3）：75-78.4chuvieco， emilio， and russell g. congalton. application of remote sensing and geographic information systems to forest fire hazard mappingj. remote sensing of environment，1989，29 （2）： 147-159.5chuvieco， emilio，

19、et al. "development of a framework for fire risk assessment using remote sensing and geographic information system technologiesj. ecological modelling，2010，221（1）： 46-58.6yan， gao， et al. "comparison of pixel-based and object-oriented image classification approachesa case study in a coal f

20、ire area， wuda， inner mongolia， chinaj. international journal of remote sensing，2006，27 （18）： 4039-4055.7roy， d. p.， et al. prototyping a global algorithm for systematic fire-affected area mapping using modis time series dataj. remote sensing of environment，2005，97（2）： 137-162.8roy， d. p.， lewis， p.

21、 e.， & justice， c. o. burned area mapping using multi-temporal moderate spatial resolution dataa bi-directional reflectance model-based expectation approachj. remote sensing of environment， 2002，83（1）：263-286.9white， joseph d.， et al. remote sensing of forest fire severity and vegetation recover

22、yj. international journal of wildland fire，1996，6 （3）： 125-136.10fraser， r. h.， z. li， and j. cihlar. hotspot and ndvi differencing synergy （hands）： a new technique for burned area mapping over boreal forestj. remote sensing of environment，2000，74 （3）： 362-376.11melgani， f.， moser， g.， & serpico

23、， s. b. unsupervised change-detection methods for remote-sensing imagesj. optical engineering， 2002，41（12）： 3288-3297.12kasischke， e. s.， & french， n. h. locating and estimating the areal extent of wildfires in alaskan boreal forests using multiple-season avhrr ndvi composite dataj. remote sensin