1、 基于高分?jǐn)?shù)據(jù)木里煤礦植被覆蓋度遙感估算 索紅燕 黃義忠 李烙摘 要： 木里煤礦是青海省最大的露天煤礦，該煤礦的大量開采影響了木里鎮(zhèn)的生態(tài)平衡，本文通過估算木里煤礦區(qū)域的植被覆蓋度驗(yàn)證該區(qū)域的生態(tài)情況，經(jīng)估算，木里煤礦區(qū)域20162018的年平均植被覆蓋度分別為0.22、0.43、0.54。從計(jì)算結(jié)果表明木里煤礦區(qū)域的植被覆蓋度在逐年增加，恢復(fù)情況較好。通過研究木里煤礦近3年的植被覆蓋情況可以反映出礦區(qū)的恢復(fù)與治理效果較好，同時(shí)為其它礦區(qū)的恢復(fù)與治理提供參考。Key：木里煤礦;像元二分法;植被覆蓋度： TP751：ADOI：10.3969/j.issn.1003-6970.2019.09.03

2、5本文著錄格式：索紅燕，黃義忠，李烙. 基于高分?jǐn)?shù)據(jù)木里煤礦植被覆蓋度遙感估算J. 軟件，2019，40（9）：153-155Remote Sensing Estimation of Vegetation Coverage in Muli Coal Mine Based on High-Fraction DataSUO Hong-yan， HUANG Yi-zhong， LI Luo（College of Land and Resources Engineering， Kunming University of Science and Technology， Kunming650051，Chi

3、na）【Abstract】： Muli Coal Mine is the largest open-pit coal mine in Qinghai Province. The ecological balance of Muli Town is affected by the large-scale mining of the coal mine. This paper validates the ecological situation of the area by estimating the vegetation coverage of Muli Coal Mine. It is es

4、timated that the annual average vegetation coverage of Muli Coal Mine from 2016 to 2018 is 0.22， 0.43 and 0.54， respectively. The calculation results show that the vegetation coverage in Muli coal mine area is increasing year by year， and the recovery situation is good. By studying the vegetation co

5、ver of Muli Coal Mine in recent three years， it can reflect that the effect of restoration and treatment of mining area is better， and provide reference for the restoration and treatment of other mining areas.【Key words】： Muli Coal Mine; Pixel dichotomy; Vegetation coverage0引言植被覆蓋度表示在植被覆蓋的區(qū)域內(nèi)植被的根莖葉經(jīng)

6、過投影與植被的比例，是評(píng)價(jià)植被生長(zhǎng)狀況的重要因素，通常用在評(píng)價(jià)海島海岸、區(qū)域環(huán)境中的生態(tài)敏感性、脆弱性評(píng)價(jià)。目前為止，植被覆蓋度的計(jì)算方法較多，傳統(tǒng)的方法為目估法1、樣方法2、樣點(diǎn)法3等，傳統(tǒng)方法與現(xiàn)代計(jì)算植被覆蓋度的方法相比較，具有誤差大，成本高，野外操作不便等缺點(diǎn)。近年來，遙感領(lǐng)域科技發(fā)展迅速，技術(shù)不斷提高，使得植被覆蓋度的計(jì)算簡(jiǎn)單易操作，目前，監(jiān)測(cè)植被覆蓋度的方法有植被指數(shù)法4-5，經(jīng)驗(yàn)?zāi)Ｐ头?-7、像元二分法8等。其中，植被指數(shù)法是通過傳感器接受到的光譜信息經(jīng)過波段計(jì)算得到植被指數(shù)圖，此方法通過ENVI軟件計(jì)算，人為干預(yù)較少，計(jì)算的結(jié)果與植被覆蓋度呈現(xiàn)較好的相關(guān)性。經(jīng)驗(yàn)?zāi)Ｐ头ㄊ峭ㄟ^外業(yè)

7、選取樣點(diǎn)計(jì)算值被覆蓋度并與遙感光譜信息建立線性回歸，以點(diǎn)代面推廣大整個(gè)研究區(qū)估算研究區(qū)植被覆蓋度，此方法不易操作，誤差較大，像元二分法是由李苗苗等人提出計(jì)算植被覆蓋度的一種方法，也是目前計(jì)算植被覆蓋度最新的方法。經(jīng)多次試驗(yàn)，此方法計(jì)算植被覆蓋度與實(shí)測(cè)結(jié)果相吻合，也是研究者研究植被覆蓋度9常用的一種方法之一。木里煤礦的開采導(dǎo)致了木里鎮(zhèn)的生態(tài)環(huán)境日益惡化，這一問題引起了國(guó)內(nèi)外的廣泛關(guān)注。國(guó)家制定了相關(guān)法律政策，木里煤礦停業(yè)整頓，未破壞的生態(tài)環(huán)境進(jìn)行保護(hù)，已經(jīng)破壞的生態(tài)環(huán)境進(jìn)行恢復(fù)治理。1研究區(qū)概況木里煤礦位于天峻縣木里鎮(zhèn)，隸屬于青海省海西藏族蒙古族自治州，木里煤礦地理坐標(biāo)為東經(jīng)990800，北緯3

8、80745。木里鎮(zhèn)總面積為2856.664平方公里，研究區(qū)面積為412.261平方千米，占木里鎮(zhèn)面積的14.43%，該區(qū)域處于青藏高原北緣，研究區(qū)海拔范圍為41004200m，屬于高寒區(qū)，受青藏高原氣流環(huán)流影響，常年氣溫溫差大，年平均溫度為4，最高氣溫為22，最低氣溫為零下35.6，該區(qū)域四季區(qū)分不明顯，夏季較短冬季長(zhǎng)，木里鎮(zhèn)總?cè)丝诩s0.11萬(wàn)，居住人口大部分為藏族，藏族占總?cè)丝诘?5%以上。2研究?jī)?nèi)容及方法2.1數(shù)據(jù)處理“高分一號(hào)”衛(wèi)星是我國(guó)發(fā)射的首顆對(duì)地觀測(cè)衛(wèi)星，該衛(wèi)星具有較高分辨率，由全色波段和多光譜組成，“高分一號(hào)”數(shù)據(jù)來源于中國(guó)資源衛(wèi)星中心10，獲取到的數(shù)據(jù)通過ENVI軟件輻射定標(biāo)、

9、大氣校正、去燥處理911等一系列的預(yù)處理之后選用紅外波段與近紅外波段進(jìn)行計(jì)算，計(jì)算完植被指數(shù)以后，通過像元二分模型計(jì)算植被覆蓋度，最后將計(jì)算的結(jié)果在ARCGIS軟件中矢量裁剪得到研究區(qū)影像。2.2像元二分法2.2.1 歸一化植被指數(shù)歸一化植被指數(shù)是反映地表植被覆蓋的指標(biāo)，他的定義為是近紅外波段與紅外波段之差與近紅外與紅外波段之和的商，也可以表示為NDVI=（rNIR rRED）/（rNIR+ rRED），它較好的反應(yīng)植被生長(zhǎng)狀況以及空間覆蓋。根據(jù)植被指數(shù)計(jì)算公式與“高分一號(hào)”衛(wèi)星特點(diǎn)選取（表1）近紅外波段為b4，紅外波段為b2，計(jì)算植被指數(shù)。2.2.2 像元二分法原理像元二分法是根據(jù)線性模型原

10、理提出的一種方法，本文在計(jì)算植被覆蓋度時(shí)采用像元二分法，該方法假定傳感器觀測(cè)到的信息植被完全覆蓋Rs與無(wú)植被覆蓋Rv通過線性加權(quán)計(jì)算像元R的值，即R=Rs+Rv （1）對(duì)于傳感器所觀察的像元信息包括兩種，一種為混合像元，該像元包含了傳感器能觀測(cè)到地物的所有信息，另一種為單一像元，它包含一種地物信息，也被稱為純像元，在實(shí)際操作中傳感器觀測(cè)到的像元都為混合像元，為了計(jì)算方便，假定像元的反射率為植被覆蓋Rveg與無(wú)植被覆蓋Rsoil兩種，混合像元由植被與土壤組成，植被完全覆蓋面積所占像元面積比例為fc，土壤覆蓋比例為1-fc，計(jì)算像元R值時(shí)，混合像元的植被部分為植被的貢獻(xiàn)量Rveg與植被覆蓋面積的乘

11、積fc與無(wú)植被覆蓋的貢獻(xiàn)量與土壤覆蓋面積的乘積之和，如下：R=Rveg*fc+Rsoil*（1fc） （2）整理公式（2）得到的植被覆蓋度公式為：Vf=（RRsoil）/（RvegRsoil） （3）經(jīng)多項(xiàng)研究表明，植被指數(shù)帶入公式（3）中，植被覆蓋度計(jì)算的誤差較小，則植被覆蓋度的公式可表達(dá)為：Vf=（NDVINDVIveg）/（NDVIveg+NDVIsoil） （4）則，NDVI表示植被指數(shù)圖，NDVIveg表示完全被植被覆蓋區(qū)域的NDVI值，NDVIsoil表示裸地區(qū)域的NDVI值。2.2.3 NDVIveg和NDVIsoil參數(shù)計(jì)算通過Reference分析，NDVIveg和NDVIs

12、oil這兩個(gè)參數(shù)取值計(jì)算各有不同，李苗苗在計(jì)算密云水庫(kù)上游的植被覆蓋度時(shí)提出像元二分法把植被覆蓋度應(yīng)該分為兩種情況，一種為植被覆蓋情況最大值接近100%，最小值接近0時(shí)NDVIveg=NDVImax，NDVIsoil=NDVImin，另一種為植被覆蓋度最大值最小值不能近似為100%與0時(shí)根據(jù)置信度值計(jì)算植被覆蓋度， 穆少杰10計(jì)算內(nèi)蒙古植被覆蓋度時(shí)同樣采用了像元二分法，植被完全覆蓋區(qū)為NDVIveg，裸土地為NDVIsoil何海燕11在基于雷達(dá)數(shù)據(jù)估算植被覆蓋度時(shí)通過設(shè)置置信度值排除水體、大于25的山坡、陰影、不透水面并確定RVIveg和RVIsoil。本文根據(jù)天峻縣的土地利用現(xiàn)狀圖12-13

13、確定植被覆蓋區(qū)域以及裸土地的分布情況選取植被指數(shù)置信度在95%與5%的值分別為NDVIveg與NDVIsoil。2.2.4 研究區(qū)植被覆蓋度估算通過上述計(jì)算分析的植被指數(shù)圖與土地利用現(xiàn)狀圖相結(jié)合，統(tǒng)計(jì)置信度95%出現(xiàn)的NDVI值與5%出現(xiàn)的NDVI值，NDVI的頻率統(tǒng)計(jì)結(jié)果如表2所示。根據(jù)表2統(tǒng)計(jì)結(jié)果，確定NDVIveg和NDVIsoil值，通過公式（4）計(jì)算植被覆蓋度圖如下圖所示。2.3結(jié)果分析通過計(jì)算，木里煤礦植被指數(shù)2016的平均植被覆蓋度分別為0.22，2017年的平均植被覆蓋度為0.43，2018年的平均植被覆蓋度為0.54。該礦區(qū)的植被在逐年恢復(fù)，2017年比2016年增長(zhǎng)了48%

14、，2018年比2017年增長(zhǎng)了20%，根據(jù)3年的植被覆蓋度分析，2018年的增長(zhǎng)速率相對(duì)于2017年較緩慢，從總體分析，這3年的植被覆蓋度呈逐漸增長(zhǎng)的趨勢(shì)。3結(jié)論研究區(qū)運(yùn)用像元二分法對(duì)“高分一號(hào)”多光譜數(shù)據(jù)估算近3年的植被覆蓋度，得出了如下結(jié)論：（1）礦區(qū)的開發(fā)嚴(yán)重影響植被的生長(zhǎng)，礦區(qū)開采區(qū)，勘探區(qū)，礦物原料區(qū)植被覆蓋為零。從礦區(qū)向外延伸，植被覆蓋在逐漸增加，充分表明，礦物的開采嚴(yán)重影響了植被的生長(zhǎng)，破壞了生態(tài)平衡。（2）從研究區(qū)3年的植被覆蓋度情況分析，植被覆蓋度逐年上升，研究表明，木里煤礦的生態(tài)環(huán)境生態(tài)環(huán)境正在逐年恢復(fù)，也說明當(dāng)?shù)卣诜e極響應(yīng)黨的號(hào)召，積極治理礦區(qū)的生態(tài)環(huán)境。（3）通過

15、“高分一號(hào)”估算植被覆蓋度信息，表明可以使用“高分一號(hào)”多光譜數(shù)據(jù)可以進(jìn)行大面積監(jiān)測(cè)植被生長(zhǎng)，也可以用在林業(yè)監(jiān)測(cè)、土壤侵蝕、土地荒漠化等領(lǐng)域。（4）NDVIveg和NDVIsoil參數(shù)的選取具有主觀性，參數(shù)的取值會(huì)影響植被覆蓋度的計(jì)算精度。Reference章文波， 符素華， 劉寶元. 目估法測(cè)量植被覆蓋度的精度分析J. 北京師范大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2001， 37（3）： 402-408.趙春玲， 李志剛， 呂海軍. 中德合作寧夏賀蘭山封山育林育草項(xiàng)目區(qū)植被覆蓋度監(jiān)測(cè)J. 寧夏農(nóng)林科技， 2000（s1）： 6-14.張光輝， 梁一民. 黃土丘陵區(qū)人工草地徑流起始時(shí)間研究J. 水土保持

16、學(xué)報(bào)， 1995（3）： 78-83.楊勝天， 劉昌明， 楊志峰等. 南水北調(diào)西線調(diào)水工程區(qū)的自然生態(tài)環(huán)境評(píng)價(jià)J. 地理學(xué)報(bào)， 2002， （1）： 11-18.ABED M A ， ZHOU S ， YAGHI A. Application of GIS and remote sensing in soil degradation assessments in the Syrian coastJ. JOURNAL OF ZHEJIANG AGRICULTURAL UNIVERSITY（AGRIC.& LIFE SCI.）， 2000.池宏康. 沙地油蒿群落覆蓋度的遙感定量化研究J. 植物生態(tài)學(xué)報(bào)， 2000， 24（4）.PUREVDORJT ，TATEISHI R ， ISHIYAMA T . Relationships between percent vegetation cover and vegetation indicesJ. International Journal of Remote Sensing， 1998， 19（18）： 3519-3535.李苗苗， 吳炳方， 顏長(zhǎng)珍. 密云水庫(kù)上游植被覆蓋度的遙感估算J. 資源科學(xué)， 2004（4）.劉鑫淼， 康朝紅， 薛樂樂.基