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:
- 800392 9422525
- 800306 9422840
- 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?
coordinate-system
proj
datum
Cavaleiro do Vento
fonte
fonte
Respostas:
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.
fonte
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/
Não sei de onde você tirou seus valores.
QGIS define EPSG: 28232 Point Noire UTM 32S como:
e o resultado de 7 parâmetros não deve estar longe disso.
EDITAR
Com o GDAL 1.10, recebo os seguintes resultados:
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:
ao transformar os graus EPSG do ponto Noire: 4262 para o Point Noire UTM resulta em:
fonte