library(tidyverse)

## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──

## ✔ ggplot2 3.3.5     ✔ purrr   0.3.4
## ✔ tibble  3.1.6     ✔ dplyr   1.0.8
## ✔ tidyr   1.2.0     ✔ stringr 1.4.0
## ✔ readr   2.1.2     ✔ forcats 0.5.1

## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()

library(tmap)
tmap_mode("view")

## tmap mode set to interactive viewing

rnet = readRDS("~/itsleeds/turas/outputs/2022-02-26/rds/kildare/rnet_work_quietest_30pc.Rds")
maynooth = tmaptools::geocode_OSM(q = "maynooth", as.sf = TRUE)
maynooth_5km = sf::st_buffer(maynooth, dist = 5000)
qtm(maynooth_5km)

rnet = rnet[maynooth_5km, ]

# colorspace::choose_palette()
pal = colorspace::sequential_hcl(n = 4, h = c(5, -162), c = c(80, NA, 65), l = c(30, 77), power = c(1, 1))
colorspace::demoplot(pal)

dal = colorspace::deutan(pal)
colorspace::demoplot(dal)

quietness_breaks = c(0, 35, 70, 85, 100)
# brks = c(0, 10, 50, 100, 500, 1000, 5000)
s = c(
  `Grey basemap` = "CartoDB.Positron",
  `Coloured basemap` = "Esri.WorldTopoMap",
  `OSM existing cycle provision (CyclOSM)` = "https://b.tile-cyclosm.openstreetmap.fr/cyclosm/{z}/{x}/{y}.png",
  `Satellite image` = "https://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}'"
)

  m2 = tm_shape(rnet) +
    tm_lines(
      legend.col.reverse = TRUE,
      col = "Quietness",
      lwd = "lwd",
      popup.vars = names(rnet)[1:4],
      palette = pal,
      scale = 15,
      breaks = quietness_breaks
    ) +
    tm_basemap(server = s) +
    tm_scale_bar()

# +
#   tm_view(bbox = tmaptools::bb(bb, 2))
m2

## Legend for line widths not available in view mode.