D'énormes différences entre les résultats de la transformation de données à 7 paramètres

12

J'essaie de transformer un point WGS84 lat / lon

11d42'32.10629" E 5d12'56.75371" S

dans trois logiciels différents (Proj4, GeoTrans et Leica GeoOffice) mais je reçois de très grandes différences entre les résultats - environ quelques centaines de mètres en X et Y! Notez que ces variations n'apparaissent qu'avec 7 paramètres et non avec 3.

Proj4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + to + proj = utm + zone = 32 + ellps = clrk80 + towgs84 = 178.3,316.7,131.5, -5.278, -6.077, -10.9792, -19.166 + sud + unités = m + no_defs


GeoTrans

Delta X 178.3 
Delta Y 316.7 
Delta Z 131.5 
Rot X -5.278 
Rot Y -6.077 
Rot Z -10.9792 
SF = -19.166 / 0.999980834 (0.000019166)

Leica GeoOffice

Capture d'écran


Les résultats sont respectivement:

  1. 800392 9422525
  2. 800306 9422840
  3. 800941 9422891

Je crois que tous les 3 packages utilisent les mêmes méthodes mathématiques pour la transformation à 7 paramètres (méthode Bursa-Wolf). Quelle pourrait être la cause de cette énorme variation?

Coureur du vent
la source
Dans Geo Office, l'ellipsoïde Clarke 1880 IGN est-il une définition personnalisée? Pouvez-vous publier ses paramètres, que ce soit ou non?
mkennedy
J'ai reçu la capture d'écran de Leica d'autres personnes et je ne sais pas encore vraiment.
WindRider

Réponses:

9

Tout d'abord, Proj4 utilise ce que EPSG appelle la version "Position Vector" de la méthode à 7 paramètres. Il est possible que GeoTrans et Leica GeoOffice utilisent l'autre version qu'EPSG a appelée "Coordinate Frame". Les deux méthodes sont équivalentes, mais les matrices de rotation sont différentes et les signes des paramètres angulaires doivent être modifiés.

Deuxièmement, merci d'avoir fourni la capture d'écran de la définition de transformation dans Leica GeoOffice. Soit la définition Proj4 ou cette définition est définie de manière incorrecte. Les paramètres de définition Ellipsoid A et Ellipsoid B doivent être commutés. Actuellement, cette transformation passe du WGS84 au Congo60. Dans Proj4, notez que l'option est + towgs84, elle est donc définie comme FROM Congo60 TO WGS84. Pour changer la direction dans la définition, changez les signes de TOUS les paramètres. Consultez également l'aide de GeoOffice et voyez si le paramètre SF veut la version pièces par million ou la valeur de facteur déjà convertie en échelle.

Je ne connais pas GeoTrans - parlez-vous du logiciel NGA? Quoi qu'il en soit, j'espère que vous pourrez faire correspondre GeoOffice et Proj4.

mkennedy
la source
1
La valeur Y (nord / latitude) pourrait être due au fait que les valeurs de Clarke 1880 ne sont pas tout à fait les mêmes, mais 60 m semblent trop grands.
mkennedy
1
@mkennedy: pouvez-vous faire les mêmes calculs que j'ai fait dans Arcgis pour voir si nous pouvons nous aligner sans Leica?
AndreJ
1
@AndreJoost, bien sûr. En utilisant les versions Esri d'EPSG 1802 et 28232, j'ai obtenu l'abscisse: 800230.139 ordonnée: 9423133.413.
mkennedy
1
J'ai essayé plusieurs variantes pour voir si je pouvais reproduire le résultat de Leica et je n'ai pas pu. Même problème que vous avez trouvé: Y off de 60+ m.
mkennedy
1
Je suis convaincu que les solutions Proj4 (et ArcGIS!) Sont correctes. La source Proj4 est ouverte. Bien que le code ArcGIS ne soit pas ouvert, je l'ai traité pendant environ 15 ans. J'ai également ajouté une transformation personnalisée à Geotrans 3.0 et obtenu le même résultat que Proj4 / ArcGIS.
mkennedy
4

Pour proj4, les signes des paramètres doivent être inversés.

Voir cette page de définition:

http://www.spatialreference.org/ref/epsg/62826405/prettywkt/

GEOGCS["Pointe Noire (deg)",
    DATUM["Congo 1960 Pointe Noire",
        SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936269,
            AUTHORITY["EPSG","7011"]],
        TOWGS84[-178.3,-316.7,-131.5,5.278,6.077,10.979,3.953271276531849],
        AUTHORITY["EPSG","6282"]],
    PRIMEM["Greenwich",0.0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.017453292519943295],
    AXIS["Geodetic latitude",NORTH],
    AXIS["Geodetic longitude",EAST],
    AUTHORITY["EPSG","62826405"]]

Je ne sais pas d'où vous tenez vos valeurs.

QGIS définit EPSG: 28232 Point Noire UTM 32S comme:

+proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-148,51,-291,0,0,0,0 +units=m +no_defs

et le résultat des 7 paramètres ne devrait pas être loin de cela.


ÉDITER

Avec GDAL 1.10, j'obtiens les résultats suivants:

4326-proj-3 parameters:
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-148,51,-291,0,0,0,0 +units=m +no_defs
800232.21   9423131.96 -1.76
4326-proj-7 parameters from spatialrefrence.org:
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-178.3,-316.7,-131.5,5.278,6.077,10.979,3.953271276531849 +units=m +no_defs
800230.13   9423133.46 91.31
4326-proj-7 parameters from proj4 datum_shift.csv (EPSG:1802):
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-178.3,-316.7,-131.5,5.278,6.077,10.979,19.166 +units=m +no_defs
800230.13   9423133.40 -5.72

Donc à moins de 3 mètres l'un de l'autre. Notez que les E et S doivent suivre les degrés sans espaces.

Et, pour mémoire: la mise à zéro de towgs84 entraîne:

800310.94   9422829.37 -109.32

lors de la transformation des degrés Point Noire EPSG: 4262 en Point Noire UTM se traduit par:

800311.21   9422892.49 0.00
AndreJ
la source
Comme vous pouvez le voir, avec les mêmes paramètres, il y a une énorme différence (> 600 m!) Avec Leica. J'ai besoin de savoir pourquoi cela se produit et comment obtenir les mêmes résultats. Serait-ce à partir de la méthode de conversion de / vers les coordonnées géocentriques utilisée dans Proj4? Il a des formules très compliquées. Peut-être que Leica utilise une méthode différente?
WindRider
1
Leica propose-t-il une conversion à 3 paramètres? Et que calculent-ils pour cela?
AndreJ