R gstat时空变异函数克立格法

Question

我正在尝试使用R包gstat中的函数variogramST来计算时空变异函数。

有12年的数据，在空间和时间上的不规则点上有20,000个数据点(没有完全网格或部分网格)。我必须使用时空包中的STIDF来处理不规则的数据集。我想要一个时间半方差图，参考点在0，90,180,270天，直到几年等。不幸的是，计算和记忆问题都会发生。当命令

samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1)

在没有进一步参数的情况下运行时，半变差函数仅考虑了非常短的时间段，这似乎没有适当地捕获固有的数据结构。

这个函数有更多的参数可供用户使用，但我不确定如何正确地参数化它们: tlag，tunit，twindow。具体地说，我想知道他们是如何互动的，以及我如何实现上面描述的目标。因此，我尝试了以下代码

samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag= ...., tunit=... , twindow= ...)

由于我的32Gbyte RAM计算机中的内存问题，以下代码结果无法工作：

samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag=90*(0:20), tunit="days")

但可能是有缺陷的，否则。此外，就计算时间而言，后一行代码似乎也不可行。

有没有人知道如何从gstat包中正确指定variogramST function，以达到所需的时间间隔？

谢谢

Answer 1

如果我理解正确的话，twindow参数应该是计算时空变异函数时要包含的观测值的数量。假设你的20k点在这12年中大致均匀分布，那么你每年大约有1600点。同样，假设我理解正确，如果您希望在时间自相关计算中包含大约两年的数据，您将这样做：

samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag=90*(0:20), tunit="days",twindow=2*1600)