word基于MODIS数据的汶川地震灾区植被覆盖度估算与动态变化分析摘要：本文以MODIS遥感数据为根底，利用像元二分模型对汶川地震灾区2008年-2013年的植被覆盖度进展监测估算，分析了植被覆盖度在这6年里的动态变化通过从整体角度分析不同植被覆盖度的面积变化和从平均植被覆盖度的角度分析平均植被覆盖度的动态变化过程，均得出植被的动态变化过程：地震前后植被受到严重破坏，低、中低植被覆盖度的面积显著增加，中高、高植被覆盖度的面积显著减少；随着植被的逐渐恢复，至2013年各植被覆盖度区域面积几乎恢复到地震前的水平同时，通过分析植被破坏与恢复前后不同植被覆盖度的面积变化得到植被的动态变化情况关键字：植被覆盖度；像元二分模型；MODIS；汶川地震引言植被，是覆盖在地面的所有植物群的总称，包括森林、灌丛、草地、农作物等等，具有截流降雨、减缓径流、防沙治沙、保持水土等功能，是生态系统的重要组成局部[1-2]植被覆盖度〔Fractional Vegetation Cover, FVC〕是指植被的枝、茎、叶垂直向下的投影面积占统计区域总面积的百分比[3]，是刻画地面植被覆盖的一个重要参数一个地区的植被覆盖情况在很多领域都要涉与，也是其研究的重要指标之一，是重要的根底数据，因此建立快速的、准确的植被覆盖度计算方法与模型在植被与其他相关领域的研究都具有十分重要的意义。
植被覆盖度的测量计算有很多的方法与模型[10,13-16]，主要是从两个大的方面来进展计算，一种是相对传统的测量方法，例如目估法、照片分析法、经验分析法等；另一种是采用建立物理的模型，通过对某地区遥感影像数据进展统计分析，最后估算出植被覆盖度的大小，例如像元分解模型法、回归模型法等归一化植被指数〔NDVI〕与植被覆盖度之间有很大的关系，根据两者之间的关系，通过计算遥感影像的NDVI值来估算植被覆盖度的大小其中，NDVI=〔NIR-R〕/〔NIR+R〕，NIR表示近红外波段的反射率，R表示红光波段的反射率[11]以省汶川地震灾区为研究区域，通过收集该地区2008年-2013年的同一时期MODIS遥感影像图，借助NEVI4.8遥感影像处理软件和ArcGIS10.2遥感影像分析软件，提取研究区的归一化植被指数〔NDVI〕，再根据NDVI与植被覆盖度之间的关系估算得到植被覆盖度大小灰度分布图根据所得到的植被覆盖度大小分布图，分析不同覆盖度在2008年-2013年之间的动态变化情况1、研究区概况汶川地区位于盆地西部边缘，阿坝藏族羌族自治州东南部，地处北纬30°45′-31°43′，东经102°51′-103°44′，东西宽84公里，南北长105公里，是平原与川西北高原的交通要道，素有阿坝州南大门之称，南低北高，沿岷江河谷，呈东北西南向，河谷山岭相对高差悬殊，地形复杂，属于典型的高山峡谷地形。
2008年5月12日14时28分，汶川地区发生里氏8.0级大地震，造成了大量的经济损失和人员的伤亡，同时由于地震的发生，对汶川地区的生态系统造成了巨大的破坏[6]，其中地面覆盖物的破坏尤其严重，地面覆盖物又以植被为主，因此汶川地区生态系统的破坏很大程度上表现为地面植被的破坏为了更好的实施灾后重建与恢复工作，需要对汶川地震灾区的植被变化进展实时监测，因此对汶川地震灾区植被覆盖度的监测估算与动态变化分析具有重要的意义图1 研究区地理位置2、研究方法本文中主要采用的数据是2008年-2013年共6年的汶川地震灾区空间分辨率为500m的MODIS遥感影像图因为需要研究地震对植被的破坏情况与地震以后的恢复变化情况，获取了研究区2008年4月份的成像影像与2009年-2013年5月-6月期间成像影像利用NEVI软件对获取得到的MODIS数据图像进展预处理，需要进展两步，首先是进展几何校正，并将投影方式转换为Lambert Conformal Conic兰伯特正形圆锥投影；其次几何校正后的图像进展图像裁剪，得到研究区的图像2.2 提取归一化植被指数NDVINDVI值是目前使用较为广泛的植被指数，其变化在一定程度上可以反映植被覆盖的变化，进而通过NDVI值得到植被覆盖度的大小。
NDVI[-1,+1],一般情况下，负值表示地表覆盖为云、水、雪等，0值表示地表为裸土，正值表示地表有植被覆盖，并且随着值越大覆盖度越大因此，通过NEVI软件对研究区MODIS影像处理得到2008年-2013年NDVI灰度分布图〔如图2所示〕本文采用像元二分模型法估算研究区植被覆盖度[7-9]像元二分模型是一种简单的遥感估算模型方法，它是假设一个像元里面的信息就只由土壤和植被两局部组成，由传感器传回的信息〔S〕可以分解为土壤局部所贡献的信息Ss和植被局部所贡献的信息Sv，即〔1〕在由植被和土壤构成的混合像元中，植被覆盖所占的比例就是该像元的植被覆盖度〔用Fc表示〕，那么与之相对应的就是土壤所占的比例对于一个由土壤和植被两局部组成的混合像元，像元中的植被覆盖面积比例即为该像元的植被覆盖度Fc，而土壤覆盖的面积比例为1-Fc那么就有〔2〕〔3〕其中，Sveg为纯植被覆盖贡献的信息，Ssoil为纯土壤覆盖的信息；Sv为混合像元中植被的贡献值，Ss为混合像元中土壤的贡献值 将式子〔2〕和式子〔3〕代入式子〔1〕可得： 〔4〕由公式〔4〕可以推出植被覆盖度的计算公式： 〔5〕其中，Ssoil和Sveg都只是参数，也就是说可以根据公式〔5〕利用遥感信息来对植被覆盖度进展估算。
由于像元二分模型形式比拟简单并且具有一定的物理意义，所以被广泛应用如ZRibi在半干旱地区用像元二分模型对雷达信号进展分解，得到了该地区的植被覆盖度[4];Qi J利用该模型完成了对美国西南部植被覆盖的时空变化分析，结果明确对遥感图像在不做大气纠正的情况下，像元二分模型同样具有较高的精度[5]总之，像元二分模型应用于估算植被动态变化方面比拟可靠基于归一化植被指数NDVI建立像元二分模型，可以获得植被覆盖度计算模型 〔6〕其中，NDVIsoil表示全是裸土覆盖的NDVI值，NDVIveg表示全是由植被覆盖度的NDVI值根据研究需要，结合苗苗等[12]在像元二分模型的根底上的模型研究，将给定置信区间中的NDVI最大值作为NDVIveg的值，将给定置信区间的NDVI最小值作为NDVIsoil的值，这样我们的计算模型转换为 〔7〕这样就建立了植被指数与植被覆盖度之间的关系，再根据这一关系进一步估算植被覆盖度将计算所得植被覆盖度划分为5个不同程度，包括低植被覆盖度〔0%-20%〕、中低植被覆盖度〔20%-40%〕、中植被覆盖度〔40%-60%〕、中高植被覆盖度〔60%-80%〕、高植被覆盖度〔80%-100%〕，再根据不同植被覆盖度的变化来分析得到汶川地区的植被的动态变化。
3、数据处理与分析通过对2008年-2013年汶川地区MODIS数据的分析计算，得到了每一年的NDVI值得灰度分布图像与其分布直方图，如图2、图3所示NDVI值在一定程度上能反映植被的生长情况在图2过a和b可以清晰的看出之间的差异，其在一定程度上也反映了地震以后植被收到了严重的破坏比拟图2中c-f的变化过程，反映出了植被的一个恢复过程，一年比一年好图3是图2中NDVI值的直方图其中b图比拟直观的反映了地震后植被大量减少，c-f反映出植被在逐渐增多，逐渐恢复a、 2008年 b、 2009年 c、 2010年d、 2011年 e、 2012年 f、 2013年图2 2008年-2013年NDVI值灰度分布图a、 2008年 b、 2009年 c、 2010年d、 2011年 e、 2012年 f、 2013年图3 2008年-2013年NDVI值灰度值直方图 植被覆盖度的计算是根据像元二分模型，得到与NDVI之间的关系〔7〕其中，将NDVI值小于的值定义为0，大于的值定义为1，进而再计算得到植被覆盖度的大小获取2008年-2013年植被覆盖度以后，再对其动态变化进展分析。
3.2.1 2008年-2013年整体植被覆盖度变化情况对2008年-2013年的植被的覆盖度大小进展计算得到以下数据表1 2008~2013年植被覆盖度像元数与比例统计表年份低中低中中高高像元总数像元数比例〔%〕像元数比例〔%〕像元数比例〔%〕像元数比例〔%〕像元数比例〔%〕2008年6504751894716271465801533614885092009年1964521098787177751603587764884862010年2137881104646527641251577514885302011年1305431262949558461721742974884392012年10937211018466857852051169204885382013年99303567916902099829163679488622表1中统计数据结果展示的是每一年不同植被覆盖度的植被占总的植被的面积比例例如，通过表1可以得到，2008年低植被覆盖度所占比例为13.30%、中低植被覆盖度所占比例为10.60%、中植被覆盖度所占比例为14.70%、中高植被覆盖度所占比例为30.00%、高植被覆盖度所占比例为31.40%。
由此就可以得到每一年不同植被覆盖度的变化情况由表1中的数据可以统计得到如下列图植被覆盖度变化情况由于2008年汶川地区发生特大地震，植被受到了严重的破坏，通过图4的统计可以对地震前后植被的变化情况进展分析a-e图各植被覆盖度所占面积比例变化显示，2008年-2009年地震前后，低植被覆盖度区域明显增加，由13.30%增加到40.20%；中低植被覆盖度区域明显增加，由10.60%增加到22.50%；中植被费盖度区域变化不明显；中高植被覆盖度区域明显减少，由30.00%减少到10.60%；高植被覆盖度区域明显减少，由31.40%减少到12.00%2009年-2010年中，各不同程度植被区域均变化不明显，说明在这一年时间里，植被恢复的速度较慢2010年-2013年中，植被覆盖呈现往地震前的覆盖水平开展且恢复的速度逐渐增加，表现为低植被覆盖度区域由最大值的43.80%恢复到20.30%；中低植被覆盖度区域由最大值的22.60%恢复到11.60%；中植被覆盖度区域的植被变化不明显；中高植被覆盖度区域由最小值的8.40%恢复到20.50%；高植被覆盖度区域由最小值的11.80%恢复到33.50%总体来说，由于地震的原因造成的植被破坏在2008年-2013年期间呈现一个向地震前开展的趋势。
a、 低值被覆盖度所占面积 b、 中低植被覆盖度所占面积c、 中植被覆盖度所占面积 d、 中高植被覆盖度所占面积e、 高植被覆盖度所占面积图4 2008-2013年不同植被覆盖度所占面积比例3.2.2 平均植被覆盖度变化分析图5 平均植被覆盖度与其变化趋势图图5中的数据展现了2008年-2013年期间汶川地区平均植被覆盖度的动态变化情况，通过分析可以发现：2008年-2009年地震前。