文档库

最新最全的文档下载
当前位置:文档库 > 计算植被覆盖度

计算植被覆盖度

用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制图。