我正在尝试使用“gstat”和“栅格”包对 R 中的降水数据执行 IDW 插值。但是,我无法使用“rgdal”包中的“writeGDAL”函数从我的插值数据创建栅格文件。
这是我使用的代码:
#Load libraries
library(sp)
library(gstat)
library(rgdal)
#Load data
data <- read.csv("20140516.csv")
#Convert to spatial points
coordinates(data) <- c("Long", "Lat")
#Create a grid
grid <- readOGR("1000_grid.shp")
#Convert CRS
proj4string(data)\<-proj4string(grid)
data <- spTransform(data, crs(grid))
#Perform IDW interpolation
idw <- gstat::idw(formula = Precipitation \~ 1, locations = data, newdata = grid)
#Write raster to file
writeGDAL(idw, "precipitation_idw.tif")
当我尝试运行代码时,出现以下错误消息:
Error in create2GDAL(dataset = dataset, drivername = drivername, type = type, : gridded(dataset) is not TRUE
任何人都可以帮我找出导致此错误的原因以及如何解决它吗?
提前感谢您的协助。