NetCDF在R中显示轨迹数据

我试图显示以下经典的NetCDF数据。这个数据有38个变量轨迹数据，我试图提取其中一个“空气温度”。虽然我能够提取数据，但不幸的是，我无法使用绘图在R中显示数据。显示应该如下图所示，我用panoply绘制了它。另一个问题是，我无法堆叠所有数据，因为轨迹都不同。是否有任何方法可以显示和堆叠一个变量的所有文件

> library(ncdf4)
> library(rasterVis)
> library(raster)
lon <- ncvar_get(ncin, "lon")
lat <- ncvar_get(ncin, "lat")
data <-lon <- ncvar_get(ncin, "air_temp_ac")          #to extract variable
> dim(data)
[1] 6639
> dput(data)
structure(c(NA, 15, 14, 13, 11, NA, 11, 11, 11, 11, 11, NA, 14, 
14, 12, NA, 12, 14, 14, 14, 16, 21, 25, 27, 26, 19, 21, 22, 22, 
21, 22, 22, 23, 24, 22, 23, 25, 23, 26, 25, 25, 23, 24, 24, 25, 
28, 28, 29, 29, 28, 28, 26, 27, 27, 29, 31, 31, 34, 36, 37, 38, 
41, 37, 37, 39, 36, 29, 33, 38, 35, 36, 36, 36, 37, 37, 36, 37, 
34, 33, 34, 38, 37, 37, 37, 37, 37, 39, 39, 39, 40, 39, 39, 40,
> plot(lon,lat,data)
Error in plot.xy(xy, type, ...) : invalid plot type

>库（ncdf4）
>图书馆（拉斯特维斯）
>图书馆（光栅）
lon r图（r）

但是，这个情节是越界的

---

这是一个使用

ggplot2

、

sf

和

rnaturalearth

的解决方案。使用您的代码，我无法从nc文件中提取变量“air\u temp\u ac”。实际变量为“air_temp_AC”（注意这是区分大小写的）

库（ncdf4）
图书馆（GG2）
图书馆（rnaturalearth）
图书馆（sf）
#>链接至GEOS 3.6.1、GDAL 2.2.3、项目4.9.3
全世界非常感谢你@Allan
> r <- raster(t(data), xmn=min(lon), xmx=max(lon), ymn=min(lat), ymx=max(lat), crs=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs+ towgs84=0,0,0"))
> r <- flip(r, direction='y')
> plot(r)