基于GIS的植被覆盖度估算

1、基于GIS的植被覆盖度估算.绪论课题研究的目的与意义植被，包括森林、灌丛、草地和农作物，既是生态系统的主要组成部分，也是 生态系统存在的基础，具有截流降雨、减缓径流、防沙治沙、保持水土等功能，联 结着土壤、大气和水分等自然过程，在陆地表面的能量交换、生物地球化学循环和 水文循环等过程中扮演着重要角色，是全球变化研究中的“指示器”1 0植被根据生态系统中水、气等的状况，调控其内部与外部的物质、能量交换。植被覆盖与 气候因子关系极为密切，研究植被覆盖变化对气候的影响是气候变化研究的主要内容之一，它影响着土壤湿度、地表温度 和地表能量与水的循环（李苗苗，2004）。植被覆盖度（vegetation

2、fractional cover,简称FC）是指植被（包括叶、茎、枝）在地面的垂直投影面积占统计区总面积的百分比 2。它是植被对地面的垂直投 影比例，对于山坡进行植被覆盖度测量时，应该采用垂直于坡面的角度。植被覆 盖度具有强烈的尺度效应，同一片植被，因被纳入统计的范围不同而表现为不同的 植被覆盖度。如一个地区的植被覆盖度很高，但平均到全国水平就大大降低了0植被覆盖度在提示地表植被分布规律，探讨植被分布影响因子，分析评价区 域生态环境，及时准确地掌握其动态变化，分析其发展趋势对维护区域生态平衡等 方面都具有重要意义。4而城市植被则是城市生态系统重要的还原组织和最重要 的元素，对于保护城市生态环境

3、具有不可忽视的作用5,如有效缓解城市“热岛效 应”，改善城市区域小气候5,7等。城市化的迅速推进，带来了多样化的生态足迹，植被覆盖度，土壤污染率，地 表侵蚀率，逐渐成为生态研究的热点，也成为环境保护的重点。借助于高速发展的 RS与GIS技术来进行植被覆盖度的估算，将是当前环境监测的必要步骤。徐州是由矿区发展起来的城市，由于长期开采矿产，导致了一系列严重的生态 问题，如塌陷地广布，植被破坏率严重，土地侵蚀率增大，等。在此背景下，研究 徐州市整体的土地覆盖情况，即是现实需要，也是未来生态城市规划的重要步骤。 研究的最终结果也会给徐州市的城市规划提供信息支持与技术保障。国内外植被覆盖度研究现状由于植

4、被覆盖度是许多学科的重要参数，为了得到准确的植被覆盖度信息，植 被覆盖度监测技术的提高，就成了多个领域发展的需要。根据检测手段，测量植被 覆盖度的方法可分为传统的地面测量和新兴的遥感测量两大类。其中，地面测量又 可以根据测量原理，分为目估法、采样法、仪器法和模型法 ；遥感测量依据对植被 光谱信息与植被覆盖度所建立的关系不同，可分为物理模型法和统计模型法。统计 模型法中使用较多的有植被指数法、回归模型法、像元分解法、分类决策树和人工 神经网络法；物理模型法中模型反演法使用最多。地面测量曾经一度是植被覆盖度监测的最主要方法。主要包括目估法、采样 法、仪器法和模型法。虽然遥感技术的发展使地面测量的主

5、导性地位有所降低， 但地面 测量依然具有其重要性，它不仅是最精确的测量方法，也为遥感测量提供 了基础标 定数据，是无可替代的。遥感技术的发展，为大范围植被覆盖信息的获取提供了一个新的发展方向。常 用 于植被覆盖信息提取的遥感数据有 NO秋叭vHRRt据、MODI嗷据、 LandsatTM与MSS数据、SPO蹴据、ATSE喷据、力行t、IEOS SARt达数据以 及AVIRIS高光谱数据等。近年来，国内外在植被遥感监测方面开展大量的研究，发展许多植被覆盖度监测方法,8 ,13,。遥感方法相较于传统地面测算使测算的外业工作极大减少，在时 效性、测算范围等方面都存在较明显优势。孙睿等 ,14,利用N

6、OA峨据，通过研究 不同时段降水量与年最大NDVI之间的相关关系分析降水对黄河流域植被覆盖度的 影响；赵彩霞等,15,通过定量分析植被覆盖度与土壤风蚀量及风蚀气候侵蚀因子3者随时间变化的相关关系，计算和比较不同类型植物防风治沙性能的动态差异、总 植被覆盖度及相应的总土壤风蚀量动态变化；RezaAmiri等,16,利用基于NDVI 值，分析植被覆盖度和土地利用/土地覆被与土地表面温度的时空变化 3者相关 性。将气候、土壤因素与植被覆盖度进行相关研究较多，而将地质条件与植被覆盖 度相结合的研究尚未有人涉及。1.3主要研究内容与方法研究区概况徐州市位于江苏省的西北部，东经 116?22 ,118?4

7、0、北纬33?43, 34?58之间。东西长约210公里，南北宽约140公里，总面积11258平方公里， 占江苏省总面积0地处苏、鲁、豫、皖四省交界，为东部沿海与中部地带上海经济区与环渤的11,海经济圈的结合部。素有“五省通衢”之称。京沪、陇海两大铁路在此交汇， 京杭大运河傍城而过贯穿徐州南北，公路四通八达，北通京津，南达沪宁，西接兰 新，东抵海滨，为全国重要水陆交通枢纽和东西、南北经济联系的重要十字路口。图1徐州市区位图1.3.2研究内容与技术路线研究内容在遥感应用领域，植被指数是一种反映地表植被信息的重要信息源之一，已广泛 用于定性和定量评价植被覆盖及其生长状况。植被指数经过 20多年的发

8、展，目前有 几十种，但常用的植被指数有：归一化植被指数(NDVI)、比值植被指数(RVI)、差 值植被指数(NVI)、土壤调节植被指数(SAVI)和修正土壤植被指数(MSAVI)等17 -18 o其中归一化植被指数是植被生长状态及植被覆盖度的最佳指示因子。本文主要利用2010年徐州市TM图像进行植被覆盖情况的研究，利用归一化植 被指数法与象元二分模型进行植被覆盖度的反演，最终得到基本植被覆盖情况与徐 州市六县的分区覆盖情况。1)归一化植被指数NDVI植被指数是利用绿色植物强吸收可见光红波段(0.6,0.7 m m)和高反射近红外波段(0.7 , 1.1 m)特点，经过变换，增强 植被信号，削弱

9、噪音组合而成。归一化植被指数(NDVI ,下式为NNDVI)是植被生 TOC o 1-5 h z 长状态及植被覆盖度最佳指示因子，19,20,被定义为近红外。波段与可见光红波 段数值之差和这两个波段数值之和比值，计算公式：NDVI = ( NIR , R)/( NIR + R) ( 1)式中，NIR为近红外波段的反射率，R为红光波段的反射率。2)像元二分模型。像元二分模型原理是假设影像上一个像元的反射率R可分为纯植被部分反射率Rv和非植被部分Rs两部分，那么，任一像元的反射率值可以表 示为由植被覆盖部分与非植被覆盖部分的线性加权的和：(2) R,R , RVS假设影像上一个像元中有植被覆盖的面

10、积比例为 fc ,即该像元的植被覆盖度, 那么非植被覆盖的面积比例为1 - f c 。如果该像元全由植被所覆盖，则所得的反 射率为R如果veg,该像元无植被覆盖，则反射率为R ,因此，混合像元的植被部分所 贡献的信息Rv可以表soil示为纯植被反射率R与像元中植被覆盖面积f c的乘积(见公式(2), 而非 植被成分所贡veg献的信息Rs可以表示为R与1 - f c的乘积(见公式(3) ) : soilR,fc*R (3) vVegR,(1-fc)*R Ssoil通过解算式(2)、式(3)和式(4),可得到计算植被覆盖度的公式，如下：fc,(R,R)/(R,R) (5) soilvegsoil其

11、中:R与R是像元二分模型的两个参数。很明显，只要求得这两个参数根据公式(5),就 vegsoil可以利用遥感信息来估算植被覆盖度。根据像元二分模型原理，我们可以将1个像元的NDVI值表示为由有植被覆盖部分地表与无植被覆盖部分地表组成的形式。因此，计算植被覆盖度的公式可表示为NDVI,NDVIsoilfc, (6) NDVI,NDVIvegsoil其中:NDVI为完全是裸土或无植被覆盖区域的 NDVI值,NDVI则代表完全由植被 所覆 soilveg盖的像元的NDVI值，即纯植被像元的NDVI值。参考李苗苗，20,、Gutman,21,等提出的估算NDVI和NDVI的方法，根据整幅影 vegso

12、il像上NDVI的灰度分布，以0.5%置信度截取NDVI的上下限阈值分别近似代表NDVI和 vegNDVL采用ARCGI球件的栅格计算器，应用式（6）计算所有影像的植被覆盖度分布soil情况。将计算得到的植被覆盖度（FC）分5级：低植被覆盖度（FC, 10%）、 较低植被覆盖度（10%?FC, 30%）、中度植被覆盖度（30% ?FC , 50%）、较高植被 覆盖度（50%?FC, 70%）和高植被覆盖度（FC?70%）.技术路线TM影像植被覆盖信息提取前期准备数据处理土地类型分类水体剔除植被指数计算植被覆盖度估算统计分析.专题图生成结果讨论与分析图2技术路线图2课题研究主要步骤数据准备采用徐

13、州市2010年的Landsat TM遥感数据，共包括了三景图像，轨道号分别为 p121r036、p121r037、p122r036,Landsat 7, TM有 7 个波段，其中 1-5, 7 波段的分辨率均为30 m; 6波段分辨率为120 m,为热红外波段，不参与此试验研究；首先利用 ERDAS8.6勺工具，依次选择 Interpreter-Utilities-Layer Stack ,将原始 tiff 数据合成多光谱图像，因为第六波段为热红外波段，所以予以舍弃。植被覆盖度遥感信息提取土地利用类型分类1）基本类型划分利用最大似然法进行监督分类，主要流程为：建立分类模 板、评价分类模板、执行

14、监督分类、精度评价。监督分类又叫训练区分类，它的前 提是已知遥感图像上样本区内地物的类别，所以选择足够数量和具有代表性的训练 区，决定了监督分类的精度。由于选取的研究区面积较大，所以在训练区数量上均 在3000以上。根据图像的分辨率，将徐州市分为六类土地类型：水体，建设用地，农田，林地，山地，未利用地。分类精度满足要求。2）分类属性重编码由于水体在植被指数图像上多显示为负值，影响植被覆盖 度的估算，所以先进行水体的剔除，这里采用的方法是掩膜法。在分类图像上进行 属性的重编码，将水体的属性设置为 0,其余地物的属性均设置为1,即是对图像 进行二值化处理。结果使水体与其他地物区分开。所得图像为 ：

15、图3二值化处理图像3）水体剔除 将二值化处理的图像与原始图像相乘，即可以将水体剔除掉，这 里利用的是ERDAS勺建模功能，所建立的模型如下所示图4剔除水体模型图5剔除水体的图像2植被指数计算利用ERDA必行植被指数的计算主要有两种方法：一种是直接利用菜单工。植被指数选择NDVI,设置相应具:Interpreter-spectral Enhancement-Indices的参数即可以自动进行计算。第二种方法是利用ERDAS勺空间建模工具进行计算。IR * VjfibleIR - VisibleI44141 Avi-WtMU/ INPUT出音注意 m ef一 lf,i| I M )141lir I

16、f 4 MlifkAM-t % in *i1.IPC7 RASTER模型如下所示:Difference I egernrimt Index、if Rlj *. * lu - WIU- MSB图6植被指数计算模型本文利用第一种方法进行计算，因为第二种方法需要已知裸地的植被覆盖情况 与最高植被覆盖情况，因为缺乏相关经验与实测数据，所以先计算出植被指数，最 后在ARCgis中进行植被植被覆盖度的估算。得到植被指数分布图如下 ：图7徐州市植被指数图基于GIS的植被覆盖度估算.3.1植被覆盖度估算根据植被指数整幅图像的灰度分布情况，截取上下限阈值，近似的估计NDVIveg =0.563 NDVI=0.0

17、078 ,采用 ARCGI漱件的栅格计算器(Spatial AnalystRaster soilCalculator),应用植被覆盖度的估算公式计算所有影像的植被覆盖度分布情 况。将计算得到的植被覆盖度(FC)分5级：低植被覆盖度(FC, 10%)、较低植 被覆盖度(10%?FC, 30%)、中度植被覆盖度(30% ?FC , 50%)、较高植被覆盖度 (50%?FC, 70%)和高植被覆盖度(FC?70%),并将分级图像彩色显示：DM乐玻学理3 LJR 可代if；CiN f果满仗痛否 味晒球押子上西舅里imur图8徐州市植被覆盖度图栅格计算得到的图像格式 为GRID格式，不具备栅格图像的灰度

18、信息，为了得 到不同植被覆盖度的面积统计值，将该图像进行属性值的重新设置低植被覆盖度的 属性值为1 ,较低植被覆盖度的属性值为 2,中度植被覆盖度的属性值为3,较高植被覆盖度的属性值为4,高植被覆盖度的属性值为5,水体与误差错分区的属性 值设置为0,这样得到按照新的属性值统计的栅格总数，根据栅格总数可以计算出每种植被覆盖情况的面积。重分类的设置如下。韦。吉5 feld;It re日书s/丁 Ctunpfr mtalng vWh to MoEWEi今3|8mpwSy匚市IE图9植被覆盖度属性重分类2.3. 2 分区植被覆盖统计利用徐州市的行政区划图与植被覆盖图进行叠加，计算出徐州市六县的植被覆

19、盖情况，这里的六县植被覆盖情况采用的平均值估计。图10植被覆盖度分区图JE建23：3 凶| jrangiu.OTLjnty二| 百Zont fled：|landlise._ma勺Value raster:|Reda55rf2CiO;J 且| loncrre NoDjfta incalcuteiti.iirisJoin tab is to zcinaf? Chart|m号an二jOutput tble；|C:Doc undents and Settingsh*_司OK I CbtiCbI图11植被覆盖度分区统计设置图3结果分析1）经过分区统计与面积统计，在所得的徐州市六区植被覆盖情况统计表（表一

20、）中显示出：徐州六县及徐州市区的植被覆盖度总体较低，没有较高及高度植被覆盖 区域。这说明徐州市的生态环境相对较差，对植被的覆盖情况不重视。新沂与丰县 的植被覆盖度是所有区域中最高的，睢宁次之，以下依次为沛县，邳州，铜山，徐 州市区，这与各个区域的经济发展水平有一定的相关性。经济发展程度高的，对城 区的开发利用程度较高，植被覆盖率就相应降低。反之，植被覆盖率较高。同时这 也反映目前徐州市的经济发展在一定程度上还是以牺牲自然资源与环境为代价的。表一徐州六县分区植被覆盖度统计睢地区 邳州 沛县 丰县 新沂 铜山 徐州宁植被覆盖0.20.2 0.22 0.28 0.28 0.17 0.07 度 7分但

21、粒被鱼反度统计叫.一,*Ml -a ijj s|手的f辅县丰县新洲铜山徐州陈宁2）从植被覆盖度面积统计表格（表二）可以看出：较低，中度植被覆盖度的面积所占份额较大，较低植被覆盖度面积占到了徐州市植被覆盖总面积的43.3%,中度植被覆盖度面积占到了徐州市植被覆盖总面积的28.3%,低度植被覆盖度面积占到了徐州市植被覆盖总面积的16.9%。高度的植被覆盖度所占比例最少，仅有 0.9%, 这与分区统计的结果有一定的相似。也进一步证明徐州地区的的植被覆盖度较低。也反映了徐州市在整顿矿区发展，发展地方经济的同时，忽略了生态环境建设。2表二各植被覆盖度面积统计（单位:hm）植被覆盖状低度较低中度较高高度况

22、160787.5413364.100843.1 面积 270254.7 8756.732 42 1务植被覆盖度面积比例统计面积折占份额4结论经过徐州市的植被覆盖情况分析可以得出结论：徐州市的植被覆盖度相对较低，没有较高及高度植被覆盖区域。这说明徐州市的生态环境相对较差，对植被的 覆盖情况不重视。新沂与丰县的植被覆盖度是所有区域中最高的，睢宁次之，徐州 市区的植被覆盖度最低，只有0.07,属于低度植被覆盖区域。这与各个区域的经 济发展水平有一定的相关性。在面积统计中，较低，中度植被覆盖度的面积所占份 额较大。高度与低度的植被覆盖度所占比例最少。如果这种情况再继续下去，将会 使徐州市的生态环境更加

23、恶劣。所以，在城市扩展，城市开发的过程中，适度的重 视环境规划，提高植被覆盖度，也有利于徐州市的长远发展。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 ofEnvironment,2002(80):76-87.3胡良军，邵明安.论水土流失研究中的植被覆盖度量指标J.西北林学院学报,2001,16(1):40-4

24、3.4祁燕王秀兰冯仲科郭祥2基于RS与GIS的北京市植被覆盖度变化研究J林业调查规划2009 34(2)1-45何云玲，张一平.城市生态环境与绿化植被相互作用研究J . 高原气 象,2004 ,23 ：2972304.6沈涛，袁春琼，刘玉安.乌鲁木齐市热岛强度分布与植被覆盖相互关系的遥 感研究J .新疆气象,2004 ,27:28230.7 刘树华，黄子琛,刘立超.半干旱区植被覆盖度对边界层气候热力影响的 数值模拟J . 气象学报，1996 ,54:3032311.,8,张春玲，余华，宫鹏，等(武汉市地表亮温与植被覆盖关系定量分析,J,(地理科学，2009, 29( 5) : 740,744,

25、9,颜长珍，吴炳方(晋陕蒙接壤区林草覆盖变化的遥感分析,J,( 地理科学，2004, 24( 4) : 465,471(,10,张云霞，李晓斌，陈云浩(草地植被覆盖度的多尺度遥感与实地测量方法综述,J,(地球科学进展，2003, 18( 1) : 85,93(,11,陈晋，陈云浩，何春阳，等(基于土地覆盖分类的植被覆盖率估算亚像元模型与应用，J,( 遥感学报，2001, 5( 6) : 416 ,422(,12, Purevdo r J S , Tateishi R , Ishiyama T , et al( Relationship s between percent vegetation cover and vegetationindices,J,( InternationalJournal of Remote Sensing1998 19 ( 18 ) :3519 ,3535(,13,万恩璞，吕宪国，王野乔，等(应用遥感信息对长白山植被覆盖的研 究,J,( 地理科学，1989,9(4) : 354,387(,14,孙睿，刘昌明，朱启疆(黄河流域植被覆盖度动态变化与降水的关系,J,( 地理学报，2001,