是否可以在R中使用gstat进行各向异性IDW估算(和交叉验证)?

时间:2019-11-15 10:10:53

标签: gstat

我正在比较R中使用gstat库的3D数据集的不同方差图模型的交叉验证(“省一分”)结果。我也想比较相似的交叉验证结果以进行逆距离估计,但是我无法看到如何在gstat中使用IDW进行各向异性估算(或交叉验证)。 gstat估计中的IDW是否仅限于各向同性估计,如果不是,可以应用什么语法?

1 个答案:

答案 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”)