R语言下载CERRA总降水数据异常：GRIB负值、NetCDF无坐标系

一、修复GRIB格式总降水的负值问题 #

1. 修正产品类型参数错误 #

你的代码中product_type=" forecast"存在两处问题：一是开头有多余空格，二是针对CERRA再分析数据，正确的product_type应为"reanalysis"（而非预报类型）。该错误可能导致下载到异常编码的预报场数据，先修正参数：

query <- r_to_py(list(
  variable= "total_precipitation",
  level_type= "surface_or_atmosphere",
  data_type= "reanalysis",
  product_type="reanalysis", # 修正此处
  year= year_n,
  month= month_n,
  day= day_n, 
  time="00:00",
  format= "grib",
  leadtime_hour = "24"
))

2. 应用GRIB编码的缩放与偏移量 #

若修正参数后仍有负值，大概率是GRIB文件采用**缩放因子（scale_factor）和偏移量（add_offset）**进行压缩编码，terra默认未自动解析应用。通过以下步骤转换为真实降水值：

nc_raster <- rast(Path)
# 提取元数据中的缩放与偏移信息
meta <- metadata(nc_raster)[[1]]
scale_factor <- as.numeric(meta[grep("scaleFactor", meta)])
add_offset <- as.numeric(meta[grep("offset", meta)])

# 转换为真实值（单位从米转为毫米）
precip_raster <- nc_raster * scale_factor + add_offset
precip_raster <- precip_raster * 1000
# 过滤异常负值
precip_raster[precip_raster < 0] <- NA

plot(precip_raster, 1)

二、修复NetCDF文件的坐标系问题 #

1. 手动指定CERRA的正确坐标系 #

CERRA再分析数据采用Lambert Conformal Conic投影，可直接使用官方对应的EPSG编码或投影参数：

# 读取NetCDF文件
nc_raster <- rast("你的NC文件路径.nc")

# 方法1：使用EPSG编码（欧洲ETRS89 LAEA投影）
crs(nc_raster) <- "EPSG:3035"

# 方法2：手动输入官方投影参数（若EPSG编码不生效）
crs(nc_raster) <- "+proj=lcc +lat_1=35 +lat_2=65 +lat_0=50 +lon_0=10 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs"

# 验证坐标系设置结果
print(crs(nc_raster))

2. 处理维度警告 #

你遇到的vobjtovarid4警告属于NetCDF文件本身的元数据缺陷（y维度未关联坐标变量），不影响坐标系的实际生效。手动指定CRS后，栅格的空间参考即可正常用于后续空间分析。

内容的提问来源于stack exchange，提问作者Macadamia Brazowy