ArcGIS

使用

Python

调用

arcpy

实现自动计算地形湿度

指数(

Topographic Wetness Index

，

TWI

)的方法

赵岩

1.

简介

TWI

是当地地形对径流流向和蓄积影响的物理指标或指标。

该指数是斜率和

上游贡献面积的函数。

TWI

的计算同时使用地理信息系统

(GIS)

和

Python (Python)

进行，

Python

是一种用于增强计算能力的编程软件。

这些指数有助于识别降雨径

流模式、潜在土壤含水量增加区域和积水区域。

地形湿润度指数

(TWI)

，

量化了地形对基本水文过程的控制

(Schillaci

等，

2015)

。

TWI

是通过细尺度地形与上梯度对地表面积的贡献相互作用，

根据以下关系得到

的

(Beven et al.

，

1979):

TWI = ln [CA/Slope]

其中，

CA

是通过网格单元排水的局部上坡集水区，

slope

是每个网格单元最

陡的向外坡度，以落差

/

距离衡量，即坡度角度的

tan

值

(Tarboton, 1997)

。

ArcGIS

通过使用工具箱可以实现地形湿度指数的计算，

但传统方法操作麻烦，

而通过

Python

调用

arcpy

可以实现自动自取，

这样我们就可以一边喝着咖啡，

一

边等待电脑帮我们完成任务了。

ArcPy

是一个

Python

站点包，可提供以实用高效的方式通过

Python

执行

地理数据分析、数据转换、数据管理和地图自动化。

2.

传统方法步骤

以

30m*30m

分辨率的图层为例。

(原文作者：

qrz110

)

。

①基于表面工具箱

Surface

计算

Slope

。