第10讲-空间插值论述.ppt

样条函数求解步骤 样条函数示例 Regularized spline with tension Thin plate spline with tension 样条函数应用 平化和连续的面 特征： 较IDW法更为平滑 受控制点分布影响较大 数据贫乏区，插值结果较实际大 张力系数一般不宜过大，在0.1-0.5之间 克里金（Kriging）插值理论简介 主要问题 合理的插值邻域控制点点数 邻域大小 除距离外的合理权值 精度评定 逐点内插 所谓逐点内插，就是以内插点为中心，确定一个邻域范围，用落在邻域范围内的采样点计算内插点的高程值 逐点计算格网点的 高程: 逐点内插法的基本步骤为： 定义内插点的邻域范围； 确定落在邻域内的采样点； 选定内插数学模型； 通过邻域内的采样点和内插计算模型计算内插点的高程。 为实现上述步骤，逐点内插法需要解决好以下几个问题： 内插函数; 邻域大小和形状; 邻域内数据点的个数; 采样点的权重; 采样点的分布; 地学统计简介 地学统计简介 法国数学家George Matheron和南非矿业工程师 D.G. Krige 区域变量理论 随机变量和确定变量之间 地理分布现象不规则，不能用平滑数学函数进行模拟 内在假定 差异的稳定性和可变形 当结构成分确定后，差异变化为同性变化，不同位置之间的差异仅为距离的函数 半方差理论 半方差理论（Semivariance） 空间相关性的度量指标 定义为所有恒定距离的成对点的方差之半 特性 随距离增加而按确定规律变化 距离=0，半方差=0； 参数 梁 值域 融核 半方差图 半方差 梁 值域 融核 距离 半方差计算步骤 半方差计算公式 公式 计算 规则分布点 不规则分布点 半方差图 半方差模型（拟和参数计算） 半方差计算：规则分布点 d d=1 d=2 d=3 半方差计算：不规则分布点 3 1 5 7 9 2 9 8 x y 半方差图与拟和方差图 d γ(d) 高斯 线性 球形 指数 圆形 克里金插值：普通克立金 克里金插值：通用克立金 其他克立金法 块克立金法 某一范围内一个变量的平均值 联合克立金法 加入一个或更多的次要变量，而次要变量与主要变量有相关关系 在降水量插值中，将高程作为附加变量会获的更好的结果 克里金插值示例 各种局部方法比较方法 交叉验证 不同方法之间的比较 方法 从控制点中删除一点 用不同方法进行删除点的估计 计算误差 误差分析方法 控制点为两部分 一部分建模，一部分检测精度 概论 Interpolation = inter + polire, means Refining by put in between. 将离散的数据点转化为连续的数据曲面 用已知点来估算其他未知点的过程 需要插值的原因 现有离散曲面的分辨率、像元大小、方向与要求不符； 现有连续曲面的数据模型与要求不一致； 现有数据不能完全覆盖所要求的区域 内插和外推 Sampled points Estimated points 空间插值基础：控制点 控制点分布 控制点密度 控制点的自相关程度 内插方法 全局方法 趋势面 回归模型 局部方法 密度估算 反距离权内插 样条函数内插技术 克里金内插方法 内插方法：全局内插 z=f( x, y) Control Points 全局内插的步骤 数学曲面函数确定 内插曲面的复杂程度 计算量 系数求解 最小二乘法 拟合精度分析 常用全局内插函数 常用全局内插函数特性分析 1st Order Trend Original surface 2nd Order Trend 3rd Order Trend 系数求解 拟合精度评定 一般认为，C在60%到70%之间拟和比较好 但由于以下原因，在空间内插中整体内插并不常用： 整体内插函数保凸性较差； 不容易得到稳定的数值解解; 多项式系数物理意义不明显; 解算速度慢且对计算机容量要求较高; 不能提供内插区域的局部地形特征. 全局内插示例：三阶趋势面 局部内插 块内逐点内插格网点高程 分块 已知点 格网点 局部内插：密度估算 Voronoi多边形 在一个多边形内，每个未知点与该多边形内的已知点最接近，而与其他已知点更远 密度估算 计算在局部范围内的点的数量 简单密度估算 核密度估算 Voronoi多边形 密度估算 简单密度估算 D = 点数/单元面积 核密度估算 已知点与核密度函数联系起来，用概率密度函数表达。 窗口范围、核 密度为窗口范围内观测点上的隆起部分的总和 核密度估算 核 带宽 密度估算示例 局部插值：反距离权插值IDW 基本思想：空