循环变量 IDW gstat R

Question

有没有办法使用 gstat 包的 IDW 插值在 R 中的 SpatialPointsDataFrame 对象内循环变量？

我的意思是这样的

for (day in list_days) {
  P.idw <- gstat::idw(day~1, P, newdata=grd, idp=2.0)
  r       <- raster(P.idw)
  plot(r)
}

编辑：因为那个结构给我错误：

Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  : 
  too many spatial dimensions: 58
In addition: Warning message:
In predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :
NAs introduced by coercion

但是如果我直接写变量名就没问题了:

P.idw <- gstat::idw(X2018.01.14~1, P, newdata=grd, idp=2.0)

Answer 1

而不是 day~1 你可以 as.formula(paste(day, "~ 1"))

循环变量 IDW gstat R

Loop variable IDW gstat R

r

gstat