基于matlab山脊线,教你如何利用水文，分析提取山脊线山谷线

原标题：教你如何利用水文，分析提取山脊线山谷线

利用水文分析提取山脊线及山谷线，山脊线相当于分水线，山谷线相当于山谷线。分水线是水流的起源点，这些栅格的水流方向只存在流出方向而不存在流入方向，所以汇流累积量为零。通过对零值的提取就可以得到山脊线。

山谷线相当于汇水线要用反地形求出，即用较大值减去DEM，DEM中山脊线就成为山谷线，山谷线变为山脊线，用求山脊线的方法求出山谷线，分别利用正反地形求交验证。

DEM进行填洼，利用水文分析求出流向流量，再提取出汇流累积量为零得值与正地形求交，即得到分水线也就是山脊线。用反地形求流向流量提取汇流累积量为零的部分与负地形求交就是山谷线。

02

操作步骤

正负地形求取

加载DEM数据，在ArcToolbox中选择Spacial Analyst Tools  NeighborhoodFocus Statistics工具，输入dem，利用11*11窗口计算平均值。设置如图2-1所示。

图2-1 焦点统计设置

在ArcToolbox中选择Spacial Analyst Tools Map Algebra Raster Calculator工具，对原始数据与焦点统计后的DEM做减法。结果如图2-2所示。

图2-2 减法计算结果

在ArcToolbox中选择Spacial Analyst Tools ReclassReclassify工具，对减法运算结果进行重分类，分级界线为0。将大于0的区域赋值为1，小于0的区域赋值为0即得到正地形；设置如图2-3所示，结果如图2-4所示。将大于0的区域赋值为0，小于0的区域赋值为1即得到负地形，设置如图2-5所示，结果如图2-6所示。

图2-3 正地形重分类设置

图2-4