tests/testthat/test_append_height.R

In gtfs2gps: Converting Transport Data from GTFS Format to GPS-Like Records

test_that("append_height", {
  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("shape804-I") %>%
    filter_single_trip() %>% 
    remove_invalid()
  
  fortaleza_gps <- gtfs2gps(gtfs, spatial_resolution = 500) %>% append_height(srtmfile)

  mysum <- sum(fortaleza_gps$height)

  expect_equal(mysum, 3314, 0.05)
  
  #fort <- fortaleza_gps[which(is.na(fortaleza_gps$height)),]
  
  #fott <- gps_as_sf(fort)
  #sf::write_sf(fott, "fort-na.shp")  
})