R中处理ERA5-Land数据遇异常：缺失值与绘图空白求助

核心错误点 #

• 维度索引错位：ERA5-Land的netCDF维度顺序通常是(time, expver, latitude, longitude)，而非你假设的(longitude, latitude, expver, time)，直接用[,,1,880]会读取到错误的切片，导致数据无效或错位。
• 经度范围不匹配：ERA5默认经度范围是0°-360°，但法国区域经度在-5°~10°区间，若未转换为-180°~180°的常规范围，创建栅格时空间范围会偏移，导致绘图空白。
• 栅格参数定义错误：手动指定经纬度范围时，若未正确读取原始数据的lon/lat变量，会导致栅格的空间参考与实际数据不匹配。

修复步骤 #

1. 先确认维度顺序 #

用ncdf4查看变量的维度结构，避免索引错误：

library(ncdf4)
library(terra)  # 推荐用terra替代raster，空间数据处理更高效
library(ggplot2)

# 读取netCDF文件
ncfile <- nc_open("你的数据文件路径.nc")
# 打印t2m变量的维度信息
print(ncfile$var$t2m$dim)

输出会明确显示维度顺序（比如time, expver, latitude, longitude），后续索引要严格对应这个顺序。

2. 正确合并expver数据集 #

ERA5中expver=1是主数据，expver=2是备用补全数据，用terra可以一键读取并合并：

# 直接读取t2m变量，自动识别维度
t2m_rast <- rast("你的数据文件路径.nc", var="t2m")
# 合并expver：主数据（expver=1）有值就用主数据，缺失的用备用数据（expver=2）填充
t2m_merged <- app(t2m_rast, function(x) {
  ifelse(is.na(x[1]), x[2], x[1])
})
# 提取第880个时间切片
t2m_slice <- t2m_merged[[880]]

3. 转换经度范围（若需要） #

检查当前经度范围，若为0-360°，转换为-180-180°：

# 查看当前经度极值
print(xmin(t2m_slice))
print(xmax(t2m_slice))

# 转换经度
if(xmax(t2m_slice) > 180) {
  t2m_slice <- rotate(t2m_slice)
}

4. 裁剪到法国区域并绘图 #

手动指定法国的经纬度范围，裁剪后绘图：

# 法国大致经纬度范围：经度-5~10，纬度41~51
france_ext <- ext(-5, 10, 41, 51)
t2m_france <- crop(t2m_slice, france_ext)

# 绘图
ggplot() +
  geom_spatraster(data = t2m_france) +
  scale_fill_viridis_c(name = "2m气温(K)") +
  coord_sf(crs = "EPSG:4326") +
  labs(title = "ERA5-Land法国区域2月均温") +
  theme_minimal()

额外提醒 #

• 确认Copernicus CDS下载页面的空间范围参数确实覆盖法国，避免下载时范围选错。
• 若坚持用raster包，需通过dim(t2m_var)查看维度大小，再对应调整索引位置，比如如果维度是(time, expver, lat, lon)，取第880个时间切片的正确索引是[880,1,,]和[880,2,,]。

内容的提问来源于stack exchange，提问作者Vo Hoang Linh Chi