Conversor de Gauss-Krüger para WGS84 e UTM

Por Que Precisa de Existir um Conversor de Coordenadas Alemão

Em quinze anos a trabalhar com dados geoespaciais, perdi a conta de quantos mapas de infraestruturas alemãs, extratos cadastrais e antigas plantas de levantamento topográfico abri só para encontrar coordenadas que não fazem sentido nenhum num mapa-múndi normal — números na casa dos milhões, concentrados à volta de valores como 4,593,000 e 5,821,000. Isso não são dados corrompidos. Isso é Gauss-Krüger, e continua por todo o lado nos arquivos geoespaciais da Alemanha, décadas depois de o país ter mudado oficialmente para outro sistema.

A minha opinião sincera, depois de lidar com isto profissionalmente: a maioria das ferramentas genéricas de "conversor de coordenadas" disponíveis online erra de uma de duas maneiras. Ou tratam o Gauss-Krüger como se fosse permutável com qualquer projeção de Mercator transversa e usam o elipsoide errado, ou aplicam uma única aproximação de transformação de datum válida para todo o país e ignoram que as redes de triangulação históricas da Alemanha Ocidental e da Alemanha Oriental nunca foram unificadas entre si antes da reunificação. Os dois erros são pequenos no código e grandes nas consequências — um erro silencioso de algumas dezenas a algumas centenas de metros, exatamente do tipo que passa despercebido até alguém tentar sobrepor o resultado num mapa real e ver a linha do lote a cair no terreno errado.

Esta ferramenta foi construída para levar isso a sério: elipsoide correto, convenção de zona correta, uma escolha regional entre a Alemanha Ocidental e a Oriental para a transformação de datum, e suporte completo para converter ficheiros GeoJSON inteiros, não apenas pontos isolados — porque, na prática, quase ninguém precisa de converter uma coordenada; precisa de converter um ficheiro cheio delas.

Dois Sistemas de Coordenadas, Uma História Dividida

O Gauss-Krüger e o WGS84/UTM não foram criados na mesma época, com o mesmo propósito, nem pelo mesmo tipo de instituição. Perceber porquê explica a maior parte das particularidades práticas com que se vai deparar.

Gauss-Krüger: Nascido de um Levantamento Topográfico do Século XIX

A história começa com Carl Friedrich Gauss, que passou grande parte do período entre 1821 e 1825 a dirigir a triangulação do Reino de Hanôver. Gauss precisava de uma forma de projetar a superfície curva da Terra sobre papel plano sem distorcer os ângulos, porque a medição de ângulos estava no centro do levantamento por triangulação. Adaptou a projeção de Mercator transversa, fazendo correr uma faixa estreita de norte a sul ao longo de um meridiano escolhido, de modo a que a distorção se mantivesse pequena dentro dessa faixa. Gauss nunca publicou uma fórmula geral completa para o método. O seu trabalho sobreviveu sobretudo como notas privadas e exemplos resolvidos ligados ao levantamento de Hanôver, e não como um procedimento público e reutilizável.

A matemática que Gauss esboçou ficou, em grande parte, por utilizar durante quase um século. Antes de 1871, a Alemanha era um mosaico de dezenas de estados independentes, e cada estado conduzia o seu próprio levantamento territorial com o seu próprio ponto de referência e o seu próprio sistema de coordenadas local. Mesmo depois da unificação, os geodesistas alemães continuavam sem uma única projeção partilhada que cobrisse todo o país. Isso mudou principalmente através do trabalho de Johann Heinrich Louis Krüger, que em 1912 publicou uma derivação rigorosa e geral da projeção de Gauss, completa com as fórmulas de expansão em série ainda hoje utilizadas para converter entre latitude/longitude e coordenadas projetadas. O artigo de Krüger transformou o esboço privado de Gauss num procedimento prático e repetível que qualquer serviço de levantamento topográfico podia aplicar. É por isso que a projeção tem os dois nomes, Gauss-Krüger, honrando a ideia original e a pessoa que a tornou utilizável.

A partir da década de 1920, os estados alemães adotaram o Gauss-Krüger como a sua projeção de levantamento padrão, dividindo o país em zonas com 3 graus de largura, cada uma centrada num meridiano a 6°E, 9°E, 12°E ou 15°E, de modo a manter a distorção dentro de limites aceitáveis para trabalhos cadastrais e de engenharia. O elipsoide de referência subjacente era o Bessel 1841, já com cerca de setenta anos quando as fórmulas de Krüger foram adotadas, mas que se ajustava suficientemente bem à rede de triangulação existente para que substituí-lo nunca compensasse a perturbação que isso causaria.

Depois veio a divisão da Alemanha após a Segunda Guerra Mundial. A Alemanha Ocidental e a Alemanha Oriental mantiveram cada uma a sua própria grelha Gauss-Krüger em funcionamento, mas as duas redes de triangulação foram ajustadas de forma independente durante mais de quatro décadas, sem qualquer coordenação entre elas. O resultado prático, ainda visível hoje: a transformação de melhor ajuste entre dados antigos de Gauss-Krüger e o WGS84 moderno difere genuinamente consoante o ponto se situe em território da antiga Alemanha Ocidental ou Oriental, um pormenor que muitas ferramentas de conversão ignoram silenciosamente. É um pequeno e concreto fragmento da história da Guerra Fria que continua incorporado nos dados cartográficos alemães.

A reunificação em 1990 não fez desaparecer o Gauss-Krüger de um dia para o outro. A Alemanha oficial só passou a adotar o ETRS89/UTM como novo sistema de referência legal de forma gradual, estado a estado, ao longo das décadas de 1990 e 2000, e alguns serviços cadastrais só terminaram a mudança bem dentro da década de 2010. Milhões de registos de levantamentos antigos, mapas de infraestruturas e parcelas cadastrais nunca foram reprojetados e continuam a existir apenas como coordenadas Gauss-Krüger. Quem digitaliza mapas alemães antigos, infraestruturas de serviços públicos ou registos fundiários históricos acaba sempre por encontrar o Gauss-Krüger mais cedo ou mais tarde, e é exatamente por isso que um conversor como este continua a justificar a sua existência.

WGS84 e UTM: Concebidos para a Era dos Satélites

O WGS84 tem uma origem muito diferente da do Gauss-Krüger. Ninguém fez o levantamento de um único reino para o produzir. Foi o Departamento de Defesa dos Estados Unidos que o construiu, precisamente para tornar possível o posicionamento por satélite. A geodesia por satélite dos primeiros tempos, ao longo das décadas de 1960 e 1970, mostrou que nenhum datum nacional existente era suficientemente preciso ou consistente em todo o planeta para suportar um único sistema de posicionamento. A rede de levantamento de cada país carregava os seus próprios pequenos erros, otimizados localmente, úteis apenas para esse país. A tarefa em causa exigia um único referencial global que se comportasse da mesma forma em toda a Terra, sem um mosaico regional por baixo.

O WGS84 foi publicado em 1984 e, desde então, foi discretamente aperfeiçoado várias vezes, a mais recente em 2013, cada revisão apertando o ajuste entre o elipsoide de referência e a forma real da Terra à medida que os dados de rastreio por satélite melhoravam. O nome e as coordenadas do dia a dia mantiveram-se estáveis ao longo de todo esse processo. É genuinamente um sistema global único, e não um conjunto de redes regionais ajustadas separadamente e depois costuradas entre si, como acontece efetivamente com o Gauss-Krüger.

O UTM, o sistema Universal Transverse Mercator, é anterior ao WGS84 em décadas. O Exército dos Estados Unidos padronizou-o na década de 1940, mas assenta na mesma ideia central que Gauss e Krüger já tinham desenvolvido: dividir o globo em zonas estreitas de Mercator transversa para manter a distorção baixa. O UTM limita-se a fazer isto a nível mundial, em zonas uniformes de 6 graus, numeradas de 1 a 60, em vez das zonas Gauss-Krüger mais estreitas, de 3 graus, usadas na Alemanha. Uma coordenada descrita como "UTM sobre WGS84" está a usar a mesma matemática de Mercator transversa que Gauss desenvolveu na década de 1820, apenas sobre um elipsoide global moderno em vez do Bessel 1841, com zonas mais largas.

A Europa adotou o ETRS89, um referencial fixo à placa europeia, funcionalmente idêntico ao WGS84 para praticamente qualquer finalidade de SIG, como seu padrão na década de 1990. As coordenadas oficiais da Alemanha hoje em dia são ETRS89/UTM: zona 32N a oeste, zona 33N a leste.

A Matemática: Como a Projeção Funciona na Realidade

Vou explicar isto da forma como explicaria a um estudante de SIG, e não da forma como costuma ser apresentado — como uma parede de letras gregas sem qualquer explicação do que cada uma significa fisicamente.

Passo 1 — Conheça o seu elipsoide

Toda a projeção parte de um modelo matemático da forma da Terra: um elipsoide, ligeiramente achatado nos polos. Dois números definem-no por completo: o semieixo maior a (o raio equatorial) e o achatamento f. Para o Bessel 1841, o elipsoide sobre o qual os dados Gauss-Krüger são construídos:

• a = 6,377,397.155 metros
• 1/f = 299.1528128 (logo, f ≈ 0.0033427731)

A partir de a e f, todas as outras grandezas nas fórmulas seguintes são derivadas, sendo a mais importante a primeira excentricidade ao quadrado, e² = 2f − f², e a segunda excentricidade ao quadrado, e'² = e²/(1−e²). Estes dois números aparecem constantemente no que se segue.

Passo 2 — Comprimento do arco de meridiano: a espinha dorsal de toda a projeção

O elemento mais importante de todos é o comprimento do arco de meridiano, M(φ) — literalmente, a distância que seria necessário percorrer sobre a superfície da Terra, a direito para norte a partir do equador, para chegar à latitude φ. Como a Terra é um elipsoide e não uma esfera, esta distância não varia linearmente com φ, e não existe para ela uma expressão fechada simples. Em vez disso, é calculada como uma série convergente em e²:

M(φ) = a · [ (1 − e²/4 − 3e⁴/64 − 5e⁶/256) · φ
            − (3e²/8 + 3e⁴/32 + 45e⁶/1024) · sin(2φ)
            + (15e⁴/256 + 45e⁶/1024) · sin(4φ)
            − (35e⁶/3072) · sin(6φ) ]

Cada termo depois do primeiro é uma correção para o achatamento elipsoidal, e cada termo sucessivo é menor do que o anterior por um fator aproximado de e² (cerca de 0.0067 para o Bessel 1841) — razão pela qual quatro termos já são suficientemente precisos, com erro inferior a uma fração de milímetro, em toda a gama de latitudes alemãs. Esta é a base do Northing: tudo o resto que se lhe soma é uma correção relativa à distância a leste ou a oeste do meridiano central em que nos encontramos.

Passo 3 — Os termos de correção este-oeste

Mais três grandezas, todas calculadas na latitude φ, determinam a correção relativa à distância em longitude ao meridiano central:

• N = a / √(1 − e²sin²φ) — o raio de curvatura no primeiro vertical (a que distância real corresponde um grau de longitude nesta latitude)
• T = tan²φ
• C = e'²cos²φ

Com l como a diferença de longitude em relação ao meridiano central, em radianos, as equações diretas completas (Snyder, 1987) são:

Easting  = N·cosφ·l + (N/6)·cos³φ·(1−T+C)·l³
                     + (N/120)·cos⁵φ·(5−18T+T²+14C−58TC)·l⁵

Northing = M(φ) + (N/2)·sinφ·cosφ·l²
                 + (N/24)·sinφ·cos³φ·(5−T+9C+4C²)·l⁴
                 + (N/720)·sinφ·cos⁵φ·(61−58T+T²+270C−330TC)·l⁶

Isto parece intimidante à primeira leitura, mas a estrutura é mais simples do que parece: é uma expansão em série de Taylor em potências de l, que corrige a aproximação simples de Terra plana (N·cosφ·l) tendo em conta como a curvatura do elipsoide distorce as coisas à medida que nos afastamos do meridiano central. Dentro das meias-zonas estreitas de 1.5 graus do Gauss-Krüger, os termos l⁵ e l⁶ já contribuem com menos de um milímetro — a série converge extremamente depressa a esta escala.

Passo 4 — A transformação de datum, antes ou depois da projeção

Tudo o que foi visto até aqui acontece inteiramente dentro de um único elipsoide. Para passar entre o Bessel 1841 (Gauss-Krüger) e o GRS80 (WGS84), é necessário um passo separado: uma transformação de Helmert de 7 parâmetros, aplicada no espaço geocêntrico 3D XYZ, e não diretamente sobre a latitude/longitude. Isto implica converter as coordenadas geográficas em XYZ geocêntrico, aplicar três translações, três rotações minúsculas e uma correção de escala, e depois converter de novo:

X' = dx + (1+s)·( X − rz·Y + ry·Z )
Y' = dy + (1+s)·( rz·X + Y − rx·Z )
Z' = dz + (1+s)·( −ry·X + rx·Y + Z )

A minha opinião como profissional da área: este é o passo em que quase toda a gente falha, não porque a fórmula seja difícil, mas porque o ignoram por completo e tentam reutilizar diretamente o easting/northing do Gauss-Krüger como se estivesse sobre o elipsoide WGS84. Nunca está. A matemática da projeção e a transformação de datum são duas operações independentes, e ambas são importantes.

Que zona e região deve escolher?

• Antiga Alemanha Ocidental: Baden-Vurtemberga, Baviera, Bremen, Hamburgo, Hesse, Baixa Saxónia, Renânia do Norte-Vestfália, Renânia-Palatinado, Sarre, Schleswig-Holstein
• Antiga Alemanha Oriental: Berlim, Brandemburgo, Meclemburgo-Pomerânia Ocidental, Saxónia, Saxónia-Anhalt, Turíngia

Não tem a certeza de qual se aplica? Converta primeiro o seu ponto com qualquer uma das regiões selecionada e depois clique em "Usar a zona GK/UTM correta para este local" — a ferramenta aplica uma regra geográfica simples (aproximadamente: a leste de 10,5°E e a norte de 50,8°N conta como antiga Alemanha Oriental) para escolher automaticamente a zona e a região corretas.

Resolvido à Mão: Convertendo o Portão de Brandemburgo

Números reais, não valores de exemplo. Esta é a cadeia de cálculo que a ferramenta efetivamente executa, mostrada passo a passo, partindo da coordenada WGS84 bem documentada do Portão de Brandemburgo.

Ponto de partida (WGS84): latitude 52.516266°N, longitude 13.377775°E.

Passo 1: Geográfico → XYZ geocêntrico (sobre GRS80/WGS84)

Conversão para coordenadas cartesianas 3D centradas na Terra:

X = 3,783,993.8378 m
Y =   899,923.5184 m
Z = 5,037,966.2582 m

Passo 2: Aplicar a transformação de Helmert inversa (parâmetros da Alemanha Ocidental)

Usando o conjunto de parâmetros da Alemanha Ocidental com sinal invertido (dx=−598.1, dy=−73.7, dz=−418.2, rx=−0.202″, ry=−0.045″, rz=+2.455″, s=−6.7 ppm) para passar do WGS84 de volta para o DHDN:

X' = 3,783,358.5749 m
Y' =   899,893.7601 m
Z' = 5,037,514.2481 m

Passo 3: Geocêntrico → geográfico (sobre o Bessel 1841)

Convertendo X', Y', Z' de novo para latitude/longitude, desta vez usando o elipsoide Bessel 1841:

DHDN latitude  = 52.51767297°
DHDN longitude = 13.37951396°

Repare que isto não é igual à coordenada WGS84 — a latitude do DHDN fica cerca de 0.0014° mais a norte e a sua longitude cerca de 0.0017° mais a leste. Esse pequeno desvio, cerca de 150-190 metros no terreno, é inteiramente o efeito acumulado de um elipsoide diferente e de um referencial histórico diferente. É exatamente por isso que não se pode simplesmente reetiquetar coordenadas WGS84 como DHDN e dar o assunto por resolvido.

Passo 4: Projetar sobre a Zona 4 de Gauss-Krüger (meridiano central 12°E)

Com φ = 52.51767297° (0.9166063 rad) e uma diferença de longitude em relação ao meridiano central de l = 1.379514° (0.0240771 rad), as grandezas intermédias resultam em:

N (radius of curvature) = 6,390,841.3354 m
T = tan²φ               = 1.700567
C = e'²cos²φ             = 0.002488
M (meridian arc length)  = 5,820,346.7101 m

Substituindo estes valores nas séries de Easting e Northing da secção de matemática acima:

Easting (before false easting)  =    93,631.7509 m
Northing                        = 5,821,241.2561 m

O false easting da Zona 4 é 4,500,000 (o "4" inicial identifica a zona), que se soma ao easting bruto calculado acima:

Rechtswert = 4,500,000 + 93,631.75 = 4,593,631.75
Hochwert   =                          5,821,241.26

Para comparação: o mesmo ponto na zona 33N do UTM

Passando a mesma coordenada original WGS84 novamente pelas fórmulas de Mercator transversa — desta vez diretamente sobre o elipsoide GRS80, sem necessidade de transformação de datum, meridiano central 15°E, fator de escala 0.9996 — obtém-se:

UTM 33N Easting  =   389,922.84 m
UTM 33N Northing = 5,819,698.02 m

Ambos os resultados descrevem exatamente o mesmo ponto físico na Terra. A diferença de quase 3.7 milhões de metros entre os dois valores de Northing nada tem a ver com distância — é puramente a diferença entre a convenção de false easting/zona do Gauss-Krüger e a do UTM, mais o pequeno desvio resultante de projetar sobre o Bessel 1841 em vez do GRS80. É exatamente este tipo de confusão que este conversor existe para eliminar.

Passo a Passo: Usar o Conversor

Passo 1

Escolha o formato de origem em "Converter de" e o formato de destino em "Converter para". Para o Gauss-Krüger, escolha também a zona correta (de 2 a 5, correspondente ao meridiano central sobre o qual os seus dados foram levantados) e se o ponto se situa na antiga Alemanha Ocidental ou Oriental.

Passo 2

Para um único ponto, mude para o separador "Coordenada Única" e digite os dois valores — Rechtswert/Hochwert, latitude/longitude, ou easting/northing, consoante o formato de origem escolhido.

Passo 3

Clique em Converter. O resultado aparece imediatamente abaixo, juntamente com a coordenada WGS84 usada como passo intermédio do cálculo, e um marcador surge no mapa para que possa verificar a localização num relance — uma seleção errada de zona feita sem querer torna-se aqui, geralmente, óbvia de imediato.

Passo 4

Mude para o separador "Ficheiro GeoJSON" para converter um ficheiro inteiro de uma só vez. Carregue-o ou cole o seu conteúdo diretamente, defina os mesmos formatos de origem/destino e clique em Converter e Pré-visualizar.

Passo 5

A forma transformada é desenhada no mapa para que possa confirmar que fica no lugar certo, e o painel indica exatamente quantos pares de coordenadas foram convertidos. Clique em "Descarregar converted.geojson" para guardar o resultado — todas as propriedades originais de cada feature são preservadas sem alterações.

Quem Realmente Precisa Deste Conversor

Sem ordem particular, a partir do que já vi, na prática, as pessoas usarem uma ferramenta como esta:

1. Empresas de serviços públicos a digitalizar mapas de redes antigos. Os desenhos de redes de gás, água e eletricidade anteriores aos anos 2000 estão, muito frequentemente, ainda em Gauss-Krüger, e modernizar um SIG de gestão de ativos implica converter décadas de geometria de tubagens e cabos sem perder precisão.
2. Topógrafos a reconciliar registos cadastrais antigos com levantamentos GNSS modernos. Um marco divisório levantado em 1975 e outro medido com um recetor RTK-GPS na semana passada precisam de ser comparados no mesmo sistema de coordenadas antes de alguém poder dizer se coincidem.
3. Arquivistas e historiadores a digitalizar mapas históricos alemães. Os mapas topográficos e cadastrais de antes da guerra e da era da Guerra Fria trazem grelhas Gauss-Krüger; converter elementos de mapas digitalizados para um formato SIG moderno exige exatamente esta transformação.
4. Consultores ambientais a trabalhar com redes de estações de monitorização mais antigas. As redes de monitorização de águas subterrâneas, solo e qualidade do ar criadas há décadas referenciam frequentemente, na sua documentação original, a localização das estações em Gauss-Krüger.
5. Programadores de software a integrar com portais de dados abertos do governo alemão. Algumas exportações de dados abertos alemãs, mais antigas ou regionais, em GeoJSON continuam a trazer coordenadas Gauss-Krüger dentro de um formato que é, tecnicamente, apenas WGS84, e um mapa web irá renderizá-las silenciosamente no local errado, normalmente algures no Golfo da Guiné, até serem convertidas.
6. Planeadores de energia eólica e solar a avaliar dados antigos de subestações e ligações à rede. Os registos de infraestruturas de rede anteriores à expansão das energias renováveis estão frequentemente ainda na projeção mais antiga.
7. Estudantes universitários de geografia e topografia a aprender geodesia. O Gauss-Krüger continua a ser ensinado nas disciplinas alemãs de Vermessungskunde precisamente pela sua importância histórica, e porque trabalhar a transformação à mão, como no exemplo resolvido acima, é a forma como o conceito realmente fica assimilado.
8. Consultores de SIG a auditar a base de dados espacial existente de um cliente quanto à consistência de datum. Bases de dados com datums misturados, em que algumas camadas continuam discretamente em DHDN e outras em ETRS89 sem que ninguém tenha identificado qual é qual, são um problema muito comum e muito dispendioso de descobrir tarde, que esta ferramenta ajuda a diagnosticar rapidamente.

Perguntas Técnicas Frequentes

Respostas mantidas deliberadamente específicas em vez de genéricas — se já sabe o suficiente para fazer estas perguntas, merece uma resposta a sério.