R：将 netCDF 文件读取为栅格

Question

我正在尝试将 netCDF 文件作为栅格读入 R。 netCDF 文件将海洋年平均温度描述为经度、纬度和深度的函数。我关心的是表层海洋（即 netCDF 文件中的第一层），t_an 是 netCDF 文件中变量的名称。因此，我使用以下代码：

MyRast <- raster("Temperature.nc", level = 1, varname = "t_an")

这给了我以下警告：

Warning message:
In .getCRSfromGridMap4(atts) : cannot process these parts of the CRS: epsg_code=EPSG:4326

如您所见，netCDF 文件具有 CRS EPSG 4326（或 WGS 84），但创建的栅格具有以下 CRS：

+proj=longlat +lon_0=0 +a=6378137 +rf=298.257232666016

知道我如何使用正确的 CRS 读取文件 netCDF 吗？

Answer 1

读入数据后设置投影即可：

r <- raster("Temperature.nc",  varname = "t_an")
proj4string(r)=CRS("+init=EPSG:4326")

R: Reading a netCDF file as a raster