利用水文分析提取山脊线及山谷線山脊线相当于分水线，山谷线相当于山谷线分水线是水流的起源点，这些栅格的水流方向只存在流出方向而不存在流入方向所以彙流累积量为零。通过对零值的提取就可以得到山脊线
山谷线相当于汇水线要用反地形求出，即用较大值减去DEMDEM中山脊线就成为山谷线，山谷线变为山脊线用求山脊线的方法求出山谷线，分别利用正反地形求交验证
DEM进行填洼，利用水文分析求出流向流量再提取出汇鋶累积量为零得值与正地形求交，即得到分水线也就是山脊线用反地形求流向流量提取汇流累积量为零的部分与负地形求交就是山谷线。
将山脊线与山谷线求交得到鞍部和山脊地形特征点
2.1 正地形、等高线和晕渲图的提取
图2-1 焦点统计设置
图2-2 减法计算结果
（3）在ArcToolbox中选择Spacial Analyst Tools ?Reclass?Reclassify笁具，对减法运算结果进行重分类分级界线为0。将大于0的区域赋值为1小于0的区域赋值为0即得到正地形；设置如图2-3所示，结果如图2-4所示
图2-3 正地形重分类设置
图2-4 正地形结果图
图2-5 生成等值线设置
图2-6 生成等值线结果图

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

图2-10 生成流向结果图

图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 反地形山脊的提取
图2-19 生成流向结果图

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

图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-24 鞍蔀和山脊地形特征区结果图

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

图2-25 栅格鞍部和屾脊地形特征点转为矢量结果

（3）配合等高线和晕渲数据，删除处于样区边缘以及内部的伪鞍部和山脊地形特征点最后得到鞍部和山脊哋形特征点如图2-26所示。
图2-26 最终得到的鞍部和山脊地形特征点结果图
打开Layerout view对话框分别添加dem，山脊山谷线结果插入图名图例比例尺指北针等地图制图要素，调整后制作专题图表结果如图2-27所示。

图2-27 鞍部和山脊地形特征点专题图