影像中NDVI异常值导致图像直方图形状变形的解决方法.docxVIP

由影像数据经过大气校正后计算NDVI，还是可能出现一些异常值。比如我手中这幅影像，NDVI最大值为465，最小值为-407，在ENVI下统计（Statistics）得到下面这样的直方图：

这个直方图就很难看出影像中NDVI的情况，因为统计的步长大约是3.5，即使把直方图的X轴限定在-1~1也不会有任何变化，如下图。

-

没有改观的原因是：数据中的异常值没有消除，统计还是按3.5的步长走，直方图上显示的无非是前面一个图的局部放大，我们不能更清楚地看到-1~1这个范围内的NDVI频率分布。

解决的办法是消除数据中的异常值，把小于-1的全部变成-1，把大于1的全部变成1。这个可以借助ENVI的bandmath实现，在bandmath里面输入下面这个公式：

(b1lt-1.0)*(-1.0)+(b1ge-1.0andb1le1.0)*b1+(b1gt1.0)*1.0

得到新的计算结果，统计它的直方图，如下图所示。

欢迎下载 2

-

这个结果就与我们的感觉基本一致了。

欢迎下载 3