地形因子LS制图过程
DEM（1:5万）
水文分析中进行填充（fill）
填充后生产河流流向图
填充后生成坡度图
90m格网DEM
非累计流量分析最大坡降计算
S=10.8sinθ+0.036 θ<5°S=16.8sinθ-0.5 5°≤θ<10° S=21.9sinθ-0.96 θ≥10° m=0.2θ<1°
m=0.3 1°≤θ<3° m=0.4 3°≤θ<5° m=0.5θ≥5°
生成坡度图生成坡长图
L=(λ/22.1)m
坡度因子S图坡长因子L图
根据上图及公示生成坡度因子图
con( [Slope] < 5,10.8 * Sin([Slope]) + 0.036,con([Slope] < 10,16.8 * Sin([Slope]) - 0.5,21.9 * Sin([Slope]) - 0.96))
计算坡度因子过程图
con( [Slope] < 5,10.8 * Sin([Slope] * 3.1415926 / 180) +
0.036,con([Slope] >= 10,21.9 * Sin([Slope] * 3.1415926 / 180) -
0.96,16.8 * Sin([Slope] * 3.1415926 / 180) - 0.5))
求得的S因子图
con([Slope] < 1,0.2,con([Slope] < 3,0.3,con([Slope] < 5,0.4,0.5)))
求得的M值图
求得流向图，为了算是否为根二
con([flowdir] == 2 | [flowdir] == 8 | [flowdir] == 32 | [flowdir] ==
128,1.414,1)
Pow(90 * [xiu] / (Cos([Slope] * 3.1415926 / 180)* 22.1),[mfac])
con( [slope] < 5,10.8 * Sin([slope] * 3.1415926 / 180) + 0.036,con([slope] >= 10,21.9 * Sin([slope] * 3.1415926 / 180) -
0.96,16.8 * Sin([slope]* 3.1415926 / 180) - 0.5))
