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

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

