使用另一个 SF 多边形剪辑 SF 多边形

Question

如果我必须为不规则形状文件（例如国家/地区）裁剪光栅，我会先执行

crop

，然后执行

mask

，如下所示：

library(raster)
r <- crop(my_raster, extent(polygon))
r <- mask(r, polygon)

是否有对多边形执行相同操作的等效操作？我尝试使用下面的方法来剪辑给定国家/地区的世界河流流域形状文件，但它只返回包含在范围内的流域，而不是形状文件边界之后的流域

  library(sf)
  cropped <- st_crop(river_basins, country_shapefile)

Answer 1

对于 sp 对象，您可以使用 gClip() 函数，如下所述：

希望这有帮助。