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")