地形湿度指数算法误差的定量评价.pdfVIP

第30 卷 第1期 地 理 科 学 进 展 Vol.30, No.1 2011 年01 月 PROGRESS IN GEOGRAPHY Jan., 2011 地形湿度指数算法误差的定量评价 1,2 1 1 包黎莉 ，秦承志 ，朱阿兴 (1. 中国科学院地理科学与资源研究所资源与环境信息系统国家重点实验室，北京 100101； 2. 中国科学院研究生院，北京 100049) 摘 要：地形湿度指数(TWI)能够定量指示地形对土壤湿度空间分布的控制，是一种应用广泛的地形属性。目前基 于栅格DEM 的TWI 计算方法结果各异，因此有必要对TWI 算法进行定量评价。对TWI 算法通常是应用实际 DEM 数据进行评价，但实际DEM 中存在的数据源误差会干扰对算法误差的评价。针对该问题，本文介绍了一种 用不含数据源误差的人造DEM 定量评价TWI 算法误差的方法。本文选择4 种代表性TWI 算法：基于经典单流向 算法的TWI 算法、基于经典多流向算法的TWI 算法、基于局域地形自适应多流向算法的TWI 算法，以及自适应多 流向算法与最大下坡相结合的TWI 算法，应用该评价方法构建了一组模拟4 种典型复合地形条件的人造表面，并 分别离散成具有不同空间分辨率(1 m 、10 m 和30 m) 的人造DEM 数据，对TWI 算法进行定量评价。评价结果显 示：在中心凸、凹坡和马鞍面地形条件下，基于多流向算法计算的TWI 优于基于单流向算法的TWI 结果；在山脊地 形条件下，基于单流向算法的结果仅次于MFD-md 结合最大下坡度算法计算的结果。随着分辨率降低，基于多流 向算法计算结果的误差逐渐增加，而基于D8 算法的计算结果随着地形条件的不同而呈现出不同的变化趋势。本 文所推荐的TWI 算法定量评价方法通过简单扩展也同样适用于其他复合地形属性算法误差定量评价。 关 键 词：数字地形分析；地形湿度指数；流向算法；栅格DEM ；人造表面；误差定量评价 1 引言 但是这种评价方式无法区分DEM 数据源误差(如 采样误差、插值误差等)和地形湿度指数算法误差 地形湿度指数(Topographic Wetness Index ， 对结果的影响，且评价结果依赖于研究区的实际地 [14] TWI)是一个重要的地形属性，被广泛应用于水文、 形条件 ；另一方面，根据实际DEM 数据无法获得 [1-3] TWI 的理论真值，难以研究各地形湿度指数算法绝 土壤和地貌等研究领域中 。TWI 被定义为单位 汇水面积(Specific Catchment Area ，SCA)与地表坡 对误差的大小和分布。 tan β [4] 本文针对TWI 算法误差定量评价问题，建立了 度正切值( )之比的自然对数 ，即： æSCA ö 一种通过构建人造DEM 模拟多种典型复合地形条 TWI = ln ç ÷ (1) ètan βø 件对TWI 算法误差进行定量评价的方法，该方法能 TWI SCA 够避免以实际DEM 评价地形湿度指数算法中存在 式中： 为地形湿度指数； 被定义为单位长 度等高线上地表水所流经的上游区域的面积，其计