Re: [問題] 兩層for迴圈問題

[andrew43 (討厭有好心推文後刪文者)]

我不是很明白 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檔   讀取多筆資料

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 60.248.222.1
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1528254529.A.9B0.html
推 AndrewShi: 謝謝andrew大~我會再好好研究的，不過我不太懂str_pad(
→ AndrewShi: 2,pad="0")這行是什麼意思??
→ andrew43: 補上前綴0。多練習自己看說明手冊 ?str_pad
推 AndrewShi: 另外想請問dt$filename裡沒有i這個變數名稱，最後在寫f
→ AndrewShi: or迴圈讀資料進來的時候應該沒辦法讀吧?!(實際操作後
→ AndrewShi: 的確出現error，沒辦法讀進來)
→ andrew43: 我的for loop裡i就是一次一次的dt$filename沒錯
→ andrew43: 你要怎麼讓for loop運作可自行調整
推 wayturn: 如果要讀多檔,可先存檔名再迴圈讀檔,關鍵字:file.list
推 AndrewShi: 好的，我再研究看看，謝謝你們~
推 AndrewShi: andrew大~我改成這樣但好像也只有讀到最後一筆資料，
→ AndrewShi: 能否請你再次指點我改如何改才能完整的把48筆資料都讀
→ AndrewShi: 進來呢??謝謝。
→ AndrewShi: https://i.imgur.com/SOGpYIo.jpg
→ andrew43: nc.res 後面要雙重中跨號
→ andrew43: 他是一個list。
→ andrew43: 還有問題可以回文並附上可重複問題的code，推文太久會
→ andrew43: 漏看。
推 AndrewShi: 了解，我解決了，忘記print也要放在for迴圈裡，謝謝and
→ AndrewShi: rew大~另外可以請教andrew大，gl指令(我有看說明)和寫f
→ AndrewShi: or迴圈的意思(重複動作)是差不多的嗎?!如果是差不多的
→ AndrewShi: ，能否順便想請問你我最一開始的code可以經過修改而達
→ AndrewShi: 到同樣的功用嗎??謝謝。
→ andrew43: 你看過說明手冊嗎？  輸入 ?gl
推 AndrewShi: 有~你們整理的指令表裡面也有，只是感覺它的功能跟for
→ AndrewShi: 迴圈很像，不太清楚它在意義上有什麼差別。
→ andrew43: gl和for風馬牛不相及。
→ andrew43: 若沒明白gl是什麼意思，看手冊推敲gl(3,2,12)是什麼結果
推 AndrewShi: 以三個分類變數，每個執行2次，列出前12筆?!
→ andrew43: 差不多了。產生一個因子，有3類，每次2個重覆，共12個
推 AndrewShi: 了解，謝謝andrew大~