《基于ndvi的色季拉山区鲁朗河流域植被覆盖度遥感估算研究》由会员分享，可在线阅读，更多相关《基于ndvi的色季拉山区鲁朗河流域植被覆盖度遥感估算研究（14页珍藏版）》请在金锄头文库上搜索。
1、 西藏大学农牧学院2013届本科毕业论文（设计）题目：基于NDVI的色季拉山区鲁朗河 流域植被覆盖度遥感估算研究姓 名： 专 业： 班 级： 学 号： 所在学院：资源与环境学院指导教师： 完成日期：2013.5.8西藏大学农牧学院本科生毕业论文（设计）原创性及知识产权声明本人郑重声明：所呈交的毕业论文是本人在导师的指导下取得的成果。对本论文研究做出重要贡献的个人和集体，及本文引用他人研究成果的均已在文中以明确方式标明。因本毕业论文引起的法律结果完全由本人承担。本毕业论文成果归西藏大学农牧学院所有。特此声明毕业论文作者签名：作者专业：作者学号：基于NDVI的色季拉山区鲁朗河流域植被覆盖度遥感估算
2、研究摘要：利用色季拉山区鲁朗河流域1999年和2011年的TM遥感影像为数据源，运用ERDAS9.2和ARCGIS10.0软件，采用归一化植被指数（NDVI）结合像元线性分解模型，测算出色季拉山鲁朗河流域2个不同时期的植被覆盖度，并利用植被覆盖等级变化的转移矩阵，探讨该区域近13年来植被覆盖等级的时空变化。关键词：植被覆盖度，NDVI，像元二分法，遥感植被覆盖度指包括乔、灌、草和农作物在内的所有植被的冠层在生长区域地面的垂直投影面积占该地区区域面积的百分比1-3，是描述植被群落及生态系统变化的重要参数。在水土保持、调节大气成分以及缓减温室气体浓度升高和维持气候稳定等方面具有不可替代的作用,而植
3、被覆盖度是衡量地表植被覆盖的重要量度。传统植被覆盖度的测量方法主要有目估法、采样法、仪器法等，但由于主观性强、野外操作不方便、尺度小和成本较高而难以在大范围内快速提取植被覆盖度2,4-6。遥感技术作为本世纪60年代发展起来的一门新术在运用遥感图像进行植被覆盖度变化研究方面得到广泛的运用。目前国内外运用遥感技术进行植被覆盖度测算而被广泛应用的方法是植被指数法，该方法由于受照相技术影响小、影像几何配准误差小和对人依赖小等特点备受研究者青睐2,7-11。植被指数经多年的发展,较常用的有:比值植被指数(RDI)、归一化植被指数(NDVI)、消除土壤影像植被指数(SADI)、垂直植被指数(PDI)等6,
4、8-10。目前,归一化植被指数(NDVI)因植被空间覆盖范围广、植物检测灵敏度高和数据具有可比性等优点成为多种植被指数中应用最多最广泛的一种2,6-11。本文即应用归一化植被指数（NDVI）结合像元线性分解模型通过对色季拉山鲁朗河流域遥感图像的处理、分析和提取，求算出色季拉山区鲁朗河流域不同时期的植被覆盖度，目的是探明应用该研究方法来确定研究区域的植被覆盖度状况及其变化情况，为相关区域的研究领域提供基础数据和技术支持，进而为进一步研究整个藏东南的植被覆盖提供相关数据和技术支持。1研究区域状况色季拉山鲁朗河流域（293210295712N、943349945217E）位于藏东南林芝地区林芝县内。
5、流于南北长约为47km，东西宽约32km，海拔24806800m,318国道纵贯流域南北。气候类型属较典型的亚高山温带半湿润气候区，冬暖夏凉、干湿季节明显。由于研究区域地处青藏高原之藏东南区域，年平均气温-0.8，最高月(7月)平均气温约9，最低月(1月)平均气温约-14。在中低海拔段为针阔混交林，林线之上为高山草甸区，海拔约4850m以上是稀疏植被群落、冰川、终年积雪区域3,12。2研究材料 研究色季拉山鲁朗河流域的遥感数据为TM影像（7波段），分别为1999-11-26的影像和2011-11-27的影像。该数据来源于中国科学院计算机网络信息中心国际科学服务平台，网址为：（http:/）。3
6、技术路线图图像降噪处理几何校正、拼接、裁剪等NDVI参数确定计算NDVI累计频率计算NDVImax和NDVIoNDVImax=NDVIgNDVImin=NDVIof=(NDVI-NDVIo)/(NDVIg-NDVIo)置信度植被覆盖度时空变化分析TM影像预处理图1植被覆盖度及NDVI的求算和分析流程4.植被覆盖度求算4.1 NDVI求算首先对图像进行预处理，也就是运用ERDAS9.2软件对TM图像进行几何校正、图层组合与裁剪等工作。本文研究选取的Landsat/TM遥感影像,共有7个波段,其中TM3(波长0.630.69gm)为红外波谱段,为叶绿素主要吸收波段；TM4(波长0.76O.90gm
7、)为近红外波谱段，对绿色植被的差异敏感，为植被通用波段。归一化植被指数NDVI的定义是:NDVI=(NIR- R)/(NIR+R)，(其中NIR 代表近红外波段, R代表红波段)要计算NDVI,就是在遥感处理软件中,计算近红外波段与红波段之差,再除以两个波段之和。利用遥感影像处理软件ERDAS9.2提取植被指数流程一般为13:(1)打开ERDAS9.2软件，打开Interpreter窗口点击Spectral Enhancement选项找到Indices并点击打开；(2)输入TM图像，并在下面Select Function中选中NDVI，进行计算；(3)生成植被指数影像文件。4.2植被覆盖度求算
8、 4.2.1像元线性分解模型像元分解模型4，6，14原理是:图像中的一个像元实际上可能由多个部分构成,每个部分对遥感传感所观测到的信息都有贡献,因此可以将遥感信息(波段或植被指数)分解,建立像元分解模型,并利用此模型估算植被覆盖度。像元二分法模型假设一个像元是由土壤与植被两部分信息组成。即通过遥感传感器所观测到的信息S，可表达为由植被所提供的信息Sz与由土壤所提供的信息St这两部分组成。即:S=Sz+St（1）；对于一个由土壤与植被两部分组成的混合像元,像元中有植被覆盖的面积比例即为该像元的植被覆盖度fg,而土壤覆盖的面积比例为1-fg。设全由植被所覆盖的纯像元所得的遥感信息为Sg。混合像元的
9、植被成分所贡献信息Sz可以表示为Sg与fg的乘积:Sz=Sgfg（2）；同理,设全由土壤所覆盖的纯像元所得的遥感信息为So，混合像元的土壤成分所提供的信息Ss可以表示为So与1-fg的乘积:Ss=(1-fg)So（3）；将式（9）与式（3）代入式（1）可得：S=fgSg+(1-fg)So（4）；对公式（4）进行变换,可得以下公式:fg=(S-So)/(Sg-So)（5）；其中SO与Sg都是参数,因而可以根据公式（5）来利用遥感信息估算植被覆盖度。此模型的参数SO与Sg则具有实际含义,即土壤与植被的纯像元所反映的遥感信息。4.2.2参数确定NDVIo值按理说是相对固定不变的,对于大多数类型的裸地
10、表面,理论上NDVIo的值应该接近零。但实际环境中由于受大气、地表湿度等条件的影响,NDVIo值也会随着变化,因此,采用一个确定的NDVIO值是不可取的14-15；同理，NDVIg代表着全植被覆盖像元的最大值，由于植被的季节变化、叶冠背景的污染（包括潮湿地面、雪、枯叶等因素）和植被类型的不同,NDVIg值的确定也存在着与NDVIo值类似的情况。考虑到遥感影像中不可避免的存在噪声、地表覆盖物及林冠背景的污染等因素，NDVIo与NDVIg并不能直接取实际计算出来的NDVI的最小值与最大值4,-7，15-16，本文对其取值时，根据图像大小和清晰程度结合NDVI值的理论变化范围，决定在NDVI频率累积
11、表上取频率为1%,99%的区间作为置信区间。4.2.3估算植被覆盖度根据归一化植被指数(NDVI)的定义：NDVI=(NIR-RED)/(NIR+RED)(6)；根据像元二分法模型,图像中每个像元的NDVI值可以看成是有植被覆盖部分的NDVIg与无植被覆盖部分的NDVIo的加权平均,由公式（4）变换可得下面的计算植被覆盖度的公式:fg=(NDVI-NDVIo)/(NDVIg-NDVIo)（7） 其中,NDVIo为裸土或无植被覆盖区域的NDVI值,即无植被像元的NDVI值；而NDVIg则代表完全被植被所覆盖的像元的NDVI值,即纯植被像元的NDVI值。5结果与分析5.1色季拉山区鲁朗河流域植被覆
12、盖度变化分析植被覆盖度的值越大表示植被覆盖度越高，本文取定植被覆盖度的取值区间为1，100%，研究参照任徳智3、陈云浩17等人的植被覆盖度等级划分标准，将A0定义为无植被覆盖；0A0.25定义为低植被覆盖；0.25A0.5定义为中植被覆盖；0.50.95定义为全植被覆盖。根据上述定义可将植被覆盖度计算结果转化为植被覆盖等级图（图2）。图2 1999年和2011年植被覆盖等级示意图由于两期TM影像均为冬季拍摄，高海拔区受积雪影响较大（如图2所示），而该区域主要为无植被覆盖区和低植被覆盖区，因此，对无植被覆盖区和低植被覆盖区不进行进行定量分析；而中、高和全植被覆盖区主要集中在低海拔区，该区受积雪影
13、响较小，因此能较为准确的反映出研究区域内中、高、全植被覆盖的变化情况。为定量分析研究区域植被覆盖度的空间格局变化，本文对1999年和2011年植被覆盖度等级的斑块数、总面积、平均斑块面积和标准差、Shannon多样性指数等景观特征参数进行了计算3,17（表1）。表1 色季拉山区鲁朗河流域各时期植被覆盖度等级的空间景观特征时 间年份植被覆盖等级面 积（Km2）斑块数（个）斑块平均面 积(Km2)斑块标准偏差多样性指数1999年No293.2111410.265.071.38L209.9427800.081.11M182.4219710.091.30H136.109490.141.32F5.001
14、440.030.10S826.6769850.122.332011年No280.6739050.072.321.40L132.87154350.010.06M169.29157500.010.07H234.0459270.041.07F9.8221570.010.02S826.69431740.020.81注：No：无植被覆盖度；L：地植被覆盖度；M：中植被覆盖度；H：高植被覆盖度；F：全植被覆盖度；S：合计；P：占有率/%；V：年变化率/%；表2同。由表1可以看出，从1999年到2011年，色季拉山区鲁朗河流域植被覆盖等级的空间景观格局有以下特点：（1）整体上景观格局表现出破碎化增大发展趋势。因考虑积雪的影响，只对研究区域中、高、全植被覆盖度斑块数进行定量比较，从1999年到2011年中、高、全植被覆盖区域呈快速破碎化的趋势。在1999年中、高、全植被覆盖区域斑块总数仅为3064块，而到了2011年中、高、全植被覆盖区
