1、1.绪论1.1 课题研究的目的与意义植被，包括森林、灌丛、草地和农作物，既是生态系统的主要组成部分，也是生态系统存在的基础，具有截流降雨、减缓径流、防沙治沙、保持水土等功能，联结着土壤、大气和水分等自然过程，在陆地表面的能量交换、生物地球化学循环和水文循环等过程中扮演着重要角色，是全球变化研究中的“指示器”1。植被根据生态系统中水、气等的状况，调控其内部与外部的物质、能量交换。植被覆盖与气候因子关系极为密切，研究植被覆盖变化对气候的影响是气候变化研究的主要内容之一，它影响着土壤湿度、地表温度和地表能量与水的循环(李苗苗，2004)。植被覆盖度(vegetation fractional cov
2、er，简称 FC)是指植被(包括叶、茎、枝)在地面的垂直投影面积占统计区总面积的百分比2。它是植被对地面的垂直投影比例，对于山坡进行植被覆盖度测量时，应该采用垂直于坡面的角度。植被覆盖度具有强烈的尺度效应，同一片植被，因被纳入统计的范围不同而表现为不同的植被覆盖度。如一个地区的植被覆盖度很高，但平均到全国水平就大大降低了3。植被覆盖度在提示地表植被分布规律, 探讨植被分布影响因子 , 分析评价区域生态环境 , 及时准确地掌握其动态变化, 分析其发展趋势对维护区域生态平衡等方面都具有重要意义。4而城市植被则是城市生态系统重要的还原组织和最重要的元素,对于保护城市生态环境具有不可忽视的作用5 ,如
3、有效缓解城市“热岛效应”,改善城市区域小气候57 等。城市化的迅速推进，带来了多样化的生态足迹，植被覆盖度，土壤污染率，地表侵蚀率，逐渐成为生态研究的热点，也成为环境保护的重点。借助于高速发展的RS与GIS技术来进行植被覆盖度的估算，将是当前环境监测的必要步骤。徐州是由矿区发展起来的城市，由于长期开采矿产，导致了一系列严重的生态问题，如塌陷地广布，植被破坏率严重，土地侵蚀率增大，等。在此背景下，研究徐州市整体的土地覆盖情况，即是现实需要，也是未来生态城市规划的重要步骤。研究的最终结果也会给徐州市的城市规划提供信息支持与技术保障。1.2 国内外植被覆盖度研究现状由于植被覆盖度是许多学科的重要参数
4、，为了得到准确的植被覆盖度信息，植被覆盖度监测技术的提高，就成了多个领域发展的需要。根据检测手段，测量植被覆盖度的方法可分为传统的地面测量和新兴的遥感测量两大类。其中，地面测量又可以根据测量原理，分为目估法、采样法、仪器法和模型法；遥感测量依据对植被光谱信息与植被覆盖度所建立的关系不同，可分为物理模型法和统计模型法。统计模型法中使用较多的有植被指数法、回归模型法、像元分解法、分类决策树和人工神经网络法；物理模型法中模型反演法使用最多。地面测量曾经一度是植被覆盖度监测的最主要方法。主要包括目估法、采样法、仪器法和模型法。虽然遥感技术的发展使地面测量的主导性地位有所降低，但地面测量依然具有其重要性
5、，它不仅是最精确的测量方法，也为遥感测量提供了基础标定数据，是无可替代的。遥感技术的发展，为大范围植被覆盖信息的获取提供了一个新的发展方向。常用于植被覆盖信息提取的遥感数据有 NOA 内叭 vHRR 数据、MODIs 数据、 LandsatTM 与MSS 数据、SPOT 数据、ATSER 数据、航片、IEOS 一 SAR 雷达数据以及 AVIRIS 高光谱数据等。近年来，国内外在植被遥感监测方面开展大量的研究，发展许多植被覆盖度监测方法8 13 。遥感方法相较于传统地面测算使测算的外业工作极大减少，在时效性、测算范围等方面都存在较明显优势。孙睿等14利用 NOAA 数据，通过研究不同时段降水量
6、与年最大 NDVI 之间的相关关系分析降水对黄河流域植被覆盖度的影响; 赵彩霞等15通过定量分析植被覆盖度与土壤风蚀量及风蚀气候侵蚀因子 3 者随时间变化的相关关系，计算和比较不同类型植物防风治沙性能的动态差异、总植被覆盖度及相应的总土壤风蚀量动态变化; RezaAmiri 等16利用基于 NDVI 值，分析植被覆盖度和土地利用/土地覆被与土地表面温度的时空变化 3 者相关性。将气候、土壤因素与植被覆盖度进行相关研究较多，而将地质条件与植被覆盖度相结合的研究尚未有人涉及。1.3 主要研究内容与方法1.3.1 研究区概况徐州市位于江苏省的西北部，东 经 11622 11840、 北 纬 3343
7、 3458之间 。 东 西 长 约 210 公 里 ， 南 北 宽 约 140 公 里 ， 总 面 积 11258 平 方 公 里 ， 占 江 苏 省 总面 积 的 11 。 地 处 苏 、 鲁 、 豫 、 皖 四 省 交 界 ， 为 东 部 沿 海 与 中 部 地 带 上 海 经 济 区 与环 渤 海 经 济 圈 的 结 合 部 。 素 有 “五 省 通 衢 ”之 称 。 京 沪 、 陇 海 两 大 铁 路 在 此 交 汇 ，京 杭 大 运 河 傍 城 而 过 贯 穿 徐 州 南 北 ， 公 路 四 通 八 达 ， 北 通 京 津 ， 南 达 沪 宁 ， 西 接 兰新 ， 东 抵 海 滨 ，
8、 为 全 国 重 要 水 陆 交 通 枢 纽 和 东 西 、 南 北 经 济 联 系 的 重 要 十 字 路 口 。图 1 徐州市区位图1.3.2 研究内容与技术路线1.3.2.1 研究内容 在遥感应用领域,植被指数是一种反映地表植被信息的重要信息源之一,已广泛用于定性和定量评价植被覆盖及其生长状况。植被指数经过20 多年的发展,目前有几十种,但常用的植被指数有:归一化植被指数(NDVI) 、比值植被指数(RVI) 、差值植被指数(NVI) 、土壤调节植被指数(SAVI) 和修正土壤植被指数(MSAVI) 等17 - 18 。其中归一化植被指数是植被生长状态及植被覆盖度的最佳指示因子。本文主要
9、利用 2010 年徐州市 TM 图像进行植被覆盖情况的研究，利用归一化植被指数法与象元二分模型进行植被覆盖度的反演，最终得到基本植被覆盖情况与徐州市六县的分区覆盖情况。1) 归一化植被指数 NDVI。植被指数是利用绿色植物强吸收可见光红波段( 0.6 0.7 m) 和高反射近红外波段 ( 0.7 1.1 m) 特点，经过变换，增强植被信号，削弱噪音组合而成。归一化植被指数( NDVI，下式为 NNDVI) 是植被生长状态及植被覆盖度最佳指示因子19,20 ，被定义为近红外。波段与可见光红波段数值之差和这两个波段数值之和比值，计算公式:NDVI = ( NIR R) /( NIR + R) (
10、1)式中，NIR 为近红外波段的反射率，R 为红光波段的反射率。2）像元二分模型。像元二分模型原理是假设影像上一个像元的反射率R 可分为纯植被部分反射率Rv 和非植被部分Rs 两部分,那么, 任一像元的反射率值可以表示为由植被覆盖部分与非植被覆盖部分的线性加权的和:(2) SV假设影像上一个像元中有植被覆盖的面积比例为f c ,即该像元的植被覆盖度, 那么非植被覆盖的面积比例为1 - f c 。如果该像元全由植被所覆盖,则所得的反射率为R veg,如果该像元无植被覆盖,则反射率为R soil ,因此,混合像元的植被部分所贡献的信息Rv 可以表示为纯植被反射率R veg与像元中植被覆盖面积f c
11、 的乘积(见公式( 2) ) , 而非植被成分所贡献的信息Rs 可以表示为R soil 与1 - f c 的乘积(见公式(3) ) :(3) Vegv*(4) soilSRf)-(通过解算式(2) 、式(3) 和式(4) , 可得到计算植被覆盖度的公式 ,如下:(5)/()(soilvegsoilfc其中:R soil与R veg是像元二分模型的两个参数。很明显 ,只要求得这两个参数根据公式 (5), 就可以利用遥感信息来估算植被覆盖度。根据像元二分模型原理,我们可以将1个像元的NDVI 值表示为由有植被覆盖部分地表与无植被覆盖部分地表组成的形式。因此,计算植被覆盖度的公式可表示为:(6)so
12、ilvegNDVIIfc其中:NDVI soil为完全是裸土或无植被覆盖区域的NDVI 值,NDVI veg则代表完全由植被所覆盖的像元的NDVI值,即纯植被像元的NDVI 值。参考李苗苗20 、Gutman21 等提出的估算NDVI veg和NDVI soil的方法，根据整幅影像上NDVI 的灰度分布，以0.5% 置信度截取NDVI 的上下限阈值分别近似代表NDVI veg和NDVI soil。采用ARCGIS软件的栅格计算器，应用式( 6) 计算所有影像的植被覆盖度分布情况。将计算得到的植被覆盖度( FC) 分5级: 低植被覆盖度( FC 10%) 、较低植被覆盖度( 10%FC 30%)
13、 、中度植被覆盖度( 30% FC 50%) 、较高植被覆盖度( 50%FC 70%) 和高植被覆盖度( FC70%) .1.3.2.2 技术路线图 2 技术路线图2 课题研究主要步骤2.1 数据准备采用徐州市2010年的Landsat TM 遥感数据,共包括了三景图像，轨道号分别为TM 影像植被覆盖信息提取前期准备数据处理土地类型分类水体剔除植被指数计算植被覆盖度估算统计分析.专题图生成结果讨论与分析p121r036、p121r037、p122r036,Landsat 7 ，TM有7个波段 ,其中1-5, 7波段的分辨率均为30 m; 6波段分辨率为120 m,为热红外波段,不参与此试验研究
14、; 首先利用ERDAS8.6的工具，依次选择Interpreter-Utilities-Layer Stack，将原始tiff数据合成多光谱图像，因为第六波段为热红外波段，所以予以舍弃。2.2 植被覆盖度遥感信息提取2.2.1 土地利用类型分类1) 基本类型划分 利用最大似然法进行监督分类，主要流程为：建立分类模板、评价分类模板、执行监督分类、精度评价。监督分类又叫训练区分类，它的前提是已知遥感图像上样本区内地物的类别，所以选择足够数量和具有代表性的训练区，决定了监督分类的精度。由于选取的研究区面积较大，所以在训练区数量上均在 3000 以上。根据图像的分辨率，将徐州市分为六类土地类型：水体，
15、建设用地，农田，林地，山地，未利用地。分类精度满足要求。2)分类属性重编码 由于水体在植被指数图像上多显示为负值，影响植被覆盖度的估算，所以先进行水体的剔除，这里采用的方法是掩膜法。在分类图像上进行属性的重编码，将水体的属性设置为 0，其余地物的属性均设置为 1，即是对图像进行二值化处理。结果使水体与其他地物区分开。所得图像为：图 3 二值化处理图像3）水体剔除 将二值化处理的图像与原始图像相乘，即可以将水体剔除掉，这里利用的是 ERDAS 的建模功能，所建立的模型如下所示:图 4 剔除水体模型图 5 剔除水体的图像2.2 2 植被指数计算利用 ERDAS 进行植被指数的计算主要有两种方法：一
16、种是直接利用菜单工具：Interpreter-spectral Enhancement-Indices。植被指数选择 NDVI，设置相应的参数即可以自动进行计算。第二种方法是利用 ERDAS 的空间建模工具进行计算。模型如下所示：图 6 植被指数计算模型 本文利用第一种方法进行计算，因为第二种方法需要已知裸地的植被覆盖情况与最高植被覆盖情况，因为缺乏相关经验与实测数据，所以先计算出植被指数，最后在ARCgis 中进行植被植被覆盖度的估算。得到植被指数分布图如下：图 7 徐州市植被指数图2.3 基于 GIS 的植被覆盖度估算2.3.1 植被覆盖度估算根据植被指数整幅图像的灰度分布情况，截取上下限
17、阈值，近似的估计NDVIveg =0.563 NDVIsoil=0.0078，采用ARCGIS软件的栅格计算器（Spatial Analyst-Raster Calculator)，应用植被覆盖度的估算公式计算所有影像的植被覆盖度分布情况。将计算得到的植被覆盖度( FC) 分5级: 低植被覆盖度( FC 10%) 、较低植被覆盖度( 10%FC 30%) 、中度植被覆盖度( 30% FC 50%) 、较高植被覆盖度( 50%FC 70%) 和高植被覆盖度( FC70%） ，并将分级图像彩色显示：图8 徐州市植被覆盖度图栅格计算得到的图像格式 为 GRID 格式，不具备栅格图像的灰度信息，为了得
18、到不同植被覆盖度的面积统计值，将该图像进行属性值的重新设置低植被覆盖度的属性值为 1 ，较低植被覆盖度的属性值为 2，中度植被覆盖度的属性值为 3，较高植被覆盖度的属性值为 4，高植被覆盖度的属性值为 5，水体与误差错分区的属性值设置为 0，这样得到按照新的属性值统计的栅格总数，根据栅格总数可以计算出每种植被覆盖情况的面积。重分类的设置如下：图 9 植被覆盖度属性重分类2.3. 2 分区植被覆盖统计利用徐州市的行政区划图与植被覆盖图进行叠加，计算出徐州市六县的植被覆盖情况，这里的六县植被覆盖情况采用的平均值估计。图10 植被覆盖度分区图图11 植被覆盖度分区统计设置图3 结果分析1）经过分区统
19、计与面积统计，在所得的徐州市六区植被覆盖情况统计表（表一）中显示出：徐州六县及徐州市区的植被覆盖度总体较低，没有较高及高度植被覆盖区域。这说明徐州市的生态环境相对较差，对植被的覆盖情况不重视。新沂与丰县的植被覆盖度是所有区域中最高的，睢宁次之，以下依次为沛县，邳州，铜山，徐州市区，这与各个区域的经济发展水平有一定的相关性。经济发展程度高的，对城区的开发利用程度较高，植被覆盖率就相应降低。反之，植被覆盖率较高。同时这也反映目前徐州市的经济发展在一定程度上还是以牺牲自然资源与环境为代价的。表一 徐州六县分区植被覆盖度统计地区 邳州 沛县 丰县 新沂 铜山 徐州睢宁植被覆盖度0.2 0.22 0.2
20、8 0.28 0.17 0.070.272)从植被覆盖度面积统计表格（表二）可以看出：较低，中度植被覆盖度的面积所占份额较大，较低植被覆盖度面积占到了徐州市植被覆盖总面积的43.3%，中度植被覆盖度面积占到了徐州市植被覆盖总面积的28.3%，低度植被覆盖度面积占到了徐州市植被覆盖总面积的16.9%。高度的植被覆盖度所占比例最少，仅有0.9%，这与分区统计的结果有一定的相似。也进一步证明徐州地区的的植被覆盖度较低。也反映了徐州市在整顿矿区发展，发展地方经济的同时，忽略了生态环境建设。表二 各植被覆盖度面积统计 （单位：hm 2)植被覆盖状况低度 较低 中度 较高 高度面积160787.52413
21、364.42270254.7100843.118756.734 结论经过徐州市的植被覆盖情况分析可以得出结论：徐州市的植被覆盖度相对较低，没有较高及高度植被覆盖区域。这说明徐州市的生态环境相对较差，对植被的覆盖情况不重视。新沂与丰县的植被覆盖度是所有区域中最高的，睢宁次之，徐州市区的植被覆盖度最低，只有0.07，属于低度植被覆盖区域。这与各个区域的经济发展水平有一定的相关性。在面积统计中，较低，中度植被覆盖度的面积所占份额较大。高度与低度的植被覆盖度所占比例最少。如果这种情况再继续下去，将会使徐州市的生态环境更加恶劣。所以，在城市扩展，城市开发的过程中，适度的重视环境规划，提高植被覆盖度，也有
22、利于徐州市的长远发展。5 参考文献1 赵英时.遥感应用分析原理与方法M.北京:科学出版社,2003.2 Gitelson A. A., Kaufman Y. J., Stark R ., et al. Novel algorithms for remote estimationof vegetation fraction J .Remote Sensing of Environment,2002(80):76-87.3 胡良军,邵明安.论水土流失研究中的植被覆盖度量指标J.西北林学院学报,2001,16(1):40-43.4祁燕王 秀兰 冯仲科 郭祥 2 基于 RS 与 GIS 的北京市植被覆
23、盖度变化研究J林业调查规划 2009 34（2）1-4 5 何云玲,张一平 . 城市生态环境与绿化植被相互作用研究J . 高原气象,2004 ,23 (3) :2972304. 6 沈涛,袁春琼,刘玉安. 乌鲁木齐市热岛强度分布与植被覆盖相互关系的遥感研究J . 新疆气象,2004 ,27 (1) :28230. 7 刘树华,黄子琛 ,刘立超 . 半干旱区植被覆盖度对边界层气候热力影响的数值模拟 J . 气象学报,1996 ,54 (3) :3032311. 8 张春玲，余华，宫鹏，等武汉市地表亮温与植被覆盖关系定量分析J 地理科学，2009，29( 5) : 740 7449 颜长珍，吴炳方
24、晋陕蒙接壤区林草覆盖变化的遥感分析J 地理科学，2004， 24( 4) : 465 47110 张云霞，李晓斌，陈云浩 草地植被覆盖度的多尺度遥感与实地测量方法综述J 地球科学进展，2003，18( 1) : 85 9311 陈晋，陈云浩，何春阳，等 基于土地覆盖分类的植被覆盖率估算亚像元模型与应用J 遥感学报， 2001，5( 6) : 416 42212 Purevdo r J S，Tateishi R， Ishiyama T，et al Relationship s betweenpercent vegetation cover and vegetation indicesJ Inte
25、rnationalJournal of Remote Sensing，1998，19 ( 18 ) : 3519 353513 万恩璞，吕宪国，王野乔，等 应用遥感信息对长白山植被覆盖的研究J 地理科学，1989，9( 4) : 354 38714 孙睿，刘昌明，朱启疆 黄河流域植被覆盖度动态变化与降水的关系J 地理学报，2001， 56( 6) : 667 67215 赵彩霞，郑大伟，何文清 植被覆盖度的时间变化及其防风蚀效应J 植物生态学报，2005， 29( 1) : 68 7316 Amiri R，Weng Q，Alimohammadi A Spetial temporal dynam
26、icsof land surface temperature in relation to fractional vegetationcover and land use /cover in the Tabriz urban area， IranJ RemoteSensing of Environment，2009， 113( 12) : 2606 261717 LYONJ G,YUAN D,LUNETTA R S ,et al.A change detection experiment usingvegetation indicesJ . Photogrammetric Engineerin
27、g and Remote Sensing ,1998 ,64 :143 - 150.18 田庆久,闵祥军. 植被指数研究进展J . 地球科学进展,1998 ,13(4) :327- 333.19 赵英时 遥感应用分析原理与方法M 北京: 科学出版社，2003: 387 398.20 李惠敏，刘洪斌，武伟 近10 年重庆市归一化植被指数变化分析J 地理科学， 2010，30( 1) : 119 12321 李苗苗 植被覆盖度的遥感估算方法研究D 中国科学院遥感应用研究所，2003: 1 11122 Gutman G G Vegetation indices from AVHRR: An update and futureprospectsJ Remote Sensing of Environment， 1991， 35( 2 3) : 121 136
