Changement de projection
Il est possible de passer d’un système de projection vers un autre avec cette commande, par exemple pour passer du Lambert zones au Lambert 93 (décret du 28 décembre 2000 concernant le RGF93). Remarquez que le système altimétrique reste inchangé (NGF 69).
La boîte de dialogue ci-contre permet le paramétrage de cette transformation.
Coordonnées dans le système de départ
Ce groupe n’est disponible que si le mode manuel est activé, sinon c’est la projection courante qui est sélectionnée.
Coordonnées dans le système d’arrivée
Permet de choisir la projection d’arrivée pour la transformation.
Mode de fonctionnement
Vous pouvez utiliser cette boîte de dialogue selon deux modes distincts.
- « Saisie manuelle des coordonnées » : vous spécifiez la projection de départ, les coordonnées à transformer (Est, Nord et Altitude) et la projection d’arrivée. L’action sur le bouton « Calculer » affiche les coordonnées transformées dans le groupe « Coordonnées dans le système d’arrivée ».
« Calcul de tous les points de la GéoBase courante » : la projection de départ est celle de la GéoBase en cours. Choisissez la projection d’arrivée et appuyez sur le bouton « Calculer ». Tous les points de la GéoBase sont transformés. L’information sur le système de projection est inscrite dans la GéoBase. Un contrôle de la latitude des points est effectué et le message ci-contre apparaît.
REMARQUE :
Vous n’avez pas accès au bouton « Calculer » si l’une des projections est « Aucune » (pas de projection).
Grilles et altitudes
Ce groupe d’options permet d’activer ou de désactiver l’utilisation des grilles planimétriques et altimétriques. Certaines transformations, pour gagner en précision, utilisent des grilles fournissant des valeurs interpolées pour les paramètres de la transformation. C’est notamment le cas pour les passages du Lambert Zones en Lambert 93.
« Tenir compte de l’altitude » : Si cette option est sélectionnée, l’altitude des points est utilisée lors des calculs. Dans ce cas, il est possible d’utiliser la grille altimétrique permettant le calcul de la hauteur au-dessus de l’ellipsoïde à partir de l’altitude du point.
« Utiliser la grille planimétrique » : Si cette case est cochée, la grille fournissant Tx, Ty et Tz est utilisée lors du passage des coordonnées cartésiennes du système1 vers les coordonnées cartésiennes du système2. Sinon ce sont les paramètres standards qui sont utilisés.
Le fait de n’utiliser aucune des grilles revient à effectuer la transformation standard dans le logiciel CIRCE de l’IGN (France).
Format de la Grille planimétrique de l’IGN
La grille de l’IGN (Gr3df97a.txt) est un fichier texte dont voici un extrait :
Les quatre premières lignes décrivent le format du fichier :
- GR3D : ce sont des codes IGN indiquant le système de départ, le système d’arrivée, les unités pour les coordonnées géographiques et le méridien origine.
- GR3D1 : longitude min.; longitude max.; latitude min.; latitude max.; pas en longitude; pas en latitude (ici 0.1° dans les deux sens).
- GR3D2 : mode d’interpolation : INTERPOLATION BILINEAIRE (cf. exemple ci-après)
- GR3D3 : codes de précision de la transformation utilisés dans l’avant-dernière colonne du tableau.
Les lignes suivantes représentent chacune un point de la grille. Les coordonnées géographiques sont lues sur l’ellipsoïde GRS80, donc relatif au système RGF93.
Il existe aussi une forme plus condensée de cette grille (Gr3df97a.mnt) qui est également lue par COVADIS.
Le but de la grille est d’obtenir pour un point connu en longitude et latitude sur l’ellipsoïde GRS80, les paramètres TX, TY et TZ de la transformation. Pour cela on utilise l’interpolation bilinéaire.
Soit l’exemple suivant :
REMARQUE :
Si le système de départ est NTF (projections Lambert zone), et le système d’arrivée RGF93 (projection Lambert 93), le programme commence par calculer les coordonnées cartésiennes approchées dans le repère lié à l’ellipsoïde GRS 80 en appliquant les paramètres standards (TX = -168, TY = -60 et TZ = 320). Ensuite, il effectue la transformation des coordonnées cartésienne vers les coordonnées géographiques, ce qui donne un point d’entrée dans la grille. Il est ensuite possible d’obtenir des paramètres TX, TY et TZ plus précis en appliquant l’interpolation bilinéaire.
Format de la Grille altimétrique de l’IGN
La grille altimétrique de l’IGN (RAF20.tac pour le continent et rac09.mnt pour la Corse) est un fichier dont voici un extrait :
-5.5 8.5 42.0 51.5 0.0333333333333 0.025
53.5019 99 53.5002 99 53.4965 99 53.4887 99 53.4753 99 53.4571 99 53.4364 99 53.4172 99 53.4013 99 53.3889 99 53.3782 99 53.3668 99 53.3526 99 53.3353 99 53.3158 99 53.2952 99 53.2730 99 53.2482 99 53.2196 99
La première ligne fournit la longitude min.(ici -5.5°), la longitude max.(ici 8.5°), la latitude min.(ici 42°), la latitude max.(ici 51.5°), le pas en longitude (ici 0.0333333333333°) et le pas en latitude (ici 0.025°).
Les lignes suivantes contiennent des couples de valeurs de l’ondulation et code de précision.
L’ondulation est la position altimétrique d’un point de l’ellipsoïde GRS 80 par rapport au Géoïde (surface de niveau zéro pour les altitudes). On a alors la relation suivante entre la hauteur ellipsoïdale he, l’altitude normale H et l’ondulation N :
H = he - N
Cette relation permet entre autres de faire du nivellement avec du GPS.
Algorithme de la conversion
Pour effectuer la conversion il est impératif de calculer les coordonnées géographiques (longitude, latitude et hauteur ellipsoïdale). S’il y a un changement d’ellipsoïde (cas d’une transformation Lambert zone <=> Lambert 93) il faut en plus convertir les coordonnées géographiques en coordonnées cartésiennes. Le programme effectue ces transformations deux fois (une fois dans chaque sens). Le diagramme suivant résume les différentes transformations effectuées.
Schéma des étapes principales lors d’une conversion :
