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

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

这个直方图就很难看出影像中NDVI的情况，因为统计的步长大约是，即使把

直方图的X轴限定在-1~1也不会有任何变化，如下图。

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

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

(b1 lt *+(b1 ge and b1 le *b1+(b1 gt *

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

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