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

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

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

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

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

图2-1 焦点统计设置

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

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

图2-5 负地形重分类设置

图2-7 填洼设置

流向计算：在ArcToolbox中选择Spacial Analyst Tools ? Hydrology ? Flow Direction工具，输入填洼后的DEM，得到流向数据，设置如图2-8所示。

图2-9 生成流向结果图

汇流累积量：在ArcToolbox中选择Spacial Analyst Tools ?Hydrology?Flow Accumulation工具，输入流向数据，得到流量数据，设置如图2-10所示，结果如图2-11所示。

图2-12 提取汇流累积量为0的值

图2-14 邻域分析设置

图2-15 0值部分求均值结果

生成等值线：在ArcToolbox中选择Spacial Analyst Tools ?Surface?Contour工具，输入原始DEM数据生成间距为40的等值线。设置如图2-16所示，结果如图2-17所示。

图2-17 生成等值线结果图

生成山体阴影：在ArcToolbox中选择Spacial Analyst Tools ?Surface?hillshade工具，输入原始DEM，其他设置按默认如图2-18所示，结果如图2-19所示。

图2-19 DEM晕渲图

在求均值后的0值汇流累积量数据上单击右键，选择属性，进行重新分级，将数据分为两级，调整临界点，以等值线图和晕渲图进行判断，属性值越接近1越可能是山脊线位置，确定分界阈值为0.5541.设置如图2-20所示。

图2-21 计算出研究区域的山脊线

04

山谷线的提取

在ArcToolbox中选择Spacial Analyst Tools ?Map Algebra?Raster Calculator工具计算反地形，公式为Abs(dem-2000)，得到反地形数据。如图2-22所示。

图2-23 生成流向结果图

汇流累积量：在ArcToolbox中选择Spacial Analyst Tools ?Hydrology?Flow Accumulation工具，输入流向数据，得到流量数据，设置如结果如图2-24所示。

图2-24 汇流累积量结果图

汇流累积量为0值得提取：在ArcToolbox中选择Spacial Analyst Tools ?Map Algebra?Raster Calculator工具，结果如图2-25所示。

图2-26 0值部分求均值结果

生成等值线：在ArcToolbox中选择Spacial Analyst Tools ?Surface?Contour工具，输入反DEM数据生成间距为40的等值线。设置如图结果如图2-27所示。

图2-28 反DEM晕渲图

在求均值后的0值汇流累积量数据上单击右键，选择属性，进行重新分级，将数据分为两级，调整临界点，以等值线图和晕渲图进行判断，属性值越接近1越可能是山脊线位置，确定分界阈值为0.65667。

将属性直接近1的的一类属性值设为1，其余的设为0。

在ArcToolbox中选择Spacial Analyst Tools ?Map Algebra?Raster Calculator工具，将重分类后的数据与负地形数据相乘，可以消除存在于正负地形区域的错误山脊线。然后将计算结果重分类，属性部位1的栅格设为NODATA，就得到了山谷线，结果如图2-29所示，浅色区域为山谷线。

图2-30 山谷线专题图

图2-31 山脊线专题图

通知：地信 地质专项群已建立，如需进群，请添加WX: dixinwang520 备注 地质返回搜狐，查看更多

相关资源：segtool.tar.gz_matlab时频脊线_小波脊线提取_脊提取_脊线_脊线…