Answers:
在sp
,SpatialPoints*
,SpatialPixels*
和SpatialGrid*
(与*
省略或替换DataFrame
)做支持多于2种空间尺寸,如OP已经完成,但是SpatialPolygons*
并SpatialLines*
没有。有了gstat
你可以做3 d块克里格与3-d块(使用block = c(10,10,10)
),但不能用于非矩形块这样做,因为OP希望。可以用时间代替第三维,这是完全可以的,但您只能使用公制ST变异函数。
library(gstat)
vignette("st")
为方差图模型提供了更多选择,但没有为块平均值提供预测(这是仅供参考,不是问题的答案)。
该问题的唯一答案是进行3D条件模拟,并在任意3D(2D多边形+时间范围)块上聚合点值。乏味,但可能;也仅沿3D路径,而不沿ST小插图中描述的路径(krigeST
不进行仿真-!)。
您确实需要对方法进行更多研究,并阅读文档以了解sp S4类对象的结构以及sp对象与相关gstat函数的交互。在sp Vignette中,详细解释了SpatialPolygons(仅多边形拓扑)对象和SpatialPolygonDataFrame(具有属性的多边形)对象之间的区别。
您要解释的不是块克里金法,并且将时间用作属性不会导致时空估计。您描述的概念方法论是完全无效的。使用多边形或多边形质心违反了Kriging关于均匀随机场,各向异性和非平稳性的假设。
这是一个使用时空包接口的时空模型上的gstat插图。我还应注意,constrainedKriging软件包可以使用非平稳均值函数和各向同性的弱平稳变异函数对任意形状的块进行块Kriging。
也就是说,要回答您的问题,可以将sp SpatialPointsDataFrame对象直接传递给gstat中的变异函数/克里格模型。在这种类型的sp对象中,属性驻留在“数据”插槽中,并且已经通过内部S4类结构附加到坐标。
# COERCE meuse DATAFRAME TO sp SpatialPointsDataFrame OBJECT
require(gstat)
data(meuse)
coordinates(meuse) <- ~ x + y
head(meuse@data)
# CREATE SEMIVARIOGRAM USING THE zinc ATTRIBUTE
# NOTE: THERE IS NO ARGUMENT FOR A "4th DIM"
v <- variogram(log(zinc) ~ 1, meuse)
plot(v, type = "l")