autokrige を使用して降雨データを補間したい automap ライブラリを使用しています。結果として得られる補間されたマップには、予測ポイント マップとは異なる次元 (行と列の数が異なる) があります。これはどのように可能ですか？以下の私のコードを見てください：
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)