Como localizar escarpas em torno de lagos usando ferramentas gratuitas
O termo "escarpas do lago localização" surge com frequência quando se precisa mapear desníveis topográficos ao redor de corpos d'água. O processo em si é direto, mas esbarra em vários detalhes práticos que só aparecem depois de gastar tempo testando. Vou explicar como funciona no dia a dia, sem enrolação.
O que é escarpas do lago localização
Refere-se basicamente à identificação e delimitação das vertentes íngremes que circundam um lago, geralmente a partir de dados de elevação. A ideia é extrair de um modelo digital de elevação (MDE) as linhas onde o gradiente supera um limite definido, gerando um polígono ou linha que representa a escarpa. Não existe uma ferramenta única chamada "Escarpas do Lago" — é mais um fluxo de trabalho que se monta com dados abertos e software de código aberto.
O fluxo prático, passo a passo
A base é sempre um MDE de boa resolução. O sensor que eu recomendo para a maioria dos casos no Brasil é o SRTM 30m, e quando o terreno é muito acidentado, vale investir no ALOS World 3D 30m. A diferença de qualidade é perceptível logo na extração dos contornos. Passo 1 — obter o dado. Baixe a área de interesse no site do Earth Observation Gateway ou diretamente no site do NASA EARHART. Para o Brasil, o Inpe também disponibiliza mosaicos nacionais. Certifique-se de que a projeção está em coordenadas geográficas (WGS84) antes de continuar.
Passo 2 — recortar o lago. Abra o QGIS. Carregue o MDE e uma shape do perímetro do lago, disponível gratuitamente no Google Earth Engine ou no dataset HydroSHEDS. Use a ferramenta "Recortar raster" para isolar apenas a área ao redor do corpo d'água. Isso reduz tempo de processamento e evita ruído nas bordas da imagem original. Passo 3 — calcular o gradiente. No processador de raster do QGIS, use a ferramenta "Cálculo do gradiente". O resultado é um raster onde cada pixel mostra a inclinação do terreno em graus ou porcentagem. Defina um limiar que faça sentido para o seu objetivo. Escarpas típicas de lago no sudeste brasileiro costumam apresentar declividade acima de 30 graus, mas isso varia conforme o tipo de solo e a geologia local.
👉 Clique no botão abaixo para saber mais sobre o assunto!
Passo 4 — converter para vetor. Aplique uma máscara binária no raster de gradiente (valores acima do limiar viram 1, abaixo viram 0). Depois use "Raster para polígono". O resultado são os trechos de terreno íngreme ao redor do lago. Finalize com dissolve para limpar pequenos fragmentos espúrios.
Problema que eu enfrentei e o que funcionou
Em um projeto para mapear escarpas ao redor de um lago na Serra da Mantiqueira, o MDE do SRTM apresentou buracos de dados (valores nulos) justamente nas áreas de transição entre a encosta e o espelho d'água. O gradiente ficou com falhas, e o contorno da escarpa apresentava lacunas grandes. Eu tentei preencher com interpolação IDW, mas o resultado distorcia a linha de crista. O workaround foi simples: substituir o trecho com buraco por dados do ASTER GDEM v3, que tinha resolução um pouco melhor naquela zona. Sobrepor os dois datasets, fazer blend com a ferramenta "Mosaic to New Raster" e recalcular o gradiente. As escarpas ficaram contínuas. Gastei cerca de 40 minutos a mais no fluxo, mas evitei refazer todo o processamento.
Pegadinhas que iniciantes ignoram
A primeira é confusing slope units. O QGIS permite calcular gradiente em degrees ou percent rise. Se você definir o limiar em percentual e comparar com um valor empírico em graus, o resultado sai completamente errado. Sempre verifique a unidade na tela de diálogo antes de rodar. A segunda é confiar cegamente no limiar fixo. Um valor de 30 graus funciona bem em terrenos cristalinos, mas em áreas de cobertura sedimentar suave, a escarpa real pode começar em 18 graus. O ideal é calibrar com uma imagem de referência, como ortofoto ou dados LiDAR, e ajustar o limiar visualmente antes de processar toda a bacia.
Quando o método simplesmente não funciona
Se o lago for plano, em área de várzea aluvionar, não haverá escarpa significativa para extrair. O gradiente vai ser baixo em toda a extensão e o resultado será um raster quase vazio. Nesse caso, o fluxo inteiro perde o sentido. Também não funcione bem em regiões com cobertura densa de vegetação onde o MDE não penetra o dossel — aí o terreno medido não é o terreno real, e as escarpas ficam deslocadas verso o alto da encosta. Alternativa para esses cenários: usar dados LiDAR aéreo ou airborne, quando disponíveis, ou recorrer a modelos de alta resolução como o RadarSat MapTEC, que oferece 12m de resolução e costuma lidar melhor com o subdosel em regiões tropicais.
Recursos e downloads úteis
O QGIS está disponível em qgis.org. Os MDEs SRTM podem ser baixados pelo Earth Exploration Tool no endereço earthexplorer.usgs.gov. O HydroSHEDS, com as bacias hidrográficas e cursos d'água, fica em hydrosheds.usgs.gov. Para lidar com buracos nos dados, o GRASS GIS, integrado ao QGIS, oferece a ferramenta r.fill.dir, que preenche depressões de forma topologicamente correta. O fluxo completo, quando bem calibrado, leva em média 25 minutos para uma área de até 500 km² em hardware comum. Não é instantâneo, mas é repetível e transparente, o que importa quando o resultado precisa ser auditado ou replicado em outra bacia.
