1、背景

          相邻两山头之间呈马鞍形的低凹部分称为鞍部。鞍部点是重要的地形控制点，它和山顶点、山谷点及山脊线、山谷线等构成地形特征点线，对地形具有很强的控制作用。因此，因此，对这些地形特征点、线的分析研究在数字地形分析中具有很重要的意义。同时，由于鞍部点的特殊地貌形态，是的鞍部点的提取方法较山顶点低谷底点更难，目前都还存在一定的技术局限性。

2、目的

           利用水文分析的方法提取地形鞍部点；

           通过多种GIS空间分析方法的应用，提高对知识的综合运用能力。

3、要求

          利用水文分析模块和空间分析模块相应功能提取样区地形鞍部点。

4、数据

          25m分辨率的DEM数据，面积约为59平方千米。

       （百度网盘链接：https://pan.baidu.com/s/1o60V6mZrPN_OSrKba8zoAA 提取码：gcm0 ）

5、算法思想

          鞍部具有独特的形态特征，可被认为是原始地形中的山脊和反地形中的山脊回合的地方，因此可通过提取正反地形的山脊线并求其交点，获得鞍部点，鞍部点的提取流程如下图所示。

                                                                            鞍部点的提取流程图

6、操作步骤

        1，正地形、等高线和晕染图的提取

                同山脊线和山谷线的提取一样，由于鞍部点的整体位置是处于山脊上的，需要提取正地形以过滤那些在负地形上的错误的点。正地形的提取方法为：先用焦点统计算出DEM的平均meandem,再使用栅格计算器计算平均meandem和dem的差值，使用重分类将其分为两级，分级界线为 0，将大于 0 的区域赋值为1（即为正地形），小于 0 的区域赋值为0，命名为zhengdixing；另一次将小于0的区域赋值为1（即为负地形），大于0的赋值为0，命名负地形。（具体操作见https://blog.csdn.net/shenshen_xu/article/details/102891772中的【6、1】）。利用【Spatial Analyst】|【表面分析】|【等值线】和【Spatial Analyst】|【表面分析】|【山体阴影】分别提取样区等高距为40m的等高线数据ctour和样区的晕染图hillshade.（具体操作见https://blog.csdn.net/shenshen_xu/article/details/102891772中的【6、2，（6）】）

                                     正地形                                                                                 负地形

                                                                                       等值线和山体阴影

        2，山脊的提取

              山脊的提取与下网页上山脊线的提取一样。分别是进行洼地填充、水流方向提取、汇流累积量的计算、汇流累积量等于0的提取。（具体操作见https://blog.csdn.net/shenshen_xu/article/details/102891772中的【6、2】）

                                                                                                   提取山脊

          3，反地形山脊的提取

               反地形山脊的提取与下网页上山谷的提取过程完全相同。分别是基于原始DEM计算出反地形DEM数据（计算中是利用原始DEM减去3000）；基于反地形DEM数据提取水流方向数据；基于水流方向数据进行汇流累积量的计算；提取汇流累积量为0的栅格。（具体操作见https://blog.csdn.net/shenshen_xu/article/details/102891772中的【6、3】）

                                                                                            提取反地形山脊

          4，鞍部点的提取

               （1）利用【Spatial Analyst 工具】|【地图代数】|【栅格计算器】工具将提取出的山脊线数据和山谷线数据相乘，结果为鞍部区。               

                                                                                   栅格计算器

                                                                                        鞍部区

               （2）利用同样的方法将鞍部区和正地形相乘，就得到了鞍部点的栅格数据anbudian。

                                                                                   鞍部点的栅格数据

               （3）重分类鞍部点栅格数据，将所有0值赋为NODATA数据，属性为一的保持不变，得到Reanbudian

                                                                               重分类

                                                                                             Reanbudian

               （4）将栅格数据Reanbudian转换成矢量数据鞍部点，配合等高线和晕染图对矢量形式的鞍部点数据进行编辑，剔除那些处于样区边缘以及内部的伪鞍部点。

                                                                                    栅格转点

                                                                                        鞍部点

注意：这里山脊和反山脊的提取只需要提取到流量等于0的地方。）