本地克里格法使用完整数据（`STFDF`）报告`krigeST`中的错误

我一直在使用gstat包中的krigeST函数，并尝试使用全套基于空间多边形的数据STFDF进行局部kriging

当我尝试以下方法时

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中，sp中带有空间多边形的某些STFDF被强制为SFIDF。其中，有一个子集步骤是SpatialPolygons类不允许的，因为其中存在重复项from@index[，1]

我想尝试不可分离的时空变异函数，所以不使用局部克里格不是一个真正的选择。这个问题有解决办法吗

---

谢谢。

克里格法通常被认为是从观测点到非观测点的插值。为什么你要把多边形传递给函数？我的数据实际上是基于区域的，所以要非常严格，我可能会考虑块时空Kriging。然而，从计算上看，这似乎要昂贵得多。也许我应该把多边形和它们的几何中心放在一起，这叫做面对点或面对面克里格法。我不建议将多边形替换为其中心点。