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

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:

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?

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.

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