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