基于自然邻点插值的数据处理方法3.pdfVIP

第 22 卷第 3 期 中 国科 学 院 研 究 生 院 学 报 Vol . 22 No . 3 2005 年 5 月 Journal of the Graduate School of the Chinese Academy of Sciences May 　2005 ( ) 文章编号 2005 基于自然邻点插值的数据处理方法 高 　洋 　张 　健 ( 中国科学院研究生院计算地球动力学实验室 , 北京 100049) (2004 年 8 月 6 日收稿 ; 2004 年 10 月 11 日收修改稿) Gao Y ,Zhang J. Data Interpolation research based on the Natural Neighbor method. J ournal of the Graduate School of the ( ) Chinese Academy of Sciences ,2005 ,22 3 :346～351 摘 　要 　利用自然邻点方法对离散的中国大陆地区大地热流观测数据和相对均匀的台湾地区 重力观测数据进行了插值计算 ,并与地学中几种常用插值方法进行了比较. 结果表明: 自然邻 点插值方法更适于处理分布不规则的观测数据. 关键词 　插值方法 , 自然邻点 ,固体地球结构数据 中图分类号 　P3 1 　引言 地球物理场本身是在时间和空间上连续分布的 ,但由于观测手段和技术上的种种限制 ,地球物理学 [ 1] 中的大多数观测数据都是任意分布的离散数据 ,例如大地热流 、地震波速 、重力位场等 . 在对这些数据 分析前一般都需进行插值 ,使不规则的观测点分布于规则的网格之中 ,用离散分布点的值来表示给定区 域内任意一点的物理量值. 通过插值方法来描述连续地球物理场的分布特征是地球科学研究中常用的 数据处理手段. 通常人们利用多项式插值将离散地球物理场连续化 ,但对于数据过分离散或测点分布比 较复杂的问题 ,简单的多项式插值方法往往不能满足要求. 因此在近几十年间产生了许多新的插值方 法 ,如地球物理学常用的 Kriging 插值法 、最小曲率插值法 、样条函数插值法 、Shepard 插值法等. 一种更先 ( ) [2 ,3 ] 进的插值方法 ———自然邻点插值方法 Natural Neighbor Interpolation ,简称 NNI ,现在广泛用于科学计 [4 ] [5 ] 算可视化 、计算流体力学 等领域. 本文研究自然邻点插值法在地球物理学中的应用 , 比较它与其他 插值方法在处理高度离散化分布的不规则节点以及描述空间尺度上剧烈变化方面的效果. 2 　方法 自然邻点插值方法是一种基于 Voronoi Diagram 的插值方法. 该方法对于处理在空间上高度离散化 分布的不规则地球物理观测数据 , 以及描述数据在空间尺度上的剧烈变化具有 良好的效果. Voronoi Diagram 是一种常用的非结构化网格 ,它所剖分的每一个网格单元称为 Voronoi Cell ( 图 1 (a) ) . Voronoi