R - 海岸线在地图图中没有正确拟合

Question

我有这个data set 它是 GPCC 1x1 降水数据：为了绘制它，我已经这样做了：

library(RNetCDF)
library(raster)
library(pracma)
library(rasterVis)
library(colorRamps)
library(rgdal)


file <- open.nc("C:/Users/Ann/Downloads/precipitacion_mensual.nc")
var <- var.get.nc(file,"precip")
ro <- var[ , ,80]

ro_2 <- pracma::circshift(ro, c(180,2))

A <- raster(t(ro_2))

extent(A) <- c(-180, 180, -90, 90)

mex <- getData("GADM", country="MEX", level=0)

B <- crop(A,mex)
plot(B,main="Tit", xlab = "lon", ylab ="lat",col=matlab.like(n = 10),zlim = c(0, 500))

plot(mex, add =TRUE)

但是当我绘制它时，它会像这样 如您所见，海岸线在地图上没有正确重叠。我怎样才能解决这个问题？

Answer 1

原来问题出在我认为的circshift 函数中，通过写入circshift(variable,c(180,2))，它将沿第二维进行180 度的循环移位。 然而，这样做是将第一个维度移动 180 度，将第二个维度移动 2 度。 我只需输入circshift(var, c(179.5,0.5))，地图就完美了。