式中:L为坡长因子;𝜆为坡长,通过汇流计算得到;m为坡长指数,是无量纲常数。
坡长𝜆计算:“流量”×像元大小;
坡长指数m计算:坡长指数m计算方式有几种,其中刘宝元老师的m取值方法使用较多(符素华 刘宝元《坡度坡长因子计算工具》),推荐采用,如下:
m=0.2,θ<0.5°;m=0.3,0.5°≤θ<1.5°;m=0.4,1.5°≤θ<3°;m=0.5,θ≥3°
arcgis具体操作:
1、加载DEM数据——系统工具箱——Spatial Analyst Tools——水文分析——填洼,输出“填洼”;
2、水文分析——流向,输入“填洼”,输出“流向”;
3、水文分析——流量,输入“流向”,输出“流量”;
4、双击流量图层(或右击图层名称,打开属性),在源中可找到像元大小数值
5、计算坡长𝜆:打开Spatial Analyst Tools——地图代数——栅格计算器,输入“流量”*0.00027777;得到坡长
6、计算坡长指数m:栅格计算器——输入Con("坡度" < 0.5,0.2,Con("坡度" < 1.5,0.3,Con("坡度" < 3,0.4,0.5)))
7、计算坡长因子L:栅格计算器——输入Power("λ" / 22.13,"m ")。