Clips an sf data frame to an area of interest polygon using
sf::st_intersection(). Handles geometry type cleanup automatically —
mixed geometry collections from intersection are filtered to the original
geometry type (e.g. polygon, linestring).
Value
An sf data frame clipped to aoi, with geometry type matching
the input. Returns an empty sf with the same columns if no features
intersect.
Details
Typical use: clip network query results (lakes, wetlands, streams) to a
watershed polygon from frs_watershed_at_measure().
Examples
if (FALSE) { # \dontrun{
# Clip wetlands to a subbasin watershed
aoi <- frs_watershed_at_measure(blk, drm, upstream_measure = urm)
wetlands <- frs_network(blk, drm, tables = list(
wetlands = "whse_basemapping.fwa_wetlands_poly"
))
wetlands_clipped <- frs_clip(wetlands, aoi)
} # }
