R：适用于R中的协同克里格的交叉协方差图不同位置的变量

如何解决R：适用于R中的协同克里格的交叉协方差图不同位置的变量

我想用两个不在同一位置采样的变量进行协同克里金法。因此，我无法使用交叉变异函数，而需要使用交叉 co 变异函数。使用软件包 gstat ，我使用函数variogram和选项covariogram=TRUE计算了交叉协方差图，但是如何使这些函数适合我的共区域化模型呢？以及如何使用交叉协方差图进行协同克里格法？

以下是使用meuse数据集的示例：

library(sp)
library(gstat)
data(meuse)
data(meuse.grid)

# subsetting the meuse dataset to have data at different locations:
Cd.data <- meuse[sample.int(nrow(meuse),0.7*nrow(meuse)),c("x","y","cadmium")]
Zn.data <- meuse[sample.int(nrow(meuse),"zinc")]
meuse.subset <- merge(Cd.data,Zn.data,by = c("x","y"),all = TRUE)

# Plot cross covariogram
g = gstat(NULL,id="Cd",formula = cadmium~1,data = meuse,locations = ~x+y)
g <- gstat(g,id = "Zn",form = zinc ~ 1,data=meuse,locations = ~x+y)
Cd.cov <- variogram(g,cutoff = 1500,width = 100,covariogram = TRUE)
plot(Cd.cov)

# fitting and cokriging (do not kNow how to do this)

使用gstat的解决方案会很好，但是如果您知道另一个可以做到这一点的软件包，将会很有帮助。

非常感谢您的时间和帮助