python实现克里金插值_python 普通克里金（Kriging）法的实现

克里金法时一种用于空间插值的地学统计方法。

克里金法用半变异测定空间要素，要素即自相关要素。

半变异公式为：

其中γ(h) 是已知点 xi 和 xj 的半变异，***h***表示这两个点之间的距离，z是属性值。

假设不存在漂移，普通克里金法重点考虑空间相关因素，并用拟合的半变异直接进行插值。

估算某测量点z值的通用方程为：

式中，z0是待估计值，zx是已知点x的值，Wx是每个已知点关联的权重，s是用于估计的已知点数目。

权重可以由一组矩阵方程得到。

此程序对半变异进行拟合时采用的时最简单的正比例函数拟合

数据为csv格式

保存格式如下：

第一行为第一个点以此类推

最后一行是待求点坐标，其中z为未知值，暂且假设为0

代码如下：

import numpy as np

from math import*

from numpy.linalg import *

h_data=np.loadtxt(open("高程点数据.csv"),delimiter=",",skiprows=0)