局部最优形态参数的RBF分块地形插值方法和实验.pdfVIP

第17卷 第3 期 Vol. 17, No.3 2015 年3 月 Mar., 2015 局部最优形态参数的RBF 分块 地形插值方法与实验 1 1* 1 1 2 吕海洋，盛业华 ，段 平 ，张思阳 ，李 佳 (1. 南京师范大学虚拟地理环境教育部重点实验室，南京210023 ； 2. 云南师范大学旅游与地理科学学院，昆明650500) 摘要：径向基函数（Radial Basis Function ，RBF ）是一种不需对数据做任何假设，能准确逼近任意维度数据的空间插 值方法。其特别适合于复杂地形的数字高程模型（Digital Elevation Model ，DEM ）插值重建，但随着已知点数量的 增加会导致插值模型求解困难或缓慢。针对这个问题，本文基于二叉树自适应递归分块原理，采用局部最优形态 参数的RBF 分块插值方法进行DEM 插值重建。首先，设定子区域最小点数阈值和相邻子区域的重叠率，自顶向 下，对研究区域进行递归分块，构建区域分块二叉树，对二叉树叶子节点区域，采用逐点交叉验证（Leave One Out Cross Validation ，LOOCV ）方法求解其最优形态参数，建立局部RBF 最优插值模型；然后，根据单元分解原理，采用 加权平均方法对相应叶子节点区域内的待插值点高程进行加权求和，自底向上递归求解，得到待插值点最终高程 值。以云南某地区DEM 进行插值实验表明，采用本文方法进行DEM 插值重建，稳定性较好，插值精度高。 关键词：形态参数；DEM ；RBF ；分块插值 DOI:10.3724/SP.J. 1047.2015.00260 1 引言 RBF 、全局多项式和局部多项式方法等）拟合得到 的空间曲面可以精确通过已知观测点，具有局部保 目前，常见的DEM 插值方法有线性内插、反距 形、准确逼近等优势。Franke 从插值参考点及权函 离权重（Inverse Distance Weighted ，IDW ）、二元样 数的选择、计算精度等方面进行实验分析，得出 [6] 条、最小二乘配置、自然邻近、三角形插值、克里金 RBF 插值方法具有更高的可靠性 。RBF 作为一种 [1-2] 维度无关的精确空间插值方法，分为全局RBF 插值 插值、径向基函数等 。李新等人将这些插值方法 [7-9] 归类为几何方法、统计方法、空间统计方法、函数方 方法和局部RBF 插值方法 。全局RBF 插值方法 法、随机模拟方法、物理模型方法和综合方法，并对 精度较高，但插值模型中存在形态参数和基函数的 [3] 选取问题，同时随着数据量增加插值模型求解困 每种方法的适用范围和优缺点进行了分析和比较 。 Li 和Heap 系统地总结和对比分析了在地学领域常 难，限制了RBF 全局插值方法的广泛运用，而DEM 用的42 个确定性插值、空间统计学插值和综合插值 插值重建一般需要大量已知高程点才能得到足够 方法的