空间插值IDW.pptx

空间插值空间插值是用已知点的数值来估算其它点的数值的过程例如：在一个没有数据记录的地点，其降水量可通过对附近气象站已知降水量记录的插值来估算出来。为什么插值为栅格？在GIS应用中主要用于估算出栅格中每个象元的值。因此空间插值是将点数据转换成面数据的一种方法，目的是使点数据也能用于空间分析和建模。空间插值方法的分类 整体拟合法 局部拟合法 确定性 随机性 确定性 随机性 趋势面（非精确） 回归（非精确） 泰森（精确）、密度估算（非精确）、反距离权重（精确）、薄板样条（精确） 克里金（精确） 空间插值主要方法 空间插值常用于将离散点的测量数据转换为连续的数据曲面，它包括内插和外推两种算法。前者是通过已知点的数据计算同一区域内其他未知点的数据，后者则是通过已知区域的数据，求未知区域的数据。主要的内插方法有：反距离加权（Inverse Distance Weighted）全局多项式（Global Polynomial Interpolation）全局多项式（Local Polynomial Interpolation）径向基函数（Radial Basis Funtions）克里格内插（ Kriging ） 空间插值的理论假设是：空间位置上越靠近的点，越可能具有相似的特征值，而距离越远的点，其特征值相似的可能性越小。空间插值方法正是依据该假设设计的，分为整体插值方法和部分插值方法两类。 整体插值：用研究区域所有采样点的数据进行全区域特征拟合，如边界内插法、趋势面分析等。 部分插值：仅仅用邻近的数据点来估计未知点的值，如最邻近点法（泰森多边形方法）、移动平均插值方法（距离倒数插值法）、样条函数插值方法、空间自协方差最佳插值方法（克里金插值）等。距离反比法(Inverse Distance) 距离反比插值方法最早由 Shepard 提出(Richard Franke,1982)提出的，并逐步得到发展。每个采样对插值结果的影响随距离增加而减弱，因此距目标点近的样点赋予的权重较大。ARCGIS——IDW控制反距离加权的参数—权重权重过高，较近点的影响较大，拟合表面更细致（不光滑）；权重过低，较远点的影响增加，拟合表面更光滑。缺省值常为 2 。控制反距离加权的参数—搜索半径搜索半径－固定 对固定型半径，搜索距离一定，所有在该半径内的样点参与计算。 可预先设定一个阈值，当给定半径内搜索到的点小于该值时可扩大搜索半径，直到达到该阈值为止。搜索半径类型－可变 设定参与计算的样点数是固定的，则搜索的半径是可变的。这样对每个插值点的搜索半径可能都不同，因为要达到规定的点数所需要搜索的区域是不一样的。控制反距离加权的参数—障碍设置可利用一线状和面状数据集来限制样点的搜索。线状数据集可作为平坦地表的悬崖或脊状障碍物：只有位于同侧的样点才符合要求。 Power = 2, search = 230Power = 2, search = 150Power = 4, search = 600Power = 2, search = 600权重系数和搜索半径的影响图示距离反比插值评价优点——简便易行；可为变量值变化很大的数据集提供一个合理的插值结果；不会出现无意义的插值结果而无法解释。不足——对权重函数的选择十分敏感；易受数据点集群的影响，结果常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式；全局最大和最小变量值都散布于数据之中。距离反比很少有预测的特点，内插得到的插值点数据在样点数据取值范围内。 空间插值接口IInterpolationOp Interface空间插值接口(IDW)IInterpolationOp.IDW MethodIDW实现主窗体的目录中添加空间插值目录两个子目录：命名为反距离插值IDW，克里金差值KrigingIDW实现新建WinForm，命名为：IDW两个下拉菜单控件，两个textbox一个ImageButton，选择图标为文件夹打开，如没有图标，将其text属性设为“浏览”二字也可两个Button：ok和closeIDW实现反距离插值IDW的单击事件下实现如下代码： IDW pIDW = new IDW(); pIDW.pMap = axMapControl1.Map; //pKriging.ShowDialog(); pIDW.Visible = false; DialogResult result = pIDW.ShowDialog(); if (result == DialogResult.Cancel) return; //ColorRampRaster(pKriging.pRasterLayer, 9); axMapControl1.AddLayer(pIDW.pRasterLayer); axMapControl1.Act