If id
is of a point object, osm_polygons
will return all
polygons containing that point. If id
is of a line or polygon object,
osm_polygons
will return all polygons which intersect the given line
or polygon.
See also
Other search:
osm_lines()
,
osm_multilines()
,
osm_multipolygons()
,
osm_points()
Examples
# Extract polygons which intersect Conway Street in London
if (FALSE) { # \dontrun{
dat <- opq ("Marylebone London") |>
add_osm_feature (key = "highway") |>
osmdata_sf ()
conway <- which (dat$osm_lines$name == "Conway Street")
id <- rownames (dat$osm_lines [conway, ])
osm_polygons (dat, id)
} # }