基于ArcGIS的批量克里金插值

前言：获取到地下水位点位数据，需要进行空间插值以获得空间连续的栅格数据。记录一下操作流程。

又前言：这里没有进行数据分布的检验，详细内容参见：

arcgis中克里金插值的操作（一） - 知乎

ArcGIS之克里金插值教学_arcgis克里金插值_树谷-胡老师的博客-CSDN博客

使用ArcGIS实现空间点的克里金（Kriging）插值详细步骤 - 程序员大本营

这里只记录用arcgis对月度数据批量克里金插值方法，便于以后使用。

准备数据集

原始数据集是整理好的每年每月地下水埋深数据（存在空值）。每个sheet为一年，每列包含1-12月。在处理时将每年的sheet单独保存为了csv文件存放在一个文件夹中。

 

 

通过ArcGIS中的模型工具，将csv文件批量导入，并转为shp文件

 

 迭代表，工作空间为存放csv文件的文件夹。导入xy坐标后要素转点，转为shp文件，命名为对应表的名称。

批量克里金插值

这里首先介绍使用的批量克里金工具箱，是从b站一个up主那免费下载的：

【空间的律动】批量插值工具箱Batch Interpolation v0.1.2使用说明_哔哩哔哩_bilibili

对于工具箱下载、应用在视频中都有介绍。这个工具箱主要是对一个表格中不同字段进行批量克里金插值。

groundwater是存放shp文件的文件夹，迭代要素为上一步生成的shp文件。同时为不同年份创建了.gdb文件（地理数据库），这里主要是因为每年1-12月字段都是相同的，所以分成不同gdb和文件夹保存，如果字段本身已经区分那么可以直接存放在一个地理数据库下面。【注：批量克里金的工具箱导出结果必须保存在gdb文件中，所以需要提前新建一个数据库】

这里创建文件夹主要是为了下一步复制栅格。

以上是批量克里金工具的界面，其中包括初始字段、结束字段、插值处理范围（可以选择边界）、插值网格大小（这里为了后续计算方便，设置为MODIS像元大小） 。

复制栅格，保存为tif文件

这里用的方法可能比较笨，因为上一步分成不同的gdb文件了，导致这里迭代不知道如何同时迭代gdb文件和其中的栅格（逐月）文件。所以直接手动一年一年进行了。。

这里需要修改第一个蓝色的圆圈，修改不同的gdb文件，以及复制栅格这个模块，将其中输出文件夹改成刚才按照对应年份生成的文件夹，然后生成逐月tif文件。 

后记：没有细看克里金插值的原理，也没有进行验证和比较。查阅资料发现在使用克里金插值之前应该检验数据是否正态分布，如果不是应该进行对数变换。通过ArcGIS里面的地统计分析也可以检查克里金插值的精度和适用性。这里简单记录后面插值流程，前面验证部分还要仔细看看原理之后再修改。