作者AndrewShi (沒有妳的我)
看板R_Language
標題Re: [問題] dimensions of z are not length(x)(-1)
時間Mon May 14 22:56:52 2018
※ 引述《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), 來自: 150.117.35.12
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1526309815.A.F20.html