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
Les résultats sont respectivement:
- 800392 9422525
- 800306 9422840
- 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?
la source
Réponses:
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.
la source
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/
Je ne sais pas d'où vous tenez vos valeurs.
QGIS définit EPSG: 28232 Point Noire UTM 32S comme:
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:
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:
lors de la transformation des degrés Point Noire EPSG: 4262 en Point Noire UTM se traduit par:
la source