R:在特定坐标下更改栅格值

Question

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

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

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

Answer 1

是的，但您需要提供一个最小的、自包含的、可重现的示例。最简单的方法可能是使用非洲多边形africa和do

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

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

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