Este es un procedimiento simple para hacer un mapa coroplético de Españo coloreado en las provincias y CCAA.
Como no tengo tus datos (i.e. la lista de provincias que participan en el convenio en cuestión) me lo voy a inventar.
Descargar los shapefiles y cargarlos como un objeto de la clase sf
En esta dirección https://www.arcgis.com/home/item.html?id=83d81d9336c745fd839465beab885ab7 se puede descargar un .zip con el shapefile necesario.
Hay que descomprimirlo y ubucar la carpeta resultante como subdirectorio del proyecto en el que estamos trabajando en R. De paso le cambié el nombre a ./provincias
para facilitarme las cosas, el nombre original es largo.
Deberían estar los siguientes archivos:
Provincias_ETRS89_30N.dbf Provincias_ETRS89_30N.sbx Provincias_ETRS89_30N.shx
Provincias_ETRS89_30N.prj Provincias_ETRS89_30N.shp
Provincias_ETRS89_30N.sbn Provincias_ETRS89_30N.shp.xml
Es necesario tener instalada la librería sf
. Tiene unas cuantas dependencias algunas a nivel sistema operativo. Puede complicarse un poco en Linux.
Carga el shapefile
library(sf)
library(tidyverse) # La vamos a usar más adelante
carto_base <- sf::read_sf("./provincias/") #Asumiendo que está en el subdirectorio provincias de getwd()
plot(carto_base) #verificacion
Si todo salió bien deberías ver 5 mapas de España. Es solo para verificar que está todo en orden.
Combinar con los datos que queremos graficar
Una de las innovaciones de sf
es que carga los objetos cartográficos en un clase dentro de R que funciona como un data.frame
con una columna especial (geometry
) que contiene una lista con toda la información de la geometría de los polígonos, proyección, etc. Lo importante es que al funcionar como data.frame toda la manipulación de datos (joins con otras bases, filtrados, etc.) se puede hacer con las herramientas usuales de manipulación de datos de R. Funciona especialmente bien con tidyverse
.
Para hacer tu mapa necesitas un data.frame
que tenga dos caractísticas:
Una columna en común con los datos geográficos en carto_base
, para poder combinarlos. Generalmente es el nombre o código numérico (mejor, menos ambiguo) de las unidades territoriales que sobre las que vas graficar. En este caso, provincias. Es mejor que el nombre de esta columna sea el mismo que la columna equivalente en los carto_base
, ya que eso simplificará el join
de los datos. En el shapefile que cargué se llama Texto
o Código
, dependiendo si usas nombres o códigos numéricos.
Otra columna con los datos a partir de los cuales se va a colorear. En tu caso podrías convenio
y tener los niveles Sí
, No
(o 0,1, lo que sea que las distinga). Supongo que ya tienes esos datos, deberías asegurarte de que se cumpla la condición 1.
Como no tengo los datos me los voy a inventar, asumiendo que Salamanca, Madrid y Guipúzcoa participan del convenio y las demas provincias no. Va a ser un poquito complicado porque tengo que simular los datos_a_graficar
Hago un df base con las tres provincias "participantes" y el Sí
para el convenio:
datos_a_graficar <- tibble(Texto = c("Salamanca", "Madrid", "Guipúzcoa"),
convenio = "Sí") #Por regla de reciclaje Sí se extiende tantas veces como el largo de Texto.
Unir los datos de interés con los datos cartográficos:
Uso el método de join porque es el más habitual si ya tienes los datos_a_graficar
left_join(carto_base, datos_a_graficar) %>% #Aquí ya quedan unidos, solo que codigo tiene el valor NA cuando no participa en el convenio y queremos que diga explíctamente No
mutate(convenio = ifelse(is.na(convenio), "No", convenio)) -> carto_unida # Aquí paso los NA a "No" usando dplyr. Nunca se pierden los datos geográficos!
Graficar el mapa básico
Ya tengo carto_unida
que incluye los datos geográficos y los datos del convenio. Primero el mapa más simple:
Nota: geom_sf
aparece en ggplot2 3.0
. Asegurarse de tener la última versión de CRAN.
ggplot
ggplot(data = carto_unida, aes(fill = convenio)) +
geom_sf() +
theme_minimal() +
labs (title = "Provincias que han suscripto el importante convenio",
caption = "Dar acá créditos por los datos y cartografía")
Graficar el mapa de la pregunta
Para hacer lo de las CCAA es un poquito más complicado, porque tu base cartográfica (las filas del objeto sf
) son provincias. Sin embargo como hay una relación jerárquica de los datos (varias provincias una CCAA) se puede hacer. De paso le hago algunos cambios estéticos al mapa.
carto_base %>%
group_by (CCAA) %>%
mutate(convenio = ifelse(any(Texto %in% datos_a_graficar$Texto), #Si alguna provincia participante está en esta CCAA marca toda la CCAA
"CCAA", "No")) %>%
mutate(convenio = ifelse(Texto %in% datos_a_graficar$Texto, #datos_a_graficar$Text es un vector, por eso funciona el %in%
"Provincia", convenio)) %>% #Marco con otro código a las provincias participantes.
mutate(convenio = factor(convenio, levels = c("Provincia", "CCAA", "No"))) %>% #Reordeno el factor, para que se vea mejor la etiqueta
filter(CCAA != "Canarias") %>% #Sólo España continental y Baleares, eliminando esta línea aparece Canarias
ggplot() +
geom_sf(aes(fill = convenio), color = "grey", size = 0.2) +
coord_sf(datum = NA) + #Elimino las coordenadas, no aportan mucho en un coroplético
theme_minimal() + #Fondo blanco, mapa simple.
theme(legend.position = "top") + # La leyenda arriba
scale_fill_manual(values = c("darkgreen", "lightgreen", "white")) + # Se pueden cambiar los colores acá.
labs(title = "Provincias y CCAA participantes del convenio",
fill = "Nivel en el que participa")
leaflet
? Hay una librería relativamente nueva que se llamasf
y facilita muchísimo hacer mapas coropléticos (coloreados de acuerdo con algún atributo de los polígonos interiores, en tu caso, formar o no formar parte del convenio). El truco desf
es que tiene los polígonos dentro de undata.frame
, por lo que puedes manipular todos los elementos no polígonos con las herramientas habituales de R y luego generar el mapa conggplot2
. Advierto: no genera un mapa interactivo javascript comoleaflet
.