Surfer7中九种插值法介绍.docVIP

等值线绘图软件Surfer 7.0 中九种插值法介绍 白世彪，陈哗，王建 (南京师范大学地理科学学院海岸与第四纪研究所，江苏 南京210097) 摘要：SURFER软件是一个功能强大的绘制等值线图及三维立体图软件包，能迅速地将离散的测量数据通过插值转换为连续的数据曲面。SURFER7.0提供的内插方法多达九种，其中每一种内插方法都有其意义及相关的参数设置。作者在本文中主要介绍了该软件中的九种内插方法并进行了应用示例。 关键词：SURFER软件；插值方法；地学应用 中图分类号：TP316 文献标识码：A 0引言 在地学的研究工作中，不可避免地要应用某种手段获取能够反映研究对象各方面特性或特征的数据，例如磁力、重力、化学勘探数据、钻探获得的岩面标高、海底或陆地的地形标高和地理要素指标等等。为了方便深人研究．要将数据处理并以等值线图的形式呈现。目前，绘制等值线的软件不少，但是，它们往往是为满足各自的需要而以小程序的形式开发出来，采用单一的插值方法，功能比较简单，适用性不广。这里介绍的SURFER是美国GOLDEN软件公司的产品，用于绘制等值线图及相应的三维图形，是科技工作者特别是地学研究人员必备的软件，应用极广。SURFER目前的最高版本为7.03，它提供的内插方法多达九种，用户可根据不同的情况选用不同的内插方法。一般可使用SURFER默认的设置进行内插生成网格文件，但要科学地选择内插方法和灵活地进行参数设置，必须要熟悉各种内插方法的基本理论知识，下面将逐一介绍这九种内插方法。 1 SURFER九种插值法介绍 1.1反距离加权插值法 反距离加权插值法(Inverse Distance to a Power)首先是由气象学家和地质工作者提出的，后来由于D．Shepard的工作被称为谢别德法(Shepard)方法“]。它的基本原理是设平面上分布一系列离散点，己知其位置坐标和属性值，为任一格网点，根据周围离散点的属性值，通过距离加权插值求P点属性值。距离加权插值法综合了泰森多边形的邻近点法和多元回归法的渐变方法的长处，它假设P点的属性值是在局部邻域内中所有数据点的距离加权平均值，可以进行确切的或者圆滑的方式插值。周围点与P点因分布位置的差异，对P(z)影响不同，我们把这种影响称为权函数，方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降。对于一个较大的方次，较近的数据点被给定一个较高的权重份额；对于一个较小的方次．权重比较均匀地分配给各数据点。计算一个格网结点时，给予一个特定数据点的权值，与指定方次的结点到观测点的距离倒数成比例。当计算一个格网结点时，配给的权重是一个分数，所有权重的总和等于1.0。当一个观测点与一个格网结点重合时，该观测点被给予一个实际为1.0的权重．所有其它观测点被给予一个几乎为0.0的权重。换言之，该结点被赋给与观测点一致的值．这就是一个准确插值。权函数主要与距离有关，有时也与方向有关．若在P点周围四个方向上均匀取点．那么可不考虑方向因素，这时： 式中，表示离散点至点的距离，为要求的待插值点。权函数，值一般取为2。 反距离加权插值法是GIS软件根据点数生成规则格网数据文件的最常见的方法，计算值易受数据点集群的影响，计算结果常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式，可在插值过程中通过动态修改搜索准则进行一定程度的改进。使用反距离加权插值法对数据进行格网化时，可以指定一个圆滑参数。圆滑参数要求大于零。对于一个特定的结点．没有哪个观测点被赋予全部的权值，即使观测点与该结点重合也是如此。圆滑参数通过修匀已被插值的格网来降低“鸭蛋”影响。 1.2克里金插值法 克里金(Kriging)插值法又称空间自协方差最佳插值法，它是以法国D．G．Krige的名字命名的一种最优内插法。克里金法广泛地应用于地下水模拟、土壤制图等领域，是一种很有用的地质统计格网化方法。它首先考虑的是空间属性在空间位置上的变异分布，确定对一个待插点值有影响的距离范围，然后用此范围内的采样点来估计待插点的属性值。该方法在数学上可对所研究的对象提供一种最佳线性无偏估计(某点处的确定值)的方法。它是考虑了信息样品的形状、大小及与待估计块段相互间的空间位置等几何特征以及品位的空间结构之后，为达到线性、无偏和最小估计方差的估计，而对每一个样品赋与一定的系数，最后进行加权平均来估计块段品位的方法。但它仍是一种光滑的内插方法。在数据点多时，其内插的结果可信度较高。 克里金法类型分常规克里金插值(常规克里金模型／克里金点模型)和块克里金插值。常规克里金插值其内插值与原始样本的容量有关，当样本数量较少的情况下，采用简单的常规克里金模型内插的结果图会出现明显的凹凸现象；块克里金插值是通过修改克里金方程以估计子块B内的平均值来克服克里金点模型的缺点，对估