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

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

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

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

思路如下:

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

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

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

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

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

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

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

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

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

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

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

05-08 08:24