ArcGIS实验-Ex18-利用水文分析方法提取山脊、山谷线

arcgis 山体范围提取摘要：1.引言2.ArcGIS 简介3.山体范围提取的方法4.提取山体范围的步骤5.结论正文：【引言】在地理信息系统(GIS) 中，山体范围提取是一项非常重要的任务。

山体是地球表面的重要地貌特征之一，对于地理学、地质学、环境科学等领域的研究都具有重要意义。

ArcGIS 是一款功能强大的GIS 软件，可以进行山体范围提取等空间分析任务。

本文将介绍如何使用ArcGIS 提取山体范围。

【ArcGIS 简介】ArcGIS 是由美国Esri 公司开发的一款GIS 软件，具有强大的空间数据处理和分析能力，被广泛应用于地理学、地质学、环境科学、城市规划等领域。

ArcGIS 包括桌面、服务器和移动端等多个平台，支持多种数据格式和空间分析算法，可以满足不同用户的需求。

【山体范围提取的方法】山体范围提取的方法有很多种，其中比较常用的方法是基于DEM(数字高程模型) 的山体提取方法。

该方法通过对DEM 进行处理，提取出山体的轮廓线和高度等信息，从而确定山体的范围。

【提取山体范围的步骤】下面是使用ArcGIS 提取山体范围的具体步骤:1.准备数据：需要一幅包含地形信息的DEM 数据，可以使用遥感影像、地形图等数据源获取。

2.填充DEM:使用ArcGIS 中的"Fill"工具，对DEM 数据进行填充，以生成完整的地形表面。

3.计算坡度：使用ArcGIS 中的"Slope"工具，计算填充后的地形表面的坡度。

4.确定山体范围：根据坡度的大小，将地形表面划分为平地、缓坡、陡坡和山峰等不同区域。

可以使用ArcGIS 中的"Raster Calculator"工具进行操作。

5.提取山体轮廓线：对山峰区域进行处理，提取出山体的轮廓线。

可以使用ArcGIS 中的"Hillshade"工具，将山峰区域转换为轮廓线。

6.聚合轮廓线：对提取出的山体轮廓线进行聚合，以生成最终的山体范围。

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

本次实验通过某区域栅格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地形分析DEM提取坡度(1)新建地图文档，加载【ArcGIS地形分析--TIN及DEM的生成，TIN的显示】经验教程中得到的DEM数据：huainan(2)在【ArcToolbox】中，执行命令[3D Analyst工具]——[栅格表面]——[坡度]，参照下图所示，指定各参数：执行后，得到坡度栅格Slope_tingri1：坡度栅格中，栅格单元的值在[0 -60] 度间变化右键点击图层[Slope_TinGrid]，执行[属性命令]，设置图层[符号系统]，重新调整坡度分级。

将类别调整为5，点[分类]按钮，用手动分级法，将中断值调整为：10，20，30，40，60。

DEM提取山顶点1.添加dem数据，制作15m和75m等高线。

2.制作阴影阴影图像【空间分析】----【表面分析】----【地表阴影】，生成地表阴影图像hillsha。

3.提取栅格数据的有效区域。

【空间分析】----【地图制图】----【栅格计算器】。

“要提取的文件名”=“huainan54”>=0，（注意：红色等号是1个=，而不是栅格计算器中的2个==）。

“back”=“huainan54”>=0，生成back文件。

5.按照等高线75m等高线15m，back，hillsha叠放。

1.提取dem数据中的最大值。

【空间分析工具】----【邻域分析】----【块统计】，生成maxpoint文件参数设置如下：Maxpoint7.提取山顶点。

【空间分析】----【地图制图】----【栅格计算器】。

输入命令：sd=（[axpoint]-[淮南54]）==0,生成sd文件。

8.山顶点栅格文件二值化。

【空间分析工具】----【重分类】----【重分类】。

生成re-sd9.生成山顶点栅格数据转换为shapefile。

【转换工具】----【栅格转换】----【栅格数据转点】，生成山顶点矢量文件。

ArcGIS利用水文分析方法提取山脊、山谷线提取方法大致可以分为以下五种：1) 基于图像处理技术的原理；2) 基于地形表面几何形态分析的原理；3) 基于地形表面流水物理模拟分析原理；4) 基于地形表面几何形态分析和流水物理模拟分析相结合的原理；5) 平面曲率与坡形组合法。

实例与练习练习1. 利用水文分析方法提取山脊、山谷线1．背景：山脊线、山谷线是地形特征线，它们对地形、地貌具有一定的控制作用。

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

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

2．目的：理解基于DEM结合水文分析的方法提取出研究区域的山脊线和山谷线的原理；掌握水流方向、汇流累积量的提取方法以及它们的提取原理；能将水文分析的方法和其它的空间分析方法相结合以解决应用问题。

3．要求：（1）利用水文分析思想和工具提取研究区域的山脊线；（2）利用水文分析思想和工具提取研究区域的山谷线。

4．数据：一幅25m分辨率的黄土地貌DEM数据，数据的区域大概有140 km2。

数据存放于…/ChP11/Ex1中，请将其拷贝到E：/ChP11/Ex1。

结果数据保存在…/ChP11/Ex1/Result 中。

5．算法思想：对于水文物理过程研究而言，由于山脊、山谷分别表示分水性与汇水性，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。

因此，对于山脊线和山谷线就可以利用水文分析的方法进行提取。

基于DEM的这种地形表面流水物理模拟分析的原理是：对于山脊线而言，由于它同时也是分水线，那么对于分水线上的那些栅格，由于分水线的性质是水流的起源点，通过地表径流模拟计算之后这些栅格的水流方向都应该只具有流出方向而不存在流入方向，也就是其栅格的汇流累积量为零。

通过对零值的汇流累积值的栅格的提取，就可以得到分水线，也就得到了山脊线；对于山谷线而言，由于其具有汇水的性质，那么对于山谷线的提取，可以利用反地形的特点，即是利用一个较大的数值减去原始的DEM数据，而得到了与原始地形完全相反的地形数据，也就是原始的DEM中的山脊变成负地形的山谷，而原始DEM中的山谷在负地形中就变成了山脊，那么，山谷线的提取就可以在负地形中利用提取山脊线的方法进行提取。

arcgis 山体范围提取【原创实用版】目录1.引言2.ArcGIS 简介3.山体范围提取方法4.山体范围提取的实际应用5.总结正文【引言】在地理信息系统 (GIS) 中，山体范围提取是一项非常重要的任务。

山体是地球表面上的一种基本地貌类型，对于地理信息系统中的空间分析、地形分析、水文分析等领域都有着重要的作用。

ArcGIS 是一款非常流行的 GIS 软件，可以用来进行山体范围提取等任务。

本文将对 ArcGIS 山体范围提取进行介绍，并对其实际应用进行探讨。

【ArcGIS 简介】ArcGIS 是由美国 Esri 公司开发的一款地理信息系统软件，具有强大的空间数据处理和分析能力，被广泛应用于地理信息系统、遥感、环境科学、城市规划等领域。

ArcGIS 包括桌面、服务器和移动端等多个平台，用户可以根据自己的需要选择不同的版本和模块。

其中，ArcGIS Desktop 是一款桌面端的软件，包括了 ArcMap、ArcCatalog、ArcToolbox 等工具，可以进行地理数据的编辑、管理、分析和可视化等操作。

【山体范围提取方法】在 ArcGIS 中，山体范围提取一般使用 DEM(数字高程模型) 数据作为输入数据。

DEM 数据是一种描述地表高程信息的数字数据，可以用来描述地形、山体等地貌类型。

下面是一种常用的山体范围提取方法:1.填充 DEM 数据，将 DEM 数据转换为填充栅格数据。

2.计算填充栅格数据的坡度和坡向。

3.根据坡度和坡向的阈值，将填充栅格数据进行分类，提取出山体范围。

4.将山体范围转换为矢量数据，以便进行后续的空间分析和可视化。

【山体范围提取的实际应用】山体范围提取在实际应用中有着广泛的应用，下面是一些例子:1.地形分析：山体范围提取可以为地形分析提供基础数据，例如计算地形坡度、坡向、地形起伏等。

2.水文分析：山体范围提取可以为水文分析提供基础数据，例如提取河流流域、水系网络等。

3.灾害评估：山体范围提取可以为灾害评估提供基础数据，例如山体滑坡、泥石流等灾害的评估。

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

本次实验通过某区域栅格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.791870117188）2）栅格计算器提取反地形DEM数据（输入公式1153.791870117188 - "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.制作立体图。

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

本次实验通过某区域栅格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•制作立体图。