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

12

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

Captura de tela


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?

Cavaleiro do Vento
fonte
No Geo Office, o elipsóide Clarke 1880 IGN é uma definição personalizada? Você pode postar seus parâmetros, seja ou não?
precisa saber é o seguinte
Recebi a captura de tela da Leica de outras pessoas e ainda não sei.
WindRider 28/09

Respostas:

9

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.

mkennedy
fonte
1
O valor Y (norte / latitude) pode ser porque os valores Clarke 1880 não são exatamente os mesmos, mas 60 m parecem muito grandes.
Mcknedy
1
@mkennedy: você pode fazer os mesmos cálculos que fiz em Arcgis para ver se podemos alinhar sem a Leica?
ANDREJ
1
@AndreJoost, com certeza. Usando as versões Esri do EPSG 1802 e 28232, obtive o leste: 800230.139 ao norte: 9423133.413.
precisa saber é o seguinte
1
Tentei várias variantes para ver se conseguia reproduzir o resultado da Leica e não consegui. O mesmo problema que você encontrou: Desativado por mais de 60 m.
Mcknedy # 26/13
1
Estou confiante de que as soluções Proj4 (e ArcGIS!) Estão corretas. A fonte do Proj4 está aberta. Embora o código do ArcGIS não esteja aberto, eu lido com ele há cerca de 15 anos. Também adicionei uma transformação personalizada ao Geotrans 3.0 e obtive o mesmo resultado do Proj4 / ArcGIS.
Mcknedy # 26/13
4

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
AndreJ
fonte
Como você pode ver, com os mesmos parâmetros, há uma enorme diferença (> 600 m!) Da Leica. Preciso saber por que isso acontece e como obter os mesmos resultados. Poderia ser do método de conversão de / para coordenadas geocêntricas usadas no Proj4? Tem fórmulas muito complicadas. Talvez a Leica use outro método?
WindRider 23/09
1
A Leica oferece conversão de 3 parâmetros? E o que eles calculam para isso?
ANDREJ