1、 Surfer 常用插值方法1 反距离加权插值法 反距离加权插值法(Inverse Distance to a Power)首先是由气象学家和地质工作者提出的，后来由于 DShepard 的工作被称为谢别德法（Shepard 方法） ，它的基本原理是设平面上分布一系列离散点，己知其位置坐标(51- ， )和属性值 a(i-1，2，)，p(z， )为任一格网点，根据周围离散点的属性值，通过距离加权插值求 P 点属性值。距离加权插值法综合了泰森多边形的邻近点法和多元回归法的渐变方法的长处，它假设 P 点的属性值是在局部邻域内中所有数据点的距离加权平均值，可以进行确切的或者圆滑的方式插值。周围点与
2、P 点因分布位置的差异，对 P( )影响不同，我们把这种影响称为权函数 m ， )，方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降对于一个较大的方次，较近的数据点 被给定一个较高的权重份额；对于一个较小的方次，权重比较均匀地分配给各数据点。计算一个格网结点时，给予一个特定数据点的权值，与指定方次的结点到观测点的距离倒数成比例。当计算一个格网结点时，配给的权重是一个分数，所有权重的总和等于 10。当一个观测点与一个格网结点重合时，该观测点被给予一个实际为 10 的权重所有其它观测点被给予一个几乎为 00 的权重。换言之，该结点被赋给与观测点一致的值这就是一个准确插值。权函数主要与距
3、离有关，有时也与方向有关若在 P 点周围四个方向上均匀取点那么可不考虑方向因素，这时 P(Z) 、 五 ， )“ 1 1 式中 ( ， )一-4 t)。+(y ) 表示由离散点( )至 P(x。)点的距离 )为要求的待插点的值。权函数w ( ，v)一 1Ed,( 一 v) “值一般取为 2。 反距离加权插值法是 GIS 软件根据点数生成规则格网数据文件的最常见的方法，计算值易受数据点集群的影响，计算结果常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式，可在插值过程中通过滤波来处理。2 克里金插值法 克里金(Kriging)插值法又称空间自协方差最佳插值法，它是以法国DGKrige 的名
4、字命名的一种最 优内插法。克里金法广泛地应用于地下水模拟、土壤制图等领域，是一种很有用的地质统计格网化方法 它 首先考虑的是空间属性在空间位置上的变异分布确定对一个待插点值有影响的距离范围，然后用此范围 内的采样点来估计待插点的属性值。该方法在数学上可对所研究的对象提供一种最佳线性无偏估计(某点 处的确定值 )的方法。它是考虑了信息样品的形状、大小及与待估计块段相互间的空间位置等几何特征以 及品位的空间结构之后，为达到线性、无偏和最小估计方差的估计，而对每一个样品赋与一定的系数，最后 进行加权平均来估计块段品位的方法。但它仍是一种光滑的内插方法 在数据点多时，其内插的结果可信 度较高 ：。克里
5、金法类型分常规克里金插值 (常规克里金模型克里金点模型)和块克里金插值。常规克里金插值 其内插值与原始样本的容量有关，当样本数量较少的情况下，采用简单的常规克里金模型内插的结果图会 出现明显的凹凸现象；块克里金插值是通过修改克里金方程以估计子块 B 内的平均值来克服克里金点模 型的缺点，对估算给定面积实验小区的平均值或对给定格网大小的规则格网进行插值比较适用。块克里金插值估算的方差结果常小于常规克里金插值，所以，生成的平滑插值表面不会发生常规克里金模型的凹凸 现象。按照空间场是否存在漂移(drift) 可将克里金插值分为普通克里金和泛克里金，其中普通克里金 (Ordinary Kriging
6、简称 OK 法)常称作局部最优线性无偏估计所谓线性是指估计值是样本值的线性组3 最小曲率法 最小曲率法(Minimum Curvature)广泛应用于地球科学 用最小曲率法生成的插值面类似于一个通过 各个数据值、具有最小弯曲量的长条形薄薄的弹性片。最小曲率法试图在尽可能严格地尊重数据的同时， 生成尽可能圆滑的曲面s。最小曲率法不是一个精确的插值法也就是说在插值的过程中不可能总是完全 4 改进谢别德法 改进谢别德法(ModifiedQuadratic Shepard)是由 Franke 及 Nielson 提出，它仍是一个与距离成反比的 加权方法。在使用反距离加权插值法时，当增加、删除或改变一个
7、点时，需要重新计算权函数 m(z， )，为 了克服反距离加权插值法的这一缺陷，改进谢别德法同样使用距离倒数加权的最小二乘方的方法，但有以 下两个方面的改进 ： (1)通过修改反距离加权插值法的权函数 wl(at“，y)一1d 扛，y)3 ，使其只能在局部范围内起作用， 以改变反距离加权插值法的全局插值性质，即它利用了局部最小二乘方法来消除或减少所生成等值线的 “鸭蛋”外观 (2)同时用节点函数 Q(一 y)来代替离散点(。 。 )的属性值，Q( ， )是一个插值于( ， )点的 二次多项式即有 Q( 。 )一( 1。 2，? ，n)。而且 Q )在点 ， )附近与函数属性值 z(x， )具 有
8、局部近似的性质。因此，如果认为距离( ， )较远的点对 Q( ， )影响不大，则可以认为在 (z ， )点 5 自然邻点插值法 自然邻点插值法(Natural Ndghbor)是 Surfer7、0 才有的网格化新方法。自然邻点插值法广泛应用于 一些研究领域中。其基本原理是对于一组泰森(Thiessen) 多边形，当在数据集中加入一个新的数据点(目 标)时，就会修改这些泰森多边形，而使用邻点的权重平均值将决定待插点的权重，待插点的权重和目标泰 森多边形成比例 。 。实际上，在这些多边形中，有一些多边形的尺寸将缩小，并且没有一个多边形的大小会 增加。同时，自然邻点插值法在数据点凸起的位置并不外推
9、等值线(如泰森多边形的轮廓线)。 6 最近邻点插值法 最近邻点插值法(Nearest Neighbor)又称泰森多边形方法，泰森多边形(Thiesen，又叫 Dirichlet 或 Voronoi 多边形)分析法是荷兰气象学家 AH_Thiessen 提出的一种分析方法。最初用于从离散分布气象 站的降雨量数据中计算平均降雨量，现在 GIS 和地理分析中经常采用泰森多边形进行快速的赋值_2。实际 上，最近邻点插值的一个隐含的假设条件是任一阿格点 p(x，y)的属性值都使用距它最近的位置点的属性 值，用每一个网格节点的最邻点值作为待的节点值 。当数据已经是均匀问隔分布，要先将数据转换为 SURFE
10、R 的网格文件，可以应用最近邻点插值法；或者在一个文件中，数据紧密完整，只有少数点没有取 值，可用最近邻点插值法来填充无值的数据点。有时需要排除网格文件中的无值数据的区域，在搜索椭圆(Search Ellipse)设置一个值，对无数据区域 赋予该阿格文件里的空白值。设置的搜索半径的大小要小于该阿格文件数据值之间的距离，所有的无数据 同格节点都被赋予空白值。在使用最近邻点插值同格化法，将一个规则间隔的 z 数据转换为一个网格7 多元回归法 多元回归(Polynomial Regression)是用来确定数据的大规模的趋势和图案 多元回归实际上不是插值 器，因为它并不预测未知的 z 值，它只是根据
11、空间的采样数据，拟台一个数学曲面，用该数学曲面来反映空间分布的变化情况，它实际上是一个趋势面分析作图程序。趋势面分析是对地质特征的空同分布进行研究和分析的一种方法，它是用某种形式的函数所代表的 曲面来逼近该地质特征的空间分布。这个函数从总体上反映了采样数据的区域性变化趋势，称为趋势面部 分；采样数据的实测值与这个函数对应值之差，称为偏差部分，它反映了，局部性的变化。这就是说，把采样 数据的实测值分解成两部分，趋势面部分和偏差部分，趋势面部分用一个函数表示，它反映采样数据的总 体变化，可以认为是由大范囝的系统性因素引起的；偏差部分反映了，局部性的变化特点，可以认为由局部 因素和随机因素引起的一
12、一。 使用多元回归法进行趋势面分析要考虑两个方面的问题 一是趋势面函数 (数学表达式)的确定；二是 拟合精度的确定 通常用的趋势面函数主要是多项式趋势面，因为多项式能够逼近任意连续函数，因此，用 多项式作趋势面能较好地反映连续变化的分布趋势，这在地质科学中常用到。一般说多项式次数越高，则 趋势面与实测数据偏差越小，但是，并不能说它就与实际情况最符合，这还要在实践中检验次数较高的趋 势面只在采样点附近效果较好，在外推和内插的效果方面不好，因而在实际应用的效果并不理想。在实际 应用中，对起伏变化比较缓和的简单采样数据配合次数较低的趋势面，就可以反映出区域背景；而变化复 杂且起伏较多的采样数据要配合
13、次数较低高的趋势面 l 。 使用多元回归法时要涉及到曲面定义和指定 xY 的最高次数设置，在曲面定义中选择所需的多项式 类型，可选用的曲面类型：简单平面(Simple planar surface)、双线性鞍(Bilinear saddle)、二次曲面8 径向基函数插值法 所谓径向基函数(Radial Basis Function)，其基函数是由单个变量的函数构成的。一个点 ， )的这种 基函数的形式往往是髓 ， )=(4)，这里的盔表示由点(z， )到第 i 个数据点的距离。径向基函数插值 法是多个数据插值方法的组合。根据生成一个圆滑曲面适应数据的能力许多人认为其中的复二次函数是 最好的方法
14、 一。所有径向基函数插值法都是准确的插值器，它们都能尽量适应你的数据。若要生成一个更 圆滑的曲面对所有这些方法都可以引入一个圆滑系数。 函数类型：最基本的函数类似于克里金中的方差图。当对于一个网格点插值时这些函数为数据点规定了一套最佳权重。 基函数类型有：倒转复二次函数(Inverse Multiquadric)： B()一 1 复对数(Multilog)：B()log(h 十尺 ) 复二次函数(Multiquadratic)：B()一2+R 自然三次样条函数(Natural Cubic Spfine) ：B()一( +R ) 薄板样条法函数(Thin Plate Spline) ：B()=(
15、 +尺 )log(h T 尺 ) 式中 h 为表示由点(z， )到第 i 个数据点的距离； R 参数是用户指定的平滑因子。其中的复二次函数(Multiquadrie)方法是由 RI Hardy 在 1971 年提出来的 它是最早提出并且应用 得最为成功的一种径向基函数插值法。现在复二次函数(Mukiquadric)方法在水文测量、大地测量、地质及 采矿、地球物理等领域都得到了，广泛应用，效果良好。在数据点数量不太大的情况下(例如，几百个数据 点)，计算也不太复杂。 薄板样条法(Thin PLate Spline)是由 RLHarder 及 RNDesmarais 在 1972年提出来的，后来得
16、到了 9 三角网线性插值法 三角网线性插值法(Triangulation with Linear Interpolation)使用最佳的 Delaunay三角形，连接数据 点间的连线形成三角形 原始数据点的连结方法是这样：所有三角形的边都不能与另外的三角形相交，其 结果构成了一张由三角形拼接起来的覆盖网格范囝的网。每一个三角形定义了一个覆盖该三角形内网格 节点的面 。三角形的倾斜和标高由定义这个三角形的三个原始数据点确定。给定三角形内的全部节点 SURFER 所采用的九种离散数据内插方法，几乎包括了目前所有的插值方法。该软件具有如此种类繁多的内插方法，这是其它同类软件所不能比拟的。一般可使用
17、SURFER 默认的设置进行内插，生成网格文件，然后再绘制所需的图件 熟悉了各种内插方法的基本理论知识，根据各种数据的不同特点，结合数据分析的目的，就能科学地选择内插方法，灵活地进行参数设置内插生成网格文件，在此基础上绘制正确的图件，供我们提取更多的地学信息。 本 文 来 自 : 博 研 石 油 论 坛 详 细 出 处 参 考 http:/ 0 所 采 用 的 九 种 离 散 数 据 内 插 方 法 ， 几 乎 包 括 了 目 前 所 有 的 插 值 方 法 。该 软 件 具 有 如 此 种 类 繁 多 的 内 插 方 法 ， 这 是 其 它 同 类 软 件 所 不 能 比 拟 的 。 一 般
18、 可 使 用SURFER 默 认 的 设 置 进 行 内 插 ， 生 成 网 格 文 件 ， 然 后 再 绘 制 所 需 的 图 件 熟 悉 了 各 种 内插 方 法 的 基 本 理 论 知 识 ， 根 据 各 种 数 据 的 不 同 特 点 ， 结 合 数 据 分 析 的 目 的 ， 就 能 科学 地 选 择 内 插 方 法 ， 灵 活 地 进 行 参 数 设 置 内 插 生 成 网 格 文 件 ， 在 此 基 础 上 绘 制 正 确的 图 件 ， 供 我 们 提 取 更 多 的 地 学 信 息 。 本 文 来 自 : 博 研 石 油 论 坛 详 细 出 处 参 考 http:/ 常 用 功 能 介 绍
