Get OSM routes geometry considering gtfs:shape_id match

Usage

osm_shapes_to_routes(
  gtfs,
  q,
  ways = FALSE,
  ways_tags = c("lanes", "psv", "bus", "way")
)

Arguments

gtfs: tidygtfs. GTFS feed.
q: osmdata::opq. Overpass query for transit network.
ways: boolean (Default False). If true, relation is disaggregated in ways.
ways_tags: character vector (Default c("lanes", "psv", "bus", "way")). List of OSM way tags to extract when ways parameter is set to true. Match is done using tidyselect::contains().

Value

A sf data.frame with the following columns:

shape_id, the shape_id attribute from shapes.txt file.
osm_id, the osm_id attribute from OSM route relation.
way_osm_id, the osm_id attribute from OSM way (if ways parameter is set to true).
*, any column that matches ways_tags parameter.
geometry, the geometrical data for the OSM route relation.

Shapes that do not have a match on OSM are ignored. If that occurs, a warning is displayed during the method execution, informing about the missing geometries.

Details

For each route, matches its trips' shapes with OSM route relations, considering the OSM gtfs:shape_id attribute.

Examples

if (FALSE) { # \dontrun{
gtfs <- GTFShift::load_feed("gtfs.zip")

q = opq("Lisbon")  |>
  add_osm_feature(key = "route", value = c("bus", "tram")) |>
  add_osm_feature(key = "network", value = "Carris", key_exact = TRUE)

shapes_geometry_osm = GTFShift::osm_shapes_to_routes(gtfs, q)
} # }