Visualizando densidades de población con R

Cuando se dispone de la densidad de población de una región en formato matricial, además de la visualización obvia de la misma pueden hacerse muchos juegos numéricos para representarla de mil formas creativas.

En este artículo vamos a sacarle jugo al dataset con la densidad de población para España en formato GeoTIFF que puede descargarse de The Humanitarian Data Exchange (he usado el fichero llamado "esp_pd_2020_1km.tif"):

Fuente: Spain population density in The Humanitarian Data Exchange

Lo primero que haremos será reproyectar los datos originales en formato longitud/latitud a una proyección cónica Lambert en km. Entonces remuestreamos al tamaño de salida deseado (Full HD de 1920px de ancho), y empezamos a jugar con los datos.

Dibujamos un mapa de España con 3 bandas horizontales de color, correspondiendo cada una al área habitada por 1/3 de la población. Con esta división la población aparenta estar muy uniformemente repartida en el eje norte/sur (clic para ver en resolución original):

A continuación recorremos filas y columnas de la matriz de densidades hasta dar con las medianas en cada uno de los dos ejes. Podemos ver que su intersección es un punto muy bien centrado en la Península, de hecho cae rozando el área urbana de Madrid (clic para ver en resolución original):

En este mapa podemos caer en la tentación de pensar que cada uno de los cuadrantes alberga 1/4 de la población del país, pero no es así. La única garantía es que a ambos lados de cada uno de los ejes habita la mitad de la población, pero los cuadrantes pueden tener diferentes poblaciones.

Sí podemos afirmar sin embargo que sea cual sea la distribución de la población, el cuadrante noroeste debe tener la misma población que el sureste, y el cuadrante noreste la misma que el suroeste. En el caso de España los primeros albergan el 22% de la población cada uno y cada uno de los segundos el 28%:

Ahora hacemos una representación que pondrá de manifiesto la extensión y baja densidad de la España vaciada. Como en el primer mapa, dividimos la población en 3 áreas cada una con la misma población pero esta vez con forma de anillos circulares centrados en el punto de cruce de las medianas anteriores (clic para ver en resolución original):

Queda claro que pese a la gran aportación poblacional de la Comunidad de Madrid, todavía necesitamos abarcar un gran área a su alrededor para lograr llegar a 1/3 de la población del país.

Terminamos con el clásico mapa que cuantifica la concentración de población en los grandes núcleos urbanos. Coloreamos la superficie total en dos tonos: uno para las zonas con mayor densidad de población sumando el 80% de la población, quedando el resto para albergar solamente el 20% restante de la población. El primer grupo solo ocupa el 4,4% del territorio (clic para ver en resolución original):

~~~

Aunque he usado algunos bucles seguramente vectorizables, las funciones matriciales de R han hecho cada ejercicio coser y cantar. Los colores y rótulos se han añadido a posteriori en Photoshop. Repositorio con el código R: GitHub.