из WGS84 (EPSG:4326) в Гауса-Крюгера 24 зона

[vlalexey]

Пробую разобраться в причинах отличия пересчета пары точек из EPSG:4326 в Гауса-Крюгера 24 зона разными инструментами:

1. Беру пару точек в epsg:4326 (территория Тюменской области)
2. Пересчитываю их в Гауса-Крюгера 1942 24 зона через Mapinfo 11 (будем считать эти результаты эталонными)
3. Эти же точки [1] пересчтываю в ГК 42 зона 24 через QGIS
получаю в 2 и 3 практически идентичные результаты.

Иду дальше:
4. Конвертирую точки [1] также в ГК 42 зона 24 через ogr2ogr, задавая итоговую СК в proj4 формате. Получаю также похожие результаты.
5. Конвертирую точки [1] также в ГК 42 зона 24 через ogr2ogr, задавая итоговую СК в WKT формате. Получаю отличие на десятки метров.
В чем ошибка? Как добиться, чтобы в [5] считалось тоже сравнимо?

Данные в виде таблицы в Google docs

Пересчет через ogr2ogr делаю так:
1. ogr2ogr.exe -f "Mapinfo File" -t_srs gk12zone24_wkt.txt out.mif in.mif
2. ogr2ogr.exe -f "Mapinfo File" -t_srs gk12zone24_зкщо4.txt out.mif in.mif

Файлы исходные и .txt с описанием СК - в архиве

[trir]

В WKT нет TOWGS84 +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1

[vlalexey]

А TOWGS84 является обязательным для задания Гауса-Крюгера 1942 зона 24?

Может можно через несколько преобразований выйти от 4326 к ней, используя WKT?

UPD:
Спасибо за наводку - изменил определение СК в WKT (добавил TOWGS84 блок c парамтерами из proj4) и добился верного результата.

Вопрос закрыт

Было:

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

PROJCS["Transverse_Mercator",GEOGCS["GCS_Krassovsky, 1942",DATUM["D_unknown",SPHEROID["krass",6378245,298.3]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",69],PARAMETER["scale_factor",1],PARAMETER["false_easting",12500000],PARAMETER["false_northing",0],UNIT["Meter",1]]

Стало:

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

PROJCS["Transverse_Mercator",GEOGCS["GCS_Krassovsky, 1942",DATUM["D_unknown",SPHEROID["krass",6378245,298.3],TOWGS84[24,-123,-94,0.02,-0.25,-0.13,1.1]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",69],PARAMETER["scale_factor",1],PARAMETER["false_easting",12500000],PARAMETER["false_northing",0],UNIT["Meter",1]]

[trir]

обязательный вообще для Крассовского!
если бы было можно, то его вообще не предумали, тем боллее, что его нет в спецефикации

[SergeyRyzhkov]

vlalexey - это хорошо что Вы разобралсиь!
Главное понимать почему так.

На всякий случай

WGS84 (4326) - географическая СК на эллипсоиде WGS84
"Гауса-Крюгера 24 зона" (как Вы ее назвали) - "это проецированная" СК на эллипсоиде Крассовского.
Важно понимать, что проекция использует тот или иной эллипсоид для конвертации географических координат в проецированные.

Соответственно какой должен быть алгоритм пересчета?

1. Из географических координат WGS84 пересчитать в географические координаты на Крассовском
2. Полученные географические координаты в п1 пересчитать (проецировать) в Гауса-Крюгера

Пункт № 1 иногда называют трансформацией датумов.
Так вот TOWGS84 как раз и определяет параметры пересчета для п1.
У Вас не задан был TOWGS84 поэтому проецирование было произведено: координата полученная на эллипсоиде WGS84 проецировалась на другой эллипсоид (Крассовского)

[Донецков]

И нужно помнить, то MapInfo далеко не эталон, более верным является не создание нового определения СК с параметрами MI (TOWGS84[24,-123,-94,0.02,-0.25,-0.13,1.1]) в том же QGIS, а наоборот добавлением пользовательских параметров в MI, но уже с параметрами по ГОСТ (что тоже не идеал, но все ближе к "истине")...
И должно быть уточнение, что это для СК42, т.к. для СК-95 уже будет по другому, хотя 24 зона есть там и там...

[vlalexey]

Спасибо за ответ и пояснения.

(Из текущих результатов у меня еще один вопрос возникает - но оформлю его лучше отдельным вопросом в ветке форума)