地形湿度指数(word文档良心出品).docVIP

对DEM填挖，Fill

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

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

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

地形湿度指数=lnAs/tanβ

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,0.0000000001,[坡度]0,[坡度])*3.1415926/180))

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

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

地形湿度指数=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,0.0000000001,[坡度]0,[坡度])*3.1415926/180))

区别在于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=0,0.000001,slope)*3.1415926/180))

河流动能指数

(Con(demacc==0,1,demacc)*900/Con((demdirection==1)|(demdirection==4)|(demdirection==16)|(demdirection==64),30,42.4264))*Tan(Con(slope=0,0.0001,slope)*3.1415926/180)

坡度增长因子

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=0,0.0001,slope)*3.1415926/180)/0.0896,1.3)