读取多个netcdf文件

Question

我正在尝试读取r中的多个nc4文件。下面是我用来执行此任务的代码;

library(ncdf4)
OSND_gpmr.df<-NULL

GPM_R.files= list.files(path,pattern='*.nc4',full.names=TRUE) 

for(i in seq_along(GPM_R.files)) {
  nc_data = nc_open(GPM_R.files[i])

  GPM_Prec<-ncvar_get(nc_data, 'IRprecipitation')

  x=dim(GPM_Prec)

  ### note start=c(42,28) are the index in image regards to real coordinates of interset
  ## R reads images from lat,long. 
  OSND_gpmr.spec =ncvar_get(nc_data, 'IRprecipitation', start = c(42,28), count = c(1,1))
  rbind(OSND_gpmr.df,data.frame(OSND_gpmr.spec))->OSND_gpmr.df 

  nc_close(nc_data)
}

但我一直都会遇到这个错误：

  

R_nc4_open出错：没有这样的文件或目录。

但是文件列表被正确识别为chr [1：1440]，如全局环境 - 值所示。

有人可以帮助我解决我做错的事吗？

Answer 1

您的工作目录可能与文件位置不同。您的 GPM_R.files 列表仅存储来自给定位置的文件名，而没有文件路径。而 nc_open() 需要具有完整路径的文件名。