ArcGIS:制作滑坡单元

一、流程图

参考论文：基于GIS的斜坡单元划分方法改进与实现-颜阁

二、流程

1.准备数据

step1.下载裁剪DEM（裁剪看往期）

step2.低通滤波过滤。Spatial Analyst工具-邻域分析-滤波器-LOW

滤波原理参考：滤波器 (Spatial Analyst)—ArcMap | 文档 (arcgis.com)

2.计算曲率

曲率 (Spatial Analyst)—ArcMap | 文档 (arcgis.com)

Spatial Analyst工具-表面分析-曲率工具

有需要的可以输出剖面曲率和平面曲率。（直接计算的曲率是平均曲率）

曲率结果

3.用曲率代替高程计算流域

step1.Spatial Analyst工具-水文分析-填洼

填洼结果

step2.Spatial Analyst工具-水文分析-流向

计算结果

step3. 流量

Step4.栅格计算器

大于5000的赋值为1（我不知道怎么确定最小径流）

细节

step5河流链接

链接后值有变化

Step6.栅格矢量化

矢量化结果

step7.流域提取

获取流域范围

step8.栅格转面

结合高程一起看，现在转为矢量面的是凹形地貌边界，

4.反转曲率重复计算一次流域

用栅格计算器反转曲率，然后再把前面步骤做一遍

(("曲率" - 曲率最大值) * -1) + 曲率最小值

流域结果如图，以凹地作为矢量边界

5.合并两次流域矢量

分的效果还行，红色是曲率计算结果，蓝色是反向曲率计算结果

数据管理工具-常规-合并

合并结果就是滑坡单元