山脊线山谷线提取实验报告材料

《山脊线山谷线提取实验报告材料》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案山脊线山谷线提取实验报告实验内容描述：山脊线和山谷线构成了地形起伏变化的分界线（骨架线），因此它对于地形地貌研究具有重要意义；另一方面，对于水文物理过程研究而言，由于山脊、山谷分别代表示分水性与汇水性，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。本次实验通过某区域栅格DEM掌握山脊线和山谷线这两个基本地形特征信息的理论及其基于DEM的提取方法与原理；同时，熟练掌握利用ArcGIS软件对这两个地形特征信息的提取方法。实验原理：1.本实验基于规则格网DEM数据使用平面曲率与坡形组合法提取山脊线和

2、山谷线，首先利用DEM数据提取地面的平面曲率及地面的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。实际应用中，由于平面曲率的提取比较繁琐，而坡向变率（SOA）在一定程度上可以很好地表征平面曲率。因此，提取过程中可以SOA代替平面曲率。2.主要用到以下理论知识：1）坡向变率：是指在提取坡向基础上，提取坡向的变化率，亦即坡向之坡度（SlopeofAspect，SOA）。它可以很好地反应等高线弯曲程度；2）反地形DEM数据：求取原始DEM数据层的最大高程值，记为H，通过公式（H-DEM），

3、得到与原来地形相反的DEM数据层，即反地形DEM数据；3）地面坡向变率SOA：地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原理，提取地表局部微小范围内坡向的最大变化情况。但是SOA在提取过程中在北面坡将会有误差产生，所以要将北坡坡向的坡向变率误差进行纠正，其公式为：SOA=(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/2其中：SOA1为原始DEM数据层坡向变率，SOA2为反地形DEM数据层坡向变率。4）焦点统计5）ArcScan自动矢量化流程图文档大全实用标准文案、文档大

4、全实用标准文案实验步骤：1.相对路径2.加载数据3.提取原始dem的坡向（利用dem数据--空间分析--表面分析--坡度工具，命名为Aspect）文档大全实用标准文案4.提取原始DEM数据的坡向变率（利用3中生成的Aspect图层--空间分析--表面分析--坡度工具，命名为SOA1）文档大全实用标准文案文档大全实用标准文案4.提取反地形DEM数据（栅格计算器--输入公式H-DEM）1）找出DEM最大高程值（右键属性---找出数据源中最大值为1153.791870117188）文档大全实用标准文案1）栅格计算器提

5、取反地形DEM数据（输入公式1153.791870117188-"dem"，命名为INdem）文档大全实用标准文案4.提取反地形DEM数据的坡向值文档大全实用标准文案4.计算反地形DEM数据的坡向变率文档大全实用标准文案8.计算进行误差纠正的地面坡向变率（栅格计算器--输入公式(("SOA1"+"SOA2")-Abs("SOA1"-"SOA2"))/2）文档大全实用标准文案8.邻域分析（原始dem--邻域分析--焦点统计focalstatistics（统计原始dem的平均值）---设置统计类型为平均值mean，

6、邻域类型为矩形(也可为圆形)，邻域大小为3*3（我发现邻域越大越模糊）(11*11)，则可得到一个邻域为3*3(11*11)的矩形的平均数据层，命名为mean文档大全实用标准文案文档大全实用标准文案文档大全实用标准文案10.计算正负地形分布区域（空间分析--地图代数--栅格计算器---输入公式为"dem"-"mean"，命名为Dvalue(差值)）文档大全实用标准文案10.利用栅格计算器提取山脊线（公式为"SOA">70&"Dvalue">0这是错的！！要加括号！！("SOA">70)&("Dvalue">0)

7、）和山谷线（("SOA">70)&("Dvalue"<0)）文档大全实用标准文案文档大全实用标准文案12.利用ArcScan自动矢量化得到山脊线山谷线的矢量图层1）在ArcCatalog中新建（方法有两种：右击文件夹--new--shapefile！或者是右击geodatabase--new--featureclass（新建要素类））山脊线图层（名称为shanjiline，类型为线）方法1：new--shapefile文档大全实用标准文案方法2：new--featureclass（但是这种方法下的线图层，在自动

8、矢量化山脊线后无法读到这个图层，所有还是选择方法1---这是因为栅格图层和矢量图层不能放在同一个geodatabase里面么？？？？？？？）文档大全实用标准文案文档大全实用标准文案1）打开开始编辑文档大全实用标准文案3）勾选扩展工具中的自动矢量化工具ArcScan文档大全实用标准文案4)在菜单栏空白处右击勾选ArcScan，打开ArcScan工具条--单击自动矢量化下的生成要素打开生成