如何选择krige插值的时空变异函数参数

Question

最近我用R进行空间 - 时间插值。我在 gstat 包中选择函数 krigeST ，它在时空数据中使用克里金法。 这是我的步骤：

1. 首先，我创建一个名为stfdf的STFDF对象：

   stfdf<-stConstruct(observation,space = list(values = 1:ncol(observation)),time = Node$Time,SpatialObj = loc,interval = T)
   

2. 然后我用空间 - 时间数据stfdf计算经验变异函数：

    var<-variogramST(values~1,data = stfdf,assumeRegular = T)
    plot(var,wireframe=T)
   

这是经验变异函数的图像： empirical variogram

1. 最后我需要将经验模型与理论模型相匹配，在这里我选择了可​​分离的模型：

   sepvgm<-vgmST("separable",space = vgm(0.02, "Exp", 8000, 0),time = vgm(0, "Exp", 2, 0.001),sill = 0.02)
   sepvgmST<-fit.stVariogram(var,sepvgm,fit.method=6,method="L-BFGS-B")
   

现在，我的问题是如何在 vgm 函数中选择参数？在上面的示例中，我将其设置为 space = vgm（0.02，＆＃34; Exp＆＃34;，8000,0）和 time = vgm（0，＆＃34; Exp＆＃ 34;，2,0.001），但它不符合经验变异函数。

有人对此有所了解吗？非常感谢！