计算植被覆盖度

用6S模型来计算植被覆盖度和地表反射率。


计算植被覆盖度
1、数据准备：第3波段和第4波段地表发射率的遥感影像
2、计算NDVI（归一化植被指数）

NDVI=（R4-R3)/(r4+r3) 在ENVI中r为b
3、计算FVC(植物覆盖度）
其中，FVC是植被覆盖度：
FVC=[(NDVI-NDVIs)/(NDVIv-NDVIs)]2


NDVIs=0.05 NDVIv=0.07 (注意：在ENVI中NDVI用b1表示）
FVC=[(b1-0.05)/(0.7-0.05)]*[(b1-0.05)/(0.7-0.05)]














EXP文件表示公式






计算植被覆盖度：计算植被覆盖度：
1，数据准备：TM3、TM4地表反射率。
2、计算NDVI(归一化植被指数)
公式：NDVI=(R4-R3)/（r4+r3）
其中r4与r3 是相应波段的地表反射率。
在envi r为b
3、计算植被覆盖度(FVC)。
FVC=[(NDVI-NDVIs)/(NDVIv-NDVIs)]2

0.05 0.07-0.05
[ndvi-0.05)/(0.07-0.05)]*[(ndvi-0.05)/(0.07-0.05)]
[(ndvi-0.05)/(0.07-0.05)]*[(ndvi-0.05)/(0.07-0.05)]ndvi在envi用b1表示



后续处理公式使1≥FVC≥0：(b1 lt 0)*0+(b1 gt 1)*1+(b1 le 1 and b1 ge 0)*b1



关于逻辑关系的说明：
ge：大于等于0
gt：大于
le：小于或等于
lt：小于
and：且
or：或
eq：等于
ne:不等于















植物覆盖度分级计算公式
b1：处理后的植被覆盖度影像，b2：掩膜影像
0-0.1：第一等级
0.1-0.3：第二等级
0.3-0.5：第三等级
0.5-0.7：第四等级
0.7-1：第五等级
植被覆盖率取值为0且位于黄石市边界范围内的像元取值为0。


(b1 ge 0 and b1 le 0.1 and b2 eq 1)*1+(b1 gt 0.1 and b1 le 0.3)*2+(b1 gt 0.3 and b1 le 0.5)*3+(b1 gt 0.5 and b1 le 0.7)*4+(b1 gt 0.7 and b1 le 1)*5+(b2 eq 0)*0

用上述公式可以得到植被覆盖度分级图。




举例：
计算植被覆盖度分级为1等级范围内的像元覆盖度影像图。【坡段计算器】
用到两个影像：b1（处理后的植物覆盖度影像FVC)和b2（植被覆盖度分级图）

(b2 eq 1)*b1+(b2 ne 1)*7.
表示分级图等于1的像元值给它赋值b1，分级图不等于1的给它赋值3（前面一部分是表示按条件删选的像元，*号是用来进行赋值的）
结果：0-0.1(第1等级）.7。




计算不同等级范围内的像元覆盖度影像图。


(b2 eq 2)*b1+(b2 ne 2)*7.
(b2 eq 3)*b1+(b2 ne 3)*7.
(b2 eq 4)*b1+(b2 ne 4)*7.
(b2 eq 5)*b1+(b2 ne 5)*7.






7、分别计算每个植被覆盖度分级等级内的植物覆盖度的平均值。

平均值得算法：将每一个等级的DN乘以Npts得到的值P求和，再将Npts值求和值Q相除得到FVC平均值
即：FVC=Q/P,其中P=DN*Npts,Q是Npts的总和。


8、利用表格制表，并将其的散点图绘制出来。



9、利用ArcGis制图。