课件_空间数据的处理空间插值.pptVIP

地理信息系统基础 Geographical Information System 空间插值 教学内容 §1 插值概论 §2 内插方法 §3 整体内插 §4 局部内插 §5 Kriging 方法 §1 插值概论 插值 用已知点来估算其他未知点的过程。 在GIS中，空间插值主要用于栅格数据，估算出栅格中每个单元的值。空间插值是将点数据转换为面数据的一种方法。 需要插值的原因: 现有数据不能完全覆盖所要求的区域 现有离散曲面的分辨率、像元大小、方向与要求不符； 现有连续曲面的数据模型与要求不一致； 内插和外推 内插和外推 内插 在已观测点的区域内估算未观测点的数据的过程。 外推 在已观测点的区域外估算未观测点的数据的过程.——预测。 空间插值基础：控制点 §2 内插方法 全局内插 在整个区域用一个数学函数来表达地形曲面。 采用全部控制点计算未知点的数据。 局部内插 将复杂的地形地貌分解成一系列的局部单元，在这些局部单元内部地形曲面具有单一的结构，由于范围的缩小和曲面形态的简化，用简单曲面就可较好的描述地形曲面。 采用样本控制点计算未知点的数据。 内插方法 整体内插： 趋势面分析 回归模型 局部方法： 泰森多边形（ Voronoi边形，边界内插） 反距离权内插 样条函数内插技术 克里金内插方法 密度估算 内插方法分类 §3 全局内插 趋势面分析 趋势面分析的步骤 数学曲面函数确定 内插曲面的复杂程度 计算量 系数求解 最小二乘法 拟合精度分析 趋势面分析内插函数 常用全局内插函数分类 常用全局内插函数 常用全局内插函数特性分析 常用全局内插函数特性分析 一阶线性平面可模拟具有单一坡度的斜坡地形表面； 二次曲面方程可表达山头、洼地区域； 三次曲面则能描述较为复杂的地形曲面。 系数求解 系数求解 拟合精度评定 （2）回归模型 建立因变量与自变量的联系 自变量选择 非空间属性、空间属性 影响程度 求解方法与趋势面类似 例：流域雪水量回归模型 swe = a + b*ROW + c*COL + d*ELV + e*PlanCur 整体内插缺点 由于以下缺点，在空间内插中整体内插并不常用： 整体内插函数保凸性较差 ：采样点的增减或移动都需要对多项式的系数作全面调整，从而采样点之间会出现难以控制的振荡现象，致使函数极不稳定，从而导致保凸性较差； 多项式系数物理意义不明显 解算速度慢且对计算机容量要求较高。 不能提供内插区域的局部地形特征； 整体内插函数保凸性较差 采样点之间会出现难以控制的振荡现象，致使函数极不稳定 整体内插优点 优点：整个区域上函数的唯一性、能得到全局光滑连续的空间曲面、充分反映宏观地形特征等。 在空间内插中，一般是与局部内插方法配合使用，例如在使用局部内插方法前，利用整体内插去掉不符合总体趋势的宏观地物特征。整体内插函数常常用来揭示整个区域内的地形宏观起伏态势。 全局内插示例：三阶趋势面等值线地图 §4 局部内插 全局内插局限 低阶多项式虽然可表达各种地形曲面，但一个地区确常常包含各种复杂的地貌形态，简单的曲面并不能很好的表达这些地形曲面。理论上任何复杂的曲面都可用多项式进行逼近，但高阶多项式的缺点，也不是理想的地形描述工具。 解决这类问题的办法 分而治之，即将复杂的地形地貌分解成一系列的局部单元，在这些局部单元内部地形曲面具有单一的结构，由于范围的缩小和曲面形态的简化，用简单曲面就可较好的描述地形曲面。 局部内插 局部内插 将地形区域按一定的方法进行分块，对每一块根据地形曲面特征单独进行曲面拟合和高程内插，每一块都可用不同的曲面进行表达。称为空间分块内插。 局部方法： (1)泰森多边形 (2)距离反比权值插值(反距离权重内插) (3)样条函数内插技术 克里金内插方法 密度估算 线性内插、 双线性内插 多项式内插 多层曲面叠加法等。 局部内插 （1）Voronoi多边形 Voronoi多边形, 泰森多边形(Thiessen,又叫Dirichlet 或Voronoi多边形) 一种极端的边界内插方法,只用最近的单个点进行区域插值. 泰森多边形按数据点位置将区域分割成子区域,每个子区域包含一个数据点,各子区域到其内数据点的距离小于任何到其它数据点的距离,并用其内数据点进行赋值.连接所有数据点的连线形成Delaunay三角形,与不规则三角网TIN具有相同的拓扑结构. 在一个多边形内，每个未知点与该多边形内的已知点最接近，而与其他已知点更远. （1）Voronoi多边形（边界内插） 泰森多边形法的基本原理 未知点的最佳值由最邻近的观测值产生。 假定任何重要的变化都发生在区域的边界上 边界内的变化则是均匀的、同质的。 Voronoi多边形 Delaunay 三角网格 Delaunay 三角网的重要性质: