maple克里金插值工具箱

1. arcgis克里金插值数据太接近怎么回事

• 首先需要一组某一地区的shp文件，在这里，我们使用之前小编制作的黄河口的一套shp文件。具体详情请参见http://jingyan..com/article/656db918e5b316e381249cd9.html（用ArcGIS在完整地图shapefile文件截取局部区域）。

  

2. 克里金插值原理

克里格方法（Kriging）又称空间局部插值法，是以变异函数理论和结构分析为基础。

在有限区域内对区域化变量进行无偏最优估计的一种方法，是地统计学的主要内容之一。 南非矿产工程师D.R.Krige（1951年）在寻找金矿时首次运用这种方法，法国著名统计学 家G.Matheron随后将该方法理论化、系统化，并命名为Kriging，即克里格方法。

克里格方法的适用范围为区域化变量存在空间相关性，即如果变异函数和结构分析的 结果表明区域化变量存在空间相关性，则可以利用克里格方法进行内插或外推；否则，是不可行的。其实质是利用区域化变量的原始数据和变异函数的结构特点，对未知样点进行 线性无偏、最优估计。

无偏是指偏差的数学期望为0，最优是指估计值与实际值之差的平 方和最小。也就是说，克里格方法是根据未知样点有限邻域内的若干已知样本点数据，在 考虑了样本点的形状、大小和空间方位，与未知样点的相互空间位置关系，以及变异函数 提供的结构信息之后，对未知样点进行的一种线性无偏最优估计。

(2)maple克里金插值工具箱扩展阅读：

应用

克里金法被广泛用于各类观测的空间插值，例如地质学中的地下水位和土壤湿度的采样；环境科学研究中的大气污染（例如臭氧）和土壤污染物的研究；以及大气科学中的近地面风场 、气温、降水等的单点观测。

克里金法在工程问题的数值试验中可作为代理模型（surrogate model）对有限的模拟结果进行插值。具体而言，若对问题全局使用确定性模拟方法（deterministic computer simulations），例如有限元方法会占用大量计算资源而无法（快速）实现时，可以仅模拟局部个别点的结果并使用克里金法插值到全局

参考资料：网络-克里金法

3. arcgis怎么在指定的研究区范围内做克里金插值，并且把各个取值范围的面积统计出来。

打开克里金插值法的界面，设置参数设置界面右下角的“环境（environment）”，“处理范围（回processing extent）”的“捕捉栅格”答选择shp文件，“栅格分析（raster analysis）”中的“掩膜”选择选择shp文件。设置好，点确定即可。

4. matlab克里金插值工具箱的使用方式，请问有没有对matlab及克里金插值熟悉的大神可以帮帮忙。

这种网上一般都有现成的代码下载，下面这个链接就可以下载
http://www.pudn.com/downloads296/sourcecode/math/detail1326372.html
另外，matlab工具箱的话，先在回File->set path中添加所要加载的工具箱（答工具箱必须复制到matlab安装目录中的toolbox中），添加之后就可以像调用系统函数一样直接调用工具箱中的函数，最好先在网上找找工具箱使用的简单例程

5. 想利用克里金(Kinging)插值法来扩充自己采集的一些数据，Kinging.m代码应该怎么实现

Bars bar = new Bars();
bar.setId(rs.getLong("id"));
bar.setName(rs.getString("name"));
bar.setType(rs.getInt("type"));
bar.setCreatorId(rs.getLong("creator_id"));
resultList.add(bar);
if (currentNum == skipEnd - 1)
break;
}

6. 克里金插值

利用克里金法插值时变异函数的确定是其关键。当区域化变量不满足二阶平稳假设存在漂移时，漂移的形式、残余（Ｒｅｓｉｄｕａｌ）变异函数参数的估计比较困难。该文提出了利用多元逐步回归法确定漂移的次数；采用矩法和最大似然法相结合估计残余变异函数参数；当区域内数据点个数比较多时，在三角网格剖分过程中一次确定三角形与其内数据点的包含关系，用于快速检索待插点邻域内的数据点，最后给出了应用实例。
【作者单位】：北京航空航天大学机械工程与自动化学院!北京100083(牛文杰;朱大培);北京航空航天大学计算机科学与工程系(陈其明)
【关键词】：泛克里金法插值;最大似然法;矩法;三角网格剖分;多元逐步回归
【分类号】：TP391.41
景资料分析、实验变异函数的图形观察来确定的，这种确定方法一般带有一定的主观性。正如文献［１１所指出的那样，对于同一个区域化变量，有些人认为满足二阶平稳假设，而另一些人则认为带有漂移，没有一个判定准则。实际应用中，漂移次数的确定可借鉴利用Ｋ阶广义多项式协方

7. 克里金插值图能转化为矢量图么

在GP工具箱有克里金插值工具，其中有一个参数就是掩膜数据，可以是矢量或者栅格，插值后的数据就是在掩膜的范围；

8. 克里金插值算法

根据项目对数据处理的要求，采用了优化的克里金插值算法，将等值线地化数据插值转换为格网数据，以便实现地化数据的三维显示（王家华等，1999）。其主要实现过程如下：

第一步，计算半变异图，用非线性最小二乘拟合半变异函数系数；

第二步，数据点进行四叉树存储；

第三步，对每一格网点搜索邻近数据点；

第四步，由待预测网格点和邻近数据点计算克里金算法中系数矩阵，及右端常数向量；

第五步，对矩阵进行LU分解，回代求解待预测点的预测值。

克里金插值算法主要包括半变异函数和邻近点搜索的计算，实现方法如下。

（1）半变异函数计算

半变异函数是地质统计学中区域化变量理论的基础。地质统计学主要完成2方面的任务：利用半变异函数生成半变异图来量化研究对象的空间结构；通过插值方法利用半变异图中拟合模型和研究对象周围的实测值来对未知值进行预测。

半变异函数是用来描述区域化变量结构性和随机性并存这一空间特征而提出的。在满足假设的条件下，随机函数z（x）和z（x+h）为某一物理参数测定值的一一对应的2组函数，h为每对数之间的距离。半变异函数γ（h）可用下式来计算：

γ（h）＝ 1/2E｛［z（x）-z（x +h）］²｝

4种基本的半变异函数模式（除了这4种基本模式以外，还有很多模式），包括：

1）线形模式（Linear Model）

浙江省农业地质环境GIS设计与实现

2）球面模式（Spherical Model）

浙江省农业地质环境GIS设计与实现

3）指数模式（Exponential Model）

浙江省农业地质环境GIS设计与实现

4）高斯模式（Gaussian Model）

浙江省农业地质环境GIS设计与实现

半变异函数γ（h）会随距离h增大而增大，并逐渐逼近一定值（C₀ +C），称为基台值（Sill）；而逼近基台值所对应的距离，称为影响范围（Range），表示空间中两位置间的距离小于影响范围时，是空间相关性的。在线性和球面模式中，影响范围等于a；在指数和高斯模式中，影响范围则分别等于3a和

。而模式于半变异函数轴的截距（C₀）成为块金系数（Nugget Effect），产生的原因主要是样本测定的误差和最小采样间距内的变异。在应用上，为探讨说明空间变异在不同方向上的差距，也可利用非等向性的变异函数模式。半变异图拟合半变异函数模式的拟合方法可采用非线性最小二乘法拟合。

（2）邻近点搜索算法

由于矩阵LU分解求解方程的算法会随着矩阵维数的增加计算量增大，所以针对大量采样数据点时不能采用全部数据进行估计，必须采用插值点的临近点数据进行计算，即采用局部数据进行克里金算法进行计算。搜索邻近点可采用四叉树结构存储总数据，以提高搜索邻近点的速度。

对于选取邻近点的数目要有所限制，因该值的大小选择会影响插值的计算结果。若太大，则内插结果过于平滑；太小，则无法反映地表的变化；距离预测点较远的实测点可能与待估样点已经不存在自相关关系，也不能参与插值计算。采取以插值点为圆心，以R为半径的圆来确定取样的范围和参加计算的实测样点数目（如果存在各向异性，则可考虑划定一椭圆作为研究区域）。为了避免方向上的偏差，将圆平均地分为4个扇区，每个扇区内实测点数目在2～5之间，这样总共参与每个待估点预测的实测点数目平均达到8个。

区域内临近点的选择，存在着两种策略。

1）以邻近点的个数为基准。通常情况下，邻近点的个数以8～12个为宜，并且个数不能少于2个。此时计算出来的图像较为光滑。

2）以邻近点的半径尺度为基准。通常情况下，选择5～10 倍栅格间距的距离为宜。此时必须定义选择邻近点的最小和最大个数，当在一定半径内查找的邻近点个数小于最小个数时，应扩大搜索半径，使之达到最小查找个数；反之在一定半径内查找的邻近点个数大于最大个数时，应缩小搜索半径，使之小于最大查找个数。通常情况下最大最小个数分别可以定为20和4。

克里金算法的优点在于它基于一些可被验证的统计假设。根据这些假设，克里金算法产生的栅格节点估计量是最佳的，所有的估计量都依赖于可获得的观测值，并且平均误差最小。克里金算法提供了方差误差分析的表达式，可以表明每一个栅格节点的估计精度。

9. 求matlab克里金插值工具箱easykrig3.0，发送到[email protected]

题发布了2天了，才求？

有。。。不传

10. matlab中DACE工具箱进行克里金插值的使用方法

您好，克里金工具箱会用了吗？