散乱分布数据集的插值及其可视化.pdfVIP

散乱分布数据集的插值及其可视化 周晓嵘林晓明周洞汝 (武汉水利电力大学计算机科学与技术系，武汉430072) [摘要】本文介绍散乱分布数据集的插值函数，重点讨沦了随机过程模型，并介 绍了绘制等值线的游动网格法及三维数据场体绘制的基本方法，最后还用实际的数 据证明了本方法的正确性．本文所提供的方法已成功地应用于流域雨量等值线的绘 制，本文方法可推广至应力、温度、压力、电位等值线的绘制． 一、序言 在水利水电工程中，经常遇到某一区域中物理量的空间分布问题，如雨量、气 压、气温、高度、电位的分布。此时若只提供表达这些物理量的数据集，则难以直 观地分析其物理含义。如果将上述数据集转换为直观可视的图象，如绘制流域雨量 等值线图．流场分布等值线图及坝内应力分布图等，就可清楚地了解其空间分布。 数据集可分为两大类，一类是数据点均在规则网格上采样，且有较大的密度：另一 类是数据点散乱分布于空间，且有较小的密度。显然后一数据集的可视化有较大的 难度。在实际的科学及工程问题中，很多数据集是属于后一类的。因此，研究散乱 分布数据集的可视化有很大的应用价值。 散乱分布数据集的可视化可有两种方法：一种方法是在空间辅设一幅较密的规 则网格．用插值方法求出网格点上的值，然后把此类问题转换为规则网格问题，这种 处理方法要建立中间网格，因此有较大的计算工作量：另一类方法是由空间散乱分布 数据点直接构造一个可描述空间分布的插值函数，由此函数直接得到其图象。由此 可见，解决空间数据集的绘制，主要应解决数据插值及插值函数的图象化的问题。本 文讨论了散乱分布数据集的插值方法及在绘制二维等值线及三维体绘制中的应用。 目前实际问题中所遇到的数据集可分为双变量及多变量两类。双变量数据集表 达的是二维平面上物理量的分布，三变量数据集表达的是三维空间物理量的分布。 pi=(xi，Yi)表示平面上的点，Fi表示此点所对应的某一物理量。如果能找到～个函 的插值函数。双变量数据集可分为规则数据集及散乱数据集两类。规则数据集通常 是科学计算中由人工建立的网格，散乱分布不规则数据集一般是由现实世界中采集 的离散点数据，如随机分布的雨量值，由测量数据绘制地形图等． 二、多维数据集的插值函数 目前，已提出多种数据集插值函数的建立方法。这些方法的基本出发点是建立 通过或逼近已知数据集中各点的光滑的插值函数。目前常用的有按距离倒数加权估 Splines)法、局部 计(Shepard)法、修正二次shepard(MoS)法、体样条(Volume 体样条法等。本文引入了多维数据集随机过程插值函数． 随机过程模型用于绘N-维空间的等值线已研究多年。这里将此模型推广到三 ——147— 维空间。假设v(x，儿z)是一个三维空间的平稳高斯随机过程，其上值均为正则分 布且未知。已知N个离散分布的数据点(x，，Y，，z．)(f’1,2，…，n)的值v． 0=i，2，…，n)。若已知随机过程的数学期望 ∥=E[V(x，Y，z)] 以及协方差矩阵s： (1) 就可以用己知条件去估计随机过程在(工¨，儿m靠+I)的值。由N个已知点可建立 一个线性估计子： 三 口。l=∑bj(vJ一∥)+∥ (2) J-l 式中bi(J=1,2，…，”)为一组待定的系数。 使用多元线性回归的标准方法，可以证明”¨的最佳无偏估计子必须满足： B=S一1C f3) 式中口=IbIb2，…，吃l，C为列向量： G=Covar(v，，Vn+1)，(f-1,2，…，n)(4) 将(3)代入(2)得：