Convertisseur Gauss-Krüger vers WGS84 et UTM

Pourquoi un convertisseur de coordonnées allemand doit exister

En quinze ans de travail avec des données géospatiales, j'ai perdu le compte du nombre de cartes de réseaux allemandes, d'extraits cadastraux et d'anciens plans de levé que j'ai ouverts pour n'y trouver que des coordonnées qui n'ont aucun sens sur une carte du monde normale — des nombres dans les millions, regroupés autour de valeurs comme 4,593,000 et 5,821,000. Ce ne sont pas des données corrompues. C'est du Gauss-Krüger, et on le trouve encore partout dans les archives géospatiales de l'Allemagne, des décennies après que le pays est officiellement passé à autre chose.

Mon opinion honnête, après avoir traité cela professionnellement : la plupart des outils génériques de « conversion de coordonnées » en ligne se trompent de l'une de ces deux façons. Soit ils traitent le Gauss-Krüger comme s'il était interchangeable avec n'importe quelle projection Mercator transverse et utilisent le mauvais ellipsoïde, soit ils appliquent une approximation unique de décalage de datum à l'échelle nationale et ignorent que les réseaux de levé historiques de l'Allemagne de l'Ouest et de l'Allemagne de l'Est n'ont jamais été unifiés entre eux avant la réunification. Les deux erreurs sont minimes dans le code et lourdes de conséquences — quelques dizaines à quelques centaines de mètres d'erreur silencieuse, exactement le genre d'erreur qui passe inaperçue jusqu'à ce que quelqu'un essaie de superposer le résultat sur une vraie carte et que la limite de la parcelle tombe dans le mauvais champ.

Cet outil est conçu pour prendre cela au sérieux : le bon ellipsoïde, la bonne convention de zone, un choix régional Ouest/Est pour le décalage de datum, et une prise en charge complète de la conversion de fichiers GeoJSON entiers, pas seulement de points uniques — car en pratique, presque personne n'a besoin de convertir une seule coordonnée ; les gens ont besoin de convertir un fichier qui en est plein.

Deux systèmes de coordonnées, une histoire divisée

Le Gauss-Krüger et le WGS84/UTM n'ont pas été construits à la même époque, dans le même but, ni par le même type d'institution. Comprendre pourquoi explique la plupart des bizarreries pratiques que vous rencontrerez.

Gauss-Krüger : né d'un levé du XIXe siècle

L'histoire commence avec Carl Friedrich Gauss, qui passa une grande partie des années 1821 à 1825 à diriger la triangulation du royaume de Hanovre. Gauss avait besoin d'un moyen de projeter la surface courbe de la Terre sur du papier plat sans déformer les angles, car la mesure des angles était au cœur du levé par triangulation. Il adapta la projection Mercator transverse, faisant courir une bande étroite du nord au sud le long d'un méridien choisi afin que la déformation reste faible à l'intérieur de cette bande. Gauss n'a jamais publié de formule générale complète pour la méthode. Son travail a survécu principalement sous forme de notes privées et d'exemples travaillés liés au levé de Hanovre, et non comme une procédure publique et réutilisable.

Les mathématiques esquissées par Gauss sont restées largement inutilisées pendant près d'un siècle. L'Allemagne était une mosaïque de dizaines d'États indépendants avant 1871, et chaque État menait son propre levé topographique avec son propre point de référence et son propre système de coordonnées local. Même après l'unification, les géodésiens allemands ne disposaient toujours d'aucune projection commune unique couvrant l'ensemble du pays. Cela changea principalement grâce au travail de Johann Heinrich Louis Krüger, qui publia en 1912 une dérivation rigoureuse et générale de la projection de Gauss, avec les formules de développement en série encore utilisées aujourd'hui pour convertir entre latitude/longitude et coordonnées projetées. L'article de Krüger transforma l'esquisse privée de Gauss en une procédure pratique et reproductible que n'importe quel bureau de levé pouvait appliquer. C'est pourquoi la projection porte les deux noms, Gauss-Krüger, honorant l'idée originale et la personne qui l'a rendue utilisable.

À partir des années 1920, les États allemands adoptèrent le Gauss-Krüger comme projection standard pour les levés, en divisant le pays en zones de 3 degrés de large, chacune centrée sur un méridien à 6°E, 9°E, 12°E ou 15°E, afin de maintenir la déformation dans des limites acceptables pour les travaux cadastraux et d'ingénierie. L'ellipsoïde de référence sous-jacent était Bessel 1841, déjà vieux d'environ soixante-dix ans au moment de l'adoption des formules de Krüger, mais il correspondait suffisamment bien au réseau de triangulation existant pour que son remplacement ne vaille jamais la perturbation occasionnée.

Vint ensuite la division de l'Allemagne après la Seconde Guerre mondiale. L'Allemagne de l'Ouest et l'Allemagne de l'Est conservèrent chacune leur propre grille Gauss-Krüger, mais les deux réseaux de triangulation furent ajustés indépendamment pendant plus de quatre décennies, sans coordination entre eux. Le résultat pratique, encore visible aujourd'hui : la transformation optimale entre les anciennes données Gauss-Krüger et le WGS84 moderne diffère véritablement selon qu'un point se situe en ancien territoire ouest-allemand ou est-allemand, un détail que de nombreux outils de conversion ignorent silencieusement. C'est un petit fragment concret de l'histoire de la guerre froide encore intégré aux données cartographiques allemandes.

La réunification de 1990 n'a pas effacé le Gauss-Krüger du jour au lendemain. L'Allemagne officielle n'est passée à l'ETRS89/UTM comme nouveau système de référence légal que progressivement, État par État, au cours des années 1990 et 2000, et certains bureaux cadastraux n'ont achevé la bascule que bien avancés dans les années 2010. Des millions d'anciens relevés de levé, de cartes de réseaux et de parcelles cadastrales n'ont jamais été reprojetés et n'existent encore que sous forme de coordonnées Gauss-Krüger. Quiconque numérise d'anciennes cartes allemandes, des infrastructures de réseaux ou des registres fonciers historiques finit tôt ou tard par se heurter au Gauss-Krüger, ce qui explique précisément pourquoi un convertisseur comme celui-ci reste utile.

WGS84 et UTM : conçus pour une ère satellitaire

Le WGS84 vient d'un tout autre horizon que le Gauss-Krüger. Personne n'a levé un seul royaume pour le produire. Le département de la Défense des États-Unis l'a construit pour rendre le positionnement par satellite possible en premier lieu. La géodésie satellitaire naissante des années 1960 et 1970 montra qu'aucun datum national existant n'était assez précis ni assez cohérent à l'échelle de la planète entière pour prendre en charge un seul système de positionnement. Le réseau de levé de chaque pays comportait ses propres petites erreurs, optimisées localement, utiles pour ce pays seul. La tâche à accomplir nécessitait un seul cadre de référence global se comportant de la même manière partout sur Terre, sans mosaïque régionale sous-jacente.

Le WGS84 a été publié en 1984 et a été discrètement affiné à plusieurs reprises depuis, la dernière fois en 2013, chaque révision resserrant l'ajustement entre l'ellipsoïde de référence et la forme réelle de la Terre à mesure que les données de poursuite satellitaire s'amélioraient. Le nom et les coordonnées du quotidien sont restés stables tout du long. C'est véritablement un système global unique, et non un ensemble de réseaux régionaux ajustés séparément puis assemblés comme l'est effectivement le Gauss-Krüger.

L'UTM, le système Universal Transverse Mercator, précède le WGS84 de plusieurs décennies. L'armée américaine l'a normalisé dans les années 1940, mais il repose sur la même idée fondamentale déjà mise au point par Gauss et Krüger : diviser le globe en zones étroites de Mercator transverse pour maintenir une faible déformation. L'UTM le fait simplement à l'échelle mondiale, en zones uniformes de 6 degrés numérotées de 1 à 60, plutôt que dans les zones Gauss-Krüger allemandes, plus étroites, de 3 degrés. Une coordonnée décrite comme « UTM sur WGS84 » applique les mêmes mathématiques de Mercator transverse que Gauss a développées dans les années 1820, simplement sur un ellipsoïde global moderne au lieu de Bessel 1841, avec des zones plus larges.

L'Europe a adopté l'ETRS89, un cadre de référence fixé à l'Europe fonctionnellement identique au WGS84 pour presque tout usage SIG, comme norme dans les années 1990. Les coordonnées officielles actuelles de l'Allemagne sont l'ETRS89/UTM : zone 32N à l'ouest, zone 33N à l'est.

Les mathématiques : comment la projection fonctionne réellement

Je vais parcourir cela comme je l'expliquerais à un étudiant en SIG, et non comme cela est habituellement présenté — un mur de lettres grecques sans aucune explication de ce que tout cela signifie physiquement.

Étape 1 — Connaître son ellipsoïde

Toute projection part d'un modèle mathématique de la forme de la Terre : un ellipsoïde, légèrement aplati aux pôles. Deux nombres le définissent complètement : le demi-grand axe a (le rayon équatorial) et l'aplatissement f. Pour Bessel 1841, l'ellipsoïde sur lequel les données Gauss-Krüger sont construites :

• a = 6,377,397.155 mètres
• 1/f = 299.1528128 (donc f ≈ 0.0033427731)

À partir de a et f, toutes les autres grandeurs des formules ci-dessous sont dérivées, la plus importante étant le carré de la première excentricité, e² = 2f − f², et le carré de la seconde excentricité, e'² = e²/(1−e²). Ces deux nombres reviennent constamment par la suite.

Étape 2 — Longueur de l'arc méridien : l'ossature de toute la projection

Le bloc de construction le plus important est la longueur de l'arc méridien, M(φ) — littéralement, la distance qu'il faudrait parcourir le long de la surface de la Terre, plein nord depuis l'équateur, pour atteindre la latitude φ. Parce que la Terre est un ellipsoïde et non une sphère, cette distance ne varie pas linéairement avec φ, et il n'existe pas d'expression simple sous forme close pour la calculer. Elle est plutôt calculée comme une série convergente en 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φ) ]

Chaque terme après le premier est une correction pour l'aplatissement ellipsoïdal, et chaque terme successif est plus petit que le précédent d'environ un facteur e² (environ 0.0067 pour Bessel 1841) — c'est pourquoi quatre termes suffisent déjà à obtenir une précision d'une fraction de millimètre sur toute la plage des latitudes allemandes. C'est le fondement du Nord : tout ce qui s'y ajoute est une correction pour la distance à laquelle vous vous trouvez à l'est ou à l'ouest du méridien central.

Étape 3 — Les termes de correction est-ouest

Trois autres grandeurs, toutes calculées à votre latitude φ, pilotent la correction pour la distance en longitude par rapport au méridien central :

• N = a / √(1 − e²sin²φ) — le rayon de courbure de la grande normale (la distance qu'un degré de longitude couvre réellement à cette latitude)
• T = tan²φ
• C = e'²cos²φ

Avec l comme différence de longitude par rapport au méridien central en radians, les équations directes complètes (Snyder, 1987) sont :

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⁶

Cela paraît intimidant à première lecture, mais la structure est plus simple qu'il n'y paraît : c'est un développement en série de Taylor en puissances de l, corrigeant l'approximation simple de Terre plate (N·cosφ·l) pour tenir compte de la façon dont la courbure de l'ellipsoïde infléchit les choses à mesure que l'on s'éloigne du méridien central. À l'intérieur des demi-zones étroites de 1.5 degré du Gauss-Krüger, les termes en l⁵ et l⁶ contribuent déjà pour moins d'un millimètre — la série converge extrêmement vite à cette échelle.

Étape 4 — La transformation de datum, avant ou après la projection

Tout ce qui précède se déroule entièrement au sein d'un seul ellipsoïde. Pour passer de Bessel 1841 (Gauss-Krüger) à GRS80 (WGS84), une étape distincte est nécessaire : une transformation de Helmert à 7 paramètres, appliquée dans l'espace géocentrique 3D XYZ plutôt que directement sur la latitude/longitude. Elle consiste à convertir les coordonnées géographiques en XYZ géocentriques, à appliquer trois translations, trois rotations minuscules et une correction d'échelle, puis à reconvertir :

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 )

Mon avis en tant que praticien : c'est l'étape que presque tout le monde rate, non pas parce que la formule est difficile, mais parce qu'on la saute entièrement et qu'on essaie de réutiliser directement l'Est/le Nord du Gauss-Krüger comme s'ils se trouvaient sur l'ellipsoïde WGS84. Ce n'est jamais le cas. Les mathématiques de projection et la transformation de datum sont deux opérations indépendantes, et toutes deux comptent.

Quelle zone et quelle région choisir ?

• Ancienne Allemagne de l'Ouest : Bade-Wurtemberg, Bavière, Brême, Hambourg, Hesse, Basse-Saxe, Rhénanie-du-Nord-Westphalie, Rhénanie-Palatinat, Sarre, Schleswig-Holstein
• Ancienne Allemagne de l'Est : Berlin, Brandebourg, Mecklembourg-Poméranie-Occidentale, Saxe, Saxe-Anhalt, Thuringe

Vous ne savez pas laquelle s'applique ? Convertissez d'abord votre point avec l'une ou l'autre région sélectionnée, puis cliquez sur « Utiliser la zone GK/UTM correcte pour cet endroit » — l'outil applique une règle géographique simple (environ : à l'est de 10,5°E et au nord de 50,8°N correspond à l'ancienne Allemagne de l'Est) pour choisir automatiquement la zone et la région correctes.

Calculé à la main : conversion de la porte de Brandebourg

Des nombres réels, pas des valeurs fictives. Voici la chaîne de calcul exacte que cet outil exécute, présentée étape par étape, à partir de la coordonnée WGS84 bien documentée de la porte de Brandebourg.

Point de départ (WGS84) : latitude 52.516266°N, longitude 13.377775°E.

Étape 1 : géographique → géocentrique XYZ (sur GRS80/WGS84)

Conversion en coordonnées cartésiennes 3D centrées sur la Terre :

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

Étape 2 : appliquer la transformation de Helmert inverse (paramètres Allemagne de l'Ouest)

En utilisant l'ensemble de paramètres inversés de l'Allemagne de l'Ouest (dx=−598.1, dy=−73.7, dz=−418.2, rx=−0.202″, ry=−0.045″, rz=+2.455″, s=−6.7 ppm) pour revenir du WGS84 vers le DHDN :

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

Étape 3 : géocentrique → géographique (sur Bessel 1841)

Reconversion de X', Y', Z' en latitude/longitude, cette fois en utilisant l'ellipsoïde Bessel 1841 :

DHDN latitude  = 52.51767297°
DHDN longitude = 13.37951396°

Remarquez que ce n'est pas la même chose que la coordonnée WGS84 — la latitude du DHDN est d'environ 0.0014° plus au nord et sa longitude d'environ 0.0017° plus à l'est. Ce minuscule décalage, environ 150-190 mètres sur le terrain, résulte entièrement de l'effet cumulé d'un ellipsoïde différent et d'un cadre de référence historique différent. C'est exactement pourquoi on ne peut pas simplement réétiqueter des coordonnées WGS84 en DHDN et considérer l'affaire close.

Étape 4 : projection sur la zone 4 du Gauss-Krüger (méridien central 12°E)

Avec φ = 52.51767297° (0.9166063 rad) et une différence de longitude par rapport au méridien central de l = 1.379514° (0.0240771 rad), les grandeurs intermédiaires donnent :

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

En intégrant ces valeurs dans les séries d'Est et de Nord de la section mathématique ci-dessus :

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

Le faux Est de la zone 4 est 4,500,000 (le « 4 » en tête marque la zone), ajouté à l'Est brut ci-dessus :

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

À titre de comparaison : le même point en UTM zone 33N

En faisant repasser la même coordonnée WGS84 d'origine par les formules de Mercator transverse — cette fois directement sur l'ellipsoïde GRS80, sans transformation de datum nécessaire, méridien central 15°E, facteur d'échelle 0.9996 — on obtient :

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

Les deux résultats décrivent exactement le même point physique sur Terre. La différence de près de 3.7 millions de mètres entre les deux valeurs de Nord n'a rien à voir avec une distance — elle correspond purement à la différence entre la convention de faux Est/zone du Gauss-Krüger et celle de l'UTM, plus le léger décalage dû à la projection sur Bessel 1841 par rapport à GRS80. C'est précisément le genre de confusion que ce convertisseur existe pour éliminer.

Étape par étape : utiliser le convertisseur

Étape 1

Choisissez votre format source sous « Convertir depuis » et votre format cible sous « Convertir vers ». Pour le Gauss-Krüger, choisissez également la zone correcte (2 à 5, correspondant au méridien central sur lequel vos données ont été levées) et si le point se situe en ancienne Allemagne de l'Ouest ou de l'Est.

Étape 2

Pour un point unique, passez à l'onglet « Coordonnée unique » et saisissez les deux valeurs — Rechtswert/Hochwert, latitude/longitude, ou est/nord, selon le format source que vous avez choisi.

Étape 3

Cliquez sur Convertir. Le résultat apparaît immédiatement en dessous, accompagné de la coordonnée WGS84 utilisée comme étape intermédiaire du calcul, et un repère se place sur la carte pour vérifier l'emplacement d'un coup d'œil — une zone mal sélectionnée par erreur de frappe devient généralement évidente immédiatement ici.

Étape 4

Passez à l'onglet « Fichier GeoJSON » pour convertir un fichier entier en une seule fois. Téléversez-le ou collez directement son contenu, définissez les mêmes formats source/cible, puis cliquez sur Convertir et prévisualiser.

Étape 5

La forme transformée s'affiche sur la carte pour que vous puissiez confirmer qu'elle se trouve au bon endroit, et le panneau indique exactement combien de paires de coordonnées ont été converties. Cliquez sur « Télécharger converted.geojson » pour enregistrer le résultat — chaque propriété d'origine de chaque entité est préservée intacte.

Qui a réellement besoin de ce convertisseur

Sans ordre particulier, d'après ce que j'ai réellement vu les gens faire avec un outil comme celui-ci :

1. Les compagnies de services publics qui numérisent d'anciennes cartes de réseaux. Les plans de réseaux de gaz, d'eau et d'électricité datant d'avant les années 2000 sont très souvent encore en Gauss-Krüger, et moderniser un SIG de gestion d'actifs implique de convertir des décennies de géométrie de canalisations et de câbles sans perdre en précision.
2. Les géomètres qui réconcilient d'anciens registres cadastraux avec des levés GNSS modernes. Une borne limite levée en 1975 et une autre mesurée avec un récepteur GPS-RTK la semaine dernière doivent être comparées dans le même système de coordonnées avant que quiconque puisse dire si elles concordent.
3. Les archivistes et historiens qui numérisent d'anciennes cartes allemandes. Les cartes topographiques et cadastrales d'avant-guerre et de l'époque de la guerre froide portent des grilles Gauss-Krüger ; convertir les entités de cartes numérisées dans un format SIG moderne exige exactement cette transformation.
4. Les consultants en environnement travaillant avec d'anciens réseaux de stations de surveillance. Les réseaux de surveillance des eaux souterraines, des sols et de la qualité de l'air mis en place il y a des décennies font souvent encore référence à l'emplacement des stations en Gauss-Krüger dans leur documentation d'origine.
5. Les développeurs logiciels qui s'intègrent aux portails de données ouvertes du gouvernement allemand. Certains exports GeoJSON de données ouvertes allemandes, plus anciens ou régionaux, contiennent encore des coordonnées Gauss-Krüger dans un format techniquement réservé au WGS84, et une carte web les affichera silencieusement au mauvais endroit, généralement quelque part dans le golfe de Guinée, jusqu'à ce qu'elles soient converties.
6. Les planificateurs de projets éoliens et solaires évaluant d'anciennes données de postes de transformation et de raccordement au réseau. Les données d'infrastructure réseau antérieures à l'essor des énergies renouvelables sont fréquemment encore dans l'ancienne projection.
7. Les étudiants universitaires en géographie et en topographie qui apprennent la géodésie. Le Gauss-Krüger est encore enseigné dans les cours allemands de Vermessungskunde précisément en raison de son importance historique, et parce que travailler la transformation à la main, comme dans l'exemple ci-dessus, est la façon dont le concept s'ancre vraiment.
8. Les consultants SIG qui vérifient la cohérence des datums dans la base de données spatiale existante d'un client. Les bases de données à datums mixtes, où certaines couches sont discrètement encore en DHDN et d'autres en ETRS89 sans que personne n'indique laquelle est laquelle, constituent un problème très courant et très coûteux à découvrir tardivement, que cet outil aide à diagnostiquer rapidement.

FAQ technique

Réponses volontairement spécifiques plutôt que génériques — si vous en savez déjà assez pour poser ces questions, vous méritez une vraie réponse.