TWI-SPI计算公式

近期进行滑坡易发性评价,需要用到TWI和SPI参数,发现ArcGIS Pro无法直接计算,需要间接计算。

在网络上检索,发现YouTube上有视频可参考。

此为记录,方便后续查看使用。

ArcGIS计算公式

需要先计算坡度(Slope)和流量(Flow Accumulation),再使用栅格计算器进行TWI和SPI计算。

计算流程图

计算公式如下:

Topographic weakness index (TWI)

1
TWI=Ln(("FlowAcc"+1)*12.5/(Tan("SlopeDegree"*0.017453)+0.001))

Stream power index (SPI)

Type 1

1
SPI="FlowAcc"*12.5* Tan("SlopeDegree"*0.017453)

Type 2

1
SPI=Ln(("FlowAcc"*12.5+0.001)*(("SlopePercent"/100)+0.001))

Tips

  1. 12.5为网格大小,需要根据实际的栅格对应的尺寸进行设置。
  2. Stream power index (SPI) 共有三种计算公式,本文给的是其中一种。

原视频

References