《SURFER常用版本中的各种插值法.docVIP

Surfer常用插值方法 1 反距离加权插值法 反距离加权插值法(Inverse Distance to a Power)首先是由气象学家和地质工作者提出的，后来由于D．Shepard的工作被称为谢别德法（Shepard方法），它的基本原理是设平面上分布一系列离散点，己知其位置坐标(51- ， )和属性值a(i-1，2，)，p(z， )为任一格网点，根据周围离散点的属性值，通过距离加权插值求P点属性值。距离加权插值法综合了泰森多边形的邻近点法和多元回归法的渐变方法的长处，它假设P点的属性值是在局部邻域内中所有数据点的距离加权平均值，可以进行确切的或者圆滑的方式插值。周围点与P点因分布位置的差异，对P( )影响不同，我们把这种影响称为权函数m ， )，方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降对于一个较大的方次，较近的数据点 被给定一个较高的权重份额；对于一个较小的方次，权重比较均匀地分配给各数据点。计算一个格网结点时，给予一个特定数据点的权值，与指定方次的结点到观测点的距离倒数成比例。当计算一个格网结点时，配给的权重是一个分数，所有权重的总和等于1．0。当一个观测点与一个格网结点重合时，该观测点被给予一个实际为1．0的权重．所有其它观测点被给予一个几乎为0．0的权重。换言之，该结点被赋给与观测点一致的值．这就是一个准确插值。权函数主要与距离有关，有时也与方向有关．若在P点周围四个方向上均匀取点．那么可不考虑方向因素，这时 P(Z) ＼、 五 [ ， )]“ ＼ 1 1 式中 ( ， )一-4 t)。+(y— ) ．表示由离散点( ． )至P(x。)点的距离 )为要求的待插点的值。权函数．w ( ，v)一1／Ed,(一v)] ．“值一般取为2。 反距离加权插值法是GIS软件根据点数生成规则格网数据文件的最常见的方法，计算值易受数据点集群的影响，计算结果常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式，可在插值过程中通过滤波来处理。 2 克里金插值法 克里金(Kriging)插值法又称空间自协方差最佳插值法，它是以法国D．G．Krige的名字命名的一种最 优内插法。克里金法广泛地应用于地下水模拟、土壤制图等领域，是一种很有用的地质统计格网化方法 它 首先考虑的是空间属性在空间位置上的变异分布．确定对一个待插点值有影响的距离范围，然后用此范围 内的采样点来估计待插点的属性值。该方法在数学上可对所研究的对象提供一种最佳线性无偏估计(某点 处的确定值)的方法。它是考虑了信息样品的形状、大小及与待估计块段相互间的空间位置等几何特征以 及品位的空间结构之后，为达到线性、无偏和最小估计方差的估计，而对每一个样品赋与一定的系数，最后 进行加权平均来估计块段品位的方法。但它仍是一种光滑的内插方法 在数据点多时，其内插的结果可信 度较高 ：。克里金法类型分常规克里金插值(常规克里金模型／克里金点模型)和块克里金插值。常规克里金插值 其内插值与原始样本的容量有关，当样本数量较少的情况下，采用简单的常规克里金模型内插的结果图会 出现明显的凹凸现象；块克里金插值是通过修改克里金方程以估计子块B内的平均值来克服克里金点模 型的缺点，对估算给定面积实验小区的平均值或对给定格网大小的规则格网进行插值比较适用。块克里金插值估算的方差结果常小于常规克里金插值，所以，生成的平滑插值表面不会发生常规克里金模型的凹凸 现象。按照空间场是否存在漂移(drift)可将克里金插值分为普通克里金和泛克里金，其中普通克里金 (Ordinary Kriging简称OK法)常称作局部最优线性无偏估计．所谓线性是指估计值是样本值的线性组 3 最小曲率法 最小曲率法(Minimum Curvature)广泛应用于地球科学 用最小曲率法生成的插值面类似于一个通过 各个数据值、具有最小弯曲量的长条形薄薄的弹性片。最小曲率法试图在尽可能严格地尊重数据的同时， 生成尽可能圆滑的曲面[s。最小曲率法不是一个精确的插值法．也就是说在插值的过程中不可能总是完全 4 改进谢别德法 改进谢别德法(ModifiedQuadratic Shepard)是由Franke及Nielson提出，它仍是一个与距离成反比的 加权方法。在使用反距离加权插值法时，当增加、删除或改变一个点时，需要重新计算权函数m(z， )，为 了克服反距离加权插值法的这一缺陷，改进谢别德法同样使用距离倒数加权的最小二乘方的方法，但有以 下两个方面的改进 ～ ： (1)通过修改反距离加权插值法的权函数wl(at，y)一1／[d 扛，y)3 ，使其只能在局部范围内起作用， 以改变反距离加权插值法的全局插值性质，即它利用了局部最小二