使用完整数据(`STFDF`)在krigeST`中报告本地克里金错误

时间:2015-10-29 15:30:13

标签: r kriging gstat

我一直在使用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]中有重复的术语。

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

感谢。

0 个答案:

没有答案