ArcGIS使用Python调用arcpy实现自动计算地形湿度指数(TWI)的方法.docx

ArcGIS 使用 Python 调用 arcpy 实现自动计算地形湿度指数（Topographic Wetness Index，TWI）的方法 赵岩 简介 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 执行地理数据分析、数据转换、数据管理和地图自动化。 传统方法步骤 以 30m*30m 分辨率的图层为例。（原文作者：qrz110）。 ①基于表面工具箱 Surface 计算 Slope。 ②单位转换：Scale_slope=Slope*pi/180。 ③填洼(Fill_dem)-Fill。 ④水流方向(FlowDir_Fill)-Flow Direation。 ⑤汇流累积量(FlowAcc_Flow)-Flow Accumulation。 ⑥单位面积的汇流量(SCA)-Raster Calculator。 ⑦地形湿度指数（TWI）-Raster Calculator：TWI=ln(SCA/Tan(Scale_slope))。 使用 Python 实现自动化计算 ①购买咖啡，烧开水。 ②打开 ArcGIS，打开 Python 窗口。 ③输入以下代码： ... TWI-filled-plus0.1constant.py ... Topographic Wetness Index ... 2013-11-12 ... Jeffrey Wolf (EEB); G. Andrew Fricker (GEOG) ... UCLA ... This script was written to be used as a tool in ArcGIS. This python script can be imported to create a TWI tool. Inputs are the workspace and the input DEM, output is the TWI layer ... We chose to fill all sinks due to some small sinks in a lidar redived DEM. We also add a small constant to the denominator to avoid dividng by zero. ... The original script was based off the arcpy script written by Prasad Pathak. ... /details.asp?dbid=16750 ... This revised script converts the terrain slope in degrees to radians ... This script also uses the default settings for the flow accumulation raster ... however different methods to calculate flow accumulation can dramtically change the results of the TWI ... import arcpy, math if name == main : ... arcpy.CheckOutExtension(Spatial) ... ... # Define workspace and set input and output files ... arcpy.env.workspace = D:/test.gdb ... inDEM = E:/DATA