1

我有一些车站数据，我正在尝试将车站数据重新网格化为网格数据。我已经尝试了该autoKrige功能，但出现错误。我想知道是否有人可以帮助我，谢谢。

测试数据可以从https://www.dropbox.com/s/68euncqtb3oijw5/testdata.rdata?dl=0下载

load("testdata.rdata",.GlobalEnv)
library(automap)
library(rgdal)

coordinates(newgrid) <- ~x+y
names(Mydata)[1:2]=c("x","y")
coordinates(Mydata) <- ~x+y

Krig=matrix(nrow = nrow(newgrid@coords),ncol = 30)
for( i in 1:30) {
  Krig[,i] = autoKrige(Mydata[,i],input_data=Mydata,new_data = newgrid)$krige_output$var1.pred
}

`Error: too many spatial dimensions: 328`
4

0 回答 0