old_warn <- getOption("warn")
options(warn = -1)

# Cargar Excel
EncuestaOrigenDestino <- read_excel("EncuestaOrigenDestino.xlsx", sheet = "Sheet4", col_types = c("date", "numeric", "text", "text", "text", "text", "text", "text", "numeric", "numeric", "text", "text", "numeric", "numeric", "numeric", "numeric", "text", "numeric", "numeric", "numeric", "numeric", "numeric", "numeric", "text", "numeric", "numeric", "numeric", "numeric"))
## New names:
## • `MUNICIPIO` -> `MUNICIPIO...7`
## • `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN` ->
##   `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN...8`
## • `MUNICIPIO` -> `MUNICIPIO...11`
## • `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN` ->
##   `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN...12`
## • `OTRO¿CUÁL?` -> `OTRO¿CUÁL?...17`
## • `OTRO¿CUÁL?` -> `OTRO¿CUÁL?...24`
options(warn = old_warn)
# Cargar Shape
comunas=shapefile("comunas_up.shp")

1. Viajes Origen

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Tabla viajes origen
viajes_origen=table(EncuestaOrigenDestino$comunaorigen)

as.numeric(names(viajes_origen))
##  [1]  0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22
orden=match(comunas$comuna,as.numeric(names(viajes_origen)))
orden
##  [1]  7  5  6  8  9 10 22 14  2  4 20 13 11 21 17 16 18 19 15 12  3 23
comunas$viajes=as.numeric(viajes_origen[orden])
comunas
## class       : SpatialPolygonsDataFrame 
## features    : 22 
## extent      : -76.59284, -76.46125, 3.331802, 3.505871  (xmin, xmax, ymin, ymax)
## crs         : +proj=longlat +datum=WGS84 +no_defs 
## variables   : 8
## names       : casos, comuna,   nombre, zona_recol,         area, perimetro, covid, viajes 
## min values  :     1,      1, Comuna 1,         NA,  2329397.941,  7983.949,     5,    348 
## max values  :    39,     22, Comuna 9,         NA, 12555929.024, 26480.361,    54,   2968
# Mapa viajes origen (leaflet)
pal <- colorNumeric(palette = "YlOrRd", domain = comunas$viajes)
l <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal(viajes),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", viajes, " viajes origen"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal, values = ~viajes_origen, opacity = 0.6)
l

2. Viajes Origen Bicicleta

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Filtrar datos viajes origen bicicleta TIPODEVEHÍCULO = 1
bicis_origen <- EncuestaOrigenDestino %>%
  filter(TIPODEVEHÍCULO == "1") %>%
  count(comunaorigen)

# Convertir nombres a numérico
bicis_origen$comunaorigen <- as.numeric(bicis_origen$comunaorigen)

# Ordenar datos basado en comunas
orden_bicis <- match(comunas$comuna, bicis_origen$comunaorigen)
bicis_origen <- bicis_origen[order(bicis_origen$comunaorigen), ]

# Añadir viajes bicicleta al dataframe comunas
comunas$origen_bicis <- as.numeric(bicis_origen$n[orden_bicis])
comunas$origen_bicis[is.na(comunas$origen_bicis)] <- 0  

# Mapa viajes origen bicicleta (leaflet)
pal_bici <- colorNumeric(palette = "Blues", domain = comunas$origen_bicis)
map_bici <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal_bici(origen_bicis),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", origen_bicis, " viajes orig bicis"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal_bici, values = ~origen_bicis, opacity = 0.6)
map_bici

3. Viajes Origen Moto

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Filtrar datos viajes origen motos TIPODEVEHÍCULO = 2
motos_origen <- EncuestaOrigenDestino %>%
  filter(TIPODEVEHÍCULO == "2") %>%
  count(comunaorigen)

# Convertir nombres a numérico
motos_origen$comunaorigen <- as.numeric(motos_origen$comunaorigen)

# Ordenar datos basado en comunas
orden_motos <- match(comunas$comuna, motos_origen$comunaorigen)
motos_origen <- motos_origen[order(motos_origen$comunaorigen), ]

# Añadir viajes moto al dataframe comunas
comunas$origen_motos <- as.numeric(motos_origen$n[orden_motos])
comunas$origen_motos[is.na(comunas$origen_motos)] <- 0  

# Mapa viajes origen motos (leaflet)
pal_moto <- colorNumeric(palette = "Reds", domain = comunas$origen_motos)
map_moto <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal_moto(origen_motos),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", origen_motos, " viajes orig motos"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal_moto, values = ~origen_motos, opacity = 0.6)
map_moto

4. Viajes Origen Carro

# Filtrar datos viajes origen carros TIPODEVEHÍCULO = 3
carros_origen <- EncuestaOrigenDestino %>%
  filter(TIPODEVEHÍCULO == "3") %>%
  count(comunaorigen)

# Convertir nombres a numérico
carros_origen$comunaorigen <- as.numeric(carros_origen$comunaorigen)

# Ordenar datos basado en comunas
orden_carros <- match(comunas$comuna, carros_origen$comunaorigen)
carros_origen <- carros_origen[order(carros_origen$comunaorigen), ]

# Añadir viajes carros al dataframe comunas
comunas$origen_carros <- as.numeric(carros_origen$n[orden_carros])
comunas$origen_carros[is.na(comunas$origen_carros)] <- 0  

# Mapa viajes origen carros (leaflet)
pal_carro <- colorNumeric(palette = "Greens", domain = comunas$origen_carros)
map_carro <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal_carro(origen_carros),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", origen_carros, " viajes orig carros"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal_carro, values = ~origen_carros, opacity = 0.6)
map_carro

5. Viajes Destino

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Tabla viajes destino
viajes_destino=table(EncuestaOrigenDestino$comunadestino)

as.numeric(names(viajes_destino))
##  [1]  0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22
orden=match(comunas$comuna,as.numeric(names(viajes_destino)))
orden
##  [1]  7  5  6  8  9 10 22 14  2  4 20 13 11 21 17 16 18 19 15 12  3 23
comunas$viajes=as.numeric(viajes_destino[orden])
comunas
## class       : SpatialPolygonsDataFrame 
## features    : 22 
## extent      : -76.59284, -76.46125, 3.331802, 3.505871  (xmin, xmax, ymin, ymax)
## crs         : +proj=longlat +datum=WGS84 +no_defs 
## variables   : 11
## names       : casos, comuna,   nombre, zona_recol,         area, perimetro, covid, viajes, origen_bicis, origen_motos, origen_carros 
## min values  :     1,      1, Comuna 1,         NA,  2329397.941,  7983.949,     5,    230,           14,          176,           123 
## max values  :    39,     22, Comuna 9,         NA, 12555929.024, 26480.361,    54,   4810,          159,         1414,          1165
# Mapa viajes destino (leaflet)
pal <- colorNumeric(palette = "YlOrRd", domain = comunas$viajes)
l <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal(viajes),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", viajes, " viajes destino"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal, values = ~viajes_destino, opacity = 0.6)
l

6. Viajes Destino Bicicleta

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Filtrar datos viajes destino bicicletas TIPODEVEHÍCULO = 1
destino_bicis <- EncuestaOrigenDestino %>%
  filter(TIPODEVEHÍCULO == "1") %>%
  count(comunadestino)

# Convertir nombres a numérico
destino_bicis$comunadestino <- as.numeric(destino_bicis$comunadestino)

# Ordenar datos basado en comunas
orden_destino_bicis <- match(comunas$comuna, destino_bicis$comunadestino)
destino_bicis <- destino_bicis[order(destino_bicis$comunadestino), ]

# Añadir viajes destino bicicleta al dataframe comunas
comunas$destino_bicis <- as.numeric(destino_bicis$n[orden_destino_bicis])
comunas$destino_bicis[is.na(comunas$destino_bicis)] <- 0  

# Mapa viajes destino bicicleta (leaflet)
pal_destino_bici <- colorNumeric(palette = "Purples", domain = comunas$destino_bicis)
map_destino_bici <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal_destino_bici(destino_bicis),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", destino_bicis, " viajes destino en bicicleta"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal_destino_bici, values = ~destino_bicis, opacity = 0.6)
map_destino_bici

7. Viajes Destino Moto

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Filtrar datos viajes destino motos TIPODEVEHÍCULO = 2
destino_motos <- EncuestaOrigenDestino %>%
  filter(TIPODEVEHÍCULO == "2") %>%
  count(comunadestino)

# Convertir nombres a numérico
destino_motos$comunadestino <- as.numeric(destino_motos$comunadestino)

# Ordenar datos basado en comunas
orden_destino_motos <- match(comunas$comuna, destino_motos$comunadestino)
destino_motos <- destino_motos[order(destino_motos$comunadestino), ]

# Añadir viajes destino motos al dataframe comunas
comunas$destino_motos <- as.numeric(destino_motos$n[orden_destino_motos])
comunas$destino_motos[is.na(comunas$destino_motos)] <- 0  

# Mapa viajes destino motos (leaflet)
pal_destino_moto <- colorNumeric(palette = "Oranges", domain = comunas$destino_motos)
map_destino_moto <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal_destino_moto(destino_motos),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", destino_motos, " viajes destino en moto"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal_destino_moto, values = ~destino_motos, opacity = 0.6)
map_destino_moto

8. Viajes Destino Carro

Pase el mouse por encima de las comunas para ver la cantidad de viajes realizados.

# Filtrar datos viajes destino carros TIPODEVEHÍCULO = 3
destino_carros <- EncuestaOrigenDestino %>%
  filter(TIPODEVEHÍCULO == "3") %>%
  count(comunadestino)

# Convertir nombres a numérico
destino_carros$comunadestino <- as.numeric(destino_carros$comunadestino)

# Ordenar datos basado en comunas
orden_destino_carros <- match(comunas$comuna, destino_carros$comunadestino)
destino_carros <- destino_carros[order(destino_carros$comunadestino), ]

# Añadir viajes destino carro al dataframe comunas
comunas$destino_carros <- as.numeric(destino_carros$n[orden_destino_carros])
comunas$destino_carros[is.na(comunas$destino_carros)] <- 0  

# Mapa viajes destino carros (leaflet)
pal_destino_carro <- colorNumeric(palette = "Greys", domain = comunas$destino_carros)
map_destino_carro <- leaflet(comunas) %>%
  addTiles() %>%
  addPolygons(color = "white", fillColor = ~pal_destino_carro(destino_carros),
              fillOpacity = 0.8,
              label = ~paste(comunas$nombre, ": ", destino_carros, " viajes destino en carro"),
              labelOptions = labelOptions(direction = "auto", noHide = F, textOnly = F)) %>%
  addLegend(pal = pal_destino_carro, values = ~destino_carros, opacity = 0.6)
map_destino_carro