山顶点的提取

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/BigSun1993/article/details/50495827

1 流程图
用DEM生成间隔为15m和75m的等高线，生成山体阴影结果图，二者构成地形晕渲图以辅助判断山顶点位置。对DEM数据进行焦点统计分析，以11*11的窗口进行处理，将生成的结果与DEM数据做差重分类后可得到栅格形式的山顶点数据。将栅格数据转为矢量后结合地形晕渲图删除不合理的山顶点，即得到山顶点的分布。流程图如图1-1所示。

图1-1 处理流程图
2、操作步骤
2.1 环境设置
（1）将DEM数据加载到ArcMap中，点击工具栏上的Geoprocessing  Enviornment，设置工作空间路径，将处理范围设置为与DEM一致，栅格分析像元值设为与DEM一致。
2.2 生成等值线
（1）在ArcToolbox中，打开Spatial Analst Tools  Surface Contour工具，输入DEM数据，设置生成等值线的间隔为15m，点击ok。具体设置如图2-1所示，结果如图2-2所示。

图2-1 生成等值线设置

图2-2 15m间隔的等值线
（2）用同样的方法，将间隔设为75m生成等值线，结果如图2-3所示。

图2-3 75m间隔的等值线
2.3 生成山体阴影
（1）在ArcToolbox中，打开Spatial Analst Tools  Surface Hillshade工具，输入DEM数据，其他设置按默认如图2-4所示，点击ok，结果如图2-5所示。

图2-4 生成山体阴影设置

图2-5 山体阴影结果图
2.4 生成山顶点
（1）在ArcToolbox中，打开Spatial Analst Tools  NeighborhoodFocal Statistics工具，输入DEM数据，设置窗口为11*11对DEM数据进行焦点统计，结果如图2-6所示。

图2-6 DEM焦点统计后结果
（2）在ArcToolbox中，打开Spatial Analst Tools  Map Alegbra Raster Calculator工具，输入如图2-7所示公式，求焦点统计后的DEM与原始值之差。

图2-7 计算焦点统计后的数据与原始数据之差
（3）在ArcToolbox中，打开Spatial Analst Tools  Reclass  Reclassify工具，输入上一步差值，将原值为0的部分设为NoData，原值为1的设为1，得到栅格山顶点图。设置如图2-8所示。

图2-8 重分类设置
（4）在ArcToolbox中，打开Convertion Tools  From Raster  Raster to Point工具，将栅格山顶点数据转为矢量格式，设置如图2-9所示。然后结合山体晕渲图删除不合理的山顶点，结果如图2-10所示。

图2-9 栅格转点设置

图2-10 山顶点分布结果图
2.5 制作专题图
将数据页面切换到制图页面，打开图中要表现的要素，插入图名、图例、比例尺、坐标格网和指北针等地图要素，调整好后输出专题图，如图2-11所示。

图2-11 山顶点分布专题图