Re: [問題] dimensions of z are not length(x)(-1)

看板R_Language作者cywhale (cywhale)時間6年前 (2018/05/15 11:23)推噓6(6推 0噓 47→)

留言53則, 2人參與討論串3/3 (看更多)

#ncfile 讀出是一個多維矩陣，你用temp[1:4]當然只會讀出4筆以下供參 #你可以針對每一次 (1:4)都畫一次，在這裡我不確定你要的操作，我取平均 #OB library(ncdf4) #library(magrittr) #library(lubridate) #library(units) infile <- "D:/Env/air/air.2m.mon.mean.nc" nx0 <- nc_open(infile) print(nx0) air <- ncvar_get(nx0, "air") latx0<- ncvar_get(nx0, "lat") lngx0<- ncvar_get(nx0, "lon") #time<- ncvar_get(nx0, "time") #tattr <- ncatt_get(nx0, "time") #get attributes #units(time) <- make_unit(tattr$units) #datex<- as.POSIXct(time) %>% as.Date(origin="1800-01-01 00:00") #:0.0") #lvl <- ncvar_get(nx0, "level") mair <- apply(air[,,1:4],c(1,2),function(x) mean(x-273.15, na.rm=T)) zt=matrix(mair,ncol=length(latx0),nrow=length(lngx0),byrow=F) image(x=lngx0-180,y=rev(latx0),z=zt[,ncol(zt):1], col = colorRamps::blue2red(128)) #rainbow(128)[50:128]) ※ 引述《AndrewShi (沒有妳的我)》之銘言： : ※ 引述《AndrewShi (沒有妳的我)》之銘言： : : [問題類型]: : : 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來) : : [軟體熟悉度]: : : 入門(寫過其他程式，只是對語法不熟悉) : : [問題敘述]: : : 這是一個全球溫度的資料(nc檔)，時間總共有467筆，我想利用R畫出前4筆的全球溫度圖 : : ，所以寫了一個簡單的while迴圈，但到了要image的時候就跑出了這個error，所以想請問大大 : : 們這個error是什麼意思，以及我該如何解決??謝謝~ : : [程式範例]: : : https://i.imgur.com/JlqAFxw.jpg

: : https://i.imgur.com/fLzCmiC.jpg

: 以下為資料來源與程式碼: : 資料來源: : https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis2.pressure.html : Air Temperature - Monthly Mean - 按圖示 -->裡面的最後一筆資料 : 程式碼: : rm(list=ls()) : library(ncdf4) : air.nc<-nc_open("C:\\Users\\TOM\\Desktop\\R(資料庫)\\air.2m.mon.mean.nc") : air.nc : "file C:\\Users\\TOM\\Desktop\\R(資料庫)\\air.2m.mon.mean.nc has 5 : dimensions:" : "lon Size:144" : "lat Size:73" : "level Size:17" : "nbnds Size:2" : "time Size:467" : "file C:\\Users\\TOM\\Desktop\\R(資料庫)\\air.2m.mon.mean.nc has 2 variables" : "short air[lon,lat,level,time] Longname:Air temperature Missval:-32767" : "short head[level,time] Longname:Missing Missval:NA" : temp.nc<-nc_open("C:\\Users\\TOM\\Desktop\\R(資料庫)\\air.2m.mon.mean.nc") : temp<-ncvar_get(temp.nc,"air") : temp[temp=="-32767"]<-NA : temp.nc$dim$lon$vals -> lon : temp.nc$dim$lat$vals -> lat : temp.nc$dim$level$vals -> lev : temp.nc$dim$nbnds$vals -> nbnds : temp.nc$dim$time$vals -> time : lat<-rev(lat) : i<-1 : while(i<5){ : print(temp[i]) : i<-i+1 : } : time="i"<-temp[i] : lon<-lon-180 : temp<-(time="i") : temp<-matrix(c("i"=1:4)) : image(lon,lat,matrix(temp,ncol=length(lat),nrow=length(lon),byrow=F)) : library(sp) : library(maptools) : gpclibPermit() : data(wrld_simpl) : plot(wrld_simpl,add=TRUE) -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.65.48 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1526354622.A.065.html

推