variogramST 如何处理距离和时间？答案

【问题标题】：How is variogramST dealing with distance and time?variogramST 如何处理距离和时间？
【发布时间】：2020-07-05 04:24:41
【问题描述】：

我是r 中spatio-temporal 分析的新手，所以我不知道我所做的是否正确。我有很多时空数据，我想绘制一个时空样本变异函数。

我的问题是我不太了解 variogramST (package: gstat) 的工作原理。 spacelag 是公里还是什么？这就是我现在所做的：

var_st <- variogramST (datast~1, data=DF, tunits="days", tlags=seq(0, 93, by= 31), 
                       boundaries=seq(0, 100, by=25), na.omit=T, assumeRegular=F)

我的第二个问题是：

它是否在精确距离 25、50、75、100 处绘制变异函数？例如，如果两个城市有距离，例如 12 是variogramST 考虑吗？或者它只是考虑精确距离为 25、50、75、100 的观察对？此外，我有每月的数据。我读到，在这种情况下，最好的选择是选择tunits="days"，而我选择tlags=seq(0, 93, by= 31)，对吗？不知道怎么处理。

非常感谢所有愿意回答我的人！

【问题讨论】：

没有参数tunits;是复制粘贴错误吗？参数tunit 只能在data 是STIDF 对象时使用。

标签： spatio-temporal r r data-visualization gstat covariogram

【解决方案1】：

空间距离单位取决于输入数据的坐标参考系，详见 sp::spDistsN1。
经验变异函数依赖于分箱 - 距离类。空间和时间上的观测按其近似空间距离分组（例如，0 公里到 25 公里；25 公里到 50 公里，...根据cutoff、width 和 boundaries）。及时，它还取决于您提供给 variogramST 的对象的类：
- 对于 STFDF，tlags 被评估为 STFDF 时隙的 id，在您的情况下，观察同时 id，具有 31、62 和 93 个 id 差异 - 很可能不是您想要的，假设DF时间槽中的每个id都指一个新的月份，你需要使用0:3。 tunit 的参数在此处无效。
- 如果数据对象是 STIDF，时间距离也会被分箱到分隔范围，其中 tlags 用作具有时间单元 tunit 的类之间的间隔。

有关将时空数据表示为 STFDF、STSDF 和 STIDF 的详细信息，请参阅 package spacetime。

HTH

【讨论】：