使用边界匹配限制参数在 R 中绘制栅格

Question

我有一张光栅图像，我想在 R 中用世界地图叠加层进行绘制。我可以同时绘制光栅图像和地图叠加层，但是为光栅图指定的 x 和 y 限制不是我指定的内容，在我不想在我的地图上显示的栅格边界之外留下了悬垂。

这是我的代码：

library(raster)
library(rworldmap)

r <- raster(ncol=500, nrow=500)
values(r) <- 1:ncell(r)
plot(r, ylim=c(-50,50), xlim=c(-100,100), col=rev(pal(100)))
world <- getMap(resolution="high")
map(world, add=TRUE, lwd=0.5)

生成这张地图的：

关于如何去除图形顶部和底部的白色非光栅悬垂物运行有什么想法吗？谢谢！！！

Answer 1

以下示例是对您的代码稍作修改。

library(raster)
library(rworldmap)
library(rgeos)

r <- raster(ncol=500, nrow=500)
values(r) <- 1:ncell(r)
r <- crop(r, extent(-100, 100, -50, 50))
world <- getMap(resolution="high")
world <- gBuffer(world, byid = T, width = 0)
world <- crop(world, r)
png(file = "C:/Example_Plot.png", height = 500, width = 880)
plot(r)
plot(world, add = T, lwd = 0.5)
dev.off()

该策略本质上是将多边形裁剪到栅格对象的范围并指定输出文件的宽度和高度。

使用边界匹配限制参数在 R 中绘制栅格

Plot raster in R with boundaries matching limits arguments

r

r-raster