Переход от одной СК к другой - наборы параметов

[updates-bot]

Обсуждение статьи "Переход от одной системы координат к другой - наборы параметров"

http://gis-lab.info/qa/datum-transform-sets.html

[Boris]

Спасибо за наборы,
но мне кажется, что в 7-ми параметрическом преобразовании надо все привести к некоторому стандарту в параметре m - либо давать его с множителем 1.0e-06, либо нет, но тогда указать, что он имеет это множитель.
А то получается, что имеется расхождение в разных записях в миллион раз.

[Максим Дубинин]

я поправил на 10-6 везде, но признаюсь честно, я до конца не понимаю, он всегда в частях на миллион? В формуле он используется вот так:

[JEY]

Во всех статьях, описывающих понятие трансформации Гельмерта, m дается, как некоторое число, умноженное на 1.0e-06. Например в статье на Википедии указано, что m всегда делится на 1,000,000.

[kmikhailov]

Хороший набор!
По работе мне нужно описать датум СК-95. Софт (Finder Schlumberger) использует "опорное" преобразование WGS-84, а параметров СК-95 в WGS-84 я нигде не мог найти. Через формулы трансформации и параметры Хелмерта ПЗ90-ВГС84, СК95-ПЗ90 рассчитал параметры прямого преобразования СК-95 в WGS-84, которые оказались один в один как в вашем наборе: 24.82 -131.21 -82.66 0 0 -0.16 -0.12*10-6.

После того, натыкаюсь на то же преобразование (СК-95 в WGS-84) в документе:
РУКОВОДСТВО пользователя по выполнению работ в системе координат 1995 года (СК-95)
ГКИНП (ГНТА)–06-278-04
ФЕДЕРАЛЬНАЯ СЛУЖБА ГЕОДЕЗИИ И КАРТОГРАФИИ РОССИИ
???
На странице 70 указанного руководства (google it или см. во вложении) приведены совершенно иные параметры преобразования:
DX = +24.653 м; DY = - 129.136 м; DZ = - 83.057 м;
RX = - 0.06696”; RY= +0.00391, RZ= - 0.12902”;
Scale = - 0.175×10 - 6.

Еще там упоминается ITRF (что это?).
Цитирую страницу 71:
"Фактические различия между WGS84 и ITRF находятся в пределах дециметра. Внутренние деформации СК-95 на всей протяженности ГГС составляют несколько дециметров. Поэтому для преобразования между СК-95 и WGS84 систему ITRF можно считать идентичной системе WGS84. Точность такого преобразования в целом по всей ГГС будет находиться в пределах нескольких дециметров и в основном определяться внутренними деформациями ГГС в СК-95."

Спрашивается, чему верить, что использовать? Судя по всему в руководстве указано преобразование СК-95 -> ITRF, правда с оговоркой, что его можно считать идентичным СК-95 в WGS-84.
С одной стороны, параметры (24.82 -131.21 -82.66 0 0 -0.16 -0.12*10-6) рассчитываются математически из известных преобразований Хелмерта (ПЗ90-ВГС84, СК95-ПЗ90), которые имеются в ГОСТ 51794. С другой стороны, официальный документ федеральной службы геодезии...

Кто-нибудь использует СК-95 -> WGS-84? Откликнитесь пожалуйста!

[old_Bibigon]

Однако, параметр +towgs84 для cs2cs для Mapinfo 1001 на деле имеет значение (24,-123,-94,0.02,-0.25,-0.13,1.1) вместо (24,-123,-94,-0.02,0.25,0.13,1.1)

Имеем MapInfo 10.1, создаём файл в проекции epsg:4284 (широта/долгота (Пулково 1942)). Вбиваем точку с координатами (47,58), сохраняем как epsg:4326.
Идём в пингвина, ogrinfo про проекцию TAB файла с СК-42 сразу говорит:

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

GEOGCS["unnamed",
    DATUM["Pulkovo_1942",
        SPHEROID["Krassovsky",6378245,298,3],
        TOWGS84[24,-123,-94,0.02,-0.25,-0.13,1.1]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0,0174532925199433]]

Ладно, достаём координаты точки для WGS84 - 46,998322999999999 58,000185999999999.
Проверяем параметры через cs2cs:

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

$ echo 47 58 | cs2cs +proj=latlong +ellps=krass +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +no_defs -f '%.10f' +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
46.9983231620	58.0001860093 -1.7521861829

$ echo 47 58 | cs2cs +proj=latlong +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,1.1 +no_defs -f '%.10f' +to +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
46.9982459642	58.0000834596 -1.7866630729

Видим что ogrinfo был прав :-/

Это какая-то багофича libproj что ей надо знак менять или не совсем верная запись MapInfo в её Internal datum parameters used in MapInfo Professional?
Версия gdal 1.7.1, libproj - 4.5.0

[old_Bibigon]

И до кучи про MapInfo, похоже в 1012 и 1014 знаки тоже попутаны.

---- Долгота / Широта (Россия ПЗ90), 1012 ----
Проекция по ogrinfo

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

GEOGCS["unnamed",
    DATUM["Russia_PZ90",
        SPHEROID["PZ90",6378136,298.257839303],
        TOWGS84[-1.08,-0.27,-0.9,-0,-0,0.16,-0.12]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

Проверка точки:

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

Mapinfo -  47.000055000000003 58.000003
ogr - 47.000054687877316 58.000003234316061

---- Долгота / Широта (Россия СК42), 1013 ----

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

GEOGCS["unnamed",
    DATUM["Russia_SK42",
        SPHEROID["PZ90",6378136,298.257839303],
        TOWGS84[23.92,-141.27,-80.9,-0,0.35,0.82,-0.12]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

Mapinfo - 46.99819 58.000211999999998
ogr - 46.998189551519104 58.000211843840418

---- Долгота / Широта (Россия СК95), 1014 ----

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

GEOGCS["unnamed",
    DATUM["Russia_SK95",
        SPHEROID["PZ90",6378136,298.257839303],
        TOWGS84[24.82,-131.21,-82.66,-0,-0,0.16,-0.12]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

Mapinfo - 46.998224 58.000208999999998
ogr - 46.998224142163863 58.00020888720158

Конвертирование проводилось ogr'ом, т.к. добавлять элипсоид ПЗ-90 в libproj было лениво. А ogr что-то похожее достал, но в целом в результате не уверен

ps. А в 1013 знаки кто-то до меня попутал, в mapinfo'шной базе 23.92,-141.27,-80.9,0,-0.35,-0.82,-0.12

[Institor]

В 1013 используется эллипсоид ПЗ90 (ка это ни парадоксально, но так написано в справке).
Странная ситуация со знаками: если использовать ГОСТ (23.92-141.27 -80.9 0 0.35 0.82 -0.12*10^-6), ничего не получается. реально работает только (23.92, -141.27, -80.9, 0, -0.35, -0.82, -0.12) (судя по картматериалу в СК-42)

[grishin77]

И до кучи про MapInfo, похоже в 1012 и 1014 знаки тоже попутаны.

---- Долгота / Широта (Россия ПЗ90), 1012 ----
Проекция по ogrinfo

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

GEOGCS["unnamed",
    DATUM["Russia_PZ90",
        SPHEROID["PZ90",6378136,298.257839303],
        TOWGS84[-1.08,-0.27,-0.9,-0,-0,0.16,-0.12]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

Проверка точки:

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

Mapinfo -  47.000055000000003 58.000003
ogr - 47.000054687877316 58.000003234316061

---- Долгота / Широта (Россия СК42), 1013 ----

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

GEOGCS["unnamed",
    DATUM["Russia_SK42",
        SPHEROID["PZ90",6378136,298.257839303],
        TOWGS84[23.92,-141.27,-80.9,-0,0.35,0.82,-0.12]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

Mapinfo - 46.99819 58.000211999999998
ogr - 46.998189551519104 58.000211843840418

---- Долгота / Широта (Россия СК95), 1014 ----

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

GEOGCS["unnamed",
    DATUM["Russia_SK95",
        SPHEROID["PZ90",6378136,298.257839303],
        TOWGS84[24.82,-131.21,-82.66,-0,-0,0.16,-0.12]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]

Mapinfo - 46.998224 58.000208999999998
ogr - 46.998224142163863 58.00020888720158

Конвертирование проводилось ogr'ом, т.к. добавлять элипсоид ПЗ-90 в libproj было лениво. А ogr что-то похожее достал, но в целом в результате не уверен

ps. А в 1013 знаки кто-то до меня попутал, в mapinfo'шной базе 23.92,-141.27,-80.9,0,-0.35,-0.82,-0.12

никто ничего не путал .ИМХО если взять.ГОСТ получается именно так как написано в mapinfo datum . Ведь переход в строке towgs84 обозначает от искомой СК к WGS8, отсюда и знаки. только не понимаю, почему в СК42 сК95 используются элипсоиды PZ-90, согласно ГОСТ вэтих СК должны использоваться элипсоид крассовского?

[grishin77]

я поправил на 10-6 везде, но признаюсь честно, я до конца не понимаю, он всегда в частях на миллион? В формуле он используется вот так:

если не сложно напишите как эти матрицы раскладываются

[kefi]

Господа специалисты, подскажите дилетанту ответ на простые вопросы :
Вообще принципиально известны или нет датумы для данной системы геодезических координат(СГК), на примере, пресловутой СК-42 ? Если нет , то - по каким причинам ?
Каков должен быть подход к выбору датума для данной СГК для данной местности на примере СК-42 ? Прав ли я в том, что датум одной СГК разный для разных местностей, если да, то - почему так ?

[Донецков]

На форуме уже неоднократно обращались к этому вопросу - например gis-lab.info/forum/viewtopic.php?f=34&t=8405 и gis-lab.info/forum/viewtopic.php?f=34&t=3469&start=30
Соответственно, по логике нужно использовать параметры из ГОСТа 2008 года..., но если работаешь на конкретной локальной территории, то для более точного пересчета следует вычислить локальные параметры и использовать их...

[kefi]

но если работаешь на конкретной локальной территории, то для более точного пересчета следует вычислить локальные параметры и использовать их...

А как ВЫЧИСЛИТЬ эти локальные параметры ? (Кстати, тот ГОСТ Р 51794-2008 дает для WGS <=> ПЗ90, а не для СК-42)

[Донецков]

В ГОСТе даны переходы WGS-ПЗ и ПЗ-СК42 и 95, соответственно найти параметры WGS-СК чистая математика путем сложения параметров с сохранением знака ...

По вычислению локальных параметров - нужны как полевые GPS измерения, должны быть известны каталожные значения координат определяемых пунктов... По самому алгоритму лучше обратиться к более знающим людям (на форуме когда-то очень активно выступал коллега из Сочи, явно имеет большой опыт в этом вопросе)...

[kefi]

Так а разве для каждой конкретной локальной территории не известны локальные параметры ( я имею ввиду - локальные datum'ы ). Т.е. я в общем не понимаю порядка вещей - что тогда дает ГОСТ , когда рассказывает , как нужно вычислять параметры перехода от одной к другой СК, ведь он же никак не регламентирует локальные регионы, где его применять ?