植被覆盖度计算与分级..

植被覆盖度计算与分级

一、植被覆盖度计算

1、数据处理

打开tm3、tm4地表反射率影像

2、计算NDVI（归一化植被指数）

Basic Tools→Band Math→输入(b4-b3)/(b4+b3)→Ok→按B3选择下面的tm3地表反射率、按B4选择下面的tm4地表反射率→Choose→2000_NDVI_gongshi→打开→Ok

3、计算FVC（植被覆盖度）

FVC= [(NDVI-NDVI S) / (NDVI V-NDVI S)]2

V,S分别是植被和裸土的NDVI值，由于缺乏详细的区域植被和土壤光谱图数据，取值NDVI V=0.70,NDVI s=0.05,

FVC= [(b1-0.05)/(0.7-0.05)]* [(b1-0.05)/(0.7-0.05)]

(此处b1指NDVI影像，FVC衡量植被生长的茂盛程度，取值0—1。)

Basic Tools→Band Math→输入[(b1-0.05)/(0.7-0.05)]* [(b1-0.05)/(0.7-0.05)]→Ok→按B1选择2000_NDVI_gongshi→Choose→2000_FVC_jisuan→打开→Ok

4、处理FVC中DN值大于1的像元

Basic Tools→Band Math→输入(b1 lt 0)*0+(b1 gt 1)*1+(b1 le 1 and b1 ge 0)*b1→Ok→按B1选择2000_FVC_jisuan→Choose→2000_FVC_chuli→打开→Ok

5、建立掩膜进行处理

打开黄石市边界矢量数据；

Vector→Open Vector File→选图（黄石市边界范围.evf)→打开

（3）以黄石市边界矢量数据建立掩膜；

Basic Tools→Masking→Bulid Mask→Display #1→Options →Import EVFS→选图（111）→Ok→Choose→命名

（2000_FVC_chuli_yanmo)→打开→Apply

（4）应用掩膜；

Basic Tools→Masking→Apply Mask→2000_FVC_chuli→Select Mask Bang→2000_FVC_chuli_yanmo→Ok→Ok→Choose→命名（2000_FVC_chuli_again)→打开→Ok

二、植被覆盖度分级

1、打开数据：2000_FVC_chuli_again、2000_FVC_chuli_yanmo

2、分级

（1）制作分级图

Basic Tools→Band Math→输入公式(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→Ok→按B1选择2000_FVC_chuli_agan、B2选择

2000_FVC_chuli_yanmo→Choose→命名为2000_FVC_chuli_fenji→打开→Ok

（2）分级

①一等级：Basic Tools→Band Math→输入公式(b2 eq 1)*b1+(b2 ne

1)*3→Ok→按B1选择2000_FVC_chuli_agan、B2选择2000_FVC_chuli_fenji→Choose→命名为2000_FVC_chuli_again_1dengji→打开→Ok

②相同方法计算2、3、4、5等级

公式分别为：

2等级：(b2 eq 2)*b1+(b2 ne 2)*7

3等级：(b2 eq 3)*b1+(b2 ne 3)*7 4等级：(b2 eq 4)*b1+(b2 ne 4)*7 5等级：(b2 eq 5)*b1+(b2 ne 5)*7 结果图：

2等级

3等级

4等级：

5等级

（3）制作分级数据

右键2000_FVC_chuli_again_1dengji→Quick Stats……→File→Save results to text flie……Choose→命名为2000_FVC_chuli_again_1dengji→打开→Ok