opendatasoft/make_map.R at main · stevecrawshaw/opendatasoft · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
p <- c("tidyverse", "here", "data.table", "lubridate", "sf", "leaflet", "leafem", "geojsonio", "glue", "tmap", "tmaptools")
library(xfun)
pkg_attach2(p)
aq_monitors <- read_sf("https://opendata.bristol.gov.uk/explore/dataset/air-quality-monitoring-sites/download/?format=geojson&timezone=Europe/London") %>%
  filter(current == "True")

wards <- read_sf("https://opendata.bristol.gov.uk/explore/dataset/wards/download/?format=geojson&timezone=Europe/London&lang=en")

tubes_sf <- read_sf("https://opendata.bristol.gov.uk/explore/dataset/no2-diffusion-tube-data/download/?format=geojson&disjunctive.location=true&q=year+%3E%3D+2019&timezone=Europe/London&lang=en")

aqma <- read_sf("https://opendata.bristol.gov.uk/explore/dataset/air-quality-management-areas/download/?format=geojson&timezone=Europe/London&lang=en")

year <- 2020

tubes_cy_sf <- tubes_sf %>%
    filter(year == {{year}})

tubes_cy_36_sf <- tubes_cy_sf %>%
    filter(conc_ugm3 >= 36)


bris_box <- st_bbox(wards)
aqma_box <- st_bbox(aqma)

 osm_bris <- readRDS("data/osm_bris.rds")

tmap_mode("view")

lp <- leaflet.providers::get_providers()
tm_basemap(server = "OpenStreetMap.Mapnik") +
    # tm_shape(osm_bris) +
    # tm_rgb() +
    tm_shape(aqma, bbox = "aqma_box") +
    tm_borders(col = "red", lwd = 1, lty = "solid") +
    tm_shape(tubes_cy_36_sf) +
    tm_symbols(size = 5, col = "conc_ugm3", scale = .5, n = 2, palette = "-RdYlGn")


#
# path <- "S:/SUSTAIN/Sustain-Common/SCCCS/write_gis_r/Mapping/os_open/local/opmplc_essh_st/OS OpenMap Local (ESRI Shape File) ST/data"
#
# roads_files <- glue("{path}/ST_Road.shp")
#
# roads <- st_read(roads_files)
#
cropbox <- c(ymax = 184000, ymin = 166000, xmin = 349000, xmax = 366000)
#
#
# bristol_roads <- roads %>% st_crop(y = cropbox)

bristol_box <- st_bbox(obj = wards_gj)


brisraster <- read_osm(bristol_box)


aqm <-   aq_monitors %>%
  as.data.frame() %>%
  separate(geo_point_2d, into = c("lat", "long"), sep = ", ") %>%
  mutate(lat = as.numeric(str_sub(lat, 3)),
         long = as.numeric(str_sub(long, 1, -2))) %>%
  mutate(color = case_when(
    instrumenttype == "Diffusion Tube" ~ "blue",
    instrumenttype == "Continuous (Reference)" ~ "red",
    instrumenttype == "Continuous (Indicative)" ~ "green"))

livenox <- aq_monitors %>%
  filter(str_detect(pollutants, "NOX")) %>%
  nrow()

leaflet(data = aqm) %>%
  setView(lng = -2.6, lat = 51.46, zoom = 12) %>%
  addTiles(urlTemplate = "http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png")  %>%
  addGeoJSON(wards_gj) %>%
  addCircleMarkers(~long,
                   ~lat,
                   popup = ~paste0(location,
                                   "<br/>",
                                   "Pollutants: ", pollutants,
                                   "<br/>",
                                   "Location class: ", locationclass),
                   label = ~ as.character(siteid),
                   color = ~ color) %>%
  addLegend(position = "bottomleft", colors = ~unique(color), labels = ~unique(instrumenttype))

#-------------------------CHLOROPLETH MAP---------------------------
dropno2 <- c("easting", "northing")
df_no2_ods <- fread("https://opendata.bristol.gov.uk/explore/dataset/no2-diffusion-tube-data/download/?format=csv&timezone=Europe/London&use_labels_for_header=false", drop = dropno2)

curr_tubes <- df_no2_ods %>%
  filter(year == max(year)) %>%
  rename(Site_ID = siteid, NO2_ugm3 = conc_ugm3) %>%
  separate(geo_point_2d, into = c("lat", "long"), sep = ",") %>%
  mutate(lat = as.numeric(lat),
         long = as.numeric(long))

getColor <- function(curr_tubes) {
  sapply(curr_tubes$NO2_ugm3, function(NO2_ugm3) {
    if(NO2_ugm3 <= 40) {
      "green"
    } else if(NO2_ugm3 <= 60) {
      "red"
    } else {
      "purple"
    } })
}

icons <- awesomeIcons(
  icon = 'circle',
  iconColor = 'black',
  library = 'fa',
  markerColor = getColor(curr_tubes)
)

#icon = 'ios-close',
#no2 <- expression("NO"[2] * "concentration in ")

curr_tubes %>%
  leaflet() %>%
  setView(lng = -2.6, lat = 51.46, zoom = 12) %>%
  addTiles(urlTemplate = "http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png")  %>%
  addAwesomeMarkers(~long,
                    ~lat,
                    popup = ~paste0(location, "<br/>",
                                    "NO2 concentration in ", year, ": ",
                                    "<b>", NO2_ugm3, "<b/>", " ugm-3"),
                    label = ~as.character(Site_ID),
                    icon = icons,
                    clusterOptions = markerClusterOptions()) %>%
  addLegend(position = "bottomleft",
            colors = c("green", "red", "purple"),
            labels = c("compliant", "breaches annual objective", "possible breach hourly objective"))