【GIS】ArcGIS中DEM异常高程点的处理

DEM异常点修正的具体步骤如下：

1. 打开 ArcMap，加载 DEM 数据。

2. 使用 Fill 工具填充坑洼。使用 ArcToolbox 窗口中的 Hydrology 工具集，选择 Fill 工具，将 DEM 中的坑洼填充。填充后可以生成无洼点 DEM，方便后续进行坡度计算和异常点修正。
   

3. 计算坡度和高程变化率。使用 Spatial Analyst 工具集中的 Surface 工具集，选择 Slope 工具和 Aspect 工具，计算 DEM 的坡度和坡向。然后使用 Raster Calculator 工具，计算 DEM 的高程变化率。计算公式为：
   ((DEM - FocalStatistics(Dem, "Rectangle", "MEAN", "NODATA")) / FocalStatistics(Dem, "Rectangle", "MEAN", "NODATA")) * 100
   其中FocalStatistics是进行邻域计算，DEM选择需要计算的 DEM栅格
   

4. 定义异常点。异常点定义为高程变化率绝对值大于一定阈值的像元点。根据实际情况选择一个合适的阈值进行定义。可以通过使用 Raster Calculator 工具，根据高程变化率和阈值的条件，生成异常点掩膜。计算公式为：
   Abs("高程变化率") > "阈值"
   或者
   Con("高程变化率" < "阈值","高程变化率" > "阈值")

5. 根据坡度修正异常点。使用 Con 工具，根据坡度和异常点判断条件，将异常点所在的像元点的高程值修正为周围像元点的高程平均值。具体地，使用 Spatial Analyst 工具集中的 Conditional 工具集，选择 Con 工具。根据坡度和掩膜的条件，使用 Con 工具生成修改后的 DEM。计算公式为：
   Con(("坡度"为合理范围) & ("异常点掩膜"不为0), FocalStatistics("DEM", "Rectangle", "MEAN"), "DEM")

6. 保存修改后的 DEM 文件。使用 Save 工具将修正后的 DEM 保存为新文件，以便后续使用。