Suite

Conversion des coordonnées des points avec un système inconnu

Conversion des coordonnées des points avec un système inconnu


J'utilise QGis depuis quelques semaines maintenant. Bien que ce site ait toujours été très utile, je ne trouve pas de solution à mon problème récent :

J'ai peu de points, et connais leurs coordonnées dans les systèmes "WGS-84" et "ISN-93". Ces points ont également des coordonnées spécifiées dans le système "Local".

Le problème est qu'il y a aussi peu de points qui n'ont que des coordonnées "Locales", et j'aimerais trouver leurs coordonnées "ISN-93".

Comment puis-je le faire?

*par exemple:

point1 ISN93 (685320.253.420614.733) Local (-3537.239.2284.487)

point2 ISN93 (685841.253,420586,536) Local (-3016,528,2248.842)

point3 ISN93 (?,?) Local (-5187.996,2584.508)

J'essayais d'utiliser une géométrie simple, mais lorsque je mesure la distance entre les points, le système ISN93 donne des longueurs légèrement différentes de Local (la différence varie entre les longueurs jusqu'à 0,01% de sa valeur), donc ça ne sert à rien.


Une approche de force brute très simple :

  1. prenez les premières coordonnées des points et convertissez-les d'ISN93 / Lambert 1993 en degrés WGS84 en utilisant cs2cs dans le shell OSGEO4W

    cs2cs +init=epsg:3057 +à +init=epsg:4326 -f "%%.8f" WGS84.txt

  2. construire un CRS omerc local sur ce point, avec les "coordonnées locales" comme faux Est et Nord

  3. convertir les coordonnées du point dans ce CRS personnalisé

    cs2cs +init=epsg:3057 +to +proj=omerc +lonc=-15.1778877053125 +lat_0=64.2376322211738 +alpha=0.01 +gamma=0 +k=1.0 +x_0=-3537.239 +y_0=2284.487 +datum=WGS84 +units= m +no_defs -f "%%.8f" local.txt

  4. calculer le décalage X et Y en mètres pour le deuxième point

  5. joue avec+alphaet+kvaleurs pour minimiser les décalages

De cette façon, j'obtiens avec un décalage de 2 mm :

+proj=omerc +lonc=-15.1778877053125 +lat_0=64.2376322211738 +alpha=2.646 +gamma=0 +k=1.00032 +x_0=-3537.239 +y_0=2284.487 +datum=WGS84 +units=m +no_defs

Je suppose que vos coordonnées sont dans l'ordre est-nord et que vous n'avez pas à vous soucier des mouvements tectoniques entre ISN93 et ​​WGS84. Sur un fond Openstreetmap, votre grille locale ressemble à ceci :