使用 R 绘制 ERA-INTERIM答案

【问题标题】：Plot ERA-INTERIM using R使用 R 绘制 ERA-INTERIM
【发布时间】：2016-02-09 09:28:53
【问题描述】：

我已经从 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)

【问题讨论】：

标签： r netcdf

【解决方案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]

【讨论】：