无法使用 writeGDAL 函数创建光栅文件

问题描述 投票:0回答:0

我正在尝试使用“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

任何人都可以帮我找出导致此错误的原因以及如何解决它吗?

提前感谢您的协助。

r raster rgdal spatial-interpolation
© www.soinside.com 2019 - 2024. All rights reserved.