1

最近我在用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)

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

  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,"Exp",8000,0) 和 time=vgm(0,"Exp",2,0.001) ，但它无法拟合经验变异函数。

有人对此有想法吗？非常感谢！

4

0 回答 0