山谷线山脊线提取

山脊线、山谷线和鞍部点的提取山脊线、山谷线和鞍部点的提取一．实习背景山脊线、山谷线是地形特征线，它们对地形、地貌具有一定的控制作用。

它们与山顶点、谷底点以及鞍部点等一起构成了地形及其起伏变化的骨架结构。

因此在数字地形分析中，山脊线和山谷线以及地形特征点等的提取和分析是很有必要的。

相邻两山头之间呈马鞍形的低凹部分称为鞍部，鞍部是两个山脊和两个山谷会合的地方。

鞍部点是重要的地形控制点，它和山顶点、山谷点以及山脊线、山谷线等构成的地形特征点线，具有对地形具有很强的控制作用。

因此，对这些地形特征点、线的分析研究在数字地形分析中具有很重要的意义。

同时，由于鞍部点的特殊地貌形态，使得鞍部点的提取方法较山顶点和山谷的提取更难，目前没有什么有效的方法来提取鞍部点，利用水文分析的方法可以来提取一些鞍部点，但是它还是具有一定局限性。

二．实习目的（1）熟练掌握基于DEM利用ArcGIS进行提取相关地形特征的方法与原理；（2）深入认识山脊线、山谷线和鞍部点3个基本地形特征；三．实习内容１.提取ｄｅｍ数据的ＳＯＡ２基于地形表面的几何形态分析方法提取山脊线山谷线3.基于DEM水文分析方法提取山脊线山谷线4.鞍部点的提取四．实习数据DEM五．实习工具Surface Analyst，model工具六．实习步骤1.提取DEM的SOA数据A．求取原始DEM数据层的最大高程值，记为H；通过Spatial Analysis 下的栅格计算器 Calculator，公式为（H－DEM），得到与原来地形相反的 DEM数据层，即反地形DEM数据；B．基于反地形 DEM数据求算坡向值；C．利用 SOA 方法求算反地形的坡向变率，记为 SOA2，由原始DEM数据求算出的坡向变率值为 SOA1；D.在 Spatial Analysis下使用栅格计算器 Calculator，公式为 SOA ＝（（[SOA1]+[SOA2]）-Abs（[SOA1]-[SOA2]））/ 2，即可求出没有误差的 DEM 的坡向变率，2.利用基于地形表面的几何形态分析方法提取山脊线山谷线（1）山脊线的提取其中在focal statistics中选择3*3的窗口，类型选择为mean,在minus中再次导入DEM，进行计算。

山脊线山谷线提取实验报告实验容描述：山脊线和山谷线构成了地形起伏变化的分界线（骨架线），因此它对于地形地貌研究具有重要意义；另一方面，对于水文物理过程研究而言，由于山脊、山谷分别代表示分水性与汇水性，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。

本次实验通过某区域栅格DEM掌握山脊线和山谷线这两个基本地形特征信息的理论及其基于DEM的提取方法与原理；同时，熟练掌握利用ArcGIS软件对这两个地形特征信息的提取方法。

实验原理：1.本实验基于规则格网DEM数据使用平面曲率与坡形组合法提取山脊线和山谷线，首先利用DEM数据提取地面的平面曲率及地面的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。

实际应用中，由于平面曲率的提取比较繁琐，而坡向变率（SOA）在一定程度上可以很好地表征平面曲率。

因此，提取过程中可以SOA代替平面曲率。

2.主要用到以下理论知识：1）坡向变率：是指在提取坡向基础上，提取坡向的变化率，亦即坡向之坡度（Slope of Aspect，SOA）。

它可以很好地反应等高线弯曲程度；2）反地形DEM数据：求取原始DEM数据层的最大高程值，记为H，通过公式（H-DEM），得到与原来地形相反的DEM数据层，即反地形DEM数据；3）地面坡向变率SOA：地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部微小围坡向的最大变化情况。

但是SOA在提取过程中在北面坡将会有误差产生，所以要将北坡坡向的坡向变率误差进行纠正，其公式为：SOA=(( [SOA1]+[ SOA2] )-Abs( [SOA1]-[ SOA2] ))/2其中：SOA1为原始DEM数据层坡向变率，SOA2为反地形DEM数据层坡向变率。

4）焦点统计5）ArcScan自动矢量化流程图、实验步骤：1.相对路径2.加载数据3.提取原始dem的坡向（利用dem数据--空间分析--表面分析--坡度工具，命名为Aspect）4.提取原始DEM数据的坡向变率（利用3中生成的Aspect图层--空间分析--表面分析--坡度工具，命名为SOA1）5.提取反地形DEM数据（栅格计算器--输入公式H-DEM）1）找出DEM最大高程值（右键属性---找出数据源中最大值为1153.8）2）栅格计算器提取反地形DEM数据（输入公式1153.8 - "dem"，命名为INdem）6.提取反地形DEM数据的坡向值7.计算反地形DEM数据的坡向变率8.计算进行误差纠正的地面坡向变率（栅格计算器--输入公式(("SOA1" + "SOA2") - Abs("SOA1" - "SOA2")) / 2）9.邻域分析（原始dem--邻域分析--焦点统计focal statistics（统计原始dem的平均值）---设置统计类型为平均值mean，邻域类型为矩形(也可为圆形)，邻域大小为3*3（我发现邻域越大越模糊）(11*11)，则可得到一个邻域为3*3(11*11)的矩形的平均数据层，命名为mean10.计算正负地形分布区域（空间分析--地图代数--栅格计算器---输入公式为"dem" - "mean"，命名为Dvalue(差值)）11.利用栅格计算器提取山脊线（公式为"SOA" > 70 & "Dvalue" > 0这是错的！！要加括号！！("SOA" > 70) &( "Dvalue" > 0)）和山谷线（("SOA" > 70) & ("Dvalue" < 0)）12.利用ArcScan自动矢量化得到山脊线山谷线的矢量图层1）在ArcCatalog中新建（方法有两种：右击文件夹--new--shapefile！或者是右击geodatabase--new--feature class（新建要素类））山脊线图层（名称为shanjiline，类型为线）方法1：new--shapefile方法2：new--feature class（但是这种方法下的线图层，在自动矢量化山脊线后无法读到这个图层，所有还是选择方法1---这是因为栅格图层和矢量图层不能放在同一个geodatabase里面么？？？？？？？）2）打开开始编辑3）勾选扩展工具中的自动矢量化工具ArcScan4)在菜单栏空白处右击勾选ArcScan，打开ArcScan工具条--单击自动矢量化下的生成要素打开生成要素对话框即可生成自动矢量化后的矢量山脊线5)用同样的方法生成矢量山谷线13.制作立体图。

山谷线、山脊线提取自动提取山脊线和山谷线arcmap 自动提取山脊线和山谷线的方法1 平面曲率与坡形组合法基于规则格网DEM是最主要的自动提取山脊线和山谷线的方法，从算法设计原理上来分，大致可以分为以下五种：1) 基于图像处理技术的原理；2) 基于地形表面几何形态分析的原理；3) 基于地形表面流水物理模拟分析原理；4) 基于地形表面几何形态分析和流水物理模拟分析相结合的原理；5) 平面曲率与坡形组合法。

平面曲率与坡形组合法提取的山脊、山谷的宽度可由选取平面曲率的大小来调节，方法简便，效果好。

该方法基本处理过程为：首先利用DEM数据提取地面的平面曲率及地面的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。

实际应用中，由于平面曲率的提取比较繁琐，而坡向变率（SOA）在一定程度上可以很好地表征平面曲率。

因此，下面的提取过程以SOA代替平面曲率。

具体提取过程为：1）激活DEM 数据，在Spatial Analysis 下使用surface 菜单下的Derive Aspect 命令，提取DEM 坡向层面，记为A；2）激活A 层面，在Spatial Analysis 下使用surface 菜单下的Derive Slope 命令，提取A 层面的坡度信息，记为SOA1；3）求取原始DEM 数据层的最大高程值，记为H；通过Spatial Analysis 下的栅格计算器Calculator，公式为（H－DEM），得到与原来地形相反的DEM 数据层，即反地形DEM 数据；4）基于反地形DEM 数据求算坡向值；5）利用SOA 方法求算反地形的坡向变率，记为SOA2；6）在Spatial Analysis 下使用栅格计算器Calculator，公式为SOA ＝（（[SOA1]+[SOA2]）-Abs （[SOA1]-[SOA2]））/ 2，即可求出没有误差的DEM 的坡向变率SOA；7）激活原始DEM 数据，在Spatial Analysis 下使用栅格邻域计算工具Neighborhood Statistics；设置Statistic type 为平均值，邻域的类型为矩形（也可以为圆），邻域的大小为275×275 MAP，则可得到一个邻域为275×275 MAP的矩形的平均值层面，记为B；8）在Spatial Analysis 下使用栅格计算器Calculator，公式为C ＝[DEM]-[B]，即可求出正负地形分布区域，9）在Spatial Analysis下使用栅格计算器Calculator，公式为D ＝[C] >0 & SOA > 70，即可求出山脊线；10）同理，在栅格计算器Calculator 中，修改公式为D ＝[C] < 0 & SOA > 70，即可求出山谷线地形特征信息提取（山谷线、山脊线）方法1：SOA法SOA法原理：山谷线和山脊线实质是平面曲率发生突变的地方，所以用SOA来近似平面曲率，提取其中变法大的就是山谷线和山脊线，其中山谷线对应的是负地形中SOA较大的值，山脊线对应的是正地形中SOA较大的值。

ArcGIS中利用水文分析提取山脊线山谷线1 流程图利用水文分析提取山脊线及山谷线，山脊线相当于分水线，山谷线相当于山谷线。

分水线是水流的起源点，这些栅格的水流方向只存在流出方向而不存在流入方向，所以汇流累积量为零。

通过对零值的提取就可以得到山脊线。

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

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

用反地形求流向流量提取汇流累积量为零的部分与负地形求交就是山谷线。

图1-1 流程图2 操作步骤2.1 正负地形求取（1）加载DEM数据，在ArcToolbox中选择Spacial Analyst Tools → Neighb orhood→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所示。

将大于0的区域赋值为0，小于0的区域赋值为1即得到负地形，设置如图2-5所示，结果如图2-6所示。

图2-3 正地形重分类设置图2-4 正地形结果图图2-5 负地形重分类设置图2-6负地形结果图2.2 山脊线的提取（1）填洼：加载DEM数据，在ArcToolbox中选择Spacial Analyst Tools ◊Hydrology◊Fill工具，输入DEM进行填洼，设置如图2-7所示。

山脊线和山谷线提取插件版本V1.0操作说明1/6山脊线和山谷线提取插件操作说明1引言1.1编写目的编写本使用说明的目的是充分叙述本脚本软件所能实现的功能及其运行环境，以便使用者了解本脚本软件的使用范围和使用方法。

1.2编写背景山脊线和山谷线的提取在很多行业，例如环境保护、水文水利、农业等行业应用很广泛。

目前主要依赖ArcGIS中的空间分析工具来完成，需要使用13个空间分析工具来完成，费时费力，同时会产生许多中间过程文件。

本插件使用Python 语言，基于ArcGIS10进行二次开发，基于DEM数字高程模型，使用1个脚本来实现山脊线和山谷线的提取。

技术问题可以在QQ群讨论：9609331152运行环境2.1软件环境win7及以上版本、ArcGIS10.0及以上版本2.2硬件环境CPU：2.4GHz以上硬盘：至少4G以上的空闲空间内存：至少4G的空闲内存显示器分辨率：1280×600以上山脊线和山谷线提取插件版本V1.0操作说明1/63运行原理本插件运行原理如下图所示：图1：山谷线和山脊线提取脚本运行过程（1）坡向分析对输入的地形图进行坡向分析（Aspect），结果记为A。

結果示例如下：（2）坡度分析Slope对第一步的结果进行坡度分析，结果记为SOA1。

（3）反地形DEM求取原始DEM数据层的最大高程值，记为H；通过栅格计算器计算，公式为H-DEM，得到与原来地形相反的DEM数据层，既反地形DEM。

结果记为rastercalc。

（4）反地形的坡向分析求反地形DEM的坡向。

结果记为Aspect_raste1.（5）反地形的坡度分析求反地形的坡向变率，记为SOA2（6）然后利用soa1和soa2求得没有误差的DEM的坡向变率（注意大小写一致）利用栅格计算器，公式：(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/2，结果记为SOA。

（7）邻域计算Block Statistics设置Statistic type（默认为MEAN平均值），领域的类型默认为矩形（也可以设置为其他类型），领域大小自行设置（可以是3×3，下图设置为275×275）。

山谷（脊）线的提取方法研究作者姓名XXX 专业测绘工程指导老师XXX摘要近年来对DEM的研究，主要是从其中提取需要的有用信息，其中山脊线与山谷线这类重要的地形信息是近几十年来国内外研究比较热门的方向。

论文首先介绍了ArcGIS关于山谷线山脊线的提取原理与步骤，然后又引入TauDEM软件，着重介绍了其关于山脊线与山谷线提取的原理，包括8方向（D8）流水方向计算以及D8有效邻域分析的原理与实现方法。

通过对负地形的生成、水流方向上洼地提取及填充、选择合适的阈值、D8水流方向的计算理论、D8有效邻域的分析与处理等的研究，提出用计算汇流量的原理提取等高线。

然后通过二值化舍去非山谷（脊）线的数据信息，即可得到清晰的山谷(脊)线。

关键词：数字高程模型，地形分析，山谷（脊）线，水流方向Study On Methods of Valley ( Ridge)Lines ExtractionAbstract：In recent years, the research based on DEM are mostly about extracting useful information from it,and extracting valley (ridge) lines is one of the most extensive direction both at home and abroad.Firstly,the article gives a detailed description of the extracting principles and steps on ArcGIS,and then describes the TauDEM software,which emphatically introduces the principle of the extraction, including eight directions (D8) water computing as well as analysis and realization methods of the D8 contribution area.This paper mainly studies the generating of the negative terrain,the extraction and filling of Depression in the direction of flow,choosing a suitable threshold to determine the accuracy of grid data and suitability,the calculation of the D8 flow, and the analysis of the D8 contribution areas,thus proposing a new principle to extract the valley (ridge) lines based on calculating the flow.Then binaries and eliminats non-data of valley (ridge) lines.Finally we can get the sharp-edged valley (ridge) lines.Key words: Digital Elevation Model , Terrain Analysis，TauDEM ，Valley (Ridge) Line , The Direction Of Flow目录第1章前言 (1)1.1 山谷（脊）线提取的目的和意义 (1)1.2 山谷（脊）线提取技术的研究现状 (2)1.3山谷（脊）线提取的主要方法分类 (4)1.4本文研究内容及方法 (5)第2章基于ArcGIS的山谷（脊）线提取方法 (6)2.1 ArcGIS软件简介 (6)2.2 ArcGIS山谷（脊）线提取的基本原理 (6)2.3 ArcGIS山谷（脊）线提取的方法步骤 (7)第3章基于TauDEM的山谷（脊）线提取方法 (17)3.1 TauDEM软件简介 (17)3.2 TauDEM山谷（脊）线提取的基本原理 (18)3.3 TauDEM山谷（脊）线提取的方法步骤 (20)第4章山谷（脊）线提取结果的对比分析 (26)4.1 TauDEM与ArcGIS山谷（脊）线提取结果对比分析 (26)4.2 TauDEM山谷（脊）线提取结果优势总结 (28)4.3 ArcGIS与TauDEM的结合处理 (29)结论 (31)致谢 (33)参考文献 (34)第1章前言1.1 山谷（脊）线提取的目的和意义山谷线和山脊线是重要的地形特征信息，对地形及其相关领域的研究有重要的意义，山谷（脊）线的提取一直都受到相关研究人员的广泛关注。

山脊线山谷线提取实验报告实验内容描述:山脊线和山谷线构成了地形起伏变化的分界线（骨架线），因此它对于地形地貌研究具有重要意义；另一方面，对于水文物理过程研究而言，由于山脊、山谷分别代表示分水性与汇水性，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。

本次实验通过某区域栅格DEM掌握山脊线和山谷线这两个基本地形特征信息的理论及其基于DEM的提取方法与原理；同时，熟练掌握利用ArcGIS软件对这两个地形特征信息的提取方法。

实验原理：1•本实验基于规则格网DEM数据使用平面曲率与坡形组合法提取山脊线和山谷线，首先利用DEM数据提取地面的平面曲率及地面的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。

实际应用中，由于平面曲率的提取比较繁琐，而坡向变率（SOA）在一定程度上可以很好地表征平面曲率。

因此，提取过程中可以SOA代替平面曲率。

2•主要用到以下理论知识：1 ）坡向变率：是指在提取坡向基础上，提取坡向的变化率，亦即坡向之坡度（Slope of Aspect，SOA）。

它可以很好地反应等高线弯曲程度；2）反地形DEM数据：求取原始DEM数据层的最大高程值，记为H,通过公式（H-DEM）, 得到与原来地形相反的DEM数据层，即反地形DEM数据；3）地面坡向变率SOA:地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部微小范围内坡向的最大变化情况。

但是SOA在提取过程中在北面坡将会有误差产生，所以要将北坡坡向的坡向变率误差进行纠正，其公式为：SOA=（（ [SOA1]+[ SOA2] ）-Abs（ [SOA1]-[ SOA2] ）2其中：SOA1为原始DEM数据层坡向变率，SOA2为反地形DEM数据层坡向变率。

4）焦点统计5）A rcScan自动矢量化流程图实验步骤：1•相对路径2•加载数据3•提取原始dem的坡向（利用dem数据--空间分析--表面分析--坡度工具，命名为9 £p曰吕I Analyst Toolslbx H 竜Conditional匡心Density1+ Dis-tance1+ W EwtrartionS. %■ G电rer^li23tian + G refund water 匡帝Hydrology l±电Interpolation 1+ 4$^- Local 1+竜Map Algebra 匡雜Matfi ld_ Multivanate l~《NeighLorhoed S ■ Overlay 11 tSv RasterCreation 田心Redact1+ ❸Solar Radiation 1 $! SurfaceConicurGonlour ListConlour with Barriers CurvatureCut FillHillshadeOb server PointsSlopeJf (Views liedAspect)X4•提取原始DEM数据的坡向变率（利用3中生成的Aspect图层--空间分析--表面分析--坡度工具，命名为S0A1）G ❺ Spatial Analyst Tools.tbxE ConditionalS 詹Density3 弧Distance3 Extract!on® S* Generalization♦J © GroundwaterS ■ Hydrology® W? InterpolationS @ Local3 Map AlgebraE 心MathS3 © Multiva riate® d Neighborhood® OverlayS 旳• Raster Creatio n3 • Reclass田& Solar Radiation曰心Surface 气Aspect 气Contour J Con tour List 弋Con*tour wrth Barriers o Curvature 气Cut Fill x - HillshadeViewshed5•提取反地形DEM数据（栅格计算器--输入公式H-DEM）1）找出DEM最大高程值（右键属性---找出数据源中最大值为1153.791870117188 ）2）栅格计算器提取反地形DEM数据（输入公式1153.791870117188 - "dem"，命名为INdem）Spatial Analyst Tools.tbxE 耐ConditiQnalB Densityl±i Di stancel+i 笔/ EjctractionFl 睦i Generalizaticn © 薦i：Groundwateir l±i 氐HydrologyFi Interpolation0 曹LoralLI 転’ Map Al get ra6•提取反地形DEM数据的坡向值■VI(■ri&7•计算反地形DEM 数据的坡向变率3回Canc& Eriwirorments... Show Help »Input rAEtsrOutput ratt&rF : \33 \e\New File Ge Q latabe . gjdb\AspectJNdem8•计算进行误差纠正的地面坡向变率(栅格计算器--输入公式(("S0A1" + "S0A2") - Abs("S0A1" -"S0A2")) / 2)9•邻域分析（原始 dem--邻域分析--焦点统计focal statistics （统计原始dem 的平均值）---设 置统计类型为平均值mean ，邻域类型为矩形（也可为圆形）,邻域大小为3*3 （我发现邻域越大越模糊）（11*11），则可得到一个邻域为 3*3（11*11）的矩形的平均数据层，命名为mean■■■ar■*jrii■H1h ■■rtryE）导Spatial Analyst Tools.tbxS ConditionalS Densityl±l & Distancei*l 心Extraction® © GeneralizationS $9 GroundwaterS • Hydrology® © Interpolation® & Local□ & Map Algebra Raster CalculatorQ 飪MathGO 抵（Ylultiyariate□心Neighborhood "弋Block Statistics 气Filter Focal Flow气Point Statisticsz ccal Statistics10.计算正负地形分布区域（空间分析 --地图代数--栅格计算器---输入公式为"dem" - "mean",命名为Dvalue （差值））EE 叵 0000□□□□ESQ i ° innnnnnTdem T'「帼酮*Output rasterI ：\33\*\Nsw File Gtod^t&base. ;db\DvalueM SL > Algebra expressionCancelEnwironiments...Show He 单 >>Conditional —*ConPickSetfJullMathAbsExp("SOA" > 70) &( "Dvalue" > 0))和山谷线(("SOA" > 70) & ("Dvalue" < 0))L*4iElt TjI I11.利用栅格计算器提取山脊线（公式为 "SOA" > 70 & "Dvalue" > 0这是错的!!要加括号！ ！ Output ras tmirfSOA' > 7Q) &{ lvalue > 0)Map Algebra expressionI:\33\e\Niw Filfe ^dbVshanjiOK Cancel Eriwronmerite...1 1Show Het »Conditional — 厶Con PickSetNull Math Abs Exp〒fSOA ・ > 70〉&CDvalue" < 0)Output rasterF:\33\e\New File Geod&t&b&se. gdb\shw.gu|OKCancel Envronments...Show Help »" Raster CalculatorJ■,Map Algebra expressionESE QESQi ° innrnrnn0]12•利用ArcScan 自动矢量化得到山脊线山谷线的矢量图层 1）在 ArcCatalog 中新建（方法有两种：右击文件夹--new--shapefile !或者是右击geodatabase--new--feature class （新建要素类））山脊线图层（名称为shanjiline ,类型为线）方法 1: new--shapefileFelde-rFile Geoddiabase zyjrA 监JS站Pers&naIDdtmbase Connection... aanewdontTo do doArcGIS Server Connectio-n... Layer...Group LayerbnewSpy1 Shapefile..1 Turn Ft• ToolbcEH d BASE ®LAS Da •Add ©Com pcNew ShapefileCrates 召 f]*vr shspefiletdset；5 Lot:ator…,]site Addre-&5 Locator..闌 XML DocumentDeleteRename Refr«5h New ItemDescnpticn...tf443B01.425 41624&0.765 MetePython TaalboK方法2：new--feature class （但是这种方法下的线图层，在自动矢量化山脊线后无法读到这个图层，所有还是选择方法1---这是因为栅格图层和矢量图层不能放在同一个geodatabase 里面么）2）打开开始编辑ArcSca n 3）勾选扩展工具中的自动矢量化工具匚jstormize Wi n d ows HelpToolbars ►Evtensiori 乩“Add-In Manager.*.C-u atomize Mode.,,Style Manager...ArcMap Options...4）在菜单栏空白处右击勾选ArcScan,打开ArcScan工具条--单击自动矢量化下的生成要素打开生成要素对话框即可生成自动矢量化后的矢量山脊线3D AnalystAdvanced EditingAnimationArcScanCOGOData Driven PagesData Frjm* TookDistributedshaqi T I 53 Vedodzalion ' ：Tl 匸 # * I RMitr Clean up ■ I Cell Selection 章0 斗耳屮J S Verrorization ▼ H:Vcrtorzatioii Settings...Sihcnv PreviewF Ge^eraTeFeatuires,,.Opticrm Raster Cleanup * Cel Selection ，w E 础卜耳十* f5）用同样的方法生成矢量山谷线sbargu VcctcrizatonVectorizatior Setting j・Show Preview■ Generate Features-・Options...• Seed -eettre Template °■W ■塞'Swrch、• g R€i3r;dire兀9壮shorjilirK一穴*13•制作立体图。

师大学H a n g z h o u N o r m a l U n i v e r s i t y《GIS分析与建模》实验报告学院班级：理学院地信141实验名称：利用水文分析方法提取山脊线和山谷线学生：文彪学生学号：2014212425授课老师：潭高提交日期：2016年12月1日目录1 实验目的 (3)2 实验步骤 (3)2.1 正负地形提取 (3)2.2 山脊线提取 (9)2.3 山谷线提取 (24)3 总结 (25)3.1 为何进行洼地填充 (25)3.2 与正地形相乘理解 (26)1实验目的●了解基于DEM水文分析方法提取山谷线和山脊线的原理，●掌握水流方向、汇流累积量提取原理及方法2实验步骤2.1正负地形提取原始dem图像2.1.1焦点统计参数如下邻域分析结果图，发现与原图相比，平滑了很多。

2.1.2栅格运算2.1.3重分类2.1.3.1正地形2.1.3.2负地形2.2山脊线提取2.2.1洼地填充洼地填充参数面板，因为是所有的洼地都填充，所以Z limit为默认。

2.2.2无洼地水流反向计算水流方向计算，输入参数为洼地填充好的dem数据，填如输出文件及其路径。

2.2.3汇流累计量计算2.2.4汇流累计量0值提取汇流累积量为0时，表示该地没有雨水聚集，所以应当为山脊。

2.2.53*3领域分析使数据光滑2.2.6Dem等值线图等值线间距为15米2.2.7山体阴影2.2.8Neiborfacc0重分级分类结果与等值线和山体阴影相比，确定最佳分类阈值。

注：可以将Neiborfacc0的0到所设阈值的部分设为no color,这样更加方便与山体阴影比较。

确定所设阈值为0.55412.2.9二值化的neiborfacc0进行重分类2.2.10与正地形相乘2.2.11栅格运算2.2.12重分类2.2.13山脊线提取结果图2.3山谷线提取山谷线提取结果图3总结3.1为何进行洼地填充来个比喻，假如你在沙地上踩两个脚印，那么，在两个脚印的中间的那条线的汇流累积量为0，进行提取的时候，就会误认为该条线为山脊，显然是不合理的，所以要进行洼地填充。