作者AndrewShi (没有你的我)
看板R_Language
标题[问题] 月平均资料
时间Mon Oct 22 17:15:18 2018
[问题类型]:
程式谘询(我想用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://webptt.com/cn.aspx?n=bbs/R_Language/M.1540199720.A.401.html
1F：→ andrew43: 你的每个loop都没有用到time这个变数，所以就只是同一件 10/22 18:27
2F：→ andrew43: 事重覆做很多次，看你设定的time有几个。 10/22 18:28
3F：→ andrew43: 建议你提供一个小nc档，不然连资料格式都不清楚没办法答 10/22 18:29
4F：→ andrew43: 不会自制nc档的话，至少把lon lat time precip都说清楚 10/22 18:30
5F：→ AndrewShi: dim(lon(经度))=1440 10/22 18:44
6F：→ AndrewShi: dim(lat(纬度))=400 10/22 18:44
7F：→ AndrewShi: dim(time)=365 10/22 18:44
8F：→ AndrewShi: (原本资料型态)dim(precip(降雨量))=1440 400 10/22 18:44
9F：→ AndrewShi: (rbind之後))dim(precip(降雨量))=1440 400 365 10/22 18:44
10F：→ AndrewShi: 以我的理解是原本的资料型态为一天一张全球地图(平面) 10/22 18:44
11F：→ AndrewShi: 上面每个网格点有一个降雨值，总共365个平面(天)，而 10/22 18:44
12F：→ AndrewShi: 我现在想要取出比如其中一个月的降雨量来做平均(全球) 10/22 18:44
13F：→ AndrewShi: ，andrew大是指这个资料资讯吗?! 10/22 18:44
14F：推 andrew43: 对的。 10/22 18:46
15F：推 andrew43: 1440和400这与时间的对应关系又如何需要知道 10/22 19:00
16F：推 andrew43: 上句不对。我想说的是结合降雨量和时间後资料是如何对 10/22 19:03
17F：→ andrew43: 应。他是array吗？ 10/22 19:03
18F：推 andrew43: 1440和400个地点有重复吗？照理说降雨量有三个维度，经 10/22 19:18
19F：→ andrew43: 、维、日。但你说得有不全然如此。老实说我还是没全明 10/22 19:18
20F：→ andrew43: 白。 10/22 19:18
21F：→ AndrewShi: 地点没有重复，结合後降雨的维度就变成像你说的有三个 10/22 19:46
22F：→ AndrewShi: 维度(经.纬.时间)，对应的关系我就不是很清楚了，只是 10/22 19:46
23F：→ AndrewShi: 合并後dim(time)=NULL了，还是我把档案放上来?! 10/22 19:46
24F：→ AndrewShi: 合并後降雨是array，时间我上面是把它设为matrix，也可 10/22 22:30
25F：→ AndrewShi: 以改成array，但我不知道合并後的时间这样设定有何差 10/22 22:30
26F：→ AndrewShi: 异。 10/22 22:30
28F：→ AndrewShi: l7SnSxJzllkN9H9QjL8vMv3wF-WSpv5?usp=sharing 10/22 22:59
29F：→ andrew43: 网址不对，缩一下 10/22 23:55
31F：→ AndrewShi: kN9H9QjL8vMv3wF-WSpv5，抱歉andrew大，我不知道怎麽缩 10/23 01:27
32F：→ AndrewShi: 网址XD 10/23 01:27