# NOT RUN {
library(gtfs2gps)
library(dplyr)
fortaleza <- read_gtfs(system.file("extdata/fortaleza.zip", package = "gtfs2gps"))
srtmfile <- system.file("extdata/fortaleza-srtm.tif", package="gtfs2gps")
subset <- fortaleza %>%
filter_week_days() %>%
filter_single_trip() %>%
filter_by_shape_id(c("shape804-I", "shape806-I"))
for_gps <- gtfs2gps(subset)
for_gps_sf_points <- gps_as_sfpoints(for_gps) # without height
for_gps <- append_height(for_gps, srtmfile)
for_gps_sf_points <- gps_as_sfpoints(for_gps) # with height
# }
Run the code above in your browser using DataLab