兰州地区植被变化的遥感分析及研究.doc

兰州地区植被变化的遥感分析及研究 摘要: 植被具有“指示器”的作用。研究特定地区植被动态变化是植被重建和生态环境恢复的前提。文利用分析了地区植被的。结果表明1986—2006年，呈现逐渐趋势；植被生长状况逐渐植被具有明显的年际和季节变化特征，也是联结土壤、大气和水分的自然纽带，在一定程度上，植被变化能在全球变化研究中充当“指示器”的作用，对植被的动态监测和预测可以从一个侧面反映气候变化的趋势。我国西北干旱、半干旱地区处于全球气候变化一、二级敏感带，对其植被动态变化进行研究处在中国的几何中心，北纬34°，东经 103°40＇Landsat5的 TM 图像作为信息处理的基础数据源。并且为了满足研究的需要整理了1986年6月、1996年8月、2001年7月、2006 年8月年4个时期轨道号为131/35 的TM影像以及矢量数据如1：400 万从省到县级的境界区划图以及重要城市和居民地、一二三级河流、公路、铁路等其他数据。 2.2数据预处理 首先用ERDAS软件按 EOSAT Fast Format B格式读取BSQ格式的TM影像，为便于以后的分析研究，将 7个波段读到一个img文件中，将生成的img 文件中投影信息改为 KRASOVSKY 1940 椭球和北京54 基准面下的105度中央经线的高斯投影，并添加东移500KM 信息项。因为ERDAS软件是按照图像的头文件中提供的信息进行读取，但图像供应商已按照高斯投影进行了投影变换，因此需要对图像的投影信息重定义。经过处理后和上面提到的矢量数据正确的套合在一起。然后用ERDAS软件对原始的TM数据进行几何校正和辐射校正。 在ERDAS IMAGINE8.7软件下Image Interpreter/Spectral Eahancement/indices从预处理好的四期研究区遥感影像数据中提取植被指数（NDVI）。ERDAS IMAGINE8.7软件提供的提取NDVI的计算模型为(band4-band3)/(band4+band3)，提取的四期NDVI如下： 1986年NDVI 1996年NDVI 2001年NDVI 2006年NDVI 3研究方法 笔者在研究中应用均值法和差值法两种方法。其中在统计某一区域的NDVI值时，采用均值法进行计算，即统计区域内所有象元的NDVI平均值,可以根据统计值来反映植被在一段时间内的总体变化情况。而差值法用于量化两个时期NDVI值的变化，即用后一时期的所有象元减去前一时期对应象元的NDVI值，用此方法可以反映出植被代际间的变化。 对于均值法，本研究需要求得不同时期NDVI的平均值。对于差值法，本研究即要获得三幅差值NDVI图。具体流程如图3-1。 图3-1差值法 在用差值法处理各期NDVI之前，首先要区分植被区与非植被区。在ERDAS IMAGINE8. 7软件下通过空间建模模块进行植被区与非植被区的区分。在建模前首先确定区分的阈值，在此取阈值为0.2，大于0.2的象元取值为1，小于0.2的象元取值为0,模型如图3-2。 图3-2 植被区与非植被区区分模型 其中中放置栅格影像，前后两个即要输入和输出的栅格影像，并分别在其中设置好路径；中需要输入处理数据所用的函数，在该处理过程中需要以下关键语句：“EITHER 1 IF (所处理数据的名称0.2) OR 0 OTHERWISE” 4. 结果及分析 4.1 植被的年际变化分析 对各年提取的NDVI值利用差值法进行处理。得到三个阶段的NDVI比较图，分别为1996年的NDVI值减1986年的NDVI值，2001年的NDVI值减1996年的NDVI值，2006年的NDVI值减2001年的NDVI值。 对生成的三个阶段的NDVI图统计各像素值的数目，进而观测分析植被的年际变化，如统计表4-1。 表4-1 三时间段改善/恶化像元个数统计 时间段 改善 保持 恶化 96年-86年 5528641025072 01年-96年 1068896408528 06年-01年 203616829856 图4-2为研究区改善/恶化象元数目的直观图比较，可以更加形象直观的反映出植被的年际变化。 图4-2研究区改善/恶化象元数目比较图 为定量的反映出兰州地区的年际植被变化，确定植被改善/恶化比来反映植被变化，即改善的象元数目/恶化的象元数目。当改善/恶化比大于1时植被状况改善，当改善/恶化比小于1时植被状况恶化，当改善/恶化比等于1时植被状况未发生改变。 据此根据上表4-1的统计数值，可得96年相对于87年NDVI的改善恶化比为：