接了一个任务,要处理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