实验四 区域植被覆盖度反演

一、实验目的要求

1、掌握波段运算与应用的基本流程以及其中各项参数的计算

2、以自己家乡的多幅Landsat8遥感影像（保证研究区域无云），经过裁剪、镶嵌和大气校正等处理后，反演近期家乡的植被覆盖度情况，并总结方法流程。

二、实验数据

1、Landsat 8 OLI_TIRS 数据

三、实验内容

1、波段运算

（1）NDVI计算

打开影像（经过FLAASH大气校正后，且经过镶嵌和裁剪的影像图）baoding-3ref—选择工具箱中spectral其下二级工具vegetation中的NDVI—选择输入文件baoding-3ref，OK—弹窗，NDVI波段红外波段为4，近红外波段为5，选择文件输出路径，命名为NDVI，OK。

（2）NDVI去除异常值及soil和vegetation值的统计

打开工具箱band algebra里的band math—输入表达式(b1 lt -1)*0+(b1 gt 1)*0+(b1 le 1 and b1 ge -1)*b1 点击加入到列表，OK—波段选择NDVI，选择文件输出位置，命名为ndvi-reg ，OK—右键单击该影像，点击快速统计，查看soil值为0.082353，vegetation值为0.537255

（3）VFC的计算

同样打开band math，输入表达式b1*0.082353,添加到列表—波段选择为NDVI—选择输出文件路径，命名为ndvi-soil；接着输入表达式b1*0.537255，添加到列表，波段选择为NDVI--选择输出文件路径，命名为ndvi-vegetation；再输入表达式（float（b1）-float（b2）/float（b3）-float（b2）），添加到列表，b1的波段选择是ndvi-reg,b2的波段选择是ndvi-soil，b3的波段选择是ndvi-vegetation，选择文件输出路径，命名为VFC，OK。

2、伪彩色增强

打开工具箱classification其下二级工具post classification中的raster color slices，选择VFC，OK—点击删除所有颜色片，再手动添加9个—从颜色改变表里找到绿色的CB-Greens,OK—从第一排到九排的颜色片范围分别设置为0-0、0.000001-0.1、0.1-0.2、0.2-0.3、0.3-0.4、0.4-0.5、0.5-0.6、0.6-0.7、0.7-0.764362，第一片颜色设为黑色，第二片颜色设为黄色，第九片设为红色，OK。