特征点、线、面提取
坡度、坡向、等高线的生成
加载30m 沿河流域栅格数据,设置色带,越红表示海拔越高,越蓝表示海拔越低
计算坡向
计算坡度(与计算坡向类似不再叙述)
计算等高线
- 使用等值线进行计算
- 设置图层属性,标注字段为contour,显示等高距
变率和曲率的计算
坡度变率
- 先对图层求取坡度
- 再对计算后的坡度图层求取坡度,即坡度的坡度(坡度变率)
坡向变率
- 先对图层求取坡向
- 再对计算后的坡向图层求取坡度,即坡向的坡度(坡向变率)
曲率计算
- 搜索曲率或工具箱-空间分析-表面分析-曲率
- 生成所需要的曲率数据
地形起伏度和地形粗糙度计算
地形起伏度
曲率的最大值与最小值的差值(高程的差值),使用空间分析-邻域分析-焦点统计工具,计算最大、最小值
- 计算最大值
- 计算最小值
- 计算起伏度(最大减最小),空间分析-地图代数-栅格计算器
地形粗糙度
地形粗糙度与坡度有关需要先计算坡度
- 计算坡度
- 使用空间分析-地图代数-栅格计算器,计算粗糙度(
1 / Cos(坡度 * 3.14 / 180)
)直接转换为弧度
- 完成计算
地表切割深度和高程变异系数计算
地表切割深度
地表切割深度:平均高程与最低高程的差值
- 计算平均值
- 计算切割深度
高程变异系数
高程变异系数是标准差与平均值的比值,反映地表宏观区域高程的相对变化
- 计算标准差
- 使用空间分析-地图代数-栅格计算器,计算变异系数
山顶点的提取
正地形求取山顶点
山顶点:区域内高程的极大值;
计算步骤:
- 取出邻域内高程最大值
- 使用计算后的最大值减去原始值
- 差值为零,即为山顶点
操作过程
- 计算最大值
- 计算差值
差值不仅有0还有大于0的部分,为0的部分为山顶点,不为0的影响分析,要去除影响 - 设置相减后的数据为0,分类差值数据
放大栅格图层会发现有,差值为0的数据有时不仅仅为一个栅格,而且还形成了面(现实中山顶点为一块平地),为了让形成面的连续栅格只显示一个栅格,需要对数据进行特殊处理
- 对计算后的差值数据进行重分类,只显示差值为0的山顶点数据
- 设置只显示原始栅格图层与差值为0的图层,黑色部分为差值为0栅格数据
- 对差值为0的栅格图层转为矢量面图层(搜索raster to polygon或栅格转面),取消简化面
- 再将矢量的面图层转化为点(feature to point 要素转点 ),在连续矢量图层内部生成点数据
- 有时会出现在相同的山顶点出现多个中心点的情况,这时需要手动处理,或者根据栅格数据的实际情况,再计算最大值时确定合适的分析区域大小(分析窗口)分析窗口过大会漏掉山顶点,分析窗口过小会出现过多山顶点,避免同一个山顶出现多个山顶点的情况发生
- 导出山顶点矢量数据
反地形求取山顶点
- 求取栅格的反地形,原始栅格乘-1,高低对调,原始高的变为低的
- 计算流向
- 计算汇(洼地),即山顶点所在的区域的栅格图层
显示原始栅格与计算后的汇图层
- 对山顶点图层进行,栅格转面(raster to polygon)-要素转点(feauter to point),计算出山顶点(含冗余点)
- 去除冗余点,计算原始栅格等高线
- 使用要素转面(feature to polygon),将等高线转为面图层,若能转化,然后通过面图层与山顶点作相交运算,相交的即为真正的山顶点,不能转化面的部分去除。
线路不闭合是不能转化为面的,所以可以过滤。
- 按位置选择
- 删除冗余山顶点。启动编辑器,打开属性表,此时选择的是相交山顶点,“切换选择”进行反选,使用delete删除选中的冗余点,保存编辑,停止编辑即可。