[問題類型]:
程式諮詢(我想用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檔 降雨

https://drive.google.com/drive/folders/1aI7SnSxJzIIkN9H9QjL8vMv3wF-WSpv5，抱歉andrew大，我不知道怎麼縮網址XDhttp://0rz.tw/JI056

網址不對，縮一下

地點沒有重複，結合後降雨的維度就變成像你說的有三個維度(經.緯.時間)，對應的關係我就不是很清楚了，只是合併後dim(time)=NULL了，還是我把檔案放上來?!合併後降雨是array，時間我上面是把它設為matrix，也可以改成array，但我不知道合併後的時間這樣設定有何差異。檔案我放在這裡，https://drive.google.com/folders/1al7SnSxJzllkN9H9QjL8vMv3wF-WSpv5?usp=sharing

對的。1440和400這與時間的對應關係又如何需要知道上句不對。我想說的是結合降雨量和時間後資料是如何對應。他是array嗎？1440和400個地點有重複嗎？照理說降雨量有三個維度，經、維、日。但你說得有不全然如此。老實說我還是沒全明白。

dim(lon(經度))=1440dim(lat(緯度))=400dim(time)=365(原本資料型態)dim(precip(降雨量))=1440 400(rbind之後))dim(precip(降雨量))=1440 400 365以我的理解是原本的資料型態為一天一張全球地圖(平面)上面每個網格點有一個降雨值，總共365個平面(天)，而我現在想要取出比如其中一個月的降雨量來做平均(全球)，andrew大是指這個資料資訊嗎?!

你的每個loop都沒有用到time這個變數，所以就只是同一件事重覆做很多次，看你設定的time有幾個。建議你提供一個小nc檔，不然連資料格式都不清楚沒辦法答不會自製nc檔的話，至少把lon lat time precip都說清楚