基于水文分析方法提取山脊线和山谷线
1.提取思路:
基于水文分析方法提取山脊线和山谷线实际上是对分水线和汇水线的提取。
对于山脊线来说,它同时也是分水线,即水流的起始点。因此通过地表径流计算后,这些栅格的水流方向应该都只具有流出方向而不存在流入方向,也就是说这些栅格的汇流累积量为0,因此通过对0值的提取也就得到了山脊线。
对山谷线的提取,可以对反地形DEM数据提取山脊线,得到的就是实际上的山谷线。
2.基础操作步骤介绍:
2.1 正负地形的提取:
选择【系统工具箱→Spatial Analyst Tools→邻域分析→焦点统计】工具,得到邻域大小为11×11(可以根据需要自行设置)的矩形的平均值数据层Mean_DEM。
选择【系统工具箱→Spatial Analyst Tools→地图代数→栅格计算器】工具,输入"DEM" - "Mean_DEM"地图代数公式,得到正负地形分布区域数据层。
选择【系统工具箱→Spatial Analyst Tools→重分类→重分类】工具,将大于0的区域赋值为1,小于0的区域赋值为0,得到正地形数据层
zhengdixing。
同样,选择【系统工具箱→Spatial Analyst Tools→重分类→重分类】工具,将大于0的区域赋值为0,小于0的区域赋值为1,得到负地形数
据层fudixing。
2.2 山脊线的提取:
选择【系统工具箱→Spatial Analyst Tools→水文分析→填洼】工
tool工具箱具,对原始DEM进行洼地填充,得到无洼地DEM数据层Fill_DEM。
具,对无洼地DEM数据进行水流方向计算,得到流向数据层FlowDir。
具,计算汇流累积量,得到数据层FlowAcc。
选择【系统工具箱→Spatial Analyst Tools→地图代数→栅格计算器】工具,输入地图表达式"FlowAcc" == 0,提取汇流累积量为0的栅格区
域,得到数据层Facc0。
查看上一步得到的数据层Facc0发现其中会有一些稀疏的点也就是说很多地方不是山脊线,因此需要对数据做平滑处理。选择【系统工具箱→Spatial Analyst Tools→邻域分析→焦点统计】工具,得到邻域大小为
3×3的矩形的平均值数据层Facc0_Mean。
具,生成原始DEM数据的等值线数据层Contour。

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。