Cómo crear un mapa de precipitaciones con QGIS (orientado a Bachillerato)

En esta entrada se explica, paso a paso, cómo elaborar un mapa de precipitaciones a partir de los datos de las estaciones automáticas de AEMET en Cantabria. El procedimiento está pensado para alumnado de Bachillerato y utiliza software libre (QGIS).

Qué necesitamos #

1. Software SIG: QGIS (gratuito, multiplataforma).
2. Datos de precipitación: descargables desde la web de AEMET (opción «exportar a Excel»).
3. Cartografía base: límites municipales, provinciales y autonómicos en formato shapefile (.shp), disponibles en el Centro de Descargas del IGN.

Preparación previa #

1. Instalar QGIS y comprobar que se ejecuta correctamente.
2. Configurar el sistema de coordenadas de referencia del proyecto a EPSG:25830 (ETRS89 / UTM zona 30N), en lugar del 4326 que aparece por defecto.
3. Preparar la hoja de cálculo de precipitaciones:
   • Eliminar columnas innecesarias.
   • Añadir dos columnas con las coordenadas UTM (X, Y) de cada estación (obtenidas previamente desde AEMET OpenData o de fuentes oficiales).
   • Exportar el archivo a formato CSV (separador de campos: punto y coma; decimal: coma, según configuración regional).

Elaboración del mapa: secuencia de trabajo #

1. Cargar la cartografía base #

Desde el menú Capa → Añadir capa → Añadir capa vectorial, cargar el shapefile de límites provinciales o autonómicos. Hacer zoom a Cantabria. Verificar que el CRS del proyecto es 25830.

2. Cargar los datos de las estaciones #

Capa → Añadir capa → Añadir capa de texto delimitado. Seleccionar el CSV, indicar las columnas de coordenadas X e Y y el CRS (25830). Las estaciones aparecerán como puntos sobre el mapa.

3. Etiquetar las estaciones #

Botón derecho sobre la capa de puntos → Propiedades → Etiquetas. Activar el etiquetado con el campo del nombre de la estación. Ajustar tamaño, color y posición según preferencia.

4. Interpolar los valores de precipitación (IDW) #

Abrir la Caja de herramientas de procesos (menú Procesos). Buscar IDW (Inverse Distance Weighting).

La interpolación Inverse Distance Weighting (IDW) estima el valor en un punto desconocido a partir de los valores conocidos, ponderados inversamente con la distancia:

$$ Z(x) = \frac{\sum_{i=1}^{n} w_i Z_i}{\sum_{i=1}^{n} w_i} \quad \text{donde} \quad w_i = \frac{1}{d_i^p} $$

donde \(p\) suele ser 2. Parámetros principales en QGIS:

• Capa de puntos de entrada.
• Campo de atributo con la precipitación.
• Extensión: la de la capa de estaciones o la de Cantabria.
• Tamaño de píxel: experimentar con valores adecuados a la escala regional (valores demasiado grandes producen un resultado tosco; demasiado pequeños generan archivos pesados).

El resultado es un ráster en escala de grises que representa la precipitación interpolada.

5. Recortar el ráster al límite de Cantabria #

Ráster → Extracción → Cortar ráster por máscara de capa. Utilizar como máscara la capa poligonal de Cantabria. Se obtiene un ráster limitado al territorio de la comunidad.

6. Aplicar una rampa de color #

Botón derecho sobre el ráster recortado → Propiedades → Simbología. Cambiar el tipo de renderizado a «Banda monocroma» o «Paleta/Valores únicos» según convenga, y aplicar una rampa de color adecuada (por ejemplo, de azules claros a azules oscuros o una rampa de precipitación estándar).

7. Generar isoyetas (curvas de nivel de precipitación) #

Ráster → Extracción → Curvas de nivel. Indicar el intervalo deseado (por ejemplo, cada 10 o 25 mm según la variabilidad de los datos). Las isoyetas se superponen al ráster y facilitan la lectura cuantitativa del mapa.

Conclusiones y posibles mejoras #

Con este flujo de trabajo se pueden elaborar no solo mapas de precipitación, sino también de temperatura u otras variables espaciales que requieran interpolación a partir de puntos.

Las posibilidades de QGIS son mucho más amplias: se puede refinar la simbología, añadir leyendas, escalas gráficas, norte, o combinar el resultado con otras capas temáticas. El objetivo de este ejercicio es, no obstante, ofrecer un procedimiento claro y reproducible para el alumnado de Bachillerato.

La interpolación IDW es solo uno de los métodos posibles. En contextos más avanzados se pueden explorar kriging, spline u otros algoritmos disponibles en la caja de herramientas de QGIS.