使用 R 绘制 ERA-INTERIM

Question

我已经从 ECMWF ERA-INTERIM (click here to download the file) 下载了温度数据 (netcdf)，但是当我尝试绘制它时，它似乎不完整。我做错了什么？

library(ncdf) 
library(raster)

UK <- getData("GADM",country="GBR",level=1)
england <- UK[UK$NAME_1 == 'England', ]

# To convert longitude from range [0,360] to [-180,180]
x <- rotate(raster(x = "netcdf-atls13-a562cefde8a29a7288fa0b8b7f9413f7-ORIAcW.nc", varname="t2m"))

plot(x, ext=extent(england))
plot(england, add=T)

Answer 1

这是一个错误。它已在光栅 2.4-24 中得到修复（可通过 R-Forge 获得）。这是针对此特定情况的解决方法：

library(raster)
r <- raster("netcdf-atls13-a562cefde8a29a7288fa0b8b7f9413f7-ORIAcW.nc", varname="t2m")
x <- rotate(r)
x[, 240] <- r[1:240,1]

使用 R 绘制 ERA-INTERIM

Plot ERA-INTERIM using R

r

netcdf