我一直在使用包中的krigeST
函数gstat
,并尝试使用全套基于 SpatialPolygons 的数据 (STFDF) 进行本地克里金法。
当我尝试以下
spatio_time_krige_sumProd = krigeST(conso~1,
some_STFDF_with_SpatialPolygons_in_sp,
newdata = some_STF_with_SpatialPolygons_in_sp,
fitted_prodSum_variogram,
nmax = 50,
stAni = 4)
Error in from@sp[from@index[, 1], ] :
SpatialPolygons selection: can't find plot order if polygons are replicated
问题似乎是在 krigeST.local 中,some_STFDF_with_SpatialPolygons_in_sp
被强制转换为SFIDF
. 在那里,有一个班级不允许的子集步骤SpatialPolygons
,因为 中有重复的术语from@index[, 1]
。
我想尝试不可分离的时空变异函数,因此不使用局部克里金法并不是一个真正的选择。有解决此问题的方法吗?
谢谢。