使用R包gstat的时空块克里金法吗?


9

在的文档中,sp我看到coordsSpatialPolygons函数的参数仅接受2列。无法创建3维多边形吗?我正在通过使用package gstat通过将时间添加为a的第三个空间维度来进行时空克里金法SpatialPointsDataFrame。我想将其扩展为以多边形为块的克里金法。有什么办法吗?

Answers:


5

spSpatialPoints*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不进行仿真-!)。


3

您确实需要对方法进行更多研究,并阅读文档以了解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")

0

时间将是第四个维度,z级别或高度将是第三个维度,锡或小数或z级别以gis为单位处理高度,时间可能由版本数据库表示。您可以将多边形映射到锡或与dem相交,或将z列与代表多边形形状点的xy列相关联。也许您还可以添加一个时间字段以及一系列值,并且值随时间变化,并以列的形式递增。当然,ESRI可能为此提供了一个工具。


我不认为Geostatistical Analyst支持时空克里金法。
Paul Hiemstra 2013年
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.