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

看板R_Language作者AndrewShi (沒有妳的我)時間6年前 (2018/05/14 22:56)推噓0(0推 0噓 0→)

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

※ 引述《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://www.ptt.cc/bbs/R_Language/M.1526309815.A.F20.html