山脊线山谷线提取实验报告.docxVIP

山脊线山谷线提取实验报告 实验内容描绘： 山脊线和山谷线组成了地形起伏变化的分界限（骨架线），所以它关于地形地貌研究拥有 重要意义；另一方面，关于水文物理过程研究而言，因为山脊、山谷分别代表示分水性与汇 水性，山脊线和山谷线的提取实质上也是分水线与汇水线的提取。 本次实验经过某地区栅格DEM掌握山脊线和山谷线这两个基当地形特点信息的理论 其鉴于DEM的提取方法与原理；同时，娴熟掌握利用ArcGIS软件对这两个地形特点信息的 提取方法。  及 实验原理： 1.本实验鉴于规则格网DEM数据使用平面曲率与坡形组合法提取山脊线和山谷线，第一利 用DEM数据提取地面的平面曲率及地面的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。实质应用中，因为平面曲率的提取比较繁琐，而坡向变率 （SOA）在必定程度上能够很好地表征平面曲率。所以，提取过程中能够SOA取代平面曲率。 2.主要用到以下理论知识： 1）坡向变率：是指在提取坡向基础上，提取坡向的变化率，亦即坡向之坡度（SlopeofAspect， SOA）。它能够很好地反响等高线曲折程度； 2）反地形DEM数据：求取原始DEM数据层的最大高程值，记为H，经过公式（H-DEM）， 获得与本来地形相反的DEM数据层，即反地形DEM数据； 3）地面坡向变率SOA：地面坡向变率在所提取的地表坡向矩阵的基础上沿袭坡度的求算原 理，提取地表局部细小范围内坡向的最大变化状况。可是SOA在提取过程中在北面坡将会 有偏差产生，所以要将北坡坡向的坡向变率偏差进行纠正，其公式为： SOA=(([SOA1]+[SOA2])-Abs([SOA1]-[SOA2]))/2 此中：SOA1为原始DEM数据层坡向变率， 4）焦点统计 5）ArcScan自动矢量化  SOA2为反地形  DEM数据层坡向变率。 流程图 、 实验步骤： 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--邻域剖析--焦点统计focalstatistics（统计原始dem的均匀值）---设 置统计种类为均匀值mean，邻域种类为矩形(也可为圆形)，邻域大小为3*3（我发现邻域越 大越模糊）(11*11)，则可获得一个邻域为3*3(11*11)的矩形的均匀数据层，命名为mean 10.计算正负地形散布地区（空间剖析--地图代数--栅格计算器---输入公式为dem-mean，命名为Dvalue(差值)） 11.利用栅格计算器提取山脊线（公式为(SOA70)(Dvalue0)）和山谷线（  SOA70Dvalue0这是错的！！要加括号！！ (SOA70)(Dvalue0)） 12.利用ArcScan自动矢量化获得山脊线山谷线的矢量图层 1）在ArcCatalog中新建（方法有两种：右击文件夹--new--shapefile！或许是右击 geodatabase--new--featureclass（新建因素类））山脊线图层（名称为shanjiline，种类为 线） 方法  1：new--shapefile 方法2：new--featureclass（可是这类方法下的线图层，在自动矢量化山脊线后没法读到这个图层，全部仍是选择方法1---这是因为栅格图层和矢量图层不可以放在同一个geodatabase里面么？？？？？？？） 2）翻开开始编写 3）勾选扩展工具中的自动矢量化工具ArcScan 4)在菜单栏空白处右击勾选ArcScan，翻开ArcScan工具条--单击自动矢量化下的生成因素翻开生成因素对话框即可生成自动矢量化后的矢量山脊线 5)用相同的方法生成矢量山谷线 13.制作立体图。用三维剖析--表面剖析--Hillshade工具鉴于原始DEM数据--做出山体暗影图层，调理各图层颜色和线宽，并用effect工具调理山体暗影图层、山脊线、山谷线图层透明度，获得最后结果。 14.制作专题图 15.导出为灰度图