[地图代数]处理DEM中的高程异常值——ArcGIS栅格计算的应用

接了一个任务，要处理DEM原始数据中的高程异常值，如图中的异常亮点。

想了一下，以前处理过建筑物附近的DEM铲平，那么高程异常值应该如何处理呢？

显然直接铲平时不太合理的，需要利用异常值周围的高程进行内插（没有高精度的DEM）。

思路如下：

1、首先建立一个面图层，手工圈出异常值的区域。

2、建立一个点图层，在异常值框选区周围采样一些点。获取DEM高程用“添加表面信息”工具。

3、根据点和面图层进行插值处理，利用反距离权重方法，哪种插值方法好需要斟酌，使用局部内插方法。

如果有高精度的DEM，这一步就不用插值了，直接裁剪（栅格计算实现）出来就可以了。

4、将面图层栅格化，这一步之前要将异常区域和DEM图层范围生成（栅格范围 (3D Analyst)）的矩形进行联合（“联合”），分类一，分类二属性值赋值如下，再生成0，1栅格图。

5、利用栅格计算：将原始DEM通过乘法运算归0，接着再运用加法运算添加高程值。至此处理完毕。

参考另外的一篇文章：http://www.cnblogs.com/yhlx125/archive/2013/03/06/2946342.html