无CRS且翻转的温度栅格无法与厄瓜多尔省界SHP匹配的技术求助

核心问题原因 #

直接读取CPC的NetCDF文件时，R的raster包可能无法自动识别地理参考信息，导致栅格范围被归一化、CRS缺失，同时NetCDF的纬度通常从北到南存储，造成地图翻转。

步骤1：正确读取并修复栅格的地理信息 #

使用brick读取NetCDF（支持多变量/时间序列），手动指定CRS和地理范围，再修正纬度方向：

library(raster)
library(ncdf4)

# 读取NetCDF，替换为你的文件路径，varname根据文件内温度变量名调整（比如CPC的日平均温通常是"tavg"）
temp_brick <- brick("2018_cpc_temp.nc", varname = "tavg")

# 提取2018年11月10日的图层：先匹配时间索引
target_date <- as.Date("2018-11-10")
layer_index <- which(getZ(temp_brick) == target_date)
temp_raster <- temp_brick[[layer_index]]

# 设置标准WGS84坐标系（CPC数据默认投影）
crs(temp_raster) <- "+proj=longlat +datum=WGS84 +no_defs"

# 重置地理范围为全球经纬度范围
extent(temp_raster) <- extent(-180, 180, -90, 90)

# 翻转纬度方向（解决地图倒置问题）
temp_raster_fixed <- flip(temp_raster, direction = "y")

运行print(temp_raster_fixed)验证输出：确认crs显示为WGS84，extent为-180, 180, -90, 90。

步骤2：匹配厄瓜多尔省界并完成裁剪/统计 #

确保省界与栅格CRS一致后，进行裁剪、掩膜和省份平均气温计算：

# 读取厄瓜多尔省界shapefile，替换为你的文件路径
ecuador_provinces <- shapefile("ecuador_provinces.shp")

# 统一坐标系：如果省界CRS与栅格不同，转换省界到栅格CRS
if (!identical(crs(ecuador_provinces), crs(temp_raster_fixed))) {
  ecuador_provinces <- spTransform(ecuador_provinces, crs(temp_raster_fixed))
}

# 裁剪栅格到厄瓜多尔范围
temp_ecuador_cropped <- crop(temp_raster_fixed, extent(ecuador_provinces))

# 掩膜保留省界内的有效温度数据
temp_ecuador_masked <- mask(temp_ecuador_cropped, ecuador_provinces)

# 计算各省份平均气温
province_avg_temp <- zonal(temp_ecuador_masked, ecuador_provinces, fun = mean)
print(province_avg_temp)

排查技巧 #

如果翻转后仍有问题，检查NetCDF的纬度维度信息：

nc <- nc_open("2018_cpc_temp.nc")
# 查看纬度的起始值和方向
print(nc$dim$lat)
nc_close(nc)

若纬度是从90°N到90°S存储，flip(direction="y")即可修正；若维度方向特殊，可尝试t()转置后再调整范围。

内容的提问来源于stack exchange，提问作者Alonso Quijano