(4)Flaash大气校正
ENVI中的FLAASH模型是基于MODTRAN4+辐射传输模型,通过参数查找表来进行大气校 正的商业化软件。实验报告中需写出为什么用进行大气校正。
在 Toolbox 中打开 FLAASH 工具/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction
报告中需要对大气校正前后同一地物的光谱曲线进行对比。
- 6 -
3、植被覆盖度估算 (1)计算NDVI
本文选取NDVI值为参数,采用像元二分模型对植被覆盖度进行反演,根据植被覆盖度的计算公式可知,要求取植被覆盖度,首先需要计算NDVI。在Envi5.1中的具体操作如下:
在Toolbox中,选择Spectral/Vegetation/NDVI,NDVI Calculation Input File面板中,选择大气校正后的图像,求算NDVI,如下:
由于大气校正后的结果有部分像元为负值,主要集中在阴影地区,这部分区域计算得到的NDVI在[-1,1]之外,为了便于后面的分析,这里统一将这部分像元进行处理,即NDVI值大于1的变为1,小于-1的变成-1。使用Bandmath工具,(写出运算公式),得到去除异常值文件。
- 7 -
(b1 lt (-1))*(-1) + (b1 gt (1))*1 + (b1 le(1) and b1 ge(-1))*b1
(2) 生成掩膜文件
该过程主要是为了计算NDVI的最大值、最小值,根据土地利用分类图(共5类,林地、农业用地、城市用地、水体与其他)制作各种土地利用类型的掩膜文件,在Envi5.1中的具体操作如下:在Toolbox中选择/Raster Management/Masking/Apply Mask,打开制作掩膜工具:
- 8 -
得到林地、农业用地、城市用地、水体与其他的掩膜文件。
应用掩膜文件
- 9 -
(3)获取阈值
计算NDVImax和NDVImin值,使用获取的掩膜文件分别对NDVI图像文件进行统计,在一定置信度范围内获取每个掩膜文件(也就是土地覆盖类型)对应的最大和最小NDVI值。
在Toolbox中,选择/Statistics/Compute Statistics,进行统计,然后在统计结果中,取一定的置信度获取最大和最小的NDVI值。林地覆盖区域的统计结果(如下图),选择NDVImin=0.3804,NDVImax=0.8667。同样的方法得到其他地物覆盖类型的NDVI阈值,其中 ,水体没有植被(水藻不属于植被),认为这部分区域的植被覆盖度为0,如下表:
- 10 -