Mapinfo и PROJ4

[Boris]

В обсуждении viewtopic.php?f=34&t=15911
Есть ответ:

Однако QGIS не MapInfo. Не забываем менять знаки параметров вращения:

Код: Выделить всё

…+towgs84=-118.754,-61.782,-93.237,2.40896,3.47502,-1.29688,6.5177…

Я его встречал уже несколько раз, но откуда берется это утверждение?
Мне нечем проверить конкретные параметры местной СК, но вот с того же сайта:

"МСК-02 зона 2", 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, 7, 58.х, 0, 1, 2300000, -х

Имеем принятый большинством де-факто расчетный вариант перехода СК-1942 -> WGS84 имени ГОСТ-2008:

Код: Выделить всё

Суммы
СК-42 — WGS-84    +23.57 −140.95 −79.8   0 −0.35 −0.79 −0.22

там же для Mapinfo (и далее по интернету)

Код: Выделить всё

"Долгота / Широта (Пулково 1942–WGS ГОСТ 51794-2008)", 1, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0

(и далее по интернету)

Код: Выделить всё

"МСК-1964 СПб ГОСТ 51794-2008", 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, 7, 30, 0, 1, 95942.17, -6552810.0

на импортном сайте то же самое - везде минусы в знаках вращения осей:

Код: Выделить всё

+proj=longlat +ellps=krass +towgs84=23.57,-140.95,-79.8,0,-0.35,-0.79,-0.22 +no_defs 

[trir]

http://geodesist.ru/forum/threads/%D0%9 ... 616/page-2
разные формулы преобразования: Helmert 7-parameter transformation и Bursa-Wolf transformation

[Игорь Белов]

откуда берется это утверждение?

Из суровой реальности, Boris, из суровой реальности. Если скормить один и тот же набор параметров разным ГИСам, получим два возможных варианта пересчёта. Познакомьтесь вот со скетчем Чака Тэйлора High-Accuracy Datum Transformations. Автор различает преобразования Гельмерта и Бурша-Вольфа.
В иных терминах это были бы "Coordinate Frame Rotation" и "Position Vector transformation". Во всяком случае, такие названия используются в именах операций EPSG. Но вот беда, EPSG — просто большой агрегатор, и его наполнители не видят разницы между значениями терминов. Они не различают операции 1032 и 1033, 1037 и 1038, 9606 и 9607. Соответственно, при конвертации параметров между разными представлениями они ничтоже сумняшеся переносят знаки as is. В итоге случается, что некоторые программы для ряда датумов используют один вид преобразования с параметрами для другого (GlobalMapper замечен в этом).

"Coordinate Frame Rotation" — вращается система координат.

"Position Vector transformation" — вращается Вселенная вокруг координатной оси.

Тэйлор принимает положительным направление вращения по правилу левого винта, поэтому для него преобразование Гельмерта соответствует вращению координатных осей, а Бурша-Вольфа — вращению вектора положения. Всякий, кому ближе правило буравчика, может с полным основанием считать иначе.

Как теперь жить на свете? Как узнать, какой вид преобразования использует данная программа? Да очень просто. Имеются источники с формулами и наборами параметров к ним (например, ГОСТ Р 51794). Есть калькулятор/Excel/компилятор/интерпретатор (выбрать по вкусу). Выберем несколько точек на глобусе и пересчитываем в другую СК, а результаты расчётов сравниваем с деятельностью ГИС-программ. По некоторым причинам мне нравится такой набор из четырёх точек:

Код: Выделить всё

-180 0
-90 0
0 0
90 0

Создаём в экзаменуемом ГИСе проекцию, например, Меркатора на анализируемом датуме и вперёд, к сияющим вершинам познания.