ArcGIS中坡度坡长等地形因子分析过程

地形因子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))