如何在R中读取.grd文件为栅格？

Question

我想读取并堆叠几个.grd文件并将其导出为R中的.nc文件。我使用以下命令

library(raster)
library(ncdf4)
library(RNetCDF)
library(rgdal)

files <-list.files(path="G:/Gridded data/", 
           pattern="GRD", all.files=FALSE, full.names=TRUE,recursive=TRUE)
s <- stack(files)
rstack <- raster(files[1])

但它出现以下错误

Error in .local(.Object, ...) : 

Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer",  : 
  Cannot create a RasterLayer object from this file.

有人可以帮忙吗？提前谢谢

Answer 1

我会尝试pattern="\\.grd$"（这意味着，文件以＆＃34; .grd＆＃34;结尾）。或者至少pattern="grd"代替pattern="GRD"。 （或使用ignore.case = TRUE）

files <- list.files(path="G:/Gridded data/",  pattern="\\.grd$", full.names=TRUE,recursive=TRUE)

如果这不是问题，那么请向我们展示files的内容，以便能够看到正在发生的事情。