空间信息技术是自20世纪60年代以来逐渐发展起来的信息技术的总称, 用于获取, 管理和分析与地理位置相关的空间信息[. 它以人类的生活环境为主要研究对象, 以遥感(Remote Sensing, RS)、全球定位系统GPS(Global Position System)和地理信息系统(Geographic Information System, GIS)技术为代表, 它不仅是客观数据的集合, 而且是从收集、处理、测量和分析到地理空间数据和信息的管理、存储、显示和发布的整个信息流程. 它具有客观获取、准确定位、灵活管理、空间分析和视觉表达的特点[.
在空间数据中, 具有不均匀位置分布的数据被称为离散数据, 在平面二维地理空间的定位中, 离散数据的坐标由不规则分布的离散样本的平面坐标实现. 高程和属性值通常用作第三维数据. 空间插值则是一种通过这些离散的空间数据计算未知空间数据的方法. 它是基于“地理学第一定律”的基本假设: 空间位置上越靠近的点, 具有相似特征值的可能性越大, 而距离远的点, 其具有相似特征值的可能性越小[. 它通常用于将离散点的测量数据转换为连续数据表面, 以便于比较其他空间现象的分布情况.
空间插值分析算法的分类方式有多种: 按插值的区域范围分类, 可以分为整体插值、局部插值、边界内插法等[; 整体插值是用研究区的所有采样点进行全区特征拟合, 在整体插值方法中, 整个区域的数值会影响单个插值点的数值, 同样单个采样点的数值的增加、减少或删除对整个区域的特征拟合都有影响[, 代表性插值方法具有趋势面分析插值方法等. 局部插值是使用相邻数据点来估计未知点的值, 首先定义邻域或搜索范围, 然后搜索属于该区域的数据点, 然后选择可以表示此有限点空间变化的数学函数, 最后通过计算为该邻域或者该区域内的未知点赋值[, 代表插值方法有样条函数插值法[、反距离权重插值法[、Kriging插值法[等. 边界内插规则假设值和属性的任何变化发生在特定区域的边界线上, 并且边界内属性的变化是均匀和同质的, 主要的插值法是泰森多边形法.
按照插值的标准分类, 可以分为确定性插值、地统计插值. 确定性插值法主要采用数学工具, 利用函数的方法来进行插值, 这种方式用来研究某区域内部的相似性, 其代表插值法有反距离加权插值法[等; 地统计插值是基于空间自相关性的, 由观测数据产生具有统计关系的曲面, 代表插值法是Kriging插值法[.
按插值的精度分类, 可以划分为精确插值、近似插值. 精确插值生成包括所有观测点的曲面, 而近似插值生成不包含所有观测点曲面[.
随着生产、研究应用的不断深入, 越来越多的插值算法被提出和不断改进. 本文将对适用性高、性能卓越、应用场景广泛的泰森多边形法、反距离权重插值法、样条函数插值法与克里金插值法进行研究综述, 并对空间插值算法的未来的研究方向进行展望.
1 泰森(Thiessen)多边形
荷兰气候学家Thiessen AH提出泰森(Thiessen)多边形法, 根据离散分布的气象站的降雨量来计算平均降雨量, 所有相邻气象站以三角形连接, 在三角形的每一边作垂直平分线, 因此气象站周围有几个垂直平分线包围的多边形[. 用某多边形内所包含的单独的气象站的降雨强度来表示该多边形区域内的降雨强度, 该多边形称为泰森多边形. 如[. 泰森多边形也被称为Voronoi图或dirichlet图[.
图 1
Fig. 1
图 1 泰森多边形示意图
泰森多边形利用离散观测点的值对该点所在的区域进行赋值, 得到的结果往往是数值的变化只发生在多边形的边界上, 而多边形内部的数值则是均匀、同质的[. 其数学表达式为:
${V_e} = {V_i}$
(1)
其中, Ve表示待插值点的距离, Vi表示i示点的离散观测值. i点必须满足如下条件:
${d_{ei}} = \min ({d_{e1}},{d_{e2}}, \cdots, {d_{en}})$
(2)
其中, dij表示点i(xi,yi) 与点j(xj,yj)间的欧几里德距离[.
泰森多边形的关键是将离散观测点合理地连接到三角网络中, 即构造Delaunay三角网络[. 构建泰森多边形步骤如
图 2
Fig. 2
图 2 泰森多边形创建流程
泰森多边形反映了离散观测点的空间控制范围或者是势力范围, 它适用于较小区域内、空间变异性不高的情况, 距离近的点比距离远的点更相似, 比较符合人的逻辑思维[. 同时, 它的实现不需要其他前提条件, 效率高, 方法简单, 但是受样本观测值的影响较大, 没有考虑空间因素、变量以及其他某些规律, 只考虑距离因素, 实际效果不是很理想.
泰森多边形插值法也在进行不断改进发展, 其中自然邻域法就是改进的一种, 它的基本原理是在插值点创建一个新的多边形, 新多边形与原始多边形的重叠比例作为观测点数值的权重, 通过这种方式计算插值点的估计值[. 另外, 刘金雅等人将泰森多边形与最小累计阻力模型相结合[, 用于估算京津冀城市群生态系统服务价值; 祁春阳等人将泰森多边形与虚拟力算法、质心算法相结合, 提出VFVP算法策略来提高无线传感网络覆盖率[; 戚远航等人提出了一种基于泰森多边形的离散