R：在某些坐标处更改栅格值

Question

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

1）我在哪里可以下载大陆的shapefile？

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

Answer 1

是的，但是您需要提供一个最小的、独立的、可重现的示例。 最简单的方法可能是使用 Africa polygons africa并做

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

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

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