matlab 山脊 提取,ArcGIS中利用水文分析提取山脊线山谷线

1 流程图

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

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

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

图1-1 流程图

2 操作步骤

2.1 正负地形求取

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

图2-1 焦点统计设置

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

图2-2 减法计算结果

(3)在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 正地形结果图

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

图2-6负地形结果图

2.2 山脊线的提取

(1)填洼：加载DEM数据，在ArcToolbox中选择Spacial Analyst Tools HydrologyFill工具，输入DEM进行填洼，设置如图2-7所示。

图2-7 填洼设置

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

图2-8 流向生成设置

图2-9 生成流向结果图

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

图2-10 生成汇流累积量设置

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

(4)汇流累积量为0值得提取：在ArcToolbox中选择Spacial Analyst Tools Map AlgebraRaster Calculator工具，输入公式如图2-12所示，结果如图2-13所示。

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

图2-13 汇流累积量为0值的结果

(5)利用邻域分析方法，对提取出来的0值结果进行3x3邻域分析，求均值，是数据光滑。设置如图2-14所示，结果如图2-15所示。

图2-14 邻域分析设置

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

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

图2-16 生成等值线设置

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

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

图2-18 生成山体阴影设置

图2-19 DEM晕渲图

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

图2-20 对均值后的数据进行分类

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

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

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

2.3 山谷线的提取

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

图2-22 反地形数据

(2)流向计算：在ArcToolbox中选择Spacial Analyst Tools HydrologyFlow Direction工具，输入反地形数据，得到流量数据，结果如图2-23所示。

图2-23 生成流向结果图

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

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

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

图2-25 汇流累积量为0值的结果

(5)利用邻域分析方法，对提取出来的0值结果进行3x3邻域分析，求均值，是数据光滑。结果如图2-26所示。

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

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

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

(7)生成山体阴影：在ArcToolbox中选择Spacial Analyst Tools Surfacehillshade工具，输入反DEM，结果如图2-28所示。

图2-28 反DEM晕渲图

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

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

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

图2-29 计算出研究区域的山谷线

2.4 制作专题图

打开Layerout view对话框，分别添加dem，山脊山谷线结果，插入图名图例比例尺指北针等地图制图要素，调整后制作专题图表，结果如图2-30,2-31所示。

图2-30 山谷线专题图

图2-31 山脊线专题图