GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли


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

Обсуждение материалов сайта: вопросы, замечания, предложения

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

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

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

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

updates-bot
Bot
 
Зарегистрирован: 04 фев 2008
 
Сообщения: 276
Репутация: 2

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

Сообщение Boris » 04 окт 2009, 00:48

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

Boris
Гуру
 
Зарегистрирован: 10 апр 2006
Откуда: Париж (на карте)
 
Сообщения: 3819
Репутация: 317
Статьи: 3
Проекты: 1

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

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

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

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8834
Репутация: 608
Статьи: 231
Проекты: 12/6

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

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

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

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

СК-95 -> WGS-84

Сообщение kmikhailov » 17 мар 2011, 13: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? Откликнитесь пожалуйста!
  • 0

Вложения
РУКск95.pdf
Страницы руководства с параметрами
(482.6 КБ) Скачиваний: 3036
kmikhailov
Новоприбывший
 
Зарегистрирован: 17 мар 2011
Откуда: Tyumen (на карте)
 
Сообщения: 1
Репутация: 0

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

Сообщение old_Bibigon » 07 май 2011, 20: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
  • 0

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

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

Сообщение old_Bibigon » 08 май 2011, 00: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
  • 0

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

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

Сообщение Institor » 27 июл 2011, 03: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)
  • 0

Аватара пользователя
Institor
Активный участник
 
Зарегистрирован: 22 апр 2011
Откуда: Хабаровск (на карте)
 
Сообщения: 240
Репутация: 27

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

Сообщение grishin77 » 17 сен 2011, 21: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, согласно ГОСТ вэтих СК должны использоваться элипсоид крассовского?
  • 0

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

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

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

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


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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

След.

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

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 4


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).