问题描述
我有超过 60 个用于特定坐标的 nc 文件,其中仅包含一个变量。我想在 csv 文件中从它们中提取变量 (tasmax)。 一个文件的代码是:
ncdata <- nc_open("data.nc")
tasmax<-ncvar_get(ncdata,"tasmax")
nc_close(ncdata)
tasmax<-as.data.frame(tasmax)
write.csv2(tasmax,file="data.csv")
我试图对所有文件进行循环,但总是出错..有什么想法吗?
解决方法
我遵循此代码并工作:
df=NULL
files<- list.files('tasmax/ssp370',pattern='*.nc',full.names=TRUE)
for(i in seq_along(files)) {
nc = nc_open(files[i])code
tasmax<-ncvar_get(nc,"tasmax")
nc_close(nc)
cbind(df,tasmax)->df
}
df<-as.data.frame(df)
write.csv2(df,file="data.csv")