Transformar os mapas de risco de enchente de Taiwan em algo clicável: os dados, os números e o que não se pode dizer
A quinta ferramenta está no ar: consulta de potencial de enchente. Aqui está o que medi e calculei ao construí-la, incluindo vários números que na página soariam como conclusões e são apenas saída de simulação.
2026-08-21 作者 William Hsu
De onde vêm os dados e como são
A Agência de Recursos Hídricos publica os “mapas de potencial de enchente” no portal de dados abertos (dataset 25766, sob a licença de dados abertos do governo). Há dois conjuntos: um pacote de 7z por município de 2018 e um conjunto nacional de agosto de 2022 organizado pelo serviço de informação espacial da agência, separado por cenário de chuva. Primeiro baixei o pacote de Keelung: dos seus 87 MB, 150 MB descompactados são dez JPGs de pôster com 14.044×9.934 pixels, e os shapefiles com que dá para calcular ocupam poucos MB. Então troquei pelo conjunto nacional: dez cenários, 140 MB no total, sem JPG, e cada polígono traz o código do município.
- Cenários 150/250/350 mm em 6 horas, 200/300/400 mm em 12 horas e 200/350/500/650 mm em 24 horas — dez ao todo
- Faixas de profundidade 0,3–0,5, 0,5–1, 1–2, 2–3 e acima de 3 metros. Abaixo de 0,3 m não é desenhado
- Coordenadas TWD97 de dois graus; de 1.204 a 2.356 polígonos por cenário, e 3,75 milhões de vértices no maior arquivo
- Versão terceira geração, construída em 12/08/2022; a agência revisa a cada cinco anos
Os limites municipais vêm à parte, do centro nacional de cartografia do Ministério do Interior (edição de março de 2025, 368 municípios), usados como denominador de área e para saber em qual município cai o ponto clicado.
Como isso vira clicável
A origem são polígonos vetoriais e o caminho óbvio seriam blocos vetoriais. Eu não tenho tippecanoe, e a consulta que as pessoas realmente querem é “em que faixa este ponto está, em cada um dos dez cenários”, o que com polígonos vira dez testes de ponto em polígono. Então fui de grade:
- Rasterizar cada cenário numa grade equidistante de 0,0004° de longitude por 0,00036° de latitude — cerca de 41×40 metros em Taiwan, a mesma ordem da grade de 40 metros do original. Mais fino seria precisão falsa.
- Cada célula guarda a faixa dos dez cenários, 10 bytes por célula, e só são guardadas as células que inundam em pelo menos um cenário: 1.276.400 células, 28 MB de SQLite.
- Também guardo grades mais grossas de 8×8 e 64×64 com a faixa máxima, para as visões afastadas.
- Os blocos são desenhados pelo servidor na hora: consulta as células de uma faixa e pinta um PNG de 256×256 em 3 a 30 milissegundos, e depois cacheia. A latitude é Mercator e a grade é equidistante, então cada linha de pixels precisa ser convertida separadamente; amostragem linear não fecha.
- A rasterização inteira roda no rasterio: dez cenários em 82 segundos.
O preço de rasterizar é a área desviar um pouco dos polígonos originais. Comparando a área dos polígonos dos dez cenários (calculada direto no plano TWD97) com a soma das células, a soma sai sempre 0,3 a 0,5% maior, no máximo 0,53%, porque a célula que encosta na borda é contada inteira. Aceitável em escala nacional; quanto menor o município, maior esse erro relativo — e a página avisa.
Números que só aparecem depois da conta
Tudo abaixo é estatística de uma simulação, não registro do que aconteceu.
Quanto de Taiwan inunda
Com 350 mm em 24 horas (mais ou menos um tufão ou uma frente de chuvas forte num dia), a área simulada inundando acima de 0,3 m é de 674,5 km², ou 1,8% do território. A 650 mm vira 2.060 km², 5,6%. Com 350 mm em 6 horas dá 1.178 km². Os mesmos 350 mm espremidos em 6 horas inundam 1,75 vez a área do caso de 24 horas — exatamente por isso a drenagem urbana teme as pancadas curtas.
Municípios com a maior fração
No cenário de 350 mm em 24 horas, as maiores frações simuladas são Lingya (38%), Qianjin (32%) e Xinxing (30%), em Kaohsiung — mas os três são pequenininhos (2 a 8 km²), o que infla a proporção. Grandes e altos ao mesmo tempo são Zhuangwei em Yilan (29% de 39,6 km²), Wujie (24%), Dongshi em Chiayi (21% de 82 km²), Liujiao, Yuanchang em Yunlin e Budai em Chiayi.
A 650 mm, Wujie chega a 67%, Zhuangwei a 65% e a cidade de Yilan a 58%. No país todo, 58 municípios passam de um quarto da área e 8 passam da metade.
Onde a inundação cresce mais rápido com a chuva
Comparando 200 mm e 650 mm em 24 horas, os maiores saltos são Wujie de 7% para 67%, a cidade de Yilan de 2,5% para 58%, Zhuangwei de 10% para 65%, Qieding em Kaohsiung de 5% para 57%, Liujiao em Chiayi de 4% para 53%, Hunei em Kaohsiung de 1% para 49% e o distrito norte da cidade de Hsinchu de 0,6% para 48%. São lugares que parecem tranquilos numa chuva forte comum e somem debaixo d'água quando o total sobe. Três municípios da planície de Lanyang aparecem juntos: o terreno é uma grande bacia.
Municípios onde a pancada curta é pior
Em 166 municípios, a área inundada simulada com 350 mm em 6 horas é mais de 1,5 vez a de 350 mm em 24 horas. As maiores diferenças são Yuli em Hualien (32,5 km² contra 5,6 km²), Wufeng em Taichung (16,9 contra 1,1), Xingang e Minxiong em Chiayi, Zhubei em Hsinchu e Houbi em Tainan. Quase todos dependem de drenagem por gravidade e não dão conta quando tudo chega de uma vez.
Estados e cidades
Ordenando pela fração com 350 mm em 24 horas: Yunlin 7,9%, condado de Chiayi 5,1%, Tainan 4,8% e Changhua 3,7% lideram. Entre as seis grandes cidades: Tainan 4,8%, Kaohsiung 2,1%, Taoyuan 2,0%, Taichung 0,8%, Nova Taipé 0,5% e Taipé 0,4%. Mas Taipé salta para 6,6% no cenário de 650 mm — com muita chuva aparece o limite de drenagem da bacia.
Há 32 municípios sem uma única célula inundada mesmo no maior cenário, todos de montanha (Wutai e Sandimen em Pingtung; Liugui, Jiaxian e Maolin em Kaohsiung; Xinyi e Lugu em Nantou…). Isso não os torna seguros: estes mapas simulam inundação, e o perigo na montanha é fluxo de detritos, que é outro mapa.
O que não se pode dizer
A parte mais demorada de construir esta ferramenta foi decidir o que a página pode dizer. O código foi a parte fácil.
- O artigo 8 do regulamento de divulgação dos dados de potencial de enchente estabelece que esses dados servem apenas para prevenção e resposta a desastres, e que qualquer controle ou restrição de uso do solo deve seguir as normas da autoridade competente. Por isso as páginas não trazem preço de imóvel, seguro nem “vale comprar”.
- O artigo 4 exige registrar que os mapas foram traçados a partir de condições de chuva de projeto, dados específicos de terreno e cálculo de um modelo hidráulico, e que devem ser usados com atenção especial. Está escrito assim.
- As premissas que a própria agência imprime no produto: modelo SOBEK, terreno de 2010–2012, uso do solo de 2008, todas as defesas e drenagens funcionando normalmente e edifícios sem obstruir o fluxo. Tudo isso vai para a página de dados, porque só “edifícios não obstruem” já invalida ler uma célula como “a minha casa”.
- “Potencial” e “vai alagar” são coisas diferentes. O mesmo ponto fica sem área marcada a 200 mm e vermelho a 650 mm; os dois estão certos, e o que muda é a premissa. Por isso um clique lista os dez cenários, mostrando como aquilo se move, em vez de entregar uma cor só.
O que ficou de fora
- Converter endereço em coordenada. Não existe serviço nacional gratuito e legal, então recorro aos 35.339 nomes de rua e lugar indexados para o mapa de acidentes mais os centroides dos 368 municípios.
- Enchentes históricas. O KML de “extensão de enchente por tufões” da agência tem só 1.231 polígonos, concentrados em alguns tufões de 2015–2016 — incompleto demais para publicar sem induzir a erro.
- Níveis de água em tempo real e sensores de enchente. Outro tipo de dado, com outra frequência de atualização; fica para depois.
A ferramenta está em lazetool.com/flood e os números por município em estatísticas municipais (ambos em chinês).