植被覆盖度遥感估算

植被覆盖度是指植被（包括叶、茎、枝）在地面的垂直投影面积占统计区总面积的百分比。容易与植被覆盖度混淆的概念是植被盖度，植被盖度是指植被冠层或叶面在地面的垂直投影面积占植被区总面积的比例。两个概念主要区别就是分母不一样。植被覆盖度常用于植被变化、生态环境研究、水土保持、气候等方面。

植被覆盖度的测量可分为地面测量和遥感估算两种方法。地面测量常用于田间尺度，遥感估算常用于区域尺度。目前已经发展了很多利用遥感测量植被覆盖度的方法，较为实用的方法是利用植被指数近似估算植被覆盖度，常用的植被指数为NDVI。

像元二分模型是一种简单实用的遥感估算模型，它假设一个像元的地表由有植被覆盖部分地表与无植被覆盖部分地表组成，而遥感传感器观测到的光谱信息也由这2个组分因子线性加权合成，各因子的权重是各自的面积在像元中所占的比率，如其中植被覆盖度可以看作是植被的权重。

下面简单介绍下在像元二分模型的基础上研究的模型（李苗苗等）：

VFC = (NDVI - NDVIsoil)/ ( NDVIveg - NDVIsoil)  （1）

其中, NDVIsoil 为完全是裸土或无植被覆盖区域的NDVI值，NDVIveg 则代表完全被植被所覆盖的像元的NDVI值，即纯植被像元的NDVI值。两个值的计算公式为：

NDVIsoil=（VFCmax*NDVImin- VFCmin*NDVImax）/( VFCmax- VFCmin)  （2）

NDVIveg=（(1-VFCmin)*NDVImax- (1-VFCmax)*NDVImin）/( VFCmax- VFCmin)  （3）

利用这个模型计算植被覆盖度的关键是计算NDVIsoil和NDVIveg。这里有两种假设：

• 当区域内可以近似取VFCmax=100%，VFCmin=0%。

公式（1）可变为：

VFC = (NDVI - NDVImin)/ ( NDVImax - NDVImin)   （4）

NDVImax 和NDVImin分别为区域内最大和最小的NDVI值。由于不可避免存在噪声，NDVImax 和NDVImin一般取一定置信度范围内的最大值与最小值，置信度的取值主要根据图像实际情况来定。

• 当区域内不能近似取VFCmax=100%，VFCmin=0%

当有实测数据的情况下，取实测数据中的植被覆盖度的最大值和最小值作为VFCmax和 VFCmin，这两个实测数据对应图像的NDVI作为NDVImax 和NDVImin。

当没有实测数据的情况下，取一定置信度范围内的NDVImax 和NDVImin。VFCmax和 VFCmin根据经验估算。

利用Landsat 8 OLI图像的像元二分模型植被覆盖度遥感估算的处理流程如下：

图： 流程图

• 数据预处理

使用的是Landsat 8 OLI的L1T级别数据，不做几何校正处理。北京市需要两景Landsat OLI数据覆盖，本流程中分别对两个图像进行辐射定标，为了简化操作流程，首先进行图像镶嵌和裁剪，然后进行大气校正处理。

注：如果两景数据成像时间不同，可以按如下流程进行：分别辐射定标、大气校正，然后进行图像镶嵌与裁剪。

• 植被覆盖度估算

在像元二分法模型中，NDVIveg代表着全植被覆盖像元的最大值，由于植被类型的影响，NDVIveg 值也会随着时间和空间而改变。因此，计算植被覆盖度时，既使同一景图像，对于NDVIsoil 和NDVIveg 值不能取固定值。一般需要土壤图和土地利用图，以及野外实测数据，根据上述模型中的两种情况分别求解。由于辅助数据的限制，本专题是利用土地覆盖图分别求解不同土地覆盖类型内的NDVImax 和NDVImin作为NDVI veg和NDVIsoil，这样就可以得到两个NDVI veg和NDVIsoil的参数图像文件。