基于DEM的地形单元多样性指数及其算法.pdfVIP

第32 卷 第1期 地 理 科 学 进 展 Vol.32, No.1 2013 年01 月 PROGRESSINGEOGRAPHY Jan., 2013 基于DEM 的地形单元多样性指数及其算法 1,2 1 1 3 田瑞云 ，王玉宽 ，傅 斌 ，刘 援 (1. 中国科学院水利部成都山地灾害与环境研究所，成都610041 ；2. 中国科学院大学，北京100049； 3. 环境保护部环境保护对外合作中心，北京100035) 摘 要：在阐述地形信息表达研究进展的基础上，提出基于DEM地形单元多样性指数的概念和算法。地形单元多 样性指数综合了高程、坡度、坡位、坡向、汇流量和水域信息等要素。其算法集成地形位置指数和地形湿度指数算 法，采用图层叠加分析，设定分类、分级指标，进行重分类组合，划分地形单元类型，利用窗口分析法计算地形单元 多样性指数。以四川省为试验区，利用精度为100m的DEM数据和水域分布数据进行模拟计算，地形单元划分为 13种典型类型，统计窗口半径设为900m，计算出的多样性指数值小于0.5的区域仅占总面积的11%，大于0.75的 区域达57%，符合实验区地形特征，并对算法进行了可行性分析和验证。结果表明，该算法提取的地形单元多样 性指数可以有效反映地表形态的多样性特征及其变化。该研究结果为进一步探讨基于DEM地表形态信息的概念 体系，以及从微观到宏观的地形信息空间分析研究创造了条件。 关 键 词：地形单元；DEM ；GIS ；地形位置指数；地形湿度指数 doi: 10.3724/SP.J. 1033.2013.00121 1 引言 地形要素特征的量化指标。 Weiss(2001)提出的地形位置指数(Topographic 地形要素是最基本的自然地形单元，影响着地 Position Index ，TPI) 可以区分出地形的坡位(如脊 表物质迁移与能量转换、地表过程以及景观格局 顶、谷底、半山坡) ，并识别出地形类别(如峡谷、平 等，其各种因子的空间分布特征是用来描述多种空 地、斜坡) ，但只综合了相对高程和坡度信息。TOP- 间变化过程的重要指标(邬伦等, 2001) 。沈泽昊等 MODEL 模型使用地形指数(Terrain Index)、亦称地 (2001)在进行生态位研究时，将地形分解为海拔、坡 形湿度指数(Topographic Wetness Index ，TWI)可以 向、坡位和坡度4 个属性指标，作为多维生态位空 有效反映地形对流域产流源面积的影响(Beven et 间的环境梯度轴。许多学者对这4 个地形因子的 al, 1979; 邓慧平等, 2002; 孔凡哲等, 2003; 任立良 空间特征及环境效应进行了研究(陈世品等, 2007; 等, 2002) ，虽然综合了汇流和坡度特征，但不能直 Ferrer-Castan et al, 2005; 李容全等, 1990; 刘新华 接反映地表的形态变化。剖面曲率和平面曲率可 等, 2001; Nevo et al, 1999; 沈泽吴等, 2001; 严岳鸿 有效反映地面的坡形与变化特征(李志林等, 2001; 等, 2011; 赵瑜等, 2010) ，并取得一定进展。例如， 闾国年等, 1998; 汤国安等, 2003; 张婷等, 2005) ，但 基于海拔提取地形起伏度可以作为水土流失评价 由于这两个指标分别在垂直与水平上描述坡形变 的地形指标（李容全等, 1990）；坡向的变化对植物 化，不能综合反映地形的空间变化特征。王雷等 的生长和分布具有显著影响（严岳鸿等, 2011 ）；坡 (2006)提