Grandes diferenças entre os resultados da transformação de dados de 7 parâmetros

Estou tentando transformar um ponto lat / lon WGS84

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

em três pacotes de software diferentes (Proj4, GeoTrans e Leica GeoOffice), mas recebo grandes diferenças entre os resultados - cerca de algumas centenas de metros em X e Y! Observe que essas variações aparecem apenas com 7 parâmetros e não com 3.

Proj4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + para + proj = utm + zone = 32 + ellps = clrk80 + towgs84 = 178.3.316.7.131.5, -5.278, -6.077, -10.9792, -19.166 + sul + unidades = 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

Os resultados são respectivamente:

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

Eu acredito que todos os 3 pacotes usam os mesmos métodos matemáticos para a transformação 7-param (método Bursa-Wolf). Qual poderia ser a causa dessa enorme variação?

Primeiro, o Proj4 usa o que o EPSG chama de versão "Position Vector" do método de 7 parâmetros. É possível que o GeoTrans e o Leica GeoOffice usem a outra versão que o EPSG chamou de "Quadro de coordenadas". Ambos os métodos são equivalentes, mas as matrizes de rotação são diferentes e os sinais dos parâmetros angulares precisam ser alterados.

Segundo, obrigado por fornecer a captura de tela da definição de transformação no Leica GeoOffice. A definição do Proj4 ou esta definição está definida incorretamente. Os parâmetros de definição Ellipsoid A e Ellipsoid B devem ser alterados. Atualmente, essa transformação é convertida de WGS84 para Congo60. No Proj4, observe que a opção é + towgs84; portanto, é definida como FROM Congo60 to WGS84. Para alterar a direção na definição, altere os sinais de TODOS os parâmetros. Verifique também a ajuda do GeoOffice e veja se o parâmetro SF deseja as partes por milhão de versão ou o valor já convertido em um fator de escala.

Não conheço a GeoTrans - você quer dizer o software NGA? De qualquer forma, espero que você consiga que o GeoOffice e o Proj4 correspondam.

Para o proj4, os sinais dos parâmetros devem ser revertidos.

Veja esta página de definição:

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"]]

Não sei de onde você tirou seus valores.

QGIS define EPSG: 28232 Point Noire UTM 32S como:

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

e o resultado de 7 parâmetros não deve estar longe disso.

EDITAR

Com o GDAL 1.10, recebo os seguintes resultados:

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

Então, a menos de 3 metros de distância. Observe que os E e S devem seguir os graus sem espaços em branco.

E, apenas para constar: definir towgs84 como zero resulta em:

800310.94   9422829.37 -109.32

ao transformar os graus EPSG do ponto Noire: 4262 para o Point Noire UTM resulta em:

800311.21   9422892.49 0.00