tests/tests_marcus/test_elevation_r5r.R

In r5r: Rapid Realistic Routing with 'R5'

# options(java.parameters = '-Xmx2G')

# library(r5r)
devtools::load_all(".")
library(data.table)
library(tidyverse)

# build transport network
data_path <- system.file("extdata/poa", package = "r5r")
r5r_core <- setup_r5(data_path = data_path, verbose = FALSE,
                     overwrite = TRUE,
                     temp_dir = FALSE,
                     elevation = "tobler")

tn <- transit_network_to_sf(r5r_core)

mapview::mapview(tn$routes, zcol = "mode")

# load origin/destination points
points <- read.csv(file.path(data_path, "poa_hexgrid.csv"))

# r5r_core$setTravelTimesBreakdown(FALSE)
ttm_r5r_walk <- travel_time_matrix(r5r_core,
                            origins = points,
                            destinations = points,
                            mode = c("WALK"),
                            max_trip_duration = 60,
                            max_walk_dist = Inf,
                            time_window = 1,
                            percentiles = c(50),
                            verbose = FALSE,
                            progress = TRUE)


ttm_r5r_bike <- travel_time_matrix(r5r_core,
                                    origins = points,
                                    destinations = points,
                                    breakdown = FALSE,
                                    mode = c("BICYCLE"),
                                    max_trip_duration = 60,
                                    max_walk_dist = Inf,
                                    time_window = 1,
                                    percentiles = c(50),
                                    verbose = FALSE,
                                    progress = TRUE)

ttm_r5r_walk$scenario <- "r5r"
ttm_r5r_walk$mode <- "walk"
ttm_r5r_bike$scenario <- "r5r"
ttm_r5r_bike$mode <- "bike"

write_csv(ttm_r5r_walk, here::here("elevation", "walk_r5r.csv"))
write_csv(ttm_r5r_bike, here::here("elevation", "bike_r5r.csv"))