R：将netCDF文件作为栅格

Question

我试图将一个netCDF文件作为栅格读入。 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")