Add a column named height to GPS data using a tif data as reference.

append_height(gps, heightfile)

Arguments

gps

A GPS data created from gtfs2gps().

heightfile

The pathname of a tif file with height data.

Value

The GPS data with a new column named height.

Examples

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)
}