空间插值算法汇总
- 格式:docx
- 大小:13.66 KB
- 文档页数:4
空间插值算法:1、距离倒数乘方法 (Inverse Distanee to a Power ) 距离倒数乘方格网化方法是一个加权平均插值法,可以进行确切的或者圆滑的方式插值。
方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降。
对于一个较大的方次,较近的数据点被给定一个较高的权重份额,对于一个较小的方次,权重比较均匀地分配给各数据点。
计算一个格网结点时给予一个特定数据点的权值与指定方次的从结点到观测点的该结点被赋予距离倒数成比例。
当计算一个格网结点时,配给的权重是一个分数,所有权重的总和等于1.0。
当一个观测点与一个格网结点重合时,该观测点被给予一个实际为1.0的权重,所有其它观测点被给予一个几乎为0.0的权重。
换言之,该结点被赋给与观测点一致的值。
这就是一个准确插值。
距离倒数法的特征之一是要在格网区域内产生围绕观测点位置的"牛眼"。
用距离倒数格网化时可以指定一个圆滑参数。
大于零的圆滑参数保证,对于一个特定的结点,没有哪个观测点被赋予全部的权值,即使观测点与该结点重合也是如此。
圆滑参数通过修匀已被插值的格网来降低"牛眼"影响。
2、克里金法 (Kriging)克里金法是一种在许多领域都很有用的地质统计格网化方法。
克里金法试图那样表示隐含在你的数据中的趋势,例如,高点会是沿一个脊连接,而不是被牛眼形等值线所孤立。
克里金法中包含了几个因子:变化图模型,漂移类型和矿块效应。
3、最小曲率法 (Minimum Curvature )最小曲率法广泛用于地球科学。
用最小曲率法生成的插值面类似于一个通过各个数据值的,具有最小弯曲量的长条形薄弹性片。
最小曲率法,试图在尽可能严格地尊重数据的同时生成尽可能圆滑的曲面。
使用最小曲率法时要涉及到两个参数:最大残差参数和最大循环次数参数来控制最小曲率的收敛标准4、多元回归法(Polynomial Regression )多元回归被用来确定你的数据的大规模的趋势和图案。
空间插值方法汇总Inverse Distance to a Power(反距离加权插值法)Kriging(克里金插值法)Minimum Curvature(最小曲率)Modified Shepard's Method(改进谢别德法)Natural Neighbor(自然邻点插值法)Nearest Neighbor(最近邻点插值法)Polynomial Regression(多元回归法)Radial Basis Function(径向基函数法)Triangulation with Linear Interpolation(线性插值三角网法)Moving Average(移动平均法)Local Polynomial(局部多项式法)1、距离倒数乘方法距离倒数乘方格网化方法是一个加权平均插值法,可以进行确切的或者圆滑的方式插值。
方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降。
对于一个较大的方次,较近的数据点被给定一个较高的权重份额,对于一个较小的方次,权重比较均匀地分配给各数据点。
计算一个格网结点时给予一个特定数据点的权值与指定方次的从结点到观测点的该结点被赋予距离倒数成比例。
当计算一个格网结点时,配给的权重是一个分数,所有权重的总和等于1.0。
当一个观测点与一个格网结点重合时,该观测点被给予一个实际为 1.0 的权重,所有其它观测点被给予一个几乎为 0.0 的权重。
换言之,该结点被赋给与观测点一致的值。
这就是一个准确插值。
距离倒数法的特征之一是要在格网区域内产生围绕观测点位置的"牛眼"。
用距离倒数格网化时可以指定一个圆滑参数。
大于零的圆滑参数保证,对于一个特定的结点,没有哪个观测点被赋予全部的权值,即使观测点与该结点重合也是如此。
圆滑参数通过修匀已被插值的格网来降低"牛眼"影响。
2、克里金法克里金法是一种在许多领域都很有用的地质统计格网化方法。
前段时间要对气象要素进行插值,翻看了多种方法,做了个PPT报告.对每个方法有简单的介绍极一些总结,不一定都是个人看法,参考了多方书面(sufer,ArcGIS应用教程)以及坛子里,百度上等搜到的资料的看后笔记,有些注了出处有些忘了.截图共享下,也不知有用没用.有错的地方请跟贴指正,谢谢啦!--------------------------------所谓空间数据插值,即通过探寻收集到的样点/样方数据的规律,外推/内插到整个研究区域为面数据的方法.即根据已知区域的数据求算待估区域值, 影响插值精度的主要因素就是插值法的选取空间数据插值方法的基本原理:任何一种空间数据插值法都是基于空间相关性的基础上进行的。
即空间位置上越靠近,则事物或现象就越相似, 空间位置越远,则越相异或者越不相关,体现了事物/现象对空间位置的依赖关系。
(/dky/nb/page/2000-3-3/2000332117262480.htm,南京师范大学地理科学学院地理信息系统专业网络课程教程)➢由于经典统计建模通常要求因变量是纯随机独立变量,而空间插值则要求插值变量具备某种程度的空间自相关性的具随机性和结构性的区域化变量。
即区域内部是随机的,与位置无关的,而在整体的空间分布上又是有一定的规律可循的,这也是不宜用简单的统计分析方法进行插值预估的原因。
从而空间统计学应用而生。
➢无论用哪种插值方法,根据统计学假设可知,样本点越多越好,而样本的分布越均匀越好。
常用的空间数据插值方法之一:趋势面分析⏹趋势面分析(Trend analyst)。
严格来说趋势面分析并不是在一种空间数据插值法。
它是根据采样点的地理坐标X,Y值与样点的属性Z值建立多元回归模型,前提假设是,Z值是独立变量且呈正态分布,其回归误差与位置无关。
⏹根据自行设置的参数可建立线性、二次…或n次多项式回归模型,从而得到不同的拟合平面,可以是平面,亦可以是曲面。
精度以最小二乘法进行验证。
一:距离加权反比法插值算法1:原理:设空间待插点为P(Xp,Yp,Zp),P点邻域内有已知散乱点Q i(x i,y i,z i),i=1,2,3….n;利用距离加权反比法对P点的属性值Zp进行插值。
其插值原理是待插点的属性值是待插点邻域内已知散乱点属性值的加权平均, 权的大小与待插点与邻域内散乱点之间的距离有关, 是距离的k(0<=k<=2)(k一般取2)次方的倒数。
其中:d i为待插点与其邻域内第i个点之间的距离。
2:克里金算法设研究区域为A, 区域化变量即欲研究的物理属性变量为{Z(x)∈A},x 表示空间位置(一维、二维或三维坐标), 在采样点x i(i=1,2,…n)处的属性值(或称为区域化变量的一次实现)为Z(x i)(i=1,2,…n),则根据普通克里金插值原理, 未采样点x0处的属性值Z(x0)估计值是n个已知采样点属性值的加权和, 即;λi为待求权系数。
假设区域化变量Z(x)在整个研究区域内满足二阶平稳假设:(1):Z(x)的数学期望存在且等于常数:E[Z(x)]=m(常数)(2):Z(x)的协方差Cov(x i,x j)存在,且只与两点之间的相对位置有关。
或满足本征假设(3)E[Z(x i)-Z(x j)]=0.(4)增量的方差存在且平稳:Var[Z(x i)-Z(x j)]= E[Z(x i)-Z(x j)]2经过无偏性要求:经推到可得:。
在无偏条件下使得估计方差达到最小,即其中:u 是拉格朗日算子。
可以得到求解权系数λi (i=1,2…n )的方程组:求出诸权系数λi (i=1,2…n )后,就可求出位采样点x 0处的属性值Z *( x 0).上述求解λi (i=1,2…n )的方程中的Cov (x i ,x j )若用变异函数表示时,其形式为:变异函数的定义为:由克里金插值所得的方差为:或。
测绘技术中的数据空间插值方法测绘技术是一门以测量、制图为基础的工程技术学科。
而数据空间插值方法则是测绘技术中的一项重要技术,它在测绘领域中起到了至关重要的作用。
本文将介绍测绘技术中的数据空间插值方法及其应用。
数据空间插值方法是指通过已知的离散点数据,通过某种方式推测出未知位置的数据。
它在测绘领域中被广泛应用于地形建模、地貌分析、图像处理等诸多方面。
常见的数据空间插值方法包括反距离加权法、克里金插值法、三角网剖分法等。
首先,反距离加权法是一种基于距离的插值方法,其原理是根据待插值点与已知点之间的距离,将距离较近的已知点的属性值加权求和,从而得到待插值点的属性值。
反距离加权法简单且易于理解,但对数据点的分布和点密度要求较高,且容易受到噪声的影响。
其次,克里金插值法是一种基于空间自相关性的插值方法,它通过测定半变异函数来建立数据点之间的关联性模型。
根据已知点之间的相互影响关系,克里金插值法可以对未知位置的属性值进行预测。
克里金插值法在测绘领域中被广泛应用于地质勘探、土地利用评价等方面,并且具有良好的预测精度。
最后,三角网剖分法是一种基于三角网格的插值方法,它将给定的离散点数据通过连接相邻点构成三角网格,从而形成一个连续的表面模型。
三角网剖分法具有较高的计算效率和较好的插值效果,广泛应用于地形建模和地表分析等领域。
除了以上介绍的几种常见的数据空间插值方法外,还有许多其他方法也被应用于测绘技术中。
例如径向基函数插值法、样条插值法、多层神经网络插值法等。
这些方法各有特点,适用于不同的场景和需求。
数据空间插值方法在测绘技术中具有重要的应用价值。
它可以通过对已知数据的合理处理,获得缺失数据的预测值,从而填补数据空间中的空缺。
这对于测绘领域中的地理数据分析、地质勘探、灾害预测等具有重要意义。
然而,数据空间插值方法在应用过程中也存在一些问题和挑战。
例如,当插值点周围的已知点密度不均匀时,插值结果可能出现不准确的情况。
空间插值算法:
1、距离倒数乘方法(Inverse Distance to a Power)距离倒数乘方格网化方法是一个加权平均插值法,可以进行确切的或者圆滑的方式插值。
方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降。
对于一个较大的方次,较近的数据点被给定一个较高的权重份额,对于一个较小的方次,权重比较均匀地分配给各数据点。
计算一个格网结点时给予一个特定数据点的权值与指定方次的从结点到观测点的该结点被赋予距离倒数成比例。
当计算一个格网结点时,配给的权重是一个分数,所有权重的总和等于1.0。
当一个观测点与一个格网结点重合时,该观测点被给予一个实际为 1.0 的权重,所有其它观测点被给予一个几乎为0.0 的权重。
换言之,该结点被赋给与观测点一致的值。
这就是一个准确插值。
距离倒数法的特征之一是要在格网区域内产生围绕观测点位置的"牛眼"。
用距离倒数格网化时可以指定一个圆滑参数。
大于零的圆滑参数保证,对于一个特定的结点,没有哪个观测点被赋予全部的权值,即使观测点与该结点重合也是如此。
圆滑参数通过修匀已被插值的格网来降低"牛眼"影响。
2、克里金法(Kriging)克里金法是一种在许多领域都很有用的地质统计格网化方法。
克里金法试图那样表示隐含在你的数据中的趋势,例如,高点会是沿一个脊连接,而不是被牛眼形等值线所孤立。
克里金法中包含了几个因子:变化图模型,漂移类型和矿块效应。
3、最小曲率法(Minimum Curvature)最小曲率法广泛用于地球科学。
用最小曲率法生成的插值面类似于一个通过各个数据值的,具有最小弯曲量的长条形薄弹性片。
最小曲率法,试图在尽可能严格地尊重数据的同时,生成尽可能圆滑的曲面。
使用最小曲率法时要涉及到两个参数:最大残差参数和最大循环次数参数来控制最小曲率的收敛标准。
4、多元回归法(Polynomial Regression)多元回归被用来确定你的数据的大规模的趋势和图案。
你可以用几个选项来确定你需要的趋势面类型。
多元
回归实际上不是插值器,因为它并不试图预测未知的Z 值。
它实际上是一
个趋势面分析作图程序。
使用多元回归法时要涉及到曲面定义和指定XY
的最高方次设置,曲面定义是选择采用的数据的多项式类型,这些类型分别是简单平面、双线性鞍、二次曲面、三次曲面和用户定义的多项式。
参数设置是指定多项式方程中X 和Y组元的最高方次。
5、径向基本函数法(Radial Basis Function)径向基本函数法是多个数据插值方法的组合。
根据适应你的数据和生成一个圆滑曲面的能力,其中的复二次函数被许多人认为是最好的方法。
所有径向基本函数法都是准确的插值器,它们都要为尊重你的数据而努力。
为了试图生成一个更圆滑的曲面,对所有这些方法你都可以引入一个圆滑系数。
你可以指定的函数类似于克里金中的变化图。
当对一个格网结点插值时,这些个函数给数据点规定了一套最佳权重。
6、谢别德法(Shepard's Method)谢别德法使用距离倒数加权的最小二乘方的方法。
因此,它与距离倒数乘方插值器相似,但它利用了局部最小二乘方来消除或减少所生成等值线的"牛眼"外观。
谢别德法可以是一个准确或圆滑插值器。
在用谢别德法作为格网化方法时要涉及到圆滑参数的设置。
圆滑参数是使谢别德法能够象一个圆滑插值器那样工作。
当你增加圆滑参数的值时,圆滑的效果越好。
7、三角网/线形插值法(Triangulation with Linear Interpolation)三角网插值器是一种严密的插值器,它的工作路线与手工绘制等值线相近。
这种方法是通过在数据点之间连线以建立起若干个三角形来工作的。
原始数据点的连结方法是这样:所有三角形的边都不能与另外的三角形相交。
其结果构成了一张覆盖格网范围的,由三角形拼接起来的网。
每一个三角形定义了一个覆盖该三角形内格网结点的面。
三角形的倾斜和标高由定义这个三角形的三个原始数据点确定。
给定三角形内的全部结点都要受到该三角形的表面的限制。
因为原始数据点被用来定义各个三角形,所以你的数据是很受到尊重的。
8.自然邻点插值法(Natural Neighbor)自然邻点插值法(NaturalNeighbor)是Surfer7.0才有的网格化新方法。
自然邻点插值法广泛应用于一些研究领域中。
其基本原理是对于一组泰森(Thiessen)多边形,当在数据集中加入一个新的数
据点(目标)时,就会修改这些泰森多边形,而使用邻点的权重平均值将决定待
插点的权重,待插点的权重和目标泰森多边形成比例[9]。
实际上,在这些多边形中,有一些多边形的尺寸将缩小,并且没有一个多边形的大小会增加。
同时,自然邻点插值法在数据点凸起的位置并不外推等值线(如泰森多边形的轮廓线)。
9.最近邻点插值法最近邻点插值法(NearestNeighbor)又称泰森多边形方法,泰森多边形(Thiesen,又叫Dirichlet或Voronoi多边形)分析法是荷兰气象学家A.H.Thiessen提出的一种分析方法。
最初用于从离散分布气象站的降雨量数
据中计算平均降雨量,现在GIS和地理分析中经常采用泰森多边形进行快速
的赋值[2]。
实际上,最近邻点插值的一个隐含的假设条件是任一网格点p(x,y)的属性值都使用距它最近的位置点的属性值,用每一个网格节点的最邻点值
作为待的节点值[3]。
当数据已经是均匀间隔分布,要先将数据转换为SURFER 的网格文件,可以应用最近邻点插值法;或者在一个文件中,数据紧密完整,只
有少数点没有取值,可用最近邻点插值法来填充无值的数据点。
有时需要排
除网格文件中的无值数据的区域,在搜索椭圆(SearchEllipse)设置一个值,对无
数据区域赋予该网格文件里的空白值。
设置的搜索半径的大小要小于该网格文件数据值之间的距离,所有的无数据网格节点都被赋予空白值。
在使用最
近邻点插值网格化法,将一个规则间隔的XYZ数据转换为一个网格文件时,可
设置网格间隔和XYZ数据的数据点之间的间距相等。
最近邻点插值网格化法没有选项,它是均质且无变化的,对均匀间隔的数据进行插值很有用,同时,它
对填充无值数据的区域很有效。
10.Moving Average(移动平均法)移动平均法是用一组最近的实际数据值来预测未来一期或几期内公司产品的需求量、公司产能等的一种常用方法。
移动平均法适用于即期预测。
当产品需求既不快速增长也不快速下降,且不存
在季节性因素时,移动平均法能有效地消除预测中的随机波动,是非常有用的。
移动平均法根据预测时使用的各元素的权重不同移动平均法是一种简单平滑预测技术,它的基本思想是:根据时间序列资料、逐项推移,依次计算包含一定项数的序时平均值,以反映长期趋势的方法。
因此,当时间序列的数值由于受周期变动和随机波动的影响,起伏较大,不易显示出事件的发展趋势时,使用移动平均法可以消除这些因素的影响,显示出事件的发展方向与趋势(即趋势线),然后依趋势线分析预测序列的长期趋势。
11.Local Polynomial(局部多项式法)
12.Modified Shepard's Method(改进谢别德法)。