我一直在尝试绘制以下网格化netcdf文件:"air.1999.nc"在以下网站上找到:
http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.html
我根据我在这里和其他地方找到的答案尝试了下面的代码,但没有运气.
library(ncdf);
temp.nc <- open.ncdf("air.1999.nc");
temp <- get.var.ncdf(temp.nc,"air");
temp.nc$dim$lon$vals -> lon
temp.nc$dim$lat$vals -> lat
lat <- rev(lat)
temp <- temp[nrow(temp):1,]
temp[temp==-32767] <- NA
temp <- t(temp)
image(lon,lat,temp)
library(maptools)
data(wrld_simpl)
plot(wrld_simpl, add = TRUE)
这段代码是从这里找到的代码修改而来的:来自netcdf文件的变量出现了翻转
有没有人对使用这些类型的netcdf文件有任何想法或经验？谢谢