普通克吕金插值

使用说明

普通克里金法是使用最普通和广泛的克里金方法。该方法假定用于插值的字段值的期望(平均值)未知且恒定。

  • 普通克里金法使用的数据应符合数据变化成正态分布的前提假设。
  • 普通克里金插值最大的特色不仅是提供一个最小估计误差的预测值,并且可明确的指出误差值的大小。
  • 普通克里金法采用两种方式来获取参与插值的采样点,进而获得相应位置点的预测值:一个是在待计算预测值位置点周围一定范围内,获取该范围内的所有采样点,通过特定的插值计算公式获得该位置点的预测值;另一个是在待计算预测值位置点周围获取一定数目的采样点,通过特定的插值计算公式获得该位置点的预测值。

打开“ExerciseData/RasterAnalysis”文件夹下的“Precipitation”数据源,其中有部分地区气象监测站点的降水量数据,我们用此数据来做示例。

功能入口

  • 空间分析 选项卡-> 栅格分析 组-> 插值分析 -> 普通克吕金
  • 工具箱 -> 栅格分析 -> 插值分析 -> 普通克吕金 。(iDesktopX)

参数说明

  • 设置插值分析的公共参数,包括源数据、插值范围、结果数据和环境设置。源数据、插值范围和结果数据等公共参数的设置请参见:公共参数说明
  • 设置样本点查找方式:支持变长查找、定长查找和块查找三种方式。关于这三种查找方式的详细介绍请参见:关于插值

    • 变长查找 :表示使用最大半径范围内的固定数目的样点值进行插值。

      • 最大半径 :输入用于变长查找的半径大小。默认值为0,表示使用最大半径查找。
      • 查找点数 :输入用于变长查找的点数目。默认点数为12。
    • 定长查找 :查找半径范围内所有的点都要参与插值运算。

      • 查找半径 :输入设定查找半径大小。默认查找半径为参与插值分析的数据集的范围的长或者宽的较大值的1/5。所有该半径范围内的采样点都要参加插值运算。
      • 最小点数 :输入用于变长查找的最少数目点。默认点数为5。当邻域中的点数小于所指定的最小值时,查找半径将不断增大,直到可以包含输入的最小点数为止。最大值为12。
    • 块查找 :根据设置的“块内最多点数”对数据集进行分块,然后使用块内的点进行插值运算。

      • 最多参与点数 :输入最多参与插值点数。默认最多参与插值点数为200。为了避免在插值时出现裂缝区,实际计算使用的插值块会在每个分块区域的基础上再均匀向外扩张,“最多参与点数”决定了块区域向外扩张的大小。一般此数值应大于设置的“块内最多点数”。
      • 块内最多点数 :输入每个块内的点的最多数量。默认块内最多点数为50。若块内点数多于此值,则继续分块;否则停止分块。“最多参与插值点数”与“块内最多点数”的参数的设置会直接影响快查找的性能。这两个值设置的越大,查找花费的时间会越久,因此建议用户在设置快查找的参数时设置比较合理的参数。
  • 其他参数 :包括半变异函数、旋转角度、平均值、基台值、自相关阈值、块金效应值等。

    • 半变异函数: 支持球函数、指数函数和高斯函数三种半变异函数。使用哪个模型需要根据数据的空间自相关性和数据现象的先验知识来决定。默认使用球函数。
    • 旋转角度: 每个查找邻域相对于水平方向逆时针旋转的角度。默认为0度。块查找不支持旋转角度设置。
    • 基台值: 半变异函数达到的顶点值,即在距离(X 轴)为0时,半变异函数与 Y 轴相交的值。默认为0。
    • 自相关阈值: 半变异函数到达基台值处的距离,即 X 轴相应的值。默认值为0。
    • 块金效应值: 在 h=0(X轴)时,半变异函数与 Y 轴相交的值。默认为0。

关于基台值、自相关阈值和块金效应三个参数的关系请参见:克吕金插值

OrdinaryKrigingResult