Message posté par : loic58471 (loic58471(a)gmail.com)
----------------------------------------
Je réponds ici à ma propre question de départ afin d'exposer la solution trouvée.
A la base je voulais simplement afficher mes données en repère local sur ma projection
ADINDAN en indiquant des correspondances de coordonnées entre différents points, en
système local puis en système national. Je me suis rendu compte que conceptuellement ça ne
fonctionnait pas car la relation n'est pas linéaire entre mes système de coordonnées
qui possèdent chacun une déformation géométrique propre à la notion de projection.
Il me fallait donc passer par le WGS84, comme l'explique bien Christian au dessus.
J'ai donc fini par créer ma propre projection et non une équivalence entre des
systèmes de coordonnées déjà projetés.
Le type de projection le plus adapté m'a paru être la projection oblique mercator qui
me permet d'orienter le nord de ma projection comme je veux, j'ai l'impression
que la plupart des autres sont naturellement orienté au nord par défaut.
Tout ensuite est assez bien expliqué sur la page web de Proj:
https://proj.org/operations/projections/omerc.html
Il me fallait donc simplement déterminer comme paramètres:
-"+lonc" et "+lat_0" les coordonnées WGS84 du centre de ma projection,
que j'ai placé au centre de mon site sur un point de coordonnées connues dans les
différents systèmes, les déformations géométriques dues à la projections seront ainsi
imperceptibles à notre échelle avant des centaines de km; mon site faisant tout au plus
20km d'ampleur, cela convenait,
- gamma et alpha dont je n'ai vraiment pas compris la raison d'être en duo,
semblant agir exactement sur la même chose c'est à dire la rotation du cylindre de
projection selon la relation:
rotation (sens horaire) = alpha - gamma
les valeurs sont en degrés.
- "x_0" et "y_0" sont les coordonnées dans mon repère local du centre
de projection dont j'ai indiqué les coordonnées wgs84 ci dessus
Cela me donne en définitive la formule proj4 suivante:
+proj=omerc +lonc=29.33132998 +lat_0=20.23684520 +gamma=0.03 +alpha=0 +x_0=5449.08
+y_0=20202.74 +ellps=WGS84
Ça fonctionne vraiment bien, le problème maintenant c'est de régler les décalages,
fixer les incohérences et erreurs de mesure topographiques mais ça c'est une autre
histoire.
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=358686#p358686
Pour y répondre : qgis_fr(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net