坡面地形因子提取
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
坡面地形因子提取
数据准备:“规则格网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.高程变异系数”进行操作。