我不是很明白 nc 档你想怎么处理,
所以这篇只简单说明一次读入多个 nc 档怎么把它们“储存在同一个 list 里面”。
假设有 2 * 24 个 nc 档案,档名叫
01-00.nc 01-01.nc 01-02.nc ... 01-23.nc 02-00.nc ... 02-23.nc
library(magrittr)
library(stringr)
library(ncdf4)
dt <-
data.frame(
day = gl(2, 24, labels = (1:2) %>% str_pad(2, pad = "0")),
hour = gl(24, 1, labels = (0:23) %>% str_pad(2, pad = "0"))
)
dt$filename <- paste0(dt$day, "-", dt$hour, ".nc")
nc.res <- vector("list", nrow(dt))
names(nc.res) <- dt$filename
for(i in dt$filename){
nc.res[[i]] <- nc_open(i)
}
则 nc.res 为所求,而 dt 可方便对照档名与日期时间。
之后对 nc.res 的呼叫就和一般 list 一样了,
例如 nc.res[[1]] 或 nc.res$`01-00.nc` 之类的。
你也可能需要用到 lapply() 方便操作 list。
※ 引述《AndrewShi (没有妳的我)》之铭言:
: [问题类型]:
: 程式咨询(我想用R 做某件事情,但是我不知道要怎么用R 写出来)
: [软件熟悉度]:
: 入门(写过其他程式,只是对语法不熟悉)
: [问题叙述]:
: 我想要读多笔名称相近的时间资料(nc档)进R,资料的时间长度为一小时一笔,而我想要
: 读2天共48笔资料进R,因此需写两层for循环(分别为日期和小时),我的程式码如下,虽然
: 过程中没有出现error,但读进来的资料日期和小时是交错的(环境叙述中的图,只有24笔
: )所以想请问该如何解决??谢谢。
: 资料名称唯一的差别在于#precip那行的"01"(日期)和"00"(小时)。
: [程式范例]:
: rm(list=ls())
: library(ncdf4)
: day=c("01","02")
: hour=c("00","01","02","03","04","05","06","07","08","09","10",11:23)
: class(day)
: class(hour)
: nc<-list()
: for(i in c(1:2)){
: for(j in c(0:23)){
: precip.file<-paste("C:\\Users\\TOM\\Desktop\\R(数据库
: )\\201505\\gsmap_mvk.201505",paste(day[i=(1:2)]),paste("."),
: paste(hour[j=(0:23)],paste("00.v6.4133.0.nc")),sep="")
: #precip<-"C:\\Users\\TOM\\Desktop\\R(数据库
: )\\20150501\\gsmap_mvk.201505"01"."00"00.v6.4133.0.nc"
: #nc[i,j]=rbind(i,j)
: }
: i<-i+1
: nc[i,j]<-nc_open(precip.file)
: }
: precip.file
: View(matrix(precip.file))
: [环境叙述]:
: 执行画面拍照截图于推文中
: [关键字]:
: 两层for循环 nc档 读取多笔资料