Переход от одной СК к другой - наборы параметов
-
- Bot
- Сообщения: 276
- Зарегистрирован: 03 фев 2008, 23:13
- Репутация: 3
Переход от одной СК к другой - наборы параметов
Обсуждение статьи "Переход от одной системы координат к другой - наборы параметров"
http://gis-lab.info/qa/datum-transform-sets.html
http://gis-lab.info/qa/datum-transform-sets.html
-
- Гуру
- Сообщения: 4231
- Зарегистрирован: 10 апр 2006, 22:34
- Репутация: -344969098
- Откуда: Париж
Re: Переход от одной СК к другой - наборы параметов
Спасибо за наборы,
но мне кажется, что в 7-ми параметрическом преобразовании надо все привести к некоторому стандарту в параметре m - либо давать его с множителем 1.0e-06, либо нет, но тогда указать, что он имеет это множитель.
А то получается, что имеется расхождение в разных записях в миллион раз.
но мне кажется, что в 7-ми параметрическом преобразовании надо все привести к некоторому стандарту в параметре m - либо давать его с множителем 1.0e-06, либо нет, но тогда указать, что он имеет это множитель.
А то получается, что имеется расхождение в разных записях в миллион раз.
- Максим Дубинин
- MindingMyOwnBusiness
- Сообщения: 9129
- Зарегистрирован: 06 окт 2003, 20:20
- Репутация: 748
- Ваше звание: NextGIS
- Откуда: Москва
- Контактная информация:
Re: Переход от одной СК к другой - наборы параметов
я поправил на 10-6 везде, но признаюсь честно, я до конца не понимаю, он всегда в частях на миллион? В формуле он используется вот так:


пристегивайтесь, турбулентность прямо по курсу
- JEY
- Активный участник
- Сообщения: 228
- Зарегистрирован: 17 июл 2008, 13:42
- Репутация: 1
Re: Переход от одной СК к другой - наборы параметов
Во всех статьях, описывающих понятие трансформации Гельмерта, m дается, как некоторое число, умноженное на 1.0e-06. Например в статье на Википедии указано, что m всегда делится на 1,000,000.
-
- Новоприбывший
- Сообщения: 1
- Зарегистрирован: 17 мар 2011, 09:53
- Репутация: 0
- Откуда: Tyumen
- Контактная информация:
СК-95 -> WGS-84
Хороший набор!
По работе мне нужно описать датум СК-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. Софт (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 КБ) 4790 скачиваний
-
- Участник
- Сообщения: 78
- Зарегистрирован: 25 янв 2010, 14:14
- Репутация: 3
Re: Переход от одной СК к другой - наборы параметов
Однако, параметр +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 сразу говорит:
Ладно, достаём координаты точки для WGS84 - 46,998322999999999 58,000185999999999.
Проверяем параметры через cs2cs:
Видим что ogrinfo был прав :-/
Это какая-то багофича libproj что ей надо знак менять или не совсем верная запись MapInfo в её Internal datum parameters used in MapInfo Professional?
Версия gdal 1.7.1, libproj - 4.5.0

Имеем 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]]
Проверяем параметры через 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
Это какая-то багофича libproj что ей надо знак менять или не совсем верная запись MapInfo в её Internal datum parameters used in MapInfo Professional?
Версия gdal 1.7.1, libproj - 4.5.0
-
- Участник
- Сообщения: 78
- Зарегистрирован: 25 янв 2010, 14:14
- Репутация: 3
Re: Переход от одной СК к другой - наборы параметов
И до кучи про MapInfo, похоже в 1012 и 1014 знаки тоже попутаны.
---- Долгота / Широта (Россия ПЗ90), 1012 ----
Проекция по ogrinfo
Проверка точки:
---- Долгота / Широта (Россия СК42), 1013 ----
---- Долгота / Широта (Россия СК95), 1014 ----
Конвертирование проводилось ogr'ом, т.к. добавлять элипсоид ПЗ-90 в libproj было лениво. А ogr что-то похожее достал, но в целом в результате не уверен 
ps. А в 1013 знаки кто-то до меня попутал, в mapinfo'шной базе 23.92,-141.27,-80.9,0,-0.35,-0.82,-0.12
---- Долгота / Широта (Россия ПЗ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
Код: Выделить всё
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
Код: Выделить всё
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

ps. А в 1013 знаки кто-то до меня попутал, в mapinfo'шной базе 23.92,-141.27,-80.9,0,-0.35,-0.82,-0.12
- Institor
- Активный участник
- Сообщения: 243
- Зарегистрирован: 22 апр 2011, 08:10
- Репутация: 27
- Откуда: Хабаровск
Re: Переход от одной СК к другой - наборы параметов
В 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)
Странная ситуация со знаками: если использовать ГОСТ (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)
-
- Новоприбывший
- Сообщения: 3
- Зарегистрирован: 17 сен 2011, 17:52
- Репутация: 0
Re: Переход от одной СК к другой - наборы параметов
никто ничего не путал .ИМХО если взять.ГОСТ получается именно так как написано в mapinfo datum . Ведь переход в строке towgs84 обозначает от искомой СК к WGS8, отсюда и знаки. только не понимаю, почему в СК42 сК95 используются элипсоиды PZ-90, согласно ГОСТ вэтих СК должны использоваться элипсоид крассовского?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]]
---- Долгота / Широта (Россия СК42), 1013 ----Код: Выделить всё
Mapinfo - 47.000055000000003 58.000003 ogr - 47.000054687877316 58.000003234316061
---- Долгота / Широта (Россия СК95), 1014 ----Код: Выделить всё
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
Конвертирование проводилось ogr'ом, т.к. добавлять элипсоид ПЗ-90 в libproj было лениво. А ogr что-то похожее достал, но в целом в результате не уверенКод: Выделить всё
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
ps. А в 1013 знаки кто-то до меня попутал, в mapinfo'шной базе 23.92,-141.27,-80.9,0,-0.35,-0.82,-0.12
-
- Новоприбывший
- Сообщения: 3
- Зарегистрирован: 17 сен 2011, 17:52
- Репутация: 0
Re: Переход от одной СК к другой - наборы параметов
если не сложно напишите как эти матрицы раскладываютсяМаксим Дубинин писал(а):я поправил на 10-6 везде, но признаюсь честно, я до конца не понимаю, он всегда в частях на миллион? В формуле он используется вот так:
-
- Интересующийся
- Сообщения: 46
- Зарегистрирован: 20 ноя 2011, 23:04
- Репутация: 0
Re: Переход от одной СК к другой - наборы параметов
Господа специалисты, подскажите дилетанту ответ на простые вопросы :
Вообще принципиально известны или нет датумы для данной системы геодезических координат(СГК), на примере, пресловутой СК-42 ? Если нет , то - по каким причинам ?
Каков должен быть подход к выбору датума для данной СГК для данной местности на примере СК-42 ? Прав ли я в том, что датум одной СГК разный для разных местностей, если да, то - почему так ?
Вообще принципиально известны или нет датумы для данной системы геодезических координат(СГК), на примере, пресловутой СК-42 ? Если нет , то - по каким причинам ?
Каков должен быть подход к выбору датума для данной СГК для данной местности на примере СК-42 ? Прав ли я в том, что датум одной СГК разный для разных местностей, если да, то - почему так ?
-
- Гуру
- Сообщения: 3058
- Зарегистрирован: 19 май 2010, 19:44
- Репутация: 189
Re: Переход от одной СК к другой - наборы параметов
На форуме уже неоднократно обращались к этому вопросу - например gis-lab.info/forum/viewtopic.php?f=34&t=8405 и gis-lab.info/forum/viewtopic.php?f=34&t=3469&start=30
Соответственно, по логике нужно использовать параметры из ГОСТа 2008 года..., но если работаешь на конкретной локальной территории, то для более точного пересчета следует вычислить локальные параметры и использовать их...
Соответственно, по логике нужно использовать параметры из ГОСТа 2008 года..., но если работаешь на конкретной локальной территории, то для более точного пересчета следует вычислить локальные параметры и использовать их...
-
- Интересующийся
- Сообщения: 46
- Зарегистрирован: 20 ноя 2011, 23:04
- Репутация: 0
Re: Переход от одной СК к другой - наборы параметов
А как ВЫЧИСЛИТЬ эти локальные параметры ? (Кстати, тот ГОСТ Р 51794-2008 дает для WGS <=> ПЗ90, а не для СК-42)Донецков писал(а): но если работаешь на конкретной локальной территории, то для более точного пересчета следует вычислить локальные параметры и использовать их...
-
- Гуру
- Сообщения: 3058
- Зарегистрирован: 19 май 2010, 19:44
- Репутация: 189
Re: Переход от одной СК к другой - наборы параметов
В ГОСТе даны переходы WGS-ПЗ и ПЗ-СК42 и 95, соответственно найти параметры WGS-СК чистая математика путем сложения параметров с сохранением знака ...
По вычислению локальных параметров - нужны как полевые GPS измерения, должны быть известны каталожные значения координат определяемых пунктов... По самому алгоритму лучше обратиться к более знающим людям (на форуме когда-то очень активно выступал коллега из Сочи, явно имеет большой опыт в этом вопросе)...
По вычислению локальных параметров - нужны как полевые GPS измерения, должны быть известны каталожные значения координат определяемых пунктов... По самому алгоритму лучше обратиться к более знающим людям (на форуме когда-то очень активно выступал коллега из Сочи, явно имеет большой опыт в этом вопросе)...
-
- Интересующийся
- Сообщения: 46
- Зарегистрирован: 20 ноя 2011, 23:04
- Репутация: 0
Re: Переход от одной СК к другой - наборы параметов
Так а разве для каждой конкретной локальной территории не известны локальные параметры ( я имею ввиду - локальные datum'ы ). Т.е. я в общем не понимаю порядка вещей - что тогда дает ГОСТ , когда рассказывает , как нужно вычислять параметры перехода от одной к другой СК, ведь он же никак не регламентирует локальные регионы, где его применять ?
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей