osm_polygons: Extract all 'osm_polygons' from an osmdata object

In osmdata: Import 'OpenStreetMap' Data as Simple Features or Spatial Objects

Extract all osm_polygons from an osmdata object

Description

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.

Usage

osm_polygons(dat, id)

Arguments

dat	An object of class osmdata
id	OSM identification of one or more objects for which polygons are to be extracted

Value

An sf Simple Features Collection of polygons

See Also

Other search: osm_lines(), osm_multilines(), osm_multipolygons(), osm_points()

Examples

## Not run: 
# Extract polygons which intersect Conway Street in London
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)

## End(Not run)

osmdata documentation built on Aug. 14, 2023, 5:08 p.m.