Recortar objeto de características simples en R

Question

¿Existe una función para recortar el objeto del mapa sf, similar a maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15)) utilizado para SpatialPolygon o SpatialLine?

Estoy considerando st_intersection() pero puede haber una forma adecuada.

Answer 1

Otra solución alternativa, para mí fue más rápido para shapefiles más grandes:

 library(sf)
library(raster)
library(rgeos)
library(ggplot2)

# Load National Forest shapefile
# https://data.fs.usda.gov/geodata/edw/edw_resources/shp/S_USA.AdministrativeForest.zip
nf.poly <- st_read("S_USA.AdministrativeForest"), "S_USA.AdministrativeForest")

crop_custom <- function(poly.sf) {
  poly.sp <- as(poly.sf, "Spatial")
  poly.sp.crop <- crop(poly.sp, extent(c(-82, -80, 35, 36)))
  st_as_sf(poly.sp.crop)
}

cropped <- crop_custom(nf.poly)