Datos de Covid por colonia disponibles en el portal de datos de la CDMX.

En el presente documento se presentan las principales bases de datos disponibles por parte del portal de datos de la CDMX a nivel colonia de casos activos, totales y por colonias prioritarias en la CDMX.

# Librerias
library(tidyverse)
library(sf)
library(leaflet)

Bases de datos

Las bases de datos son las siguientes:

• Mapa de todas las colonias de la CDMX.

• Base histórica de casos a nivel colonia, con datos al 15 de Agosto del 2020 (fecha más reciente encontrada).

• Mapa de Municipios de la CDMX.

• Base de datos de las colonias de acción prioritaria de la CDMX.

# Bases de datos
shp <- st_read("01_Datos/coloniascdmx.geojson", quiet = TRUE) %>% 
  rename(clave_colonia = cve_col)

bd <- read_csv("https://archivo.datos.cdmx.gob.mx/sinave-colonias/200815_base_historica.csv")   
mpios <- st_read("https://raw.githubusercontent.com/JuveCampos/Shapes_Resiliencia_CDMX_CIDE/master/geojsons/Division%20Politica/mpios2.geojson", quiet = TRUE) %>% filter(CVE_ENT == "09")

# Fuente: https://datos.cdmx.gob.mx/explore/dataset/colonias-de-atencion-prioritaria-covid-kioscos/information/?sort=-id
col_prioritarias <- st_read("01_Datos/colonias-de-atencion-prioritaria-covid-kioscos/colonias-de-atencion-prioritaria-covid-kioscos.shp", quiet = TRUE)

# Juntamos la información 
mapa <- merge(shp, bd, by = "clave_colonia", all.x = TRUE) %>% 
  mutate(total_200815 = ifelse(is.na(total_200815), 
                               yes = 0,
                               no = total_200815))

sin_datos <- mapa %>% 
  filter(total_200815 == 0)

con_datos <- mapa %>% 
  filter(total_200815 != 0)

# Sacamos el mapa ----

label = lapply(paste0(con_datos$clave_colonia, ": ", con_datos$nombre, "<br>", 
                               con_datos$alcaldia.x, ", ", con_datos$entidad), htmltools::HTML)

popup = paste0("<b>Colonia: </b>", con_datos$clave_colonia, ": ", con_datos$nombre, "<br>", 
               "<b>Datos al 01 de Julio del 2020: </b>", con_datos$total_200815, " casos.")

# Paleta de colores ----
pal <- colorNumeric(palette = "magma", rev = TRUE,
                    domain = con_datos$total_200815)

# Codigo del mapa ----
leaflet(con_datos) %>% 
  addProviderTiles("CartoDB.Positron") %>% 
  addPolygons(fillColor = pal(con_datos$total_200815), 
              label = label, 
              popup = popup,
              fillOpacity = 1, 
              color = "black", 
              opacity = 0.1, 
              weight = 0.2) %>% 
  addPolygons(data = sin_datos, fill = NA, 
              color = "black", weight = 0.2) %>% 
  addPolygons(data = mpios, 
              color = "green", 
              fill = NA,
              weight = 0.7) %>% 
  addLegend(pal = pal, 
            values =  mapa$total_200815, 
            position = "bottomright", 
            labFormat = labelFormat(suffix = " casos"),
            opacity = 1,
            title = "<b>Casos históricos<br>registrados por<br>colonia<br>en la<br>CDMX")

# Mapas de colonias de atención prioritaria: 
col_prioritarias$casos_activ

##   [1]  9 15 19  4 12 12 26 20 27  4  8  7  6 11  4  6  6  4 10 14 16 13  3
##  [24] 17 21  8 11  8  6  9 12 19 22 18 10  8  7 14  8 16  9 12 18 15 20 10
##  [47] 14 22 23 26  0 13 14 16  7 10  9  5  6 10 10 19 11  6  9 12 13 28  4
##  [70]  8  9  7 16  9  4  5  6  7 15  8  8 15 16  6 11 15  7 30 13  5 17 12
##  [93] 15  8 14 15  4  6 19 23 18  4  4 12 10  6 19  3  6 10  8  6  9 10 20
## [116]  7 10 15  4 18  9  8 10 13 15 10 11  8  8 19 12 23 29  7  7  7 14 11
## [139] 19 17  7 11  9  8  5 10 30 29 14 13  7 30  2  5  7 12 12 17

pal <- colorNumeric(palette = "magma", rev = TRUE,
                    domain = c(0, max(col_prioritarias$casos_activ, na.rm = T)))

label = lapply(paste0("<b>Colonia: </b>", col_prioritarias$cve_col, "-", col_prioritarias$nombre, "<br>"), htmltools::HTML)

popup = paste0("<b>Colonia: </b>", col_prioritarias$cve_col, "-", col_prioritarias$nombre, "<br>",
               "<b>Población: </b>", col_prioritarias$poblacion, "<br>",
               "<b>Casos Activos: </b>", col_prioritarias$casos_activ, "<br>",
               "<b>Tasa de Positividad: </b>", col_prioritarias$tasa_positi, "<br>",
               "<b>Tasa de Casos: </b>", col_prioritarias$tasa_casos_, "<br>", 
               "<b style = 'color:gray; size:10px;'>Datos al 23 de Agosto de 2020 (última modificación)"
               )

leaflet(col_prioritarias) %>%
  addProviderTiles("CartoDB.Positron") %>% 
  addPolygons(color = "black", 
              fillColor = pal(col_prioritarias$casos_activ), 
              opacity = 1, 
              fillOpacity = 0.9, 
              weight = 0.2,
              label = label, 
              popup = popup) %>% 
  addPolygons(data = mapa, 
              fill = NA,
              color = "blue", 
              weight = 0.3) %>% 
  addPolygons(data = mpios, 
              color = "black", 
              fill = NA,
              weight = 0.7) %>% 
  addLegend(pal = pal, 
            values =  col_prioritarias$casos_activ, 
            position = "bottomright", 
            labFormat = labelFormat(suffix = " casos"),
            opacity = 1,
            title = "<b>Casos activos<br>registrados en<br>colonias de atención<br>prioritaria, CDMX")