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

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

Ситуашия следующая:
Имею точки в местной системе координат. С помощью 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 метров, во втором больше и явно не правильные.

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

 

Что по вашему мнению есть местный датум, если в него можно конвертировать координаты?

Разница между чем и чем? Что в качестве "истины"?

А в контроллере откуда параметры взялись?

Что такое rx, ry, rz и ы, В каких единицах?

 

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

 

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

 

Догадываюсь, но хотелось услышать от топистартера об источнике. Коме того, меня зело смущает сама цифирь.

 

Всем добрый день.
Спасибо за вопросы, попробую ответить на них, если получится конечно.
1) местный датум (или местная ИГД), я могу ошибаться в терминах, в топкон выглядит это так.

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

 

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

 

Эти специалисты достойны быть если не академиками Азербайджанской академии наук, то членами-корреспондентами уж точно. (Или, мягко выражаясь, просто членами?) Мой стёб основан на одной хохме. Когда-то давным-давно в базе EPSG появились 7 параметров, для преобразования из СК-42 в WGS 84, которые были определены Азербайджанской академией наук. Там углы вращения составляли тысячи секунд дуги чего не может быть в принципе. В худшем случае угол меньше 6" (чел стоит на Гринвичском меридиане)

Именно малость углов позволяет заменять синус угла самим углом, а косинус – единицей в полноценной матрице вращения Гельмерта (т.е. использовать преобразование Бурши-Вольфа). И только ПО от ляйки всегда работало с полноценной матрицей Гельмерта (Хотя с этим возникает свой затык, все определяют углы из упрощённой модели, т.е. с какой-то ошибкой, а Leica использует их как полноценные углы)

Так как вы уже внесли ЭТО недоразумение в TT (Topcon Tools), то попробуйте с этими параметрами обратный пересчёт координат. Насколько координаты будут отличаться от первоначальных?

 

Добрый день.
Эту же точку с "чистого листа" завел в виде lat/lon WGS84, после обратного перобразования получил первоначальную координату без каких либо изменений.

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

В параметрах ИГД стоит указание DX при этом единицы измерения - метры


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

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

 

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

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

Это разные DX. Это то, что в параметрах Гельмерта обозначено как +x +y +z

 

https://proj.org/operations/transformations/helmert.html

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

 

Сам расчет для Мурманской области.
Расчеты для topcon делали спеиалисты из геоприбора, если не ошибаюсь

 

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

Вот эти три определения эквивалентны, попробуйте с +towgs84

Code:

+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

 

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

 

Spoiler

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

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

--- Сообщения объединены, Dec 16, 2022, Оригинальное время сообщения: Dec 16, 2022 ---

Поясните пожалуйста эту загадочную для меня фразу. Что есть "местная ИГД"?

 

Меня и эта картинка

сразу смутила. Но ввиду малости эффекта оставил вопрос на потом.

--- Сообщения объединены, Dec 16, 2022, Оригинальное время сообщения: Dec 16, 2022 ---

Maloy_ksb, надеюсь что вы понимаете, что вопросы к вам не с целью зачмырить, а разобраться.

 

Там масштаб почти -54 ppm, это перебор