最近我用R进行空间 - 时间插值 . 我在包 gstat 中选择函数krigeST，它在时空数据中使用克里金法 . 这是我的步骤：
- 首先，我创建一个名为stfdf的STFDF对象：
stfdf<-stConstruct(observation,space = list(values = 1:ncol(observation)),time = Node$Time,SpatialObj = loc,interval = T)
- 然后我用空间 - 时间数据stfdf计算经验变异函数：
var<-variogramST(values~1,data = stfdf,assumeRegular = T)
plot(var,wireframe=T)
这是经验变异函数的图像：empirical variogram
- 最后我需要将经验模型与理论模型相匹配，在这里我选择了可分离的模型：
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，"Exp"，8000,0）和time = vgm（0，"Exp"，2,0.001），但它无法拟合经验变异函数 .
有没有人对此有所了解？非常感谢！