pct-international.R
In pct: Propensity to Cycle Tool

## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## -----------------------------------------------------------------------------
knitr::opts_chunk$set(eval = FALSE)

## ----libraries, message=FALSE-------------------------------------------------
#  # devtools::install_github("ATFutures/geoplumber")
#  # require("geojsonsf")
#  library(pct)

## -----------------------------------------------------------------------------
#  head(santiago_od)

## -----------------------------------------------------------------------------
#  sf:::print.sf(santiago_zones)
#  plot(santiago_zones)

## ---- warning=FALSE-----------------------------------------------------------
#  desire_lines = stplanr::od2line(flow = santiago_od, zones = santiago_zones)

## ---- out.width="100%"--------------------------------------------------------
#  plot(santiago_zones$geometry)
#  plot(santiago_lines["pcycle"], lwd = santiago_lines$n / 3, add = TRUE)
#  # gj = geojsonsf::sf_geojson(santiago_lines)
#  # path = file.path(tempdir(), "dl.geojson")
#  # write(gj, path)
#  # html_map = geoplumber::gp_map(path, browse_map = FALSE)
#  # htmltools::includeHTML(html_map)

## -----------------------------------------------------------------------------
#  desire_lines$hilliness = 0

## -----------------------------------------------------------------------------
#  desire_lines$distance = as.numeric(sf::st_length(desire_lines))

## -----------------------------------------------------------------------------
#  desire_lines$godutch_pcycle = uptake_pct_godutch(distance = desire_lines$distance, gradient = 0)

## -----------------------------------------------------------------------------
#  cor(x = desire_lines$pcycle, y = desire_lines$godutch_pcycle)
#  plot(x = desire_lines$pcycle, y = desire_lines$godutch_pcycle)
#  plot(x = desire_lines$distance, y = desire_lines$godutch_pcycle, ylim = c(0, 1))

## ---- out.width="50%", fig.show='hold'----------------------------------------
#  library(leaflet)
#  leaflet(width = "100%") %>%
#    addTiles() %>%
#    addPolylines(data = desire_lines, weight = desire_lines$pcycle * 5)
#  leaflet(width = "100%") %>%
#    addTiles() %>%
#    addPolylines(data = desire_lines, weight = desire_lines$godutch_pcycle * 5)

## ---- eval=FALSE--------------------------------------------------------------
#  santiago_routes_cs = stplanr::line2route(desire_lines)
#  # > 10 % out of 200 distances calculated
#  # > 20 % out of 200 distances calculated
#  # > 30 % out of 200 distances calculated
#  # > 40 % out of 200 distances calculated
#  # > 50 % out of 200 distances calculated
#  # > 60 % out of 200 distances calculated
#  # > 70 % out of 200 distances calculated
#  # > 80 % out of 200 distances calculated
#  # > 90 % out of 200 distances calculated
#  # > 100 % out of 200 distances calculated
#  # > Warning message:
#  # > In value[[3L]](cond) : Fail for line number 32

## -----------------------------------------------------------------------------
#  leaflet() %>%
#    addTiles() %>%
#    addPolylines(data = santiago_routes_cs[32, ])

## -----------------------------------------------------------------------------
#  routes = sf::st_sf(
#    cbind(sf::st_drop_geometry(santiago_routes_cs),
#    sf::st_drop_geometry(desire_lines)),
#    geometry = santiago_routes_cs$geometry
#  )

## -----------------------------------------------------------------------------
#  routes$godutch_slc = round(routes$godutch_pcycle * routes$all)
#  rnet = stplanr::overline2(routes, "godutch_slc")
#  plot(rnet, lwd = rnet$godutch_slc / mean(rnet$godutch_slc))
#  # library(tmap)
#  # tmap_mode("view")
#  #   tm_shape(rnet) +
#  #   tm_lines(lwd = "godutch_slc", scale = 9)

## ---- eval=FALSE--------------------------------------------------------------
#  route_segments_1_5 = route(l = desire_lines[1:5, ], route_fun = cyclestreets::journey)
#  mapview::mapview(route_segments_1_5)