Convertir los mapas de riesgo de inundación de Taiwán en algo que se pueda pulsar: los datos, los números y lo que no se puede decir
La quinta herramienta está en marcha: consulta de potencial de inundación. Esto es lo que medí y calculé al construirla, incluidos varios números que en la página se leerían como conclusiones y solo son la salida de una simulación.
2026-08-21 作者 William Hsu
De dónde vienen los datos y qué forma tienen
La Agencia de Recursos Hídricos publica los «mapas de potencial de inundación» en el portal de datos abiertos (dataset 25766, con licencia de datos abiertos del gobierno). Hay dos conjuntos: un paquete de 7z por municipio de 2018 y un conjunto nacional de agosto de 2022 preparado por su servicio de información espacial, separado por escenario de lluvia. Primero descargué el paquete de Keelung: de sus 87 MB, 150 MB descomprimidos son diez JPG tipo póster de 14.044×9.934 píxeles, y los shapefiles con los que de verdad se puede calcular ocupan unos pocos MB. Así que cambié al conjunto nacional: diez escenarios, 140 MB en total, sin JPG, y cada polígono lleva su código de municipio.
- Escenarios 150/250/350 mm en 6 horas, 200/300/400 mm en 12 horas y 200/350/500/650 mm en 24 horas: diez en total
- Bandas de profundidad 0,3–0,5, 0,5–1, 1–2, 2–3 y más de 3 metros. Por debajo de 0,3 m no se dibuja
- Coordenadas TWD97 de dos grados; entre 1.204 y 2.356 polígonos por escenario, y 3,75 millones de vértices en el archivo mayor
- Versión tercera generación, construida el 12-08-2022; la agencia la revisa cada cinco años
Los límites municipales vienen aparte del centro nacional de topografía del Ministerio del Interior (edición de marzo de 2025, 368 municipios), y sirven de denominador de superficie y para saber en qué municipio cae un punto pulsado.
Cómo se vuelve pulsable
El origen son polígonos vectoriales y lo evidente serían teselas vectoriales. No tengo tippecanoe, y la consulta que la gente quiere es «en qué banda está este punto en cada uno de los diez escenarios», que con polígonos son diez pruebas de punto en polígono. Así que fui a una rejilla:
- Rasterizar cada escenario sobre una rejilla equidistante de 0,0004° de longitud por 0,00036° de latitud: unos 41×40 metros en Taiwán, el mismo orden que la rejilla de 40 metros del original. Más fino sería precisión falsa.
- Cada celda guarda la banda de los diez escenarios, 10 bytes por celda, y solo se guardan las celdas que se inundan en al menos un escenario: 1.276.400 celdas, 28 MB de SQLite.
- También se guardan rejillas más gruesas de 8×8 y 64×64 tomando la banda máxima, para las vistas alejadas.
- Las teselas las dibuja el servidor al vuelo: consulta las celdas de un rango y pinta un PNG de 256×256 en 3 a 30 milisegundos, y luego se cachea. La latitud es Mercator y la rejilla es equidistante, así que cada fila de píxeles hay que convertirla por separado; un muestreo lineal no cuadra.
- Toda la rasterización corre con rasterio: diez escenarios en 82 segundos.
El precio de rasterizar es que las áreas se desvían un poco de los polígonos originales. Comparando el área de los polígonos de los diez escenarios (calculada directamente en el plano TWD97) con la suma de celdas, la suma sale siempre entre un 0,3 y un 0,5% más alta, como mucho un 0,53%, porque una celda que toca el borde se cuenta entera. Aceptable a escala nacional; cuanto más pequeño es el municipio, mayor es ese error relativo, y la página lo advierte.
Números que solo aparecen al calcularlos
Todo lo que sigue es estadística de una simulación, no el registro de algo ocurrido.
Cuánto de Taiwán se inunda
Con 350 mm en 24 horas (aproximadamente un tifón o un frente de lluvias intenso en un día), la superficie que la simulación inunda por encima de 0,3 m es de 674,5 km², el 1,8% del territorio. A 650 mm pasa a 2.060 km², el 5,6%. Con 350 mm en 6 horas son 1.178 km². Los mismos 350 mm apretados en 6 horas inundan 1,75 veces la superficie del caso de 24 horas: justo por eso el drenaje urbano teme los chaparrones cortos.
Municipios con mayor proporción
En el escenario de 350 mm en 24 horas, las mayores proporciones simuladas son Lingya (38%), Qianjin (32%) y Xinxing (30%), en Kaohsiung, aunque los tres son diminutos (2 a 8 km²) y eso infla el porcentaje. Grandes y altos a la vez son Zhuangwei en Yilan (29% de 39,6 km²), Wujie (24%), Dongshi en Chiayi (21% de 82 km²), Liujiao, Yuanchang en Yunlin y Budai en Chiayi.
A 650 mm, Wujie llega al 67%, Zhuangwei al 65% y la ciudad de Yilan al 58%. En todo el país, 58 municipios superan la cuarta parte de su superficie y 8 superan la mitad.
Dónde crece más rápido al subir la lluvia
Comparando 200 mm y 650 mm en 24 horas, los mayores saltos son Wujie del 7% al 67%, la ciudad de Yilan del 2,5% al 58%, Zhuangwei del 10% al 65%, Qieding en Kaohsiung del 5% al 57%, Liujiao en Chiayi del 4% al 53%, Hunei en Kaohsiung del 1% al 49% y el distrito norte de la ciudad de Hsinchu del 0,6% al 48%. Son lugares que parecen no tener problemas con una lluvia fuerte normal y se inundan enteros cuando el total sube. Tres municipios de la llanura de Lanyang salen juntos: el terreno es un cuenco grande.
Municipios donde el chaparrón corto es peor
En 166 municipios, la superficie inundada simulada con 350 mm en 6 horas supera en más de 1,5 veces la de 350 mm en 24 horas. Las mayores diferencias son Yuli en Hualien (32,5 km² frente a 5,6 km²), Wufeng en Taichung (16,9 frente a 1,1), Xingang y Minxiong en Chiayi, Zhubei en Hsinchu y Houbi en Tainan. Casi todos son sitios cuyo drenaje depende de la gravedad y no da abasto cuando llega de golpe.
Provincias y ciudades
Ordenando por proporción con 350 mm en 24 horas: Yunlin 7,9%, condado de Chiayi 5,1%, Tainan 4,8% y Changhua 3,7% encabezan la lista. Entre las seis grandes ciudades: Tainan 4,8%, Kaohsiung 2,1%, Taoyuan 2,0%, Taichung 0,8%, Nuevo Taipéi 0,5% y Taipéi 0,4%. Pero Taipéi salta al 6,6% en el escenario de 650 mm: con mucha lluvia aparece el límite de drenaje de la cuenca.
Hay 32 municipios sin una sola celda inundada ni en el escenario mayor, y todos son de montaña (Wutai y Sandimen en Pingtung; Liugui, Jiaxian y Maolin en Kaohsiung; Xinyi y Lugu en Nantou…). Eso no los hace seguros: estos mapas simulan inundación, y el peligro de montaña son los flujos de derrubios, que es otro mapa.
Lo que no se puede decir
La parte más lenta de construir esta herramienta fue decidir qué puede decir la página. La parte de código fue la más sencilla.
- El artículo 8 del reglamento de publicación de datos de potencial de inundación establece que estos datos son solo para prevención y respuesta ante desastres, y que cualquier control o limitación del uso del suelo debe seguir la normativa de la autoridad competente. Por eso las páginas no hablan de precios de vivienda, seguros ni de si comprar.
- El artículo 4 obliga a indicar que los mapas se han trazado a partir de condiciones de lluvia de diseño, datos concretos del terreno y el cálculo de un modelo hidráulico, y que deben usarse con especial cautela. Está escrito tal cual.
- Los supuestos que la propia agencia imprime en su producto: modelo SOBEK, terreno de 2010–2012, usos del suelo de 2008, todas las defensas y drenajes funcionando con normalidad y edificios que no obstruyen el flujo. Todo eso está en la página de datos, porque solo con «los edificios no obstruyen» ya deja de tener sentido leer una celda como «mi casa».
- «Potencial» y «se va a inundar» no son lo mismo. El mismo punto no tiene color a 200 mm y es rojo a 650 mm; ambas cosas son correctas y lo que cambia es el supuesto. Por eso al pulsar se listan los diez escenarios, para ver cómo se mueve, en vez de entregar un solo color.
Lo que no se hizo
- Convertir direcciones en coordenadas. No hay ningún servicio nacional gratuito y legal, así que se recurre a los 35.339 nombres de calles y lugares indexados para el mapa de accidentes más los centroides de los 368 municipios.
- Inundaciones históricas. El KML de «extensión de inundación por tifones» de la agencia tiene solo 1.231 polígonos, concentrados en unos pocos tifones de 2015–2016: demasiado incompleto para publicarlo sin inducir a error.
- Niveles de agua en tiempo real y sensores de inundación. Otro tipo de dato con otra frecuencia de actualización; queda para más adelante.
La herramienta está en lazetool.com/flood y las cifras por municipio en estadísticas municipales (ambas en chino).