如何在R中读取.grd栅格文件并堆叠导出为.nc文件？

我来帮你搞定这个读取GRD文件的问题～你遇到的错误是因为R没法从目标文件创建RasterLayer，咱们一步步排查解决：

错误原因排查 #

首先先理清可能导致问题的几个关键点：

• GRD配套文件缺失：R原生的GRD格式需要成对存在.grd（头文件）和.gri（数据文件），如果其中一个丢失，raster包就无法识别这个文件。
• 文件名匹配不准确：你用pattern="GRD"会匹配所有文件名里包含GRD的文件，可能混进了非GRD格式的文件，导致读取失败。
• 单个文件损坏/路径问题：你测试raster(files[1])出错，说明第一个文件本身就无法被读取，大概率是路径错了、文件损坏，或者不是标准的R GRD格式。

分步解决步骤 #

1. 先验证单个文件的可读性 #

先手动指定一个确定的GRD文件路径，测试能不能正常读取：

library(raster)
# 替换成你实际的单个GRD文件完整路径
test_raster <- raster("G:/Gridded data/具体文件名.grd")
print(test_raster)

如果这一步报错，优先检查：

• 该文件的同目录下是否有同名的.gri文件；
• 文件是否损坏（可以尝试重新获取或生成该GRD文件）；
• 路径是否正确，Windows路径里记得用/或者转义的\\。

2. 修正文件列表的获取 #

调整list.files的参数，精确匹配.grd后缀的文件，避免混入无关文件：

files <- list.files(
  path = "G:/Gridded data/",
  pattern = "\\.grd$",  # 精确匹配以.grd结尾的文件
  all.files = FALSE,
  full.names = TRUE,
  recursive = TRUE,
  ignore.case = TRUE  # 如果文件名有大写GRD，加上这个参数兼容大小写
)
# 打印文件列表，确认都是正确的GRD文件
print(files)

3. 正确堆叠并导出为NetCDF文件 #

当确认所有文件都能正常读取后，执行堆叠和导出操作：

# 堆叠所有GRD文件
s <- stack(files)
# 导出为NetCDF文件，format指定为CDF即可
writeRaster(
  s,
  filename = "G:/Gridded data/stacked_output.nc",
  format = "CDF",
  overwrite = TRUE,  # 如果目标文件已存在，允许覆盖
  varname = "your_variable",  # 自定义变量名
  longname = "Your Variable Description"  # 变量的详细描述
)

注意：writeRaster导出NetCDF时，不需要额外加载ncdf4或RNetCDF包，raster包本身就能处理。

额外提示 #

如果你的GRD文件不是R原生格式（比如其他软件生成的GRD），可以尝试用rgdal包转换读取：

library(rgdal)
test_raster <- raster(readGDAL("G:/Gridded data/你的文件.grd"))

内容的提问来源于stack exchange，提问作者UseR10085