第六讲空间插值概要.ppt

二、空间插值方法 6.5 kriging插值方法 式中：Z0为待插入点的值； Zx为已知点的值 Wx为每个点的权重值 计算Wi，按采样点数据的变异函数图的统计分析原理来计算 二、空间插值方法 6.5 kriging权重计算方法 h11 r(h11) 式中：r(hij)是已知点i和j之间的变异值，r(hi0)是已知点i和待估点之间的变异值， r(hij)和r(hi0)可以通过变异函数求得； Wi为每个已知点的权重，u是拉格朗日系数，是为了将估计误差降低到最小； 根据上式，通过解算联立方程可以计算出Wi 二、空间插值方法—待续 习题 构建泰森多边形 Kriging估值 块金 变程 基台 1.0 110(m) 3.1 谢 谢 ！ 第六讲 空间插值 金毅 河南理工大学资源环境学院 454003 ＱＱ电话E_mail:jinyi2005@hpu.edu.cn 基本内容 空间插值：定义及应用 空间插值方法及特征 泰森多边形（ Voronoi ）及不规则三角网（TIN） 距离反比加权法（IDW） 地质统计学（Geostatistics) 利用样条曲线优化插值结果 插值精度评估 三参数插值方法（体数据或者动态演化特征） 为何进行插值？ 如基于空间离散点，剖面数据和等高线等来构建连续的栅格数据 生成更密、更光滑的TIN或者栅格单元。 1. 2D离散点转化为连续面，如地表、地层界面 为何进行插值？ 1. 2D离散点向规则栅格数据的转化 给定的离散样品数据 规则格网数据 为何进行插值？ 2. 3D离散点转化为连续体数据，土壤属性、金属分布 高程[m] 土样空间构型 利用3D空间样品数据转化为3D连续的体数据 砂岩 黏土 为何进行插值？ 3. 动态采样的2D、3D数据转化为动态演化的时空连续模型 不同时期地表采 样点高程的变化 恢复的三维形貌 地下水污染扩散模型 一、空间插值的基本原理 基于离散样品的空间和属性信息，根据预定的精度要求或确定的规则，推演出时空演化模型（函数）。 利用该函数，借助或完全脱离样品数据，估算未采样位置的目标属性（属性和空间信息）的过程。 利用插值结果，估算插值的精度和误差。 1. 空间插值的定义 一、空间插值的基本原理 1. 空间插值的定义-二维解释 一、空间插值的基本原理 2. 地质现象的多解性 r rj 需要额外条件！！！ 二、空间插值方法 1. 插值条件—spatial interpolation conditions 根据不同的条件假设，可以定义如下几类插值方法 定域条件(Locality Condition)插值方法： 样品点只对其周围一定范围内的数据有影响 即一给定位置的数值同其领域相似—相近者相似 地质统计学方法(Geostatistical Condition) 任何地学现象都是某一随机函数的一个实现 光滑插值方法： 插值函数在确保通过或接近采样点的同时尽量光滑 二、空间插值方法 2. 局部插值方法(Local interpolation) 在局部插值的过程中，待估点属性值得获取只需要全局样本的一部分子集。 泰森多边形插值： 基于不规则三角网的插值 距离反比加权插值 二、空间插值方法 2.1 泰森多边形和Delaunay三角剖分 Delaunay三角剖分：将样品点集投影到某一平面，根据这些点构建三角形，任何三角形的外接圆内无样本点。 二、空间插值方法 2.1 泰森多边形和Delaunay三角剖分 最近邻点法又叫泰森多边形方法。它采用一种极端的边界内插方法—只用最近的单个点进行区域插值（区域赋值）。 泰森多边形按数据点位置将区域分割成子区域，每个子区域包含一个数据点，各子区域到其内数据点的距离小于任何到其它数据点的距离，并用其内数据点进行赋值。 二、空间插值方法 2.1 泰森多边形评价 特征：得到的结果图变化只发生在边界上，在边界内都是均质的和无变化的； 适用于较小的区域内，变量空间变异性也不很明显的情况； 优点是不需其他前提条件，方法简单，效率高； 缺点是受样本点的影响较大，只考虑距离因素，对其他空间因素和变量所固有的某些规律没有过多地考虑。实际应用中，效果常不十分理想。 二、空间插值方法 3. 局部邻域算数平均法 算术平均值方法以区域内所有测值的平均值来估计插值点的变量值(Creutin, 1982) ； 算术平均值的算法比较简单，容易实现。但只考虑算术平均，根本没有顾及其他的空间因素，这也是其一个致命的弱点，因而在实际应用中效果不理想。 二、空间插值方法 4. 距离反比加权法 距离反比插值方法最早由 Shepard 提出(Richard Franke,1982)提出的，并逐步得到发展； 每个采样点对插值结果的影响随距离增加而减弱，因此