地理信息系统课程GIS空间插值要点.pptVIP

空间插值Spatial Interpolation 空间插值的概念 空间插值的类型 空间插值的方法 空间插值概念 空间插值——空间插值常用于将离散点的测量数据转换为连续的数据曲面，以便与其它空间现象的分布模式进行比较，它包括了空间内插和外推两种算法。空间内插算法：通过已知点的数据推求同一区域未知点数据。空间外推算法：通过已知区域的数据，推求其它区域数据。 例如：在一个没有数据记录的地点，其降水量可通过对附近气象站已知降水量记录的插值来估算出来。 进行空间插值要有两个基本条件：已知点和插值方法 控制点是已知数值的点，也称为已知点、样本点或观测点。 控制点提供了为空间插值建立插值方法的必要数据。 空间插值的一个基本假设是估算点的数值受到邻近控制点的影响比较远控制点的影响更大。 空间插值方法的分类 整体插值、局部插值和边界内插法； 确定性插值和地统计插值； 精确插值和近似插值。 1、整体插值、局部插值和边界内插法 整体插值 整体插值：用研究区所有采样点数据进行全区特征拟合。 整个区域的数据都会影响单个插值点，单个数据点变量值的增加、减少或者删除，都对整个区域有影响。 典型例子是：全局趋势面分析 某些区域化地质变量，如地层面的深度、地层厚度、地层水矿化度等，均可认为其分布在三维空间的某个曲面G 上。若已知G，则可根据它来研究这些地质变量在区域上的分布规律和局部特征. 实际的工作中无法得到准确的G ，但却可以根据已知的观测数据：Mi (xi , yi , zi) (i=1,2,…,n)，构造(拟合)一个近似于G 的数学曲面L。一般把这个拟合的曲面L称为趋势面。 趋势面分析就是在空间中已知点Mi (xi , yi , zi) 的控制下，拟合一个连续的数学曲面，并以此研究地质变量在区域上和局部范围内变化规律的一种统计方法。 拟合的数学曲面叫做趋势面。多项式和付立叶级数是趋势面分析常用的数学模型。 一、多项式曲面的一般形式 1、当数据为一维时 1）线性回归: 2）非线性：二次或高次多项式： 式中 z —变量； x , y —观测点的地理坐标。 若多项式中自变量的最高次数为k，则称这种多项式为k次多项式。多项式曲面的形态将随着k的增大而趋向复杂。 2、数据是二维的二元二次或高次多项式 图 1-3次多项式趋势面示意图 最小二乘法 对于变量间的相关关系，我们可以根据大量的统计资料，找出它们在数量变化方面的规律，这种统计规律所揭示的关系就是回归关系,所表示的数学方程就是回归方程。 图中的直线可表示为 根据上式，在确定α、β的情况下，给定一个x值，我们就能够得到一个确定的y值，然而根据上式得到的y值与实际的y值存在一个误差 如果我们以ｕ表示误差，则方程变为： α、β为回归系数 ｕ为随机误差项 使直线与各散点的距离的平方和最小 RSS= = 根据最小化的一阶条件，将式分别对x,y求偏导，并令其为零，即可求得α, β 一阶线性平面可模拟具有单一坡度的斜坡地形表面； 二次曲面方程可表达山头、洼地区域； 三次曲面则能描述较为复杂的地形曲面。 2 调和趋势面分析 调和趋势面分析的数学模型是正弦和余弦函数的组合，它具有波动的特征。 地层界面的波状起伏、沉积旋回、地球磁场的变化等都表现出不严格的周期性重复。因此,调和趋势面有时更有利于把具有波状特征的地质变量的趋势部分和异常分开,进而研究地质变量的波动特征。 多项式与调和趋势面的比较 由于以下缺点，在空间内插中整体内插并不常用： 采样点的增减或移动都需要对多项式的系数作全面调整，从而采样点之间会出现难以控制的振荡现象，致使函数极不稳定； 解算速度慢且对计算机容量要求较高。 不能提供内插区域的局部地形特征； 优点：整个区域上函数的唯一性、能得到全局光滑连续的空间曲面、分反映宏观地形特征等。 在空间内插中，一般是与局部内插方法配合使用，例如在使用局部内插方法前，利用整体内插去掉不符合总体趋势的宏观地物特征。整体内插函数常常用来揭示整个区域内的地形宏观起伏态势。 全局多项式插值法适用于在数据集中创建平滑表面及标识长期趋势。然而，在地球科学中，除了长期趋势之外，感兴趣的变量通常还具有短程变化。当数据集显示出短程变化时，局部多项式插值法地图可捕获这种变化。 将区域按一定的方法进行分块，对每一块根据地形曲面特征单独进行曲面拟合和高程内插，每一块都可用不同的曲面进行表达，称