[問題] 月平均資料

[AndrewShi (沒有妳的我)]

[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)

[軟體熟悉度]:
入門(寫過其他程式，只是對語法不熟悉)

[問題敘述]:
各位大大好，
我放入的這筆資料是2010年全球每天的降雨(量)資料，現在我想把每日的降雨量計算成月
平均.年平均降雨量，下面我所想到的迴圈是可以畫得出圖來，但畫出來感覺不太正確，所以想請
教大大們我的迴圈是否有問題，能否給我一些提點，謝謝。
p.s:原本的資料型態中降雨值的維度只包含經度和緯度(2維)，所以我用rbind把時間的維
度也併到降雨值裡。

[程式範例]:
rm(list=ls())
library(ncdf4)

TRMM_data <- "C:\\Users\\TOM\\Desktop\\R(資料庫)\\TRMM資料\\trmm_2010.nc"
obs <- nc_open(TRMM_data)
print(obs)

lon <- ncvar_get(obs, "lon")
lat <- ncvar_get(obs, "lat")
time <- ncvar_get(obs, "time")
precip <- ncvar_get(obs,"r")
time <- matrix(seq(as.Date("2010-01-01"), as.Date("2010-12-31"),1))
rbind(dim(time),precip[[3]])

time <- c()
for(time in seq_along(1:31)){
  mean(precip)
}

time <- c()
for(time in seq_along(1:365)){
  mean(precip)
}

lon <- lon-180
#lat <- rev(lat)
precip <- precip[,,time]
library(RColorBrewer)
image(lon,lat,precip,col=rev(brewer.pal(10,"RdBu")))
library(maptools)
gpclibPermit()
data(wrld_simpl)
plot(wrld_simpl,add=TRUE)

[環境敘述]:

[關鍵字]:
月平均   nc檔   降雨

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 118.169.229.203
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1540199720.A.401.html
→ andrew43: 你的每個loop都沒有用到time這個變數，所以就只是同一件
→ andrew43: 事重覆做很多次，看你設定的time有幾個。
→ andrew43: 建議你提供一個小nc檔，不然連資料格式都不清楚沒辦法答
→ andrew43: 不會自製nc檔的話，至少把lon lat time precip都說清楚
→ AndrewShi: dim(lon(經度))=1440
→ AndrewShi: dim(lat(緯度))=400
→ AndrewShi: dim(time)=365
→ AndrewShi: (原本資料型態)dim(precip(降雨量))=1440 400
→ AndrewShi: (rbind之後))dim(precip(降雨量))=1440 400 365
→ AndrewShi: 以我的理解是原本的資料型態為一天一張全球地圖(平面)
→ AndrewShi: 上面每個網格點有一個降雨值，總共365個平面(天)，而
→ AndrewShi: 我現在想要取出比如其中一個月的降雨量來做平均(全球)
→ AndrewShi: ，andrew大是指這個資料資訊嗎?!
推 andrew43: 對的。
推 andrew43: 1440和400這與時間的對應關係又如何需要知道
推 andrew43: 上句不對。我想說的是結合降雨量和時間後資料是如何對
→ andrew43: 應。他是array嗎？
推 andrew43: 1440和400個地點有重複嗎？照理說降雨量有三個維度，經
→ andrew43: 、維、日。但你說得有不全然如此。老實說我還是沒全明
→ andrew43: 白。
→ AndrewShi: 地點沒有重複，結合後降雨的維度就變成像你說的有三個
→ AndrewShi: 維度(經.緯.時間)，對應的關係我就不是很清楚了，只是
→ AndrewShi: 合併後dim(time)=NULL了，還是我把檔案放上來?!
→ AndrewShi: 合併後降雨是array，時間我上面是把它設為matrix，也可
→ AndrewShi: 以改成array，但我不知道合併後的時間這樣設定有何差
→ AndrewShi: 異。
→ AndrewShi: 檔案我放在這裡，https://drive.google.com/folders/1a
→ AndrewShi: l7SnSxJzllkN9H9QjL8vMv3wF-WSpv5?usp=sharing
→ andrew43: 網址不對，縮一下
→ AndrewShi: https://drive.google.com/drive/folders/1aI7SnSxJzII
→ AndrewShi: kN9H9QjL8vMv3wF-WSpv5，抱歉andrew大，我不知道怎麼縮
→ AndrewShi: 網址XD
→ AndrewShi: http://0rz.tw/JI056