OrdinaryKriging 克里金插值

x,y,d = 一维数组（散点数据）
x0 = np.linspace(x.min, x.max, 150)
y0 = np.linspace(y.min, y.max, 100)
x1, y1 = np.meshgrid(x0, y0)

for i in ['linear', 'power', 'gaussian', 'spherical', 'exponential', 'hole-effect']:
    Krin = OrdinaryKriging(x, y, d, variogram_model=i)
    data, ssl = Krin.execute('grid', x0, y0)
    plt.contourf(x1, y1, np.array(data), cmap=cmap)
    plt.title(i)
    plt.show()