Skip to contents

2. Filter GTFS feeds

Source: vignettes/filter.Rmd
filter.Rmd

Introduction

GTFS feeds do not have a defined scope regarding its coverage of the transportation system. Some can be bounded to one agency, whereas others can aggregate several modes in the same city, or even national wise.

From the simpler to the most complex feeds, some analysis require to narrow the perspective. GTFShift provides some to help in this process.

Filters

Filter by agency

This example uses a national GTFS feed for long distance rail in Germany, retrieved from https://gtfs.de/en/feeds/de_fv/.

# Load GTFS
gtfs = load_feed("https://download.gtfs.de/germany/fv_free/latest.zip")
summary(gtfs)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, attributions, calendar, calendar_dates, feed_info, stops
#> agencies     DB Fernverkehr AG, ZSSK, HZZP ... 11 more
#> service      from 2026-05-16 to 2026-06-15
#> uses         stop_times (no frequencies)
#> # routes       96
#> # trips      6209
#> # stop_ids   1255
#> # stop_names  753
#> # shapes     1755

Multimodal feeds aggregate several agencies. GTFShift::filter_by_agency(), allows to filter by the agency id, name, or both. You can filter by the id and/or name of the agency.

# Filter by agency id
gtfs_5 = GTFShift::filter_by_agency(gtfs, id = 5)
summary(gtfs_5)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, attributions, calendar, calendar_dates, feed_info, stops
#> agencies     DB Fernverkehr AG, ZSSK, HZZP ... 11 more
#> service      from 2026-05-16 to 2026-06-15
#> uses         stop_times (no frequencies)
#> # routes      2
#> # trips      68
#> # stop_ids    8
#> # stop_names  7
#> # shapes      7

# Filter by agency name
gtfs_sncf = GTFShift::filter_by_agency(gtfs, name = "SNCF")
summary(gtfs_sncf)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, attributions, calendar, calendar_dates, feed_info, stops
#> agencies     DB Fernverkehr AG, ZSSK, HZZP ... 11 more
#> service      from 2026-05-16 to 2026-06-15
#> uses         stop_times (no frequencies)
#> # routes      2
#> # trips      68
#> # stop_ids    8
#> # stop_names  7
#> # shapes      7

Original GTFS 

shape_agency = gtfs$trips |>
  left_join(gtfs$routes, by = "route_id") |>
  left_join(gtfs$agency, by = "agency_id") |>
  select(shape_id, agency_id, agency_name) |>
  distinct() |> 
  mutate(name = paste(agency_id, agency_name, sep = ": "))

shapes_sf = tidytransit::shapes_as_sf(gtfs$shapes) |>
  left_join(shape_agency, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "name", legend = TRUE, layer.name="Agency")

GTFS filtered for agency id = 5 

shape_agency = gtfs_5$trips |>
  left_join(gtfs_5$routes, by = "route_id") |>
  left_join(gtfs_5$agency, by = "agency_id") |>
  select(shape_id, agency_id, agency_name) |>
  distinct() |> 
  mutate(name = paste(agency_id, agency_name, sep = ": "))

shapes_sf = tidytransit::shapes_as_sf(gtfs_5$shapes) |>
  left_join(shape_agency, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "name", legend = TRUE, layer.name = "Agency")

GTFS filtered for agency name = SNCF 

shape_agency = gtfs_sncf$trips |>
  left_join(gtfs_sncf$routes, by = "route_id") |>
  left_join(gtfs_sncf$agency, by = "agency_id") |>
  select(shape_id, agency_id, agency_name) |>
  distinct()|> 
  mutate(name = paste(agency_id, agency_name, sep = ": "))

shapes_sf = tidytransit::shapes_as_sf(gtfs_sncf$shapes) |>
  left_join(shape_agency, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "name", legend = TRUE, layer.name="Agency")

Filter by mode

This example uses a metropolitan wide GTFS feed for Los Angeles rail services, retrieved from https://developer.metro.net/gtfs-schedule-data/.

# Load GTFS
gtfs = GTFShift::load_feed("https://gitlab.com/LACMTA/gtfs_rail/raw/master/gtfs_rail.zip", create_transfers=FALSE)
summary(gtfs)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, fare_attributes, fare_rules, shapes, calendar, calendar_dates, feed_info, stops
#> agency       Metro - Los Angeles
#> service      from 2026-05-21 to 2026-06-04
#> uses         stop_times (no frequencies)
#> # routes        6
#> # trips      7129
#> # stop_ids    463
#> # stop_names  357
#> # shapes       18

Multimodal feeds aggregate several modes. GTFShift::filter_by_mode() allows to restrict the feed to some modes only. Refer to routes.txt ‘route_type’ parameter on GTFS documentation for more details on the modes id that should be used as parameters.

# Filter by mode tram
gtfs_tram = GTFShift::filter_by_modes(gtfs, modes = list(0))
summary(gtfs_tram)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, fare_attributes, fare_rules, shapes, calendar, calendar_dates, feed_info, stops
#> agency       Metro - Los Angeles
#> service      from 2026-05-21 to 2026-06-04
#> uses         stop_times (no frequencies)
#> # routes        4
#> # trips      5147
#> # stop_ids     95
#> # stop_names   95
#> # shapes        8

Original GTFS 

shape_route = gtfs$trips |>
  left_join(gtfs$routes, by = "route_id") |>
  select(shape_id, route_id, route_type) |>
  distinct() |> 
  mutate(route_type = as.character(route_type))

shapes_sf = tidytransit::shapes_as_sf(gtfs$shapes) |>
  left_join(shape_route, by = "shape_id") |>
  filter(!is.na(route_type))
mapview::mapview(shapes_sf, zcol = "route_type", legend = TRUE, layer.name="Route type")

GTFS filtered for mode tram (0) 

shape_route = gtfs_tram$trips |>
  left_join(gtfs_tram$routes, by = "route_id") |>
  select(shape_id, route_id, route_type) |>
  distinct() |> 
  mutate(route_type = as.character(route_type))

shapes_sf = tidytransit::shapes_as_sf(gtfs_tram$shapes) |>
  left_join(shape_route, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "route_type", legend = TRUE, layer.name="Route type")

Filter by route_name

This article uses a GTFS feed from the library GTFS database for Portugal as an example. Refer to the vignette(“download”) for more details.

# Get GTFS from library GTFS database for Portugal
data = read.csv(system.file("extdata", "gtfs_sources_pt.csv", package = "GTFShift"))
gtfs = GTFShift::load_feed(data[data$ID=="faro",]$URL, create_transfers=FALSE)
summary(gtfs)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, fare_attributes, fare_rules, shapes, calendar, calendar_dates, feed_info, stops
#> agency       Próximo - Transportes Urbanos de Faro
#> service      from 2024-01-01 to 2026-12-31
#> uses         stop_times (no frequencies)
#> # routes      14
#> # trips      745
#> # stop_ids   189
#> # stop_names 127
#> # shapes      29

GTFS feeds aggregate several routes. GTFShift::filter_by_route_name(), allows to filter the feed for specific routes, given a partial or total match with the short or the long name.

# Filter by short_name with exact match
gtfs_1 = GTFShift::filter_by_route_name(
  gtfs,
  values = list("1"),
  short_name = TRUE,
  exact_match = TRUE
)
summary(gtfs_1)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, fare_attributes, fare_rules, shapes, calendar, calendar_dates, feed_info, stops
#> agency       Próximo - Transportes Urbanos de Faro
#> service      from 2024-01-02 to 2026-12-31
#> uses         stop_times (no frequencies)
#> # routes      1
#> # trips      82
#> # stop_ids   23
#> # stop_names 19
#> # shapes      1

# Filter by long_name with partial match
gtfs_terminal = GTFShift::filter_by_route_name(
  gtfs,
  values = list("Terminal", "Rodoviário"),
  short_name = FALSE,
  exact_match = FALSE
)
summary(gtfs_terminal)
#> tidygtfs object
#> files        agency, routes, stop_times, trips, fare_attributes, fare_rules, shapes, calendar, calendar_dates, feed_info, stops
#> agency       Próximo - Transportes Urbanos de Faro
#> service      from 2024-01-01 to 2026-12-31
#> uses         stop_times (no frequencies)
#> # routes       6
#> # trips      406
#> # stop_ids   104
#> # stop_names  73
#> # shapes      16

Original GTFS 

route_name = gtfs$trips |>
  left_join(gtfs$routes, by = "route_id") |>
  select(shape_id, route_id, route_short_name, route_long_name) |>
  distinct() |> 
  mutate(name = paste(route_short_name, route_long_name, sep = ": "))

shapes_sf = tidytransit::shapes_as_sf(gtfs$shapes) |>
  left_join(route_name, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "name", legend = TRUE, layer.name="Route name")

GTFS filtered for route short name 1 (exact match) 

route_name = gtfs_1$trips |>
  left_join(gtfs_1$routes, by = "route_id") |>
  select(shape_id, route_id, route_short_name, route_long_name) |>
  distinct() |> 
  mutate(name = paste(route_short_name, route_long_name, sep = ": "))

shapes_sf = tidytransit::shapes_as_sf(gtfs_1$shapes) |>
  left_join(route_name, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "name", legend = TRUE, layer.name="Route name")

GTFS filtered for route long name Terminal, Rodoviário (partial match) 

route_name = gtfs_terminal$trips |>
  left_join(gtfs_terminal$routes, by = "route_id") |>
  select(shape_id, route_id, route_short_name, route_long_name) |>
  distinct() |> 
  mutate(name = paste(route_short_name, route_long_name, sep = ": "))

shapes_sf = tidytransit::shapes_as_sf(gtfs_terminal$shapes) |>
  left_join(route_name, by = "shape_id")
mapview::mapview(shapes_sf, zcol = "name", legend = TRUE, layer.name="Route name")