基于像元二分模型的植被覆盖度反演--像元二分法

一、理论部分

1、植被覆盖度
是指植被（包括叶、茎、枝）在地面的垂直投影面积占统计区总面积的百分比。常用于植被变化、生态环境研究、水土保持、气候等方面。
2、像元二分法模型
像元二分模型是一种简单实用的遥感估算模型，它假设一个像元的地表由有植被覆盖部分地表与无植被覆盖部分地表组成，而遥感传感器观测到的光谱信息也由这2个组分因子线性加权合成，各因子的权重是各自的面积在像元中所占的比率，其中植被覆盖度可以看作是植被的权重。
VFC=(S-Ssoil)/(Sveg-Ssoil) S为遥感信息（即像元值）

3、改进的像元二分法模型

1）改进的像元二分法—遥感信息选择为NDVI
VFC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)
2）
在这里插入图片描述
注：由于本研究区域植被覆盖度不是很密，所以使用NDVI来做，但对于其他研究区域有可能会只用别的，不一定使用NDVI。
4、流程图

二、实践部分

1、数据预处理
本文所用数据为Landsat8，北京地区数据
2、NDVI计算
1）打开数据，加载大气校正后的数据
2）打开Toolbox–搜索NDVI计算工具
在这里插入图片描述

3）查看结果，由于大气校正后，植被、水体部分有负值。
右击–密度分割（New Raster Color Slice）

查看异常值存在地方：删除所有原来的值，重新添加

异常值结果：

4）统计最大值与最小值
打开土地利用分类图
File --Open–土地利用分类图
在这里插入图片描述
5) Toolbox–Classification–Class Statistics

注：只要选择分类的5种，背景不需要

注：查看直方图，一般为单峰，服从正态分布，双峰表示地物特殊。
5）去除异常值
波段运算
-1>b1<1
6）再次打开
Toolbox–Classification–Class Statistics
7）计算参数公式**（选择2%与92%的取值）**
模型参数：
在这里插入图片描述

结合计算公式
(b1-b2)/(b3-b2)
其中：b1–NDVI（去除异常值）；b2–NDVI_soil；b3–NDVI_veg
波段计算

注：计算NDVI与视频学习的不同，还在解决中，但方法思路一致
结果（直接加载提供的数据）：

对结果进行处理，同上：
有波段运算：0>b1<1 b1为植被覆盖度数据
得：
在这里插入图片描述

（密度分割可以使用）
可以划分等级，制作专题图等。