空间插值.ppt

空间插值 贺欣欣 空间插值是用已知点的数值来估算其它点的数值的过程 例如：在一个没有数据记录的地点，其降水量可通过对附近气象站已知降水量记录的插值来估算出来。 在GIS应用中主要用于估算出栅格中每个象元的值。因此空间插值是将点数据转换成面数据的一种方法，目的是使点数据也能用于空间分析和建模。 为什么插值为栅格？ 插值之所以可称为一种可行的方案，是因为我们假设，空间分布对象都是空间相关的，也就是说，彼此接近的对象往往具有相似的特征。例如，如果街道的一边正在下雨，那么便可以充分肯定地预测街道的另一边也在下雨。而对于是否整个城镇都在下雨，则无法确定，同时也无法确定邻县的天气情况。 通过上述类比很容易发现，接近采样点的点值相对于距离采样点较远的点，与采样点相似的可能性更大。这是插值的基础。点插值的典型应用是通过一组采样测量值创建一个高程面。 一、空间插值的要素 进行空间插值要有两个基本条件：已知点和插值方法 1 控制点 控制点是已知数值的点，也称为已知点、样本点或观测点。 控制点提供了为空间插值建立插值方法的必要数据。 空间插值的一个基本假设是估算点的数值受到邻近控制点的影响比较远控制点的影响更大。 二、空间插值的类型 空间插值有多种分类方法 第一，它可以分为全局和局部拟合法。 全局插值法利用现有的每个已知点来估算未知点的值。 而局部插值法则是用已知点的样本来估算位置点的值。 这两种方法的区别就是用于估算的控制点数目不一样 其次空间插值方法还分为精确和非精确插值法。 对某个数值已知的点，精确插值法在该点位置的估算值与该点已知值相同。换句话说，精确插值所生成的面通过所有控制点。 相反，非精确插值或叫做近似插值，估算的点值与该点已知值不同。 第三，空间插值方法可以分成确定性和随机性两种。 确定性插值方法不提供预测值的误差检验。 随机性插值方法则用估计变异提供预测误差的评价 空间插值方法的分类 趋势面模型 作为一种非精确的插值方法，趋势面分析用多项式方程拟合已知点，并用于估算其它值的点。线性或一阶趋势面用如下方程： zx,y=b0+b1x+b2y 式中：属性值y是x和y的函数。系数b由控制点估算 在概念上，趋势插值法类似于取一张纸将其插入各凸起点之间（凸起到一定高度）。下图展示的是从平缓山丘采集而来的一组高程采样点。使用的纸张为洋红色 图一五个已知值的气象站，这五个值又围绕着未知值的0号站。下表显示各个站点的x、y坐标（以象元大小为2000m的栅格的行列测量）及其已知值 实例 本例说明如何用方程（或线性趋势面）对未知点的0号站点进行插值。最小二乘法通常用于计算方程中的待定系数b0、b1和b2。 因此，第一步是建立如下三个正规方程。 ∑z=b0n+b1 ∑x+b2 ∑y ∑xz=b0 ∑x+b1 ∑x2+b2 ∑xy ∑yz=b0 ∑y+b1 ∑xy+b2 ∑y2 以上方程可以改写为如下矩阵： n ∑x ∑y b0 ∑z ∑x ∑x2 ∑xy × b1 = ∑xz ∑y ∑xy ∑y2 b2 ∑yz 用已知五个点的值，我们就能计算出统计值并将其带入方程： 5 377 318 b0 67.270 29007 23862 × b1 = 5043.650 318 23862 20714 b2 4445.800 将左边第一个矩阵的逆矩阵与右边的矩阵相乘，便能算出系数b -10.094 0.020 0.347 0号站点的未知值可用这些系数由下式估算 Z0=-10.094+(0.020)(69)+(0.347)(67) =14.535 大部分自然现象的分布通常比有一阶趋势面生成的倾斜面更复杂。因而需要更高阶的趋势面模型来拟合更复杂的表面 例如，包含山和谷的三阶面。该模型基于下方程 zx,y=b0+b1x+b2y+b3x2+b4xy+b5y2+b6x3+