我目前正在R中与澳大拉西亚地区范围内的环境数据合作,对物种分布进行建模。我使用raster
包将这些环境文件转换为栅格进行分析,并将它们堆叠在一起。
所有这些光栅都以格林威治子午线为中心。但是,出于可视化目的,我希望能够将这些光栅绘制在太平洋上。
我试图修改CRS,但它没有改变绘图输出。
有人能建议一种方法来重新定位子午线180中的光栅吗?
非常感谢
此代码是我使用过的代码的简单版本。
xy <- matrix(rnorm(400),20,20)
rast <- raster(xy)
extent(rast) <- c(-180, 180, -60, -3.197442e-14)
projection(rast) <- CRS("+proj=longlat +datum=WGS84")
plot(rast)