Estoy tratando de transformar un punto lat / lon WGS84
11d42'32.10629" E 5d12'56.75371" S
en tres paquetes de software diferentes (Proj4, GeoTrans y Leica GeoOffice), pero recibo grandes diferencias entre los resultados: ¡unos pocos cientos de metros en X e Y! Tenga en cuenta que estas variaciones solo aparecen con 7 parámetros, no con 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 + sur + 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
Los resultados son respectivamente:
- 800392 9422525
- 800306 9422840
- 800941 9422891
Creo que los 3 paquetes usan los mismos métodos matemáticos para la transformación de 7 parámetros (método Bursa-Wolf). ¿Cuál podría ser la causa de esta gran variación?
coordinate-system
proj
datum
Jinete del viento
fuente
fuente
Respuestas:
Primero, Proj4 usa lo que EPSG llama la versión "Vector de posición" del método de 7 parámetros. Es posible que GeoTrans y Leica GeoOffice usen la otra versión que EPSG llamó "Marco de coordenadas". Ambos métodos son equivalentes, pero las matrices de rotación son diferentes y los signos de los parámetros angulares deben cambiarse.
En segundo lugar, gracias por proporcionar la captura de pantalla de la definición de transformación en Leica GeoOffice. La definición de Proj4 o esta definición se define incorrectamente. Los parámetros de definición Elipsoide A y Elipsoide B deben cambiarse. Actualmente, esta transformación se convierte de WGS84 a Congo60. En Proj4, tenga en cuenta que la opción es + towgs84, por lo que se define como FROM Congo60 TO WGS84. Para cambiar la dirección en la definición, cambie los signos de TODOS los parámetros. Consulte también la ayuda de GeoOffice y vea si el parámetro SF quiere las partes por millón de versiones o el valor de factor de escala ya convertido.
No sé sobre GeoTrans, ¿te refieres al software NGA? De todos modos, es de esperar que GeoOffice y Proj4 coincidan.
fuente
Para proj4, los signos de los parámetros deben invertirse.
Ver esta página de definición:
http://www.spatialreference.org/ref/epsg/62826405/prettywkt/
No estoy seguro de dónde obtuviste tus valores.
QGIS define EPSG: 28232 Point Noire UTM 32S como:
y el resultado de los 7 parámetros no debería estar muy lejos de eso.
EDITAR
Con GDAL 1.10, obtengo los siguientes resultados:
Por lo tanto, a menos de 3 metros de distancia. Tenga en cuenta que E y S deben seguir los grados sin espacios en blanco.
Y, solo para el registro: establecer towgs84 en cero da como resultado:
mientras se transforma de Point Noire grados EPSG: 4262 a Point Noire UTM resulta en:
fuente