自然邻点插值算法及其在二维不规则数据网格化中的应用.pdfVIP

第33卷第3期 物探化探计算技术 2011年5月 文章编号：100l—1749(2011)03珈29l—05 自然邻点插值算法及其在二维不规则数据 网格化中的应用 张 伟，覃庆炎，简兴祥 (成都理工大学信息工程学院，四川成都610059) 摘要：由于观测数据的不足。对不规则散点数据的二维插值在地球物理中应用非常广泛。这 里研究并实现了自然邻点二维插值算法，并以一个大地电磁测深反演数据为例，验证了该方法的 可行性与效果。实际资料试算与成图结果表明，该方法具有精度高、插值效果好、速度快以及便 于模拟地形数据等优点。 关键词：自然邻点；二维插值；网格化 中图分类号：P63l 文献标识码：A NNI所具有的优良特性，使其成为地球物理网 O前言 格化处理最为合适的一种插值算法之一。作者在 对NNI算法原理研究的基础上，实现了该插值算 地球物理场在空间上是连续分布的，但由于 法，并将其应用到大地电磁反演数据的实时网格化 数据采集技术和工区客观条件上的种种限制，地球 成图中。实验结果表明，该算法对于横向尺度、纵 物理中所得到的数据都是不规则的散点数据。在 向尺度。都有比较大的二维数据体网格化插值计 对地球物理数据成图解释前，需要对不规则分布的 算，具有速度快、精度高、运算时间短等优势，可以 散点数据进行网格化处理，通过插值方法来描述连 达到实时网格化成图的要求。 续地球物理场的分布特征，是地球科学研究中常用 的数据处理手段。 l 自然邻点插值算法原理 目前，普遍使用的网格化方法主要有最小曲率 Cunrature 1．1 自然邻点 插值法(Minimum Splines)、陆百ng插值 法、双线性插值法，以及贝塞尔插值等uJ。自然邻 自然邻点插值方法是一种基于Dela叽ay三角 111terl)0lation，NNI)， 点插值算法(Natu融Neighbor 也被称为AreaSte址ng插值法，该方法广泛应用在 地球物理建模，表面重构，科学计算可视化以及计 的各个三角形中，对每一边画出其垂直平分线，即 算流体力学等领域。该方法主要有以下三个特 性‘2l： 图是一种常用的非结构化网格，它可将数据点所在 的平面，划分成数个多边形，每一个多边形只包含 (1)原始点在经过NNI插值后，仍然保持不变。 (2)插值过程是局部的，待插点的值只受其周 一个数据点。下页图l(b)是一个包含六个数据点 围的点(自然邻点)的影响。 的v0咖i图，该图所剖分的每一个网格单元，被 (3)除了原始点以外，插值函数的一阶导数处 称为Vomnoi 处连续。 对应的VomnoiceU具有共同边的邻近节点称为自 国地质调查局地质调童工作项目(1212010914049) 收稿日期：2010—10一19改回日期：201l—03—17 物探化探计算技术 33卷 然邻点(Nat删融NeigIlbors，NN)，如图1(b)中，单 元节点p4的自然邻点就是与其具有共同边界的节 点pl、p2、p3、几、p6。