uk-cities.R
In pct: Propensity to Cycle Tool

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

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

## ----setup, message=FALSE, warning=FALSE--------------------------------------
#  library(pct)
#  library(dplyr)
#  library(sf)
#  library(tmap)

## -----------------------------------------------------------------------------
#  region_name = "wales"
#  zones_all = get_pct_zones(region_name)

## ----include=TRUE-------------------------------------------------------------
#  unique(zones_all$lad_name)

## -----------------------------------------------------------------------------
#  la_name = "Cardiff"

## ----include=TRUE, fig.align='center', fig.width = 7, fig.height = 6----------
#  zones = zones_all %>%
#    filter(lad_name == la_name)
#  plot(zones["bicycle"])

## -----------------------------------------------------------------------------
#  rnet_all = pct::get_pct_rnet(region_name)
#  rnet = rnet_all[zones, ]

## ----include=TRUE, fig.align='center', fig.width = 7, fig.height = 6----------
#  plot(zones$geometry)
#  plot(rnet["dutch_slc"], add = TRUE)

## -----------------------------------------------------------------------------
#  pct_zones_rnet = function(la_name, region_name = "devon") {
#    zones_all = pct::get_pct_zones(region_name)
#    zones = zones_all %>%
#    filter(lad_name == la_name)
#    plot(zones["bicycle"])
#    rnet_all = pct::get_pct_rnet(region_name)
#    rnet = rnet_all[zones, ]
#    plot(zones$geometry)
#    plot(rnet["dutch_slc"], add = TRUE)
#    list(zones = zones, rnet = rnet)
#  }
#  
#  plymouth_results = pct_zones_rnet(la_name = "Plymouth")
#  exeter_results = pct_zones_rnet(la_name = "Exeter")

## -----------------------------------------------------------------------------
#  tmap_mode("view")
#  b = c(0, 100, 200, 500, 1000)
#  m1 = tm_shape(plymouth_results$rnet) +
#    tm_lines("dutch_slc", breaks = b, palette = "viridis", lwd = 2) +
#    tm_scale_bar()
#  m2 = tm_shape(exeter_results$rnet) +
#    tm_lines("dutch_slc", breaks = b, palette = "viridis", lwd = 2) +
#    tm_scale_bar()
#  tmap_arrange(m1, m2, ncol = 2)