水文分析和地形分析、提取山谷和鞍部点全解.doc

Arcgis地形分析 DEM提取坡度 新建地图文档，加载【ArcGIS地形分析--TIN及DEM的生成，TIN的显示】经验教程中得到的DEM数据 (2) 在【ArcToolbox】中，执行命令[3D Analyst工具]——[栅格表面]——[坡度]， 参照下图所示，指定各参数： 执行后，得到坡度栅格Slope_tingri1：坡度栅格中，栅格单元的值在[0 -] 度间变化 右键点击图层[Slope_TinGrid]，执行[属性命令]，设置图层[符号系统]，重新调整坡度分级。将类别调整为5，点[分类]按钮，用手动分级法，将中断值调整为：，，，，。 1.添加dem数据，制作15m和75m等高线。 2.制作阴影阴影图像【空间分析】----【表面分析】----【地表阴影】，生成地表阴影图像hillsha。 3.提取栅格数据的有效区域。【空间分析】----【地图制图】----【栅格计算器】。 “要提取的文件名”=“huainan54”=0，（注意：红色等号是1个=，而不是栅格计算器中的2个==）。“back”=“huainan54”=0，生成back文件。 5.按照等高线75m等高线15m，back，hillsha叠放。 提取dem数据中的最大值。【空间分析工具】----【邻域分析】----【块统计】，生成maxpoint文件 参数设置如下： Maxpoint 7.提取山顶点。【空间分析】----【地图制图】----【栅格计算器】。 输入命令：sd=（[axpoint]-[淮南54]）==0,生成sd文件。 8.山顶点栅格文件二值化。【空间分析工具】----【重分类】----【重分类】。生成re-sd 9.生成山顶点栅格数据转换为shapefile。【转换工具】----【栅格转换】----【栅格数据转点】，生成山顶点矢量文件。 ArcGIS利用水文分析方法提取山脊、山谷线 提取方法大致可以分为以下五种：1) 基于图像处理技术的原理；2) 基于地形表面几何形态分析的原理；3) 基于地形表面流水物理模拟分析原理；4) 基于地形表面几何形态分析和流水物理模拟分析相结合的原理；5) 平面曲率与坡形组合法。 平面曲率与坡形组合法提取的山脊、山谷的宽度可由选取平面曲率的大小来调节，方法简便，效果好。该方法基本处理过程为：首先利用 DEM数据提取地面的平面曲率及地面的正负地形，取正地形上平面曲率的大值即为山脊，负地形上平面曲率的大值为山谷。实际应用中，由于平面曲率的提取比较繁琐，而坡向变率（SOA）在一定程度上可以很好地表征平面曲率。因此，下面的提取过程以 SOA代替平面曲率。 1. Spatial Analysis 下使用surface 菜单下的Aspect，提取 DEM 坡向aspect-huainan aspect-huainan 2. 在 Spatial Analysis 下使用 surface 菜单下的 Slope 命令，提取A 层面的坡向变率，记为 SOA1。 SOA1 3. Spatial Analysis 下的map algebra，栅格计算器Calculator，得到与原来地形相反的 DEM 数据层，即反地形 DEM 数据，记为 IN。公式为：IN= IN=最大值-huinan_54 IN 4. Aspect-IN 5. 利用 SOA 方法求算反地形的坡向变率，记为 SOA2。 soa2 6. 在Spatial Analysis 下map algebra使用栅格计算器 Calculator，求出没有误差的 DEM 的坡向变率SOA。公式为SOA ＝ （（[SOA1]+[SOA2]）-Abs（[SOA1]-[SOA2]））/ 2。 Soa 7. 在 Spatial Analysis 下使用栅格邻域计算工具 Neighborhood，focal Statistics，输入原始dem数据huainan54，设置 Statistic type 为平均值，生成栅格图像为mean。 Mean 8. 在 Spatial Analysis 下map algebra使用栅格计算器 Calculator，求出正负地形分布区域dvalue。公式为 davlue＝[huainan54]-[mean]。 Dvalue 9. Spatial Analysis下map algebra使用栅格计算器 Calculator，求出山脊线shanjixian。公式为 shanjixian＝（[dvalue]0） （SOA 0）。 10. Spatial Analysi