01流程图

利用水文分析提取山脊线及山谷线，山脊线相当于分水线，山谷线相当于山谷线。分水线是水流的起源点，这些栅格的水流方向只存在流出方向而不存在流入方向，所以汇流累积量为零。通过对零值的提取就可以得到山脊线。

山谷线相当于汇水线要用反地形求出，即用较大值减去DEM，DEM中山脊线就成为山谷线，山谷线变为山脊线，用求山脊线的方法求出山谷线，分别利用正反地形求交验证。

DEM进行填洼，利用水文分析求出流向流量，再提取出汇流累积量为零得值与正地形求交，即得到分水线也就是山脊线。用反地形求流向流量提取汇流累积量为零的部分与负地形求交就是山谷线。

将山脊线与山谷线求交得到鞍部点。

图1-1 流程图

02操作步骤

2.1 正地形、等高线和晕渲图的提取

(1)加载DEM数据，在ArcToolbox中选择Spacial Analyst Tools → Neighborhood→Focus Statistics工具，输入dem，利用11*11窗口计算平均值。设置如图2-1所示。

图2-1 焦点统计设置

(2)在ArcToolbox中选择Spacial Analyst Tools →Map Algebra →Raster Calculator工具，对原始数据与焦点统计后的DEM做减法。结果如图2-2所示。

图2-2 减法计算结果

(3)在ArcToolbox中选择Spacial Analyst Tools →Reclass→Reclassify工具，对减法运算结果进行重分类，分级界线为0。将大于0的区域赋值为1，小于0的区域赋值为0即得到正地形；设置如图2-3所示，结果如图2-4所示。

图2-3 正地形重分类设置

图2-4 正地形结果图

(4)生成等值线：在ArcToolbox中选择Spacial Analyst Tools →Surface→Contour工具，输入原始DEM数据生成间距为40的等值线。设置如图2-5所示，结果如图2-6所示。

图2-5 生成等值线设置

图2-6 生成等值线结果图

(5)生成山体阴影：在ArcToolbox中选择Spacial Analyst Tools →Surface→hillshade工具，输入原始DEM，其他设置按默认如图2-7所示，结果如图2-8所示。

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

图2-8 DEM晕渲图

2.2 山脊线的提取

(1)填洼：加载DEM数据，在ArcToolbox中选择Spacial Analyst Tools →Hydrology→Fill工具，输入DEM进行填洼.

(2)流向计算：在ArcToolbox中选择Spacial Analyst Tools → Hydrology → Flow Direction工具，输入填洼后的DEM，得到流向数据，设置如图2-9所示，结果如图2-10所示。

图2-9 流向生成设置

图2-10 生成流向结果图

(3)汇流累积量：在ArcToolbox中选择Spacial Analyst Tools →Hydrology→Flow Accumulation工具，输入流向数据，得到流量数据，结果如图2-11所示。

图2-11 汇流累积量结果图

(4)汇流累积量为0值得提取：在ArcToolbox中选择Spacial Analyst Tools →Map Algebra→Raster Calculator工具，输入公式如图2-12所示，结果如图2-13所示。

图2-12 提取汇流累积量为0的值

图2-13 汇流累积量为0值的结果

(5)利用邻域分析方法，对提取出来的0值结果进行3x3邻域分析，求均值，是数据光滑。设置如图2-14所示，结果如图2-15所示。

图2-14 邻域分析设置

图2-15 0值部分求均值结果

(6)在求均值后的0值汇流累积量数据上单击右键，选择属性，进行重新分级，将数据分为两级，调整临界点，以等值线图和晕渲图进行判断，属性值越接近1越可能是山脊线位置，确定分界阈值为0.5541.设置如图2-16所示。

图2-16 对均值后的数据进行分类

(7)将属性直接近1的的一类属性值设为1，其余的设为0。

(8)在ArcToolbox中选择Spacial Analyst Tools →Map Algebra→Raster Calculator工具，将重分类后的数据与正地形数据相乘，可以消除存在于负地形区域的错误山脊线。然后将计算结果重分类，属性部位1的栅格设为NODATA，就得到了山脊线，结果如图2-17所示。

图2-17 计算出研究区域的山脊线

2.3 反地形山脊的提取

(1)在ArcToolbox中选择Spacial Analyst Tools →Map Algebra→Raster Calculator工具计算反地形，公式为Abs(dem-2000)，得到反地形数据。如图2-18所示。

图2-18反地形数据

(2)流向计算：在ArcToolbox中选择Spacial Analyst Tools →Hydrology→Flow Direction工具，输入反地形数据，得到流量数据，结果如图2-19所示。

图2-19 生成流向结果图

(3)汇流累积量：在ArcToolbox中选择Spacial Analyst Tools →Hydrology→Flow Accumulation工具，输入流向数据，得到流量数据，设置如结果如图2-20所示。

图2-20 汇流累积量结果图

(4)汇流累积量为0值得提取：在ArcToolbox中选择Spacial Analyst Tools →Map Algebra→Raster Calculator工具，结果如图2-21所示。

图2-21 汇流累积量为0值的结果

(5)利用邻域分析方法，对提取出来的0值结果进行3x3邻域分析，求均值，是数据光滑。结果如图2-22所示。

图2-22 0值部分求均值结果

(6)在求均值后的0值汇流累积量数据上单击右键，选择属性，进行重新分级，将数据分为两级，调整临界点，以等值线图和晕渲图进行判断，属性值越接近1越可能是山脊线位置，确定分界阈值为0.65667。

(7)将属性直接近1的的一类属性值设为1，其余的设为0。

(8)在ArcToolbox中选择Spacial Analyst Tools →Map Algebra→Raster Calculator工具，将重分类后的数据与负地形数据相乘，可以消除存在于正负地形区域的错误山脊线。然后将计算结果重分类，属性部位1的栅格设为NODATA，就得到了山谷线，结果如图2-23所示。

图2-23 计算出研究区域的山谷线

2.4 鞍部点的提取

(1)在ArcToolbox中选择Spacial Analyst Tools →Map Algebra→Raster Calculator工具，将山脊线数据与山谷线数据相乘，得到结果为鞍部区。结果如图2-24所示。

图2-24 鞍部区结果图

(2)鞍部区与正地形数据相乘得到鞍部点的栅格数据，重分类栅格数据所有0值赋为NoData，属性为1的值保持不变，得到重分类栅格鞍部点。在ArcToolbox中选择Convertion Tools →From Raster→Raster to Point工具，得到矢量的点。结果如图2-25所示。

图2-25 栅格鞍部点转为矢量结果

(3)配合等高线和晕渲数据，删除处于样区边缘以及内部的伪鞍部点，最后得到鞍部点如图2-26所示。

图2-26 最终得到的鞍部点结果图

2.5 制作专题图

打开Layerout view对话框，分别添加dem，山脊山谷线结果，插入图名图例比例尺指北针等地图制图要素，调整后制作专题图表，结果如图2-27所示。

图2-27 鞍部点专题图

文章转载于：CSDN，作者：BigSun1993

版权归原作者所有，如有侵权请告知删除

- END -

ArcGIS中使用DEM山体阴影制作精美地形图ArcGIS实现由DEM创建明暗等高线7.72GB ARCGIS  ERDAS  ENVI干货教程网盘下载ArcGIS Python资源收集帖在ArcGIS中利用等高线生成DEM技巧教学