我一直在使用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]
中有重复的术语。
我想尝试不可分离的时空变异函数,所以不使用局部克里金法并不是一个真正的选择。这个问题有解决方法吗?
感谢。