Leaflet showcase

Ejemplo de leaflet

Cargo paquetes y data.

library(dplyr)
library(sf)
library(lubridate)
library(here)
'%!in%' <- function(x,y)!('%in%'(x,y))
library(leaflet)
library(RColorBrewer)

load(here::here("data", "df_gps_152.RData"))

Tomo las coordenadas de una de las cabeceras.

cabecera_lat_long <- c(long = -58.51257, lat = -34.51122)

cabecera <- data.frame(long = cabecera_lat_long["long"],
                       lat = cabecera_lat_long["lat"]) %>% 
    st_as_sf(coords = c("long", "lat"),
             crs = 4326) %>% 
    st_transform(crs = 5347)

Le lavo la cara a los datos y armo las columnas necesarias.

df_gps_152_staged <- df_gps_152 %>%
    as.data.frame() %>% 
    filter(interno %!in% c("3", "99", "117", "72", "5", "26", "122")) %>%
    filter(!(latitude == 0 & longitude == 0)) %>% 
    distinct(interno, date_time, .keep_all = TRUE)

df_gps_152_5347 <- df_gps_152_staged %>%
    st_as_sf(coords = c("longitude", "latitude"),
             crs = 4326,
             remove = FALSE) %>%
    st_transform(crs = 5347)

df_gps_152_5347 <- cbind(df_gps_152_5347,
                         df_gps_152_5347 %>% st_coordinates())

df_gps_152_5347 <- df_gps_152_5347 %>% 
    mutate(dist_cabecera = st_distance(geometry, cabecera) %>% 
               round() %>% 
               as.numeric())

# armo columnas extras
df_gps_152_5347 <- df_gps_152_5347 %>% 
    group_by(interno) %>% 
    arrange(interno, date_time) %>% 
    mutate(se_aleja_de_cabecera = dist_cabecera > lag(dist_cabecera),
           en_cabecera = dist_cabecera <= 400,
           delta_x_m = dist_cabecera - (dist_cabecera %>% lag()),
           delta_t_s = difftime(date_time,
                                date_time %>% lag(),
                                units = "secs") %>% as.numeric()) %>% 
    mutate(v_kmph = (delta_x_m/delta_t_s*3.6) %>% abs() %>% round(2)) %>% 
    mutate(detenido = v_kmph <= 2)

A mapear ! Yupi.

Primero: obtengo subset de datos a mapear.

data <- df_gps_152_5347 %>%
    filter(interno == "95") %>%
    arrange(interno, date_time) %>% 
    mutate(rn = row_number()) %>% 
    mutate(popup_text = paste0("v: ", v_kmph, " km/h <br/> delta_t_s: ", delta_t_s, " s <br/> delta_x_m: ", delta_x_m, " m <br/> date_time: ", date_time, "."))

Lo segundo: ¡paletas!

pal <- colorFactor(
    palette = topo.colors(2),
    domain = data$se_aleja_de_cabecera)
pal2 <- colorFactor(
    palette = rainbow(2),
    domain = data$en_cabecera)
pal3 <- colorNumeric(
    palette = "viridis",
    domain = data$v_kmph)
pal4 <- colorFactor(
    palette = RColorBrewer::brewer.pal(3, "Set1"),
    domain = data$detenido)

Y finalmente: el mapa.

leaflet(data = data %>%
            dplyr::filter(rn < 60),
        width = "100%") %>%  
    addTiles() %>% 
    addMarkers(lng = cabecera_lat_long["long"],
               lat = cabecera_lat_long["lat"]) %>%
    addCircleMarkers(lng = ~longitude,
                     lat = ~latitude,
                     radius = 5,
                     color = ~pal(se_aleja_de_cabecera),
                     fillOpacity = 0.7,
                     label = ~(rn %>% as.character()),
                     labelOptions = labelOptions(noHide = T, opacity = 0.8),
                     group = "Alejamiento cabecera") %>% 
    addLegend("bottomright", pal = pal, values = ~se_aleja_de_cabecera, title = "Alejamiento cabecera", opacity = 1) %>% 
    addCircleMarkers(lng = ~longitude,
                     lat = ~latitude,
                     radius = 5,
                     color = ~pal2(en_cabecera),
                     fillOpacity = 0.7,
                     label = ~(rn %>% as.character()),
                     labelOptions = labelOptions(noHide = T, opacity = 0.8),
                     group = "En cabecera") %>%
    addLegend("bottomright", pal = pal2, values = ~en_cabecera, title = "En cabecera", opacity = 1) %>% 
    addCircleMarkers(lng = ~longitude,
                     lat = ~latitude,
                     radius = ~v_kmph,
                     color = ~pal3(v_kmph),
                     fillOpacity = 0.7,
                     label = ~(rn %>% as.character()),
                     labelOptions = labelOptions(noHide = T, opacity = 0.8),
                     popup = ~popup_text,
                     group = "Velocidad") %>% 
    addLegend("bottomright", pal = pal3, values = ~v_kmph, title = "Velocidad", opacity = 1) %>% 
    addCircleMarkers(lng = ~longitude,
                     lat = ~latitude,
                     radius = 5,
                     color = ~pal4(detenido),
                     fillOpacity = 1,
                     label = ~(rn %>% as.character()),
                     labelOptions = labelOptions(noHide = T, opacity = 0.8),
                     group = "Detenido") %>%
    addLegend("bottomright", pal = pal4, values = ~detenido, title = "Detenido", opacity = 1) %>%
    addLayersControl(overlayGroups = c("Alejamiento cabecera", "En cabecera", "Velocidad", "Detenido"),
                     options = layersControlOptions(collapsed = FALSE),
                     position = "bottomleft") %>%
    hideGroup(c("En cabecera", "Velocidad", "Detenido")) %>%
    addMeasure(primaryLengthUnit = "meters", position = "bottomleft") %>%
    leaflet.extras::addFullscreenControl()

Nota del autor: Se apreciará el no-juicio estético de un mapa que ha sido creado para analizar y testear herramientas.

Nota del autor 2: “En mi PC funcionaba”. En realidad se veía más grande. Leaflet arreglame el full-screen.

Anyways: Carpe díem.