get_trip_geometry: Get trip geometry

In gtfstools: General Transit Feed Specification (GTFS) Editing and Analysing Tools

Description

Returns the geometry of each specified trip_id, based either on the shapes or the stop_times file (or both).

Usage

get_trip_geometry(
  gtfs,
  trip_id = NULL,
  file = c("shapes", "stop_times"),
  crs = 4326
)

Arguments

gtfs	A GTFS object as created by read_gtfs.
trip_id	A string vector including the trip_ids to have their geometries generated. If NULL (the default), the function generates geometries for every trip_id in the GTFS.
file	The file from which geometries should be generated. By default uses both shapes and stop_times.
crs	The CRS of the resulting object. Defaults to 4326 (WGS 84).

Value

A LINESTRING sf.

Details

The geometry generation works differently for the two files. In the case of shapes, the shape as described in the text file is converted to an sf object. For stop_times the geometry is the result of linking subsequent stops along a straight line (stops' coordinates are retrieved from the stops file). Thus, the resolution of the geometry when generated with shapes tends to be much higher than when created with stop_times.

Examples

data_path <- system.file("extdata/spo_gtfs.zip", package = "gtfstools")

gtfs <- read_gtfs(data_path)

trip_geometry <- get_trip_geometry(gtfs)
head(trip_geometry)

trip_ids <- c("CPTM L07-0", "2002-10-0")
trip_geometry <- get_trip_geometry(gtfs, trip_id = trip_ids)
trip_geometry
plot(trip_geometry["origin_file"])

gtfstools documentation built on Feb. 23, 2021, 9:06 a.m.