Конвертация координат из МСК в WGS с помощью pyProj

Тема в разделе "Общие вопросы", создана пользователем Maloy_ksb, 13 дек 2022.

  1. Добрый день всем.
    Зашёл в тупик со своими изысканиями...
    Прошу помочь знающих людей.

    Ситуашия следующая:
    Имею точки в местной системе координат. С помощью Topcon Tools могу их конвертировать как в местный датум так и в WGS84. Получаю в итоге именно те координаты, которые должны быть (сравнивал с контроллером Topcon)
    Проекция TransverseMercator, ИГД с преобразованием по Гельмерту.

    В pyProj получаю на выходе правильные координаты в местном датуме, но не могу перевести в WGS84, разница порядка 12 метров.

    Эта строка дает верные координаты в местном датуме
    "+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=" + lon_0 + " +k=1 +x_0=" + x_0 + " +y_0=" + y_0 + " +ellps=WGS84 +step +proj=push +v_3 +step +proj=cart +ellps=WGS84 +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1"

    но добавив преобразование по гельмерту, координыт далеки от ожидаемых
    "+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=" + lon_0 + " +k=1 +x_0=" + x_0 + " +y_0=" + y_0 + " +ellps=WGS84 +step +proj=push +v_3 +step +proj=cart +ellps=WGS84 +step +proj=helmert +x=152.2213 +y=176.626 +z=277.808 +rx=218.862377 +ry=147.524412 +rz=658.798277 +s=-53.702102 +convention=coordinate_frame +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1"

    пробовал оба тива конвертации coordinate_frame и position_vector
    в первом случае расстояние между ожидаемыми и полученными координатами составляет 12 метров, во втором больше и явно не правильные.

    Прошу знатоков помочь в вопросе, я новичек в этом, да и профиль не мой. Я по части ИТ специализируюсь
     
  2. stout

    stout Форумчанин

    Что по вашему мнению есть местный датум, если в него можно конвертировать координаты?
    Разница между чем и чем? Что в качестве "истины"?
    А в контроллере откуда параметры взялись?

    Что такое rx, ry, rz и ы, В каких единицах?
     
    zvezdochiot нравится это.
  3. zvezdochiot

    zvezdochiot Форумчанин

    А откуда эта строка? И почему эллипсоид WGS84? Это что за "местная СК"?
     
    stout нравится это.
  4. cr2

    cr2 Форумчанин

    Это так теперь параметры +towgs84 полагается писать.
     
  5. stout

    stout Форумчанин

    Догадываюсь, но хотелось услышать от топистартера об источнике. Коме того, меня зело смущает сама цифирь.
     
  6. Всем добрый день.
    Спасибо за вопросы, попробую ответить на них, если получится конечно.
    1) местный датум (или местная ИГД), я могу ошибаться в терминах, в топкон выглядит это так.
    upload_2022-12-14_8-23-0.png
    Grid - координаты в МСК, Datum - кооринаты в местной ИГД (которые я в pyProj получаю без проблем), WGS84 - координаты в WGS84
    2) в качестве истины - координаты перобразованные с помошью topcon tools. Их я сравнивал с прибором, который используют маркшейдера для работы, все четко.
    3) контроллер был настроен специалистами, которые предоставили мне информацию о преобразованиях, эти параметры я уже внес в topcon tools.
    4) rx, ry, rz, в секундах.
    5) пытался использовать +towgs84 результат аналогичный
    6) элипс именно wgs84 потому, что расчеты преобразований из МСК в WGS84 были построены на нем
     
  7. zvezdochiot

    zvezdochiot Форумчанин

    Это что за МСК такая? В каком месте Луны она находится?
     
  8. stout

    stout Форумчанин

    Эти специалисты достойны быть если не академиками Азербайджанской академии наук, то членами-корреспондентами уж точно. (Или, мягко выражаясь, просто членами?) Мой стёб основан на одной хохме. Когда-то давным-давно в базе EPSG появились 7 параметров, для преобразования из СК-42 в WGS 84, которые были определены Азербайджанской академией наук. Там углы вращения составляли тысячи секунд дуги чего не может быть в принципе. В худшем случае угол меньше 6" (чел стоит на Гринвичском меридиане)
    [​IMG]
    Именно малость углов позволяет заменять синус угла самим углом, а косинус – единицей в полноценной матрице вращения Гельмерта (т.е. использовать преобразование Бурши-Вольфа). И только ПО от ляйки всегда работало с полноценной матрицей Гельмерта (Хотя с этим возникает свой затык, все определяют углы из упрощённой модели, т.е. с какой-то ошибкой, а Leica использует их как полноценные углы)

    Так как вы уже внесли ЭТО недоразумение в TT (Topcon Tools), то попробуйте с этими параметрами обратный пересчёт координат. Насколько координаты будут отличаться от первоначальных?
     
    Trade_geo и zvezdochiot нравится это.
  9. Добрый день.
    Эту же точку с "чистого листа" завел в виде lat/lon WGS84, после обратного перобразования получил первоначальную координату без каких либо изменений.

    вчера разбирался, смутило вот что:
    В преобразование гельмерта есть параметры +x +y +z и +dx +dy +dz
    первые измеряются в метрах, вторые в метрах в год

    В параметрах ИГД стоит указание DX при этом единицы измерения - метры
    upload_2022-12-15_10-6-43.png

    но это все- равно не так влияет.
    Можете правильно подсказать, правильную последовательность проекций и преобразований.

    сейчас первым идет проекция по tmerc, далее преобразование из lat/lon в X,Y далее преобразование гельмерта, после переход обратно к lat/lon и перевод радиан в градусы с реверсом осей.
    это вчера полностью разбирал строчку шагов преобразования пошагово.
     
  10. stout

    stout Форумчанин

    У на какой-то испорченный телефон получается. Для того чтобы помочь вам, неплохо бы посмотреть на скриншоты ваших пошаговых действий. Начиная с того, как вы описали свою местную СК.

    Можно ссылку на описание формата? Опять таки только для того чтобы не было разночтений.
    Почти уверен, что они действительно в м/год.

    Это разные DX. Это то, что в параметрах Гельмерта обозначено как +x +y +z
     
  11. https://proj.org/operations/transformations/helmert.html

    Согласен, что телефон не очень. Ввиду ограничения конфиденциальность не могу показать описание местной системы координат.(
    транскормация из МСК в местную ИГД проходит без проблем, на выходе то же, что в Topcon.
    Преобразования между ИГД никак не идут(
     
  12. X-Y-H

    X-Y-H Администратор Команда форума

    Maloy_ksb, город то какой? кто вам мог такое посчитать
     
    zvezdochiot нравится это.
  13. Сам расчет для Мурманской области.
    Расчеты для topcon делали спеиалисты из геоприбора, если не ошибаюсь
     
  14. X-Y-H

    X-Y-H Администратор Команда форума

    Maloy_ksb, для всей области? я не верю в такие цифры
     
  15. cr2

    cr2 Форумчанин

    Кто создатель местной системы координат ? Случайно не МСК-51 ?

    Вот эти три определения эквивалентны, попробуйте с +towgs84
    Код:
    +proj=helmert +x=152.2213 +y=176.626 +z=277.808 +rx=218.862377 +ry=147.524412 +rz=658.798277 +s=-53.702102 +convention=coordinate_frame
    
    +proj=helmert +x=152.2213 +y=176.626 +z=277.808 +rx=-218.862377 +ry=-147.524412 +rz=-658.798277 +s=-53.702102 +convention=position_vector
    
    +towgs84=152.2213,176.626,277.808,-218.862377,-147.524412,-658.798277,-53.702102
    
     
  16. andrew_klikunov

    andrew_klikunov Форумчанин

    Я, конечно, всё понимаю, но такие чумачечие развороты.... Это надо сильно постараться, либо оси крутятся несколько сотен раз, как юла. У меня, при вычислении параметров, меньше 2 сек. всегда значения были, да и то не всегда.
     
    stout нравится это.
  17. stout

    stout Форумчанин

    Оффтоп

    Даже для "вздыбленного" относительно эллипсоида Красовского эллипсоида Бесселя получается чуть больше 3"
    Даже для "вздыбленного" относительно эллипсоида Красовского эллипсоида Бесселя получается
    Посмотреть вложение 186498


    --- Сообщения объединены, 16 дек 2022, Оригинальное время сообщения: 16 дек 2022 ---
    Поясните пожалуйста эту загадочную для меня фразу. Что есть "местная ИГД"?
     
    andrew_klikunov и zvezdochiot нравится это.
  18. X-Y-H

    X-Y-H Администратор Команда форума

    Оффтоп
    Даже с нормальными высотами такого разворота не будет
     
    andrew_klikunov нравится это.
  19. stout

    stout Форумчанин

    Меня и эта картинка
    [​IMG]
    сразу смутила. Но ввиду малости эффекта оставил вопрос на потом.
    --- Сообщения объединены, 16 дек 2022, Оригинальное время сообщения: 16 дек 2022 ---
    Maloy_ksb, надеюсь что вы понимаете, что вопросы к вам не с целью зачмырить, а разобраться.
     
  20. andrew_klikunov

    andrew_klikunov Форумчанин

    Там масштаб почти -54 ppm, это перебор ::blink.gif::
     
  1. Этот сайт использует файлы cookie. Продолжая пользоваться данным сайтом, Вы соглашаетесь на использование нами Ваших файлов cookie.
    Скрыть объявление
  1. Этот сайт использует файлы cookie. Продолжая пользоваться данным сайтом, Вы соглашаетесь на использование нами Ваших файлов cookie.
    Скрыть объявление