基于ARCGIS计算RUSLE中LS因子值

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

另外还有很多计算方法

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值

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
对结果重新设置样式
​

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