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

Обсуждение материалов сайта: вопросы, замечания, предложения
updates-bot
Bot
Сообщения: 276
Зарегистрирован: 03 фев 2008, 23:13
Репутация: 3

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

Сообщение updates-bot » 16 авг 2009, 05:22

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

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

Boris
Гуру
Сообщения: 3927
Зарегистрирован: 10 апр 2006, 22:34
Статьи: 3
Проекты: 1
Репутация: 349
Откуда: Париж

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

Сообщение Boris » 03 окт 2009, 23:48

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

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 8896
Зарегистрирован: 06 окт 2003, 20:20
Статьи: 231
Проекты: 12/6
Репутация: 640
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

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

Сообщение Максим Дубинин » 04 окт 2009, 05:48

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

Аватара пользователя
JEY
Активный участник
Сообщения: 228
Зарегистрирован: 17 июл 2008, 13:42
Статьи: 1
Репутация: 1

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

Сообщение JEY » 05 окт 2009, 13:23

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

kmikhailov
Новоприбывший
Сообщения: 1
Зарегистрирован: 17 мар 2011, 09:53
Репутация: 0
Откуда: Tyumen
Контактная информация:

СК-95 -> WGS-84

Сообщение kmikhailov » 17 мар 2011, 12:11

Хороший набор!
По работе мне нужно описать датум СК-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? Откликнитесь пожалуйста!
Вложения
РУКск95.pdf
Страницы руководства с параметрами
(482.6 КБ) 3242 скачивания

old_Bibigon
Участник
Сообщения: 78
Зарегистрирован: 25 янв 2010, 14:14
Проекты: 2
Репутация: 3

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

Сообщение old_Bibigon » 07 май 2011, 19:49

Однако, параметр +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
Участник
Сообщения: 78
Зарегистрирован: 25 янв 2010, 14:14
Проекты: 2
Репутация: 3

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

Сообщение old_Bibigon » 07 май 2011, 23:04

И до кучи про 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
Активный участник
Сообщения: 240
Зарегистрирован: 22 апр 2011, 08:10
Репутация: 27
Откуда: Хабаровск

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

Сообщение Institor » 27 июл 2011, 02:47

В 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
Новоприбывший
Сообщения: 3
Зарегистрирован: 17 сен 2011, 17:52
Репутация: 0

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

Сообщение grishin77 » 17 сен 2011, 20:59

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

grishin77
Новоприбывший
Сообщения: 3
Зарегистрирован: 17 сен 2011, 17:52
Репутация: 0

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

Сообщение grishin77 » 17 сен 2011, 21:10

Максим Дубинин писал(а):я поправил на 10-6 везде, но признаюсь честно, я до конца не понимаю, он всегда в частях на миллион? В формуле он используется вот так:
Изображение
если не сложно напишите как эти матрицы раскладываются

kefi
Интересующийся
Сообщения: 46
Зарегистрирован: 20 ноя 2011, 23:04
Репутация: 0

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

Сообщение kefi » 04 дек 2011, 12:12

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

Донецков
Гуру
Сообщения: 3051
Зарегистрирован: 19 май 2010, 19:44
Репутация: 187

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

Сообщение Донецков » 04 дек 2011, 13:28

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

kefi
Интересующийся
Сообщения: 46
Зарегистрирован: 20 ноя 2011, 23:04
Репутация: 0

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

Сообщение kefi » 06 дек 2011, 17:26

Донецков писал(а): но если работаешь на конкретной локальной территории, то для более точного пересчета следует вычислить локальные параметры и использовать их...
А как ВЫЧИСЛИТЬ эти локальные параметры ? (Кстати, тот ГОСТ Р 51794-2008 дает для WGS <=> ПЗ90, а не для СК-42)

Донецков
Гуру
Сообщения: 3051
Зарегистрирован: 19 май 2010, 19:44
Репутация: 187

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

Сообщение Донецков » 06 дек 2011, 21:51

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

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

kefi
Интересующийся
Сообщения: 46
Зарегистрирован: 20 ноя 2011, 23:04
Репутация: 0

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

Сообщение kefi » 06 дек 2011, 22:27

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

Ответить

Вернуться в «Материалы сайта»