克里金插值是一种基于变异函数理论和结构分析的线性、无偏插值,在地理信息和气象学方面有着广泛的应用,如进行空间高程的插值、气温的插值,克里金插值的原理和过程可以参考(
http://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/3d-analyst-toolbox/how-kriging-works.htm
) ,下面以插值高程为例演示ArcGIS克里金插值的详细过程(不足或错误之处还望不吝指出)
一 、打开ArcMap,导入离散点坐标
原始数据如下图所示
打开ArcMap,导入上图离散点数据并显示(不清楚如何导入的可以百度),如下图所示
二、打开克里金插值对话框
右键菜单栏空白处打开地学统计分析Geostatical Analyst
.打开地学统计向导
选择“地学统计方法”中的克里金法,
注意右侧值域的选择,此处为对高程Z插值,故选择高程Z
!
三、克里金的详细设置
克里金又分为简单克里金、普通克里金和泛克里金(常用语不连续区域)等,对于连续的区域,选择普通克里金即可,如下图所示
点击下一步
点击下一步进入变异函数选择的部分,如下图所示
在窗口右侧模型中选择合适的变异函数模型,常见的有球函数,指数函数等,使得窗口左边的曲线与蓝点实现最佳拟合。
四、结果分析
插值结果如下图所示,渲染图表示等高区域分布,符合实际情况。如果想要获取离散点的克里金插值,可更改窗口右侧的X、Y坐标值,Z的预测值即可更新显示。
五、精度分析
克里金插值相比于其他插值方法的一重要优点就是可以给出插值点的插值精度,如下图窗口左侧所示,使用内部检核方法,对插值精度进行分析;窗口右侧通过插值均方根、标准均方根等给出总体插值精度。如下图所示:
克里金插值是一种基于变异函数理论和结构分析的线性、无偏插值,在地理信息和气象学方面有着广泛的应用,如进行空间高程的插值、气温的插值,克里金插值的原理和过程可以参考(http://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/3d-analyst-toolbox/how-kriging-works.htm) ,下面以插值高程为例演示ArcGIS克里金插值的详细...
h为样本
点
之间的距离;n为由h分开的成对样本
点
的数量;z为
点
的属性值(高程或其他属性值)。
计算半方差时步骤如下:
(1)求所有样本
点
之间的距离,共有n(n-1)/2个不同的距离;
(2)对所有距离从小到大排序并分为n组,计算
(3)计算n组距离每组的平均距离,将平均距离代入半方差公式中,计算出每组距离所对应的实验变差值;
(4)绘制距离-半方差散
点
图,求出拟合系数,即对应的参数。
2.在不同距离的半方差值都计算出来后,
今天是
克里金插值
最后一讲,主要内容为
插值
步骤,和每种方法在创建表面过程中需要考虑的重要参数。
克里金插值
创建预测表面的核心参数有:数据变换、趋势、半变异函数/协方差模型和搜索领域(均在之前的章节进行了详细讲解)。
由于每种方法步骤均为加载数据源和属性字段,设置对应的方法参数,通过执行下一步,最终生成所需
插值
结果,因此,本文以普通
克里金
的预测表面为例进行演示,其他方法以此为参考。
步骤一:打开地统计...
详细步骤
:
一、准备数据
在excel中准备需要进行
克里金插值
的内容,应该包括至少三列:X坐标,Y坐标,数据。在每一列的列首做好标记,方便导入ArcMap后对应。注意存储时存为xls文件。
二、将数据导如导ArcMap
点
击下面这个标记,把数据文件添加到ArcMap。
如果是初次
使用
的化,可能会不知道进入添加数据对话框之后如何操作。这时应该先
点
击红框标记的图标,连接到文件夹。
学过
空间
插值
的人都知道和反距离
插值
(IDW)和
克里金插值
, 本文讲简单介绍基本
克里金插值
的原理,以及在
Arcgis
中
实现
的详细过程。由于IDW操作和
克里金
很相似,并且最常用的是
克里金
,因此实操部分给了
克里金
的。读者可以根据如下教程摸索IDW是完全可以的呢。
文章目录一、反距离
插值
(IDW)二、
克里金插值
法三、
Arcgis
实际操作3.1
插值
实操 - 以
克里金
为例四、出图
一、反距离
插值
(IDW)...
克里金
法是通过一组具有 z 值的分散
点
生成估计表面的高级地统计过程。与
插值
工具集中的其他
插值
方法不同,选择用于生成输出表面的最佳估算方法之前,有效
使用
克里金
法工具涉及 z 值表示的现象的
空间
行为的交互研究。
什么是
克里金
法?
IDW(反距离加权法)和样条函数法插 值工具被称为确定性
插值
方法,因为这些方法直接基于周围的测量值或确定生成表面的平滑度的指定数学公式。第二类
插值
方法由地统计方法(如...
IDW反距离权重
插值
介绍
反距离权重 (IDW)
插值
:彼此距离较近的事物要比彼此距离较远的事物更相似。当为任何未测量的位置预测值时,反距离权重法会采用预测位置周围的测量值。与距离预测位置较远的测量值相比,距离预测位置最近的测量值对预测值的影响更大。反距离权重法假定每个测量
点
都有一种局部影响,而这种影响会随着距离的增大而减小。由于这种方法为距离预测位置最近的
点
分配的权重较大,而权重却作为距离的函数而减小,因此称之为反距离权重法。
反距离权重计算过程:
ArcGIS
实现
.
1. 准备好需要进行
插值
的多个
点
数据图层,并确保它们都具有相同的投影和坐标系。
2. 在ArcMap中打开一个新的地图文档,然后将所有需要进行
插值
的
点
数据图层添加到地图中。
3.
点
击“Geoprocessing”菜单栏下的“ArcToolbox”按钮,打开ArcToolbox窗口。
4. 在ArcToolbox窗口中,展开“Spatial Analyst Tools”文件夹,然后双击“Interpolation”文件夹。
5. 选择“
Kri
gin
g”工具,双击打开“
Kri
gin
g”对话框。
6. 在“
Kri
gin
g”对话框中,设置所需的参数,如
克里金插值
方法、输出栅格分辨率、输出栅格范围等。
7.
点
击“Environments”按钮,设置环境变量,如输出坐标系、输出文件格式等。
8.
点
击“Batch”按钮,选择需要进行批量
克里金插值
的
点
数据图层,然后
点
击“OK”按钮。
9.
点
击“OK”按钮开始批量
克里金插值
,等待处理完成。
注意事项:在批量
克里金插值
之前,需要确保每个
点
数据图层都具有相同的字段结构和数据类型,否则会出现
插值
结果不一致的情况。另外,批量
克里金插值
过程中需要耗费大量计算资源,因此需要根据计算机配置和数据量合理设置
插值
参数,以免导致计算时间过长或计算失败。