[问题类型]:
程式咨询(我想用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档 降雨