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