基于遥感影像的北京市植被覆盖率反演解析.docVIP

基于ETM遥感影像的北京市植被覆盖度反演 1 绪论 1.1 研究区概况 北京中心位于北纬39度54分20秒，东经116度25分29秒。位于华北平原西北边缘。毗邻渤海湾，上靠辽东半岛，下临山东半岛。北京与天津相邻，并与天津一起被河北省环绕。西部是太行山山脉余脉的西山，北部是燕山山脉的军都山，两山在南口关沟相交，形成一个向东南展开的半圆形大山弯，人们称之为“北京弯”，它所围绕的小平原即为北京小平原。诚如古人所言：“幽州之地，左环沧海，右拥太行，北枕居庸，南襟河济，诚天府之国”。全市平均海拔43.5米。北京平原的海拔高度在20～60米，山地一般海拔1,000～1,500米。 北京市国土面积16410.54平方公里，市区面积12187平方公里，建成区面积1386平方公里 植被覆盖度的测量可分为地面测量和遥感估算两种方法。地面测量常用于田间尺度，遥感估算常用于区域尺度。目前已经发展了很多利用遥感测量植被覆盖度的方法，较为实用的方法是利用植被指数近似估算植被覆盖度，常用的植被指数为NDVI。 2.2 象元二分法 像元二分模型是一种简单实用的遥感估算模型，它假设一个像元的地表由有植被覆盖部分地表与无植被覆盖部分地表组成，而遥感传感器观测到的光谱信息也由这2个组分因子线性加权合成，各因子的权重是各自的面积在像元中所占的比率，如其中植被覆盖度可以看作是植被的权重。 下面简单介绍下在像元二分模型的基础上研究的模型： ………………………………………………(2-1) 其中, NDVIsoil?为完全是裸土或无植被覆盖区域的NDVI值，NDVIveg?则代表完全被植被所覆盖的像元的NDVI值，即纯植被像元的NDVI值。两个值的计算公式为： ………………………(2-3) ……………(2-4) 利用这个模型计算植被覆盖度的关键是计算NDVIsoil和NDVIveg。这里有两种假设： 1）VFCmax=100%，VFCmin=0%。公式（2-1）可变为： ………………………………………………(2-4) NDVImax和NDVImin分别为区域内最大和最小的NDVI值。由于不可避免存在噪声，NDVImax和NDVImin一般取一定置信度范围内的最大值与最小值，置信度的取值主要根据图像实际情况来定。 2）VFCmax=100%，VFCmin=0% 当有实测数据的情况下，取实测数据中的植被覆盖度的最大值和最小值作为VFCmax和?VFCmin，这两个实测数据对应图像的NDVI作为NDVImax?和NDVImin。 2.3 植被指数NDVI 本研究选择NDVI估算植被覆盖度，尽管NDVI增强了近红外与红色通道反射率的对比度，增强了低值部分，导致对少植被区域较高、对多植被取区域较低的敏感度，但有其自身的优点： NDVI是植被生长状况的最佳因子。 经过比值处理的NDVI可以消除与太阳高度角、地形、云等和大气条件有关的辐射条件变化的影响，使在研究区植被覆盖度提取中的误差一定程度上地减小。 北京市的几种典型的地面覆盖类型在NDVI图像上区分明显，植被突出。 鉴于以上的优点，本研究用NDVI提取北京市的植被覆盖度。NDVI，即近红外波段NIR与可见光红波段RED数值之差与这两个波段数值之和的比值即： NDVI=……………………………………………(2-5) 式中：NIR为近红外波段数值；RED为可见光波段数值。 2.4 主要技术流程 图2-1 技术路线 3 数据的预处理 3.1 数据的打开 本研究中图像为美国landsat5卫星所拍摄的TM5数据，包含有7个波段，格式为TIF。带有地理坐标。因本研究所使用的软件为ENVI5.0,在其中可以直接打开，跳过了以往的需要进行波段的叠合过程，具体过程为：Open As→Landsat→Geo TIFF With Metadata 打开影像如下： 图3-1 图像（01） 图3-2 图像（02） 3.2 影像的辐射定标 随着遥感技术的迅速发展, 遥感应用逐渐从定性走向定 量, 辐射定标是实现遥感信息定量化的基本前提。辐射定标是指建立遥感传感器的数字量化输出值DN与其所对应视场中辐射亮度值之间的定量关系。…………………………………………………(3-1) 式中：L指定标后的反射率值；单位（mW/(cm2*Sr*)） gain指增益系数；之间，多分布在高度在5大气校正的目的是消除大气和光照等因素对地物反射的影响，获得地物反射率和辐射率、地表温度等真实物理模型参数，用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响，消除大气分子和气溶胶散射的影响。FLAASH?可以处理任何高光谱数据、卫星数据和航空数据(860nm/1135nm)，这些数据是由HyMAP、AVIRIS、CASI、?HYDICE、HYPERION(