Add a column named height to GPS data using a tif data as reference.
append_height(gps, heightfile)
A GPS data created from gtfs2gps().
The pathname of a tif file with height data.
The GPS data with a new column named height.
if (FALSE) {
# this example takes more than 10s to run
library(magrittr)
fortaleza <- system.file("extdata/fortaleza.zip", package = "gtfs2gps")
srtmfile <- system.file("extdata/fortaleza-srtm.tif", package = "gtfs2gps")
gtfs <- read_gtfs(fortaleza) %>%
gtfstools::filter_by_shape_id("shape836-I") %>%
filter_single_trip()
fortaleza_gps <- gtfs2gps(gtfs, spatial_resolution = 500) %>% append_height(srtmfile)
}