【ArcGIS】关于剔除DEM中的异常值(上)

时间:2024-04-04 07:24:50

前提:得到一幅DEM图像,通过查看属性可知它的地理坐标系为GCS_WGS_1984,加载到arcmap中发现该DEM图像存在几块很明显的异常值区域,加载进Global_Mapper中则可以看出异常区域更加明显。

任务:使用SRTM1 DEM数据替换该DEM图像里存在的异常值。

思路:裁剪出原DEM中存在的异常区域,用裁剪出异常区域继续裁剪/掩膜SRTM1,将SRTM1填补到异常区域,这样做的话填补的边界一定会很明显,所以后期是不是要进行高斯平滑,处理一下边界数据?由于现在还没有处理到后面的流程,所以记录就暂时先这样了。如果有变化或者新的思路后期再继续修改吧。

1、新建一个面图层并赋予它和已有图层一致的地理坐标系【数据管理-投影与变换-定义投影】;

2、绘制面,圈出异常区域【工具条中的编辑器】;

3、提取原栅格范围【3D Analyst-转换-由栅格转出-栅格范围】;

4、异常面域域栅格范围联合形成一个整体面域【分析工具-叠加分析-联合】;

5、打开属性表,新建两个字段“fenlei”&“fenlei2”,并进行如下赋值;

【ArcGIS】关于剔除DEM中的异常值(上)

6、整体面域转成栅格图像【转换-转为栅格-面转栅格】;注意像元大小要与原DEM一致,选择字段时分别选择"fenlei"和“fenlei2”

7、栅格图像二值化成功,蓝色区域值为1,绿色区域值为0;

【ArcGIS】关于剔除DEM中的异常值(上)

8、原DEM与二值化的栅格相乘,得到异常区域值变为0,其余区域保持不变的DEM;黑色区域值为0;

【ArcGIS】关于剔除DEM中的异常值(上)

(未完待续,希望明天能够成功完成这个任务)