Katso lisää: http://data.markuskainu.fi/opendata/kaupunkipyorat/
# 1. Ladataan aineisto ja valitaan päiväksi kesäkuun 1.
library(dplyr)
library(sf)
library(readr)
stations <- read_csv2("http://data.markuskainu.fi/opendata/kaupunkipyorat/tellingit_2017-2020.csv.gz")
kesakuu1 <- stations %>%
# valitaan kesäkuun 1. päivä kultakin vuodelta
filter(grepl("06/01", time))
# Lasketaan telineiden määrät vuosittain
kesakuu1 %>%
count(year)
## # A tibble: 4 x 2
## year n
## <dbl> <int>
## 1 2017 147
## 2 2018 220
## 3 2019 356
## 4 2020 450
# poistetaan testit / remove tests
kesakuu1 %>%
filter(x <= 30, y >= 55) %>%
sf::st_as_sf(coords = c(3,4)) -> spatdat
library(leaflet)
leaflet() %>%
leaflet::addProviderTiles(providers$CartoDB.Positron) %>%
addCircleMarkers(data = spatdat[spatdat$year == "2017",],
color = "blue",
label = ~paste0(name," (", id, ")"),
group = "2017") %>%
addCircleMarkers(data = spatdat[spatdat$year == "2018",],
color = "red",
label = ~paste0(name," (", id, ")"),
group = "2018") %>%
addCircleMarkers(data = spatdat[spatdat$year == "2019",],
color = "orange",
label = ~paste0(name," (", id, ")"),
group = "2019") %>%
addCircleMarkers(data = spatdat[spatdat$year == "2020",],
color = "purple",
label = ~paste0(name," (", id, ")"),
group = "2020") %>%
# Layers control
addLayersControl(
baseGroups = c("2017", "2018", "2019", "2020"),
options = layersControlOptions(collapsed = FALSE)
)