坡面地形因子提取

坡面地形因子提取

数据准备：“规则格网DEM、散点DEM、TIN、等高线”文件夹中的“dem-grid”数据。

一、坡度√

操作按P32中的“1.坡度”进行。

二、坡向√

操作按P34中的“2.坡向”进行。

三、变率（坡度变率、坡向变率）

操作按P35中的“1.变率”进行。

三、四、曲率（综合曲率、平面曲率、剖面曲率）√

操作按P38中的“2.曲率”进行。

四、坡形√

1、打开arctoolbox，spatial analyst工具—邻域分析—焦点统计，弹出如下窗口：

单击确定，得到邻域均值计算结果图。

2、spatial analyst工具—地图代数—栅格计算器，弹出如下窗口：

输入公式"dem-grid" - "mean"，点确定，得到下图。

3、对该图进行分类，双击得到的坡向图打开其属性对话框，符号系统—已分类—分类，将其分为两类，＞0和＜0，＞0的是凸形坡，＜0的是凹形坡。

五、坡长

由于坡长是指在地面上沿水流方向逆流而上，到其水流起点间的最大地面距离在水平面上的投影长度。所以对于坡长提取的操作将放于水文分析的章节中进行。

六、坡位

A\沟间地

1、对坡度进行重分类。坡度大于25°的标记为0，小于25°的标记为1. Spatial Analyst 工具—重分类—重分类，弹出窗口，设置可参考下图：

单击确定，得到坡度重分类计算结果图。

2、进行邻域分析，分割正负地形。spatial analyst工具—邻域分析—焦点统计，弹出窗口，设置参考下图：

点击确定，得到邻域平均值计算结果图，如下：

3、栅格计算正负地形，原始DEM减去求平均值后的DEM。spatial analyst工具—地图代数—栅格计算器，弹出窗口，设置可参考下图：

点击确定，得到正负地形结果图，如下：

1、对正负地形进行重分类，大于0的为1，小于0 的为0。spatial analyst工具—重分类—

—重分类，设置参考下图：

点击确定，得到地形重分类结果图，如下：

5、进行栅格计算，判定候选沟间地。spatial analyst工具—地图代数—栅格计算器，弹出窗口，设置可参考下图：

当计算结果的栅格值为1时，说明该栅格满足坡度小于25°，且位于正地形，就判定其为候选沟间地，结果图如下：

6、在正负地形的分割中，会有把较宽沟谷谷底的某些区域判定为正地形，所以需要沟谷图层来将候选沟间地中的沟谷部分剔除。具体的操作步骤为：将“河流缓冲区”图层添加进来，然后spatial analyst工具—地图代数—栅格计算器，弹出窗口，设置可参考下图：

点击确定，计算结果如下：

7、对上步所得的结果图进行重分类。spatial analyst工具—重分类—重分类，弹出窗口，设置参考下图，直接将-1和0的新值改为0即可。

点击确定，得到最终沟间地分布图：

B\沟底地

1、对坡度进行重分类，大于20°的标记为0，小于20°的标记为1.

2、栅格计算。将沟谷图层乘以坡度分级图层。计算结果中，栅格值为1的满足坡度小于20，

且位于沟谷缓冲区内，即为沟底地。

C、沟坡地。

除去沟间地与沟底地以为的区域即为沟坡地。思考计算沟坡地的方法并通过软件操作实现。

七、地形起伏度

1、计算分析窗口内的最大值。打开arctoolbox，spatial analyst工具—邻域分析—焦点统

计，弹出窗口，设置参考下图：

单击确定，得到邻域最大值计算结果图。

分析窗口范

围是7*7的

统计类型选择的是计算分析窗口内的最大值（MAXIMUM）

2、计算分析窗口内的最小值。spatial analyst工具—邻域分析—焦点统计，弹出窗口，

设置参考下图：

分析窗口为

7*7

此处选择的是最小

值MINIMUM

3、计算地形起伏度。spatial analyst工具—地图代数—栅格计算器，弹出窗口，设置可

参考下图：

此处的公式根据自己前面

操作时的命名来进行点击确定，得到地形起伏度计算结果图。

八、地形粗糙度

根据书上p52的“2.地形粗糙度”进行操作。

九、地表切割深度

根据书上p52的“3.地表切割深度”进行操作。

十、高程变异系数

根据书上p53的“4.高程变异系数”进行操作。