PTT
Submit
Submit
選擇語言
正體中文
简体中文
PTT
R_Language
[問題] 月平均資料
作者:
AndrewShi
(沒有妳的我)
2018-10-22 17:15:18
[問題類型]:
程式諮詢(我想用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檔 降雨
作者:
AndrewShi
(沒有妳的我)
2018-10-23 01:23:00
https://drive.google.com/drive/folders/1aI7SnSxJzII
kN9H9QjL8vMv3wF-WSpv5,抱歉andrew大,我不知道怎麼縮網址XD
http://0rz.tw/JI056
作者:
andrew43
(討厭有好心推文後刪文者)
2018-10-22 23:55:00
網址不對,縮一下
作者:
AndrewShi
(沒有妳的我)
2018-10-22 19:46:00
地點沒有重複,結合後降雨的維度就變成像你說的有三個維度(經.緯.時間),對應的關係我就不是很清楚了,只是合併後dim(time)=NULL了,還是我把檔案放上來?!合併後降雨是array,時間我上面是把它設為matrix,也可以改成array,但我不知道合併後的時間這樣設定有何差異。檔案我放在這裡,
https://drive.google.com/folders/1a
l7SnSxJzllkN9H9QjL8vMv3wF-WSpv5?usp=sharing
作者:
andrew43
(討厭有好心推文後刪文者)
2018-10-22 18:46:00
對的。1440和400這與時間的對應關係又如何需要知道上句不對。我想說的是結合降雨量和時間後資料是如何對應。他是array嗎?1440和400個地點有重複嗎?照理說降雨量有三個維度,經、維、日。但你說得有不全然如此。老實說我還是沒全明白。
作者:
AndrewShi
(沒有妳的我)
2018-10-22 18:44:00
dim(lon(經度))=1440dim(lat(緯度))=400dim(time)=365(原本資料型態)dim(precip(降雨量))=1440 400(rbind之後))dim(precip(降雨量))=1440 400 365以我的理解是原本的資料型態為一天一張全球地圖(平面)上面每個網格點有一個降雨值,總共365個平面(天),而我現在想要取出比如其中一個月的降雨量來做平均(全球),andrew大是指這個資料資訊嗎?!
作者:
andrew43
(討厭有好心推文後刪文者)
2018-10-22 18:27:00
你的每個loop都沒有用到time這個變數,所以就只是同一件事重覆做很多次,看你設定的time有幾個。建議你提供一個小nc檔,不然連資料格式都不清楚沒辦法答不會自製nc檔的話,至少把lon lat time precip都說清楚
繼續閱讀
Re: [問題] 請問怎麼把rawdata分組對應?? (資料清洗)
andrew43
[問題] 請問怎麼把rawdata分組對應?? (資料清洗)
ronald736
[問題] notepad++編輯R script
k97231
[問題] 同義字轉成標準字(中文文字處理)
ardodo
[問題] 從Rademacher distribution找limitation
AmigoSafin
Re: [問題] 利用R下載外匯資料
TreeMan
[問題] 利用R下載外匯資料
ericchin0404
[問題] Bernoulli distribution驗證weak large n
AmigoSafin
[問題] 用curve command 將一個pdf畫圖
AmigoSafin
[問題] 如何將function放入迴圈? (求救)
ronald736
Links
booklink
Contact Us: admin [ a t ] ucptt.com