地形湿度指数

1. 对DEM 填挖，Fill

2. 对填挖后的DEM 求水流方向，Flow Direction

3. 对第二步生成的水流方向求汇流累积量，Flow Accumulation

4. 用以下公式使用栅格计算器计算地形湿度指数

2. 地形湿度指数=lnAs/tanβ

3. As 是垂直于水流方向的特定汇流面积，流经地表某点的单位等高线长度上的汇流面积

地形湿度指数 = Ln(([汇流累积量] * 900 / con([水流方向] == 1 , 30 , [水流方向] == 4 , 30 , [水流方向] == 16 , 30 , [水流方向] == 64 , 30 , [水流方向] == 2 , 30 * Sqrt(2) , [水流方向] == 8 , 30 * Sqrt(2) , [水流方向] == 32 , 30 * Sqrt(2) ,

[水流方向] == 128 , 30 * Sqrt(2)))/ Tan(con([坡度] 0 , [坡度]) * 3.1415926 / 180))

4. 最后的那个公式我觉得有点问题

5. 我觉得正确的应该是这样的

6. 地形湿度指数 = Ln((con([汇流累积量] == 0 , 1 , [汇流累积量]) * 900 / con([水

流方向] == 1 , 30 , [水流方向] == 4 , 30 , [水流方向] == 16 , 30 , [水流方向] == 64 , 30 , [水流方向] == 2 , 30 * Sqrt(2) , [水流方向] == 8 , 30 * Sqrt(2) ,

[水流方向] == 32 , 30 * Sqrt(2) , [水流方向] == 128 , 30 * Sqrt(2)))/ Tan(con([坡度] 0 , [坡度]) * 3.1415926 / 180))

7. 区别在于con([汇流累积量] == 0 , 1 , [汇流累积量]) 这里的意思是判断汇流累积

量是否为0，如果为0，则令其为1，因为你这里 [汇流累积量]*900其实就是算栅格汇流面积，对于汇流累积量为0的栅格其汇流面积应该至少为其本身的大小即栅格单元大小，在你这里就是1*900=900

地形湿度指数

Ln((Con("demacc" == 0,1,"demacc") * 900 / Con(("demdirection" ==1) |

("demdirection" == 4) | ("demdirection" == 16) | ("demdirection" ==

64),30,42.4264))/Tan(Con("slope"

河流动能指数

(Con("demacc" == 0,1,"demacc") * 900 / Con(("demdirection" ==1) | ("demdirection" == 4) | ("demdirection" == 16) | ("demdirection" ==

64),30,42.4264))*Tan(Con("slope"

坡度增长因子

Power(Con("demacc" == 0,1,"demacc") * 900 / Con(("demdirection" ==1) | ("demdirection" == 4) | ("demdirection" == 16) | ("demdirection" ==

64),30,42.4264)/22.13,0.6)*Power(Sin(Con("slope"

3.1415926 / 180)/0.0896 ,1.3)