在R中绘图时netCDF的经度/经度问题

Question

我试图使用netCDF4软件包从R中的a global climate model读取/绘制每月的海底温度数据，但遇到了经度数据不能正确显示的障碍。我在Panoply中检查了数据，它在那里绘制得很好，但当我在R中绘制它时，经度是不正确的：

​

​

我已经尝试过旋转和转置数据，但都无济于事，所以如果能有一个解决方案，我将非常感激。我使用的代码如下所示。谢谢。

nc_data <- nc_open("tob_Omon_MPI-ESM1-2-HR_ssp126_r1i1p1f1_gn_201501-201912.nc")

# Save the print(nc) dump to a text file
{
  sink('tob_Omon_MPI-ESM1-2-HR_ssp126_r1i1p1f1_gn_201501-201912_metadata.txt')
  print(nc_data)
  sink()
}

lon <- ncvar_get(nc_data, "longitude", verbose = F)
lon[lon > 180] <- lon[lon > 180] - 360
range(lon)
lat <- ncvar_get(nc_data, "latitude", verbose = F)
range(lat)
t <- ncvar_get(nc_data, "time")

tob.array <- ncvar_get(nc_data, "tob") # store the data in a 3-dimensional array
dim(tob.array) 

fillvalue <- ncatt_get(nc_data, "tob", "_FillValue")
fillvalue


tob.array[tob.array == fillvalue$value] <- NA # change NA values to the standard 'NA'

tob.2015 <- tob.array[,,1:12] # subset of 2015 data

# convert 2015 to raster brick
brick_2015 <- brick(tob.2015, xmn=min(lon), xmx=max(lon), ymn=min(lat), ymx=max(lat), crs=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs+ towgs84=0,0,0"))
plot(brick_2015) # inspect 
brick_2015 <- (t(brick_2015)) # transpose 
plot(brick_2015) # inspect
mean.2015 <- calc(brick_2015, mean) # save as one raster object
plot(mean.2015)

Answer 1

根据NetCDF文件中的元数据，这是一个很难处理的网格。在R中解决这些问题通常是不值得的。在这些情况下，我自己的偏好是使用气候数据运算符重新生成常规的latlon网格：

cdo remapbil,r360x180 tob_Omon_MPI-ESM1-2-HR_ssp126_r1i1p1f1_gn_201501-201912.nc out.nc