基于像元二分模型的植被覆盖度反演以北京市为例.docx

1、基于像元二分模型的植被覆盖度反演以北京市为例基于像元二分模型的植被覆盖度反演-以北京市为例 王玲（西北大学 城市与环境学院，陕西 西安 710127）摘要：采用遥感技术监测植被覆盖度具有重要意义。本文以北京市为例，基于2013年的Landsat8 OLI影像，选取NDVI值为参数，采用像元二分模型对植被覆盖度进行反演，最终反演的结果与实际情况符合，说明采用该方法反演植被覆盖度可行。关键词：植被覆盖度、像元二分模型、NDVI、植被指数引言植被覆盖度（Vegetation fractional cover，简称fc）是指植被（包括叶、茎、枝）在地面的垂直投影面积占统计区总面积的百分比，即植土比。通

2、常林冠称郁闭度，灌草等植被称覆盖度1。它是衡量地表植被覆盖的一个最重要的指标，被覆盖度及其变化是区域生态系统环境变化的重要指示，对水文、生态、全球变化等都具有重要意义2。根据监测手段, 测量植被覆盖度的方法可分为地面测量和遥感测量两大类，测量常用于田间尺度，遥感估算常用于区域尺度。目前已经发展了很多利用遥感测量植被覆盖度的方法，较为实用的地面测量单的方法就是目估法，缺点主要是主观性太强。客观的测量方法有样点法、样方法、样带法等，借助于采样仪器的测量方法，空间定量计、移动光量计等。这些方法虽然提高了测量精度，但野外操作不便，并且成本较高, 难以在大范围内快速提取植被覆盖度。而采用遥感技术为监测大

3、面积区域的植被覆盖度，甚至全球的植被覆盖度监测提供了可能3。目前已有许多利用遥感技术测量植被覆盖度的方法，其中应用最广泛的方法是利用植被指数近似估算植被覆盖度，常用的植被指数为NDVI。一、数据源 本文选取两景覆盖北京市的Landsat8 OLI影像、土地覆盖类型图以及北京行政边界矢量数据为数据源。其中，土地覆盖类型图是作为掩膜文件使用，其目的是为了便于植被覆盖度的估算；北京行政边界矢量数据是为了将两景镶嵌好的影像数据进行裁剪使用的，其目的是裁剪出北京市行政区内的范围。另外，Landsat8 OLI影像是从地理空间数据云网站上下载得到的，其成像时间为2013年10月份。与Landsat7的ET

4、M+成像仪相比，OLI成像仪获取的遥感图像辐射分辨率达到12比特，图像的几何精度和数据的信噪比也更高。OLI成像仪包括9个短波谱段（波段1波段9），幅宽185km，其中全色波段地面分辨率为15m，其他谱段地面分辨率为30m4。表1 Landsat8 OLI陆地成像仪波段参数 二、研究方法 本文反演植被覆盖度所采用的是像元二分模型方法，像元二分模型是一种简单实用的遥感估算模型，它假设一个像元的地表由有植被覆盖部分地表(SV)与无植被覆盖部分地表(SS)组成，而遥感传感器观测到的光谱信息(S)也由这2个组分因子线性加权合成，各因子的权重是各自的面积在像元中所占的比率，如其中植被覆盖度可以看作是植被

5、的权重。因此，像元二分模型的原理如下：遥感传感器观测到的光谱信息(S)由有植被覆盖部分地表(SV)与无植被覆盖部分地表(SS)组成，可得出： S = SV + SS 公式1假设一个像元中有植被覆盖的面积比例为fc , 即该像元的植被覆盖度, 则裸土覆盖的面积比例为1 -fc ，如果全由植被所覆盖的纯像元所得的遥感信息为Sveg , 则混合像元的植被部分所贡献的信息Sv可以表示为Sveg与fc的乘积： Sv =fcSveg 公式2那么， Ss =(1 -fc )Ssoil 公式3将公式2与公式3代入到公式1中，可得到：S =fc Sveg +(1 -fc)S soil 公式4对公式4进行变换,

6、可得以下计算植被覆盖度的公式：fc =(S -Ssoil) (Sveg -Ssoil ) 公式5其中Ssoil 为纯土壤像元的信息, Sveg 为纯植被像元的信息, 因而可以根据公式5利用遥感信息来估算植被覆盖度。将归一化植被指数(NDVI)代入公式5可以被近似为:fc =(NDVI - NDVIsoil ) (NDVIveg -NDVIsoil) 公式6其中, NDVIsoil 为裸土或无植被覆盖区域的NDVI值, 即无植被像元的NDVI 值;而NDVIveg 则代表完全被植被所覆盖的像元的NDVI 值, 即纯植被像元的NDVI 值2。当区域内可以近似取VFCmax=100%，VFCmin=

7、0%，VFC = (NDVI -NDVImin)/ ( NDVImax -NDVImin)，NDVImax 和NDVImin分别为区域内最大和最小的NDVI值。由于不可避免存在噪声，NDVImax 和NDVImin一般取一定置信度范围内的最大值与最小值，置信度的取值主要根据图像实际情况来定；当区域内不能近似取VFCmax=100%，VFCmin=0%，当有实测数据的情况下，取实测数据中的植被覆盖度的最大值和最小值作为VFCmax和VFCmin，这两个实测数据对应图像的NDVI作为NDVImax 和NDVImin。当没有实测数据的情况下，取一定置信度范围内的NDVImax 和NDVImin。VF

8、Cmax和VFCmin根据经验估算。三、数据处理1、数据预处理本文使用的Landsat8 OLI为L1T级别数据，不需做几何校正处理。而北京市需要两景Landsat OLI数据覆盖，因此首先要进行图像镶嵌和裁剪，然后进行大气校正等预处理过程。（1）辐射定标辐射定标是将传感器记录的电压或数字值转换成绝对辐射亮度的过程。目的是消除传感器本身所产生的误差,由于传感器在不断的运行中光学器件性能逐渐退化,因此定标的系数也随之不同,这些定标系数也在不断的更改,在用户获得数据的时候,这些定标系数也在影像的头文件中同时提供给用户。遥感数据辐射定标就是将传感器得到的灰度值转换为星上的辐射亮度值或星上反射率,即表

9、观辐射度或表观反射率。辐射定标主要校正由传感器的灵敏度带来的辐射误差5。其目的是为FLAASH大气校正准备数据：定标符合单位要求的辐射量数据、转换数据储存顺序等。该处理过程在Envi5.2中实现，具体操作：在ENVIToolbox中，选择Toolbox/Radiometric Correction/ Radiometric Calibration，选择*_MultiSpectral多光谱组（7个波段），打开辐射定标工具，对两景影像分别做辐射定标。（2）影像镶嵌因本文所使用的影像数据源是两景Landsat OLI影像，因此需进行影像镶嵌，镶嵌的目的是将不同的影像文件无缝地拼接成一幅完整的包含研究

10、区域的影像。该处理过程在Envi5.2中实现，具体操作：在Toolbox中，选择/Mosaicking/Seamless Mosaic，打开无缝镶嵌工具，然后进行相关参数设置，如下所示：（3）影像裁剪因本文所使用的影像数据包含了北京市行政区划以外的部分地区，因此需进行影像裁剪，以将研究区裁剪出来，并且减小了数据量，加快了数据处理速度，本文使用北京行政边界矢量裁剪图像。过程在Envi5.2中的具体操作如下： 在Toolbox中，选择/Regions of Interest/Subset Data from ROIs，打开裁剪工具：影像裁剪结果如下所示：（4）Flaash大气校正电磁波在大气中的传

11、输和遥感器观测过程中受光照条件以及大气作用等的影响,只有小部分(在0.85um波段80%,在0.45um波段50%)太阳辐射能反射到遥感器,导致遥感器的测量值与地物实际的光谱辐射率不一样。辐射损失主要发生在大气吸收和散射过程，因此地表参数的遥感定量反演研究中,必须纠正目标辐射的不确定性信息6。ENVI中的FLAASH模型是基于MODTRAN4+辐射传输模型，通过参数查找表来进行大气校正的商业化软件。FLAASH大气校正模块支持多种传感器数据，其光谱处理范围0.4m-2.5m，可以有效地去除水蒸气/气溶胶散射效应，同时该方法基于图像像素级的校正，能够解决目标像元和邻近象元的“邻近效应”问题，校正

12、结果精度高，简单易行7。然后，对大气校正前后同一地物的光谱曲线进行对比，这里以植被为例，光谱曲线如下图所示：校正前的植被光谱曲线 校正后的植被光谱曲线2、植被覆盖度估算（1）计算NDVI 本文选取NDVI值为参数，采用像元二分模型对植被覆盖度进行反演，根据植被覆盖度的计算公式可知，要求取植被覆盖度，首先需要计算NDVI。在Envi5.2中的具体操作如下：在Toolbox中，选择Spectral/Vegetation/NDVI，NDVI Calculation Input File面板中，选择LC8_rad_beijing_ref.dat 图像，求算NDVI，如下：NDVI求算结果如下：由于大气

13、校正后的结果有部分像元为负值，主要集中在阴影地区，这部分区域计算得到的NDVI在-1,1之外，为了便于后面的分析，我们这里统一将这部分像元进行处理，即NDVI值大于1的变为1，小于-1的变成-1。在Bandmath中的表达式为：-1b1b11.0，其中b1为植被覆盖度 ；-NaN 可以用掩膜进行处理，即在Build Mask中用 -NaN生成掩膜。去掉异常值之后，并对其进行分类显示，最终得到的植被覆盖度图如下：（5）结果验证目前业内植被覆盖度的验证方法主要是通过野外调查结果进行验证，如下为一种方法：以与地面垂直的角度用数码相机拍摄采样点的地面照片，使用GPS定位获得采样点的经纬度坐标。为获得准

14、确的植被覆盖度实测数据减少像片边缘变形误差，将数码相机得到的数字图像截取长、宽各三分之二的中心地带，采用非监督分类法为10类，并将分类结果分为植被、非植被两类，以此来计算出植被覆盖度。为了保证验证时能够正确定位，野外一般选取3*3个像元大小的样方，即90m*90m，并在样方中均匀拍摄多张照片。取从照片计算的平均值作为样方的植被覆盖，取样方中心点所在3*3像元的植被覆盖度平均值作为对应的遥感估算值进行验证。三、 结论 本文通过归一化植被指数（NDVI）像元二分模型来估算北京市的植被覆盖度，从估算结果来看，北京市市中心的植被覆盖度相比远离市中心地区的植被覆盖度相对较小，尤其是西南地区人口密集、城市

15、化水平高，植被覆盖度相对较小。但是，对于市内绿化程度较高的地区，植被覆盖度相对较高。因此，从实际情况来看，本文中对植被覆盖度的估算结果与实际情况较为符合，说明采用NDVI像元二分模型估算植被覆盖度的效果较好。参考文献1 吴云,曾源,赵炎.基于MODIS数据的海河流域植被覆盖度估算及动态变化分析J.资源科学,2010,32(7):1417-1424.2陈巧,永富. QuickBird遥感数据监测植被覆盖度的研究J.林业学研究,2005,8（4）：375-380.3苗正红，刘志明.基于MOdis NDVI的吉林省植被覆盖度动态遥感监测J.遥感技术与应用，2010，25（3）：387-393.4 初庆伟，张洪群，吴业炜.Landsat-8卫星数据应用探讨J.遥感信息，2013，28（4）：110-114.5周婷婷.遥感影像辐射校正研究与应用D.福建师范大学,2010年6月.6钱金波.祁连山狼毒盖度地面测量与卫星遥感估算D.兰州大学,2009年5月.7任唯敏.植被指数计算区域植被覆盖度的适用性研究D.西北大学,2012年6月.