Question

我正在使用自动化库，我想使用autokrige来插入降雨数据。生成的插值映射具有与我的预测点映射不同的维度（不同的行数和列数）。这怎么可能？请参阅下面的代码：

library(automap)
library(raster)
library(rgdal)

r <- read.table("c:\\Active\\2013-001_Qfever\\model\\data\\output_kriging\\raindata_old1.txt", header = TRUE, skip = 2, sep="")
coordinates(r) = ~x+y

predmap <- readGDAL(fname="c:\\Active\\2013-001_Qfever\\model\\NLclone_scalar.map")
kriging_result <- autoKrige(rain~1, r, predmap)

plot(kriging_result)

xyz <- as.data.frame(kriging_result$krige_output)                                                 
head(xyz)
p = xyz[,c('x','y','var1.pred')]

output <- rasterFromXYZ(p, res=c(250,250), crs="+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs", digits=5)

writeRaster(output, "c:\\Active\\2013-001_Qfever\\model\\data\\output_kriging\\test.asc", format="ascii")
plot(output)

Answer 1

在代码中替换以下行，它将起作用 -

p = xyz[,c('x1','x2','var1.pred')]

autokriging导致不同的维度

1 个答案: