• Главная
• Вопросы и ответы

Сравнение преобразований координат в разном ПО

Как считают используя одни и те же параметры разные программные пакеты.

Тесты сгруппированы по программному продукту, каждый тест, кроме названия продукта включает также порядковый номер показанный римской цифрой, например ArcGIS I, ArcGIS II и т. д.

Тестируемая исходная точка: 50 гр. с.ш., 50 гр. в.д., исходная система координат: WGS84 (shapefile, покрытие Arcinfo).

Результаты приводятся в формате: долгота широта (что-то третье у cs2cs?). Если результатом является также векторный набор данных, он также делается доступным для скачки. Атрибутивная таблица векторного набора данных должна содержать два поля: Lat-new и Long-new с рассчитанными значениями широты и долготы.

Проверяются 3-х параметрические преобразования.

• Arcview Projection Utility
• ArcGIS
• Arcinfo Workstation
• ERDAS
• GeoCalc
• ogr2ogr
• cs2cs
• Trackmaker
• Формулы пересчета

Arcview Projection Utility

Arcview Projection Utility I

Результат: 49.9999631298 50.0014646551 (загрузить)
Комментарий: Параметры: 28, -130, -95, Geocentric translation.

ArcGIS

ArcGIS I

Pulkovo_1942_To_WGS_1984

Результат: 49.99996312980 50.0014646515 (загрузить)
Комментарий: Абсолютно идентичен Arcview Projection Utility. Настройка трансформации: 28, -130, -95, Geocentric translation.
Источник: Supported geographic (datum) transformations (geographic_transformations.pdf)

ArcGIS II

Pulkovo_1942_To_WGS_1984_13

Результат: 49.9998716889 50.0013449067 (загрузить)
Комментарий: Настройка трансформации: 27, -135, -84.5, 0, 0, 0.554, 0.2263, Position vector.
Источник информации о преобразовании: Supported geographic (datum) transformations (geographic_transformations.pdf)

Arcinfo Workstation

Arcinfo Workstation I

Projection GEOGRAPHIC
Units DD
datum user_defined 0 0 0 0 0 0 0
Parameters 6378137.00 6356752.31
output
projection geogr
units dd
Datum user_defined 28 -130 -95 0 0 0 0
Parameters 6378245.0 6356863.019

Результат: 49.9999618530 50.0014648438

Комментарий:

ERDAS IMAGINE

ERDAS IMAGINE

 "Pulkovo 1942-1 (Russia)" 28 -130 -95 0 0 0 0

Результат: 49.9999631385638 50.0014646841

Комментарий: Совпадение с Arcview Projection Utility c точностью до 7 знака.

Источник информации о преобразовании: %ERDAS%/etc/spheroid.tab

ogr2ogr

ogr2ogr I

ogr2ogr -s_srs "+proj=latlong +ellps=WGS84 +datum=WGS84" -t_srs "+proj=latlong +ellps=krass +towgs84=23.92,-141.27,-80.9,0,0,0,0" -f "ESRI Shapefile" pul-ogr-I.shp wgs84-5050.shp

Результат: 49.9998041685 50.0015220974 (загрузить)
Комментарий: полностью аналогичен cs2cs II, что и должно быть. Убедившись что результат у ogr и cs2cs один и тот же (что логично, так как они используют одну и туже библиотеку трансформации), используем дальше cs2cs, чтобы не генерировать лишних shape-файлов.

cs2cs

cs2cs I

cs2cs +proj=latlong +ellps=WGS84 +to +proj=latlong +ellps=krass +towgs84=23.92,-141.27,-80.9,0,0,0,0 -f '%.10f'

Результат: 50 50 0
Комментарий: При переходе между системами координат обязательно требуется указание датума, например +datum=WGS84

cs2cs II

cs2cs +proj=latlong +ellps=WGS84 +datum=WGS84 +to +proj=latlong +ellps=krass +towgs84=23.92,-141.27,-80.9,0,0,0,0 -f '%.10f'

Результат: 49.9998041685 50.0015220974 12.0732616797

cs2cs III

cs2cs +proj=latlong +ellps=WGS84 +datum=WGS84 +to +proj=latlong +ellps=krass +towgs84=28,-130,-95,0,0,0,0 -f '%.10f'

Результат: 49.9999631698 50.0014646551 15.639
Комментарий: Результат равен результату трансформации с помощью Arcview Projection Utility с точностью до 7 знака.

GeoCalc

О GeoCalc

Параметры: dX=-28, dY=130, dZ=95, остальные 0
Результат: 49.9999631301 50.0014646551 (15.6339)
Комментарий: Результат равен cs2cs III и Arcview Projection Utility I с точностью до 7 знака (обсуждение). Инвертация параметров сдвига.

Trackmaker

Результат: 49.9999631388026 50.0014646841461
Комментарий: Результат равен cs2cs III и Arcview Projection Utility I с точностью до 7 знака.

Формулы пересчета

О формулах

Параметры: dX=28, dY=-130, dZ=-95, остальные 0
Результат: 49.999963139 50.001464672
Комментарий: Результат равен cs2cs III и Arcview Projection Utility I с точностью до 7 знака.

Обсудить в форуме Комментариев — 5

Последнее обновление: September 09 2021

Дата создания: 02.12.2009
Автор(ы): Максим Дубинин