ANUDEM----专业化数字高程模型插值软件-Read.PDFVIP

ANUDEM 专业化数字高程模型插值软件 1* 2 2 3 1 杨勤科 ，Tim R. McVicar ，李领涛 ，Tom G. Van Niel ，张彩霞 (1. 中国科学院水利部水土保持研究所；陕西，杨凌，西农路 26 号，712100; 2. CSIRO Land and Water, GPO Box 1666, Canberra, 2601, ACT, Australia; 3. CSIRO Land and Water, Private Bag No. 5, Wembley, 6913,WA) 摘要：基于已有地形图的数字化，通过内插方法建立 DEM ，是迄今为止最为成熟和简单易行的 DEM 建立 方法。ANUDEM 是利用等高线、高程点和河流等数据插值建立 DEM 的专业化软件。该软件能够自动诊断 数据错误，客观地确定适应于插值基础数据精度的分辨率，建立能正确反映地面水文地貌特征的 DEM ，因 而在全球范围内已得到广泛应用，其算法被集成于 ARCGIS 之中。本文根据有关资料和我们的初步应用对 ANUDEM 做一简单介绍，以期改进我国DEM 建立方法，为流域水文、土壤侵蚀等研究提供基础数据支持。 关键词：ANUDEM ；DEM ；插值算法；地形强化算法 1 引言 数字高程模型 （digital elevation model, DEM ）是利用一组有序数据阵列对地表高程空间变化的定 量表现[1] 。DEM在地貌定量分析、流域水文和土壤侵蚀模拟分析、植物分布、遥感图像辅助分类等方 面具有重要用途[1,2,3] 。生成DEM 的方法，根据基础数据不同分为两类，一是利用遥感数据提取高程生 成DEM[4] ，二是根据已有测绘资料经过插值生成DEM 。迄今为止最为成熟和可行的方法是将现有地 形图数字化，通过插值生成DEM[1,2] 。就后者而言，国内最常用的方法是利用多种高程信息 （等高线、 高程点，骨架点），通过构建不规则三角网 （Triangle Irregular Network, TIN ）建立DEM[5-7]，该方法的 缺点是如果不加入新的高程数据而仅利用等高线和高程点，生成DEM将不能很好的反映河流 （沟道） 和流域边界等重要的水文地貌特征[1] 。国外较为流行的方法是Hutchinson教授于 20 世纪 80 年代提出 的DEM插值算法[8]及其基础上开发的专业化DEM生产软件ANUDEM[9,10] 。ANUDEM 以其简便易行、 能够建立正确反映水文地貌关系的DEM而在全球范围内得到广泛应用[11-14] ，其算法已被ARC/INFO地 理信息系统软件集成（TOPOGRID ）。但是，国内对Hutchinson算法极少应用，仅有的文献也因为对 [15] (1) DEM这个专业术语翻译的不甚准确 或介绍不甚详细 而尚未引起重视。本文力图根据有关文献和我 们在研究中的初步应用[16]，对ANUDEM做一简要介绍，以期改进我国DEM建立方法，为流域水文、 土壤侵蚀等相关研究提供更有力的支持。 2 ANUDEM 的算法 ANUDEM 的算法主要包括四个方面，即插值方法（Interpolation Algorithm ）、数据平滑（roughness Penalty Algorithms ）、地形强化 （drainage Enforcement Algorithm ）和局地适应性处理 （locally adaptive strategies ）。 2.1 插值算法 已有的插值方法可粗略划分为两类，即整体插值法和局部插值法。整体插值方法生成的表面是 连续的，但一般算法插值计算所需时间与观测点数量的 3 次方成正比，所以当数据点较多时对计算 机的计算性能要求很高以至难以负担。局部方法一般是将区域划分为若干小块，为每小块拟合一函 数并使这些函数在边界处连续，这种方法的主要问题是对插值区域的划分不好把握。 * 通讯作者：杨勤科，男，43 岁，主要从事基于遥感和地理信息系统技术的区域水土流失评价研究。email: qkyang@ 资助项目：中国科学院知识创新重要方向项目 :黄土高原水土保持的区域环境效应研究，合同号：KZCX3-SW-421 ；中澳合作项目黄土高 原植被恢复重建设水环境效应研究，合同号 LWR