作者:
andrew43 (討厭有好心推文後刪文者)
2018-06-06 11:08:45我不是很明白 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檔 讀取多筆資料
作者: wayturn (岳) 2018-06-06 19:33:00
如果要讀多檔,可先存檔名再迴圈讀檔,關鍵字:file.list