Próbuję przekształcić punkt WGS84 lat / lon
11d42'32.10629" E 5d12'56.75371" S
w trzech różnych pakietach oprogramowania (Proj4, GeoTrans i Leica GeoOffice), ale otrzymuję bardzo duże różnice między wynikami - około kilkuset metrów w X i Y! Zauważ, że te odmiany pojawiają się tylko z 7 parametrami, a nie z 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 + południe + jednostki = 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
Wyniki są odpowiednio:
- 800392 9422525
- 800306 9422840
- 800941 9422891
Wierzę, że wszystkie 3 pakiety używają tych samych metod matematycznych dla transformacji 7-param (metoda Bursa-Wolfa). Co może być przyczyną tej ogromnej różnorodności?
coordinate-system
proj
datum
Powietrzny jeźdźca
źródło
źródło
Odpowiedzi:
Po pierwsze, Proj4 korzysta z tego, co EPSG nazywa „wersją wektora pozycji” metody 7-parametrowej. Możliwe, że GeoTrans i Leica GeoOffice używają innej wersji, którą EPSG nazywa „Ramką współrzędnych”. Obie metody są równoważne, ale macierze obrotu są różne i znaki parametrów kątowych muszą zostać zmienione.
Po drugie, dziękuję za udostępnienie zrzutu ekranu definicji transformacji w Leica GeoOffice. Definicja Proj4 lub ta definicja jest niepoprawnie zdefiniowana. Parametry definiujące Ellipsoid A i Ellipsoid B powinny zostać przełączone. Obecnie transformacja ta przekształca się z WGS84 w Congo60. Zauważ, że w Proj4 opcja jest + towgs84, więc jest zdefiniowana jako FROM Congo60 DO WGS84. Aby zmienić kierunek w definicji, zmień znaki WSZYSTKICH parametrów. Sprawdź także pomoc dla GeoOffice i sprawdź, czy parametr SF chce części na milion wersji, czy też już przekonwertowanej na współczynnik skali.
Nie wiem o GeoTrans - masz na myśli oprogramowanie NGA? W każdym razie, mam nadzieję, że uda ci się dopasować GeoOffice i Proj4.
źródło
W przypadku proj4 znaki parametrów muszą zostać odwrócone.
Zobacz tę stronę definicji:
http://www.spatialreference.org/ref/epsg/62826405/prettywkt/
Nie jestem pewien, skąd masz swoje wartości.
QGIS definiuje EPSG: 28232 Point Noire UTM 32S jako:
a wynik 7 parametrów nie powinien być od tego daleki.
EDYTOWAĆ
Dzięki GDAL 1.10 otrzymuję następujące wyniki:
Czyli mniej niż 3 metry od siebie. Zauważ, że E i S muszą podążać za stopniami bez odstępów.
I dla przypomnienia: ustawienie towgs84 na zero powoduje:
podczas przekształcania ze stopni Point Noire EPSG: 4262 na Point Noire UTM powoduje:
źródło