Questões Sobre Coordenadas Geográficas

O que realmente acontece quando você trabalha com coordenadas geográficas no dia a dia

A gente começa a mexer com latitude e longitude achando que é só colocar os números num arquivo CSV e pronto. A realidade é bem diferente. Já vi projeto inteiro de mapeamento travar porque alguém salvou graus decimais como strings em vez de float, e o sistema de visualização interpretou 45.123456 como texto vazio. Isso não é erro bobo, é o tipo de coisa que custa horas pra diagnosticar. Questões sobre coordenadas geográficas aparecem sempre que você tenta cruzar dados de fontes diferentes. O INPE usa SAD69, a IBGE usa SIRGAS2000, e o Google Maps trabalha com WGS84. Na prática, isso significa que o mesmo ponto na superfície da Terra pode ter coordenadas ligeiramente diferentes dependendo do sistema de referência que você escolhe. A diferença costuma ficar entre 1 e 2 metros, o que parece pouco até você precisar fazer um recorte topográfico preciso ou um zoneamento ambiental.

Diferenças práticas entre sistemas de referência

WGS84 é o padrão que todo mundo conhece porque é o que o GPS civil usa. SIRGAS2000 é o sistema oficial do Brasil desde 2015, adotado pelo IBGE e pela cartografia nacional. SAD69 é o antigo sistema sul-americano que ainda aparece em mapas físicos antigos e documentos de registro imobiliário. A conversão entre eles não é trivial. O problema real é que a transformação não é uma simples troca de parâmetros. Quando eu precisei migrar uma base de imóveis rurais de Cuiabá para SIRGAS2000, usei o projeto PROJ.4 com os parâmetros ntl_grid (null) e a diferença saltou para cerca de 1,8 metros no sentido norte-sul. Se você usar apenas uma transformação padrão de 7 parâmetros ( Helmert ), o erro pode chegar a 3 metros em regiões como o Centro-Oeste. A solução foi aplicar o grid de deformação do IBGE, que leva em conta o movimento das placas tectônicas na América do Sul.

Como converter sem errar

Se você tá começando agora e precisa só plotar pontos num mapa, WGS84 resolve. Use EPSG:4326. Mas se o trabalho exige precisão legal ou técnica, aí tem que tomar cuidado. O comando básico com ogr2ogr seria algo como transformar de SAD69 para SIRGAS2000 aplicando o grid NADCON_BR, mas isso exige baixar os arquivos de grid do site do IBGE primeiro. Aqui vai um exemplo prático de linha de comando que eu uso:

ogr2ogr -t_srs EPSG:4674 -s_srs EPSG:4618 entrada.shp saida.shp -tr 0.0001 O código 4618 é SAD69 geodésico, o 4674 é SIRGAS2000 geodésico. O parâmetro -tr define a resolução de transformação, e 0.0001 graus é aproximadamente 11 metros na superfície. Para obras de engenharia, eu costumo usar 0.00001 que dá cerca de 1 metro de precisão. O processo completo de conversão de uma shapefile com 50 mil polígonos leva em média 3 a 5 minutos numa máquina padrão, dependendo da complexidade geométrica.

Armazenamento e formatação

Formato decimal graus é o mais comum. Graus minutos segundos aparece em documentos cartoriais e mapas impressos antigos. O problema é que alguns sistemas aceitam vírgula, outros ponto, e isso gera confusão. Eu já vi arquivo onde a latitude vinha como 23,5 e o sistema interpretou como 23 graus e 5 minutos, quando na verdade era 23,5 graus decimais. Para evitar isso, padronize o uso de ponto decimal e inclua sempre a descrição do sistema de referência. Um header mínimo deveria ter: fonte dos dados, datum, sistema de projeção, data de coleta e margem de erro estimada. Sem essas informações, qualquer pessoa que receber o arquivo vai precisar gastar tempo investigando, e na maioria das vezes não consegue.

Erros comuns que todo mundo comete

Inverter latitude e longitude. Simples assim. O padrão é (latitude, longitude), mas muitos sistemas de GIS aceitam (x, y) que corresponde a (longitude, latitude). A APIs de mapas modernos como Leaflet e Mapbox usam [longitude, latitude] na maioria dos formatos, enquanto o padrão OGC mantém (latitude, longitude). Essa contradição causa confusão constante. Outro erro frequente é usar coordenadas de projeção como se fossem geográficas. Se você extraiu dados num sistema UTM e tentou visualizar direto num mapa mundial, os pontos vão aparecer deslocados. O UTM divide o mundo em 60 fusos, e cada fuso tem seu próprio sistema de coordenadas. Um mesmo ponto pode ter valores completamente diferentes em fusos adjacentes.

Quando não usar coordenadas geográficas tradicionais

Para trabalhos de alta precisão como monitoramento de falhas geológicas, plataformas GNSS RTK chegam a 1 centímetro de erro. Para agricultura de precisão, drones com RTK dão resultados confiáveis abaixo de 5 centímetros. Se o seu projeto exige essa precisão, coordenadas WGS84 genéricas não são suficientes e você precisa de dados diferenciados. Já para coisas simples como marcar um ponto num aplicativo ou fazer uma busca por endereço, WGS84 com 6 casas decimais é mais do que suficiente. Cada casa decimal adiciona aproximadamente 11 metros de precisão, então 6 casas dão cerca de 11 metros, o que é adequado para a maioria das aplicações de consumo.

Recursos úteis

O site do INPE oferece conversores online gratuitos. O PROJ.4 é a biblioteca open source mais robusta pra transformações. O QGIS permite visualizar e converter camadas com diferentes sistemas de referência sem precisar escrever código. Pra quem trabalha com Python, a biblioteca pyproj faz as transformações de forma transparente. Uma dica prática: antes de processar uma base grande, converta um ponto de teste e compare com dados oficiais do IBGE ou do INPE. Se a diferença for maior que 2 metros, revise os parâmetros de transformação antes de continuar.