我正在比较R中使用gstat库的3D数据集的不同方差图模型的交叉验证(“省一分”)结果。我也想比较相似的交叉验证结果以进行逆距离估计,但是我无法看到如何在gstat中使用IDW进行各向异性估算(或交叉验证)。 gstat估计中的IDW是否仅限于各向同性估计,如果不是,可以应用什么语法?
答案 0 :(得分:0)
假设无法在gstat中进行各向异性IDW估计,我将数据位置转换为各向同性,然后使用krige.cv进行了交叉验证。
从一个带有X,Y,Z,HM值的csv文件开始(X,Y和Z转换为使各向同性的点-在这种情况下,XY旋转15度,然后X坐标乘以3.333,Z值乘以30:“ hm_assays_iso.csv”)。我使用的R会话如下:
>库(gstat)
>库(sp)
> hm_iso <-read.csv(“ hm_assays_iso.csv”)
>坐标(hm_iso)<-c(“ X”,“ Y”,“ Z”)
> hm_idw.cv <-krige.cv(HM〜1,location = hm_iso,set = list(idp = 3))
> write.csv(hm_idw.cv,file =“ hm_cv_idw3.csv”)