Я пытаюсь преобразовать точку широты / долготы WGS84
11d42'32.10629" E 5d12'56.75371" S
в трех разных программных пакетах (Proj4, GeoTrans и Leica GeoOffice), но я получаю очень большие различия между результатами - около нескольких сотен метров по X и Y! Обратите внимание, что эти вариации появляются только с 7 параметрами, а не с 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 + юг + юниты = m + no_defs
Геотранс
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
Результаты соответственно:
- 800392 9422525
- 800306 9422840
- 800941 9422891
Я полагаю, что все 3 пакета используют одни и те же математические методы для преобразования с 7 параметрами (метод Бурса-Вольфа). Что может быть причиной этого огромного изменения?
coordinate-system
proj
datum
WindRider
источник
источник
Ответы:
Во-первых, Proj4 использует то, что EPSG называет «позиционным вектором» версии метода с 7 параметрами. Возможно, что GeoTrans и Leica GeoOffice используют другую версию, которую EPSG назвала «Координатный фрейм». Оба метода эквивалентны, но матрицы вращения различны, и знаки угловых параметров должны быть изменены.
Во-вторых, спасибо за предоставленный снимок экрана определения преобразования в Leica GeoOffice. Либо определение Proj4, либо это определение определено неправильно. Определяющие параметры Эллипсоид A и Эллипсоид B должны быть переключены. В настоящее время это преобразование преобразуется из WGS84 в Congo60. В Proj4, обратите внимание, что опция + towgs84, поэтому она определена как ОТ Congo60 TO WGS84. Чтобы изменить направление в определении, измените знаки ВСЕХ параметров. Также проверьте справку для GeoOffice и посмотрите, хочет ли параметр SF количество частей на миллион в версии или значение, уже преобразованное в масштабный коэффициент.
Я не знаю о GeoTrans - вы имеете в виду программное обеспечение NGA? В любом случае, мы надеемся, что вы сможете настроить GeoOffice и Proj4.
источник
Для proj4 знаки параметров необходимо поменять местами.
Смотрите эту страницу определения:
http://www.spatialreference.org/ref/epsg/62826405/prettywkt/
Я не уверен, откуда вы взяли свои ценности.
QGIS определяет EPSG: 28232 Point Noire UTM 32S как:
и результат 7 параметров не должен быть далеко от этого.
РЕДАКТИРОВАТЬ
С GDAL 1.10 я получаю следующие результаты:
Так что на расстоянии менее 3 метров. Обратите внимание, что E и S должны следовать градусам без пробелов.
И, просто для записи: установка towgs84 в ноль приводит к:
при преобразовании степеней Point Noire EPSG: 4262 в Point Noire UTM получается:
источник