基于ARCGIS计算RUSLE中LS因子值

LS因子针对不同的研究区,计算方法不同,本文参考别人的文献,选择了计算方法:

基于ARCGIS计算RUSLE中LS因子值_第1张图片
另外还有很多计算方法

1、首先对DEM数据进行填洼处理:arcgis中fill工具。
之后计算坡度:坡度输出单位直接设置为degree就行。但在接下来的三角函数运算中,度要转换为弧度。

2、使用栅格计算器,计算S因子:
Con(“slope.tif” <= 5 ,10.8 * Sin(“slope.tif” * 0.01745) + 0.03 , Con(“slope.tif” <=10 ,16.8 * Sin(“slope.tif” * 0.01745) - 0.50 , Con(“slope.tif” <= 25,20.204 * Sin(“slope.tif” *0.01745 ) -1.2404,29.585 * Sin(“slope.tif” *0.01745 ) -5.6079 ) ) )

​这里乘以0.01745的目的是将角度转换为弧度,进而计算sin值
基于ARCGIS计算RUSLE中LS因子值_第2张图片

3、进行流向分析:
在流向分析和流量分析时,使用arcmap速度很慢,尤其是数据量很大的时候,推荐使用正版的arcgis pro。申请一个试用版也许。
流向分析也是只设置出入和输出就行。
4、流量分析,也是只设置输入输出就好。
5、Lambda的计算:lanbda使用流量栅格分辨率,我这里是使用的带投影的。分辨率大约30m。
6、计算L
我这里把5和6合并在一起计算了。
Con(“slope.tif” <=1 ,Power(“flowacc.tif” * 29.4481461287987 /22.13 ,0.2) , Con(“flowacc.tif” <=3 , Power(“flowacc.tif”
29.4481461287987 /22.13 ,0.3) , Con(“slope.tif” <=5 , Power(“flowacc.tif”* 29.4481461287987 /22.13 ,0.4), Power(“flowacc.tif”* 29.4481461287987 /22.13 ,0.5)) ) )


也可以先使用con函数计算M,再使用power(flowacc,m)进行计算

7、计算LS=L*S
对结果重新设置样式

基于ARCGIS计算RUSLE中LS因子值_第3张图片

在运算过程中,尤其是数据量很大的情况下,尽量放在根目录下面,我不知道为啥目录很深,也会运算不出来结果。
对于结果,要根据不同人的研究,去进行进一步的验证与纠正。

你可能感兴趣的:(学习笔记,学习)