R：在特定坐标处更改栅格值

Question

我 运行 R 中的物种分布模型，想为没有岛屿的非洲大陆创建变量栅格。我只能找到非洲及其岛屿的 shapefile，而不仅仅是来自大陆。

1)哪里可以下载大陆的shapefile？

2) 如果没有 shapefile，我想从我的栅格中手动删除这些岛屿。有没有办法做到这一点，f.e。将某些坐标之间的部分栅格设置为 NA?

Answer 1

是的。这是一个最小的、独立的、可重现的例子。最简单的方法可能是使用非洲多边形 africa 并执行

library(raster)
afr <- aggregate(africa)
v <- disaggregate(afr)
a <- area(v)
afnois <- v[which.max(a), ]

然后在 mask 中使用它从栅格中移除岛屿

您还可以使用 raster::drawPoly 创建多边形并将其用于遮罩。