Добро пожаловать!

Войдите или зарегистрируйтесь сейчас!

Войти

Помогите перевести из МГГТ в МСК50

Тема в разделе "Исходные данные", создана пользователем v320, 30 окт 2013.

  1. v320

    Регистрация:
    30 окт 2013
    Сообщения:
    15
    Симпатии:
    0
    Здравствуйте! В первый раз столкнулся с пересчетом координат. Помогите пожалуйста перевести из МГГГТ в МСК50 два участка!
    Химкинсикй район. Так же интересует процесс пересчета. Хотелось бы самому научится,чтобы в последствии не дергать людей. Заранее спасибо!

    Х У
    1,32805.783,-14746.352
    2,32822.982,-14769.748
    3,32783.443,-14793.657
    4,32779.375,-14782.325
    5,32811.170,-14790.415
    6,32791.409,-14769.006
    7,32843.839,-14769.380
    8,32792.151,-14752.962
    9,32793.851,-14789.547
    10,32826.392,-14777.401
    11,32825.780,-14775.853
    12,32817.332,-14762.340

    Х У
    1,29234.557,-855.953
    2,29202.286,-900.448
    3,29153.809,-862.631
    4,29217.986,-844.069
     
    #1
  2. В.Шуфотинский

    В.Шуфотинский Модератор
    Команда форума Форумчанин

    Регистрация:
    10 дек 2008
    Сообщения:
    17.358
    Симпатии:
    4.976
    Не надо в трёх темах размещать одинаковые сообщения.
     
    #2
  3. texkarta

    Форумчанин

    Регистрация:
    10 ноя 2010
    Сообщения:
    94
    Симпатии:
    10

    Пожалуйста!
    1 492143.1212173843.570
    2 492160.6522173820.425
    3 492121.4552173795.953
    4 492117.2262173807.225
    5 492149.1352173799.588
    6 492129.0702173820.716
    7 492181.5042173821.081
    8 492129.5842173836.766
    9 492131.8042173800.211
    10 492164.1712173812.814
    11 492163.5372173814.355
    12 492154.8972173827.744
    1 488374.5762187682.604
    2 488342.9392187637.657
    3 488293.9282187674.776
    4 488357.8372187694.257
     
    #3
  4. v320

    Регистрация:
    30 окт 2013
    Сообщения:
    15
    Симпатии:
    0
    Спасибо!
     
    #4
  5. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная
    Уважаемый texkarta, не могли бы вы пояснить, как получили эти результаты?
    Дело вот в чём. На основе документа Положение о пространственной местной системе координат города Москвы (ПМСК Москвы) написал маленькую утилитку. Но результаты расходятся с вашими на 2-4 см.
    Данных в открытой части этого документа вполне достаточно для преобразования из МГГГТ в МСК50 и обратно. Заявлена точность на уровне 0.006×√2.
    Вот и хочу понять, в чём причина.
    v320.png
    З.Ы. По хорошему, тут надо бы и высоту учитывать, если речь идёт о миллиметрах.
    Неучёт высоты 200 м даёт ошибку 3-4 мм в плане.
     
    #5
  6. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    21.851
    Симпатии:
    7.105
    Адрес:
    Россия
    stout, а что выдают другие программы?
     
    #6
  7. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная

    fromPinnacle_v320.jpg
     
    #7
  8. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    21.851
    Симпатии:
    7.105
    Адрес:
    Россия
    Параметры эллипсоида у вас разные могут быть?
     
    #8
  9. texkarta

    Форумчанин

    Регистрация:
    10 ноя 2010
    Сообщения:
    94
    Симпатии:
    10

    Используем программы перевода,разработанные в одном из предприятий ГУГК
    Переводил МГГТ в 63 год,а потом в мск50.
    По поводу учета высоты,а в каких формулах учитывается высота?
    По контролю перевода:
    Космические корабли запускать не собираемся,а для всех выполняемых работ точности достаточно,не известно сколько ошибок
    можно притащить в измерениях.
     
    #9
  10. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная
    Специально задавал в Pinnacle с таким же сжатием, что и в документе. В Blue Marble GC третий результат. ::biggrin24.gif::
    Очевидно, что где-то путаю знаки. Это самое поганое, когда результат не сходится на несколько миллиметров.::sad24.gif::
    Pic_01.png
    Точка P (h≠0) и её проекция (по нормали) на первом эллипсоиде P'(h=0) дают разные проекции на втором эллипсоиде.
     
    #10
    Geoshaman и Lex K-G нравится это.
  11. X-Y-H

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

    Регистрация:
    18 май 2007
    Сообщения:
    21.851
    Симпатии:
    7.105
    Адрес:
    Россия
    stout, параметры эллипсоида Бесселя и МГГТ можно в личку? Послезавтра посмотрю.
     
    #11
  12. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная
    Так все величины взял из документа о котором написал в #5 посте. В том то и прикол, что для перехода из системы «МСК Москвы» в систему МСК-50 ничего, кроме открытой части документа, не надо. В координатных калькуляторах перевычисления осуществляются через WGS 84. Можно отождествить ПМСК с WGS 84 и использовать параметры перехода из системы «МСК Москвы» в ПМСК и из ПМСК в систему МСК-50. Т.е. ПМСК выступает в качестве промежуточной с.к., координаты в которой определять не надо.
    Единственное, что может и должно смущать, так это кривость МСК-50. Т.е. локально там может быть всё хорошо, но вот во всей зоне ответственности ПМСК — это вряд ли.
    Не верю я пункту 7.2
    "7.2 Параметры в таблице 1 определены по эмпирическим данным. Среднеквадратическая ошибка преобразования плановых координат между«ПМСК Москвы» и МСК-50 составляет 0.006 м. "

    Случайно не знаете, кто писал эти программы?
     
    #12
  13. Дядя Вова

    Форумчанин

    Регистрация:
    29 дек 2009
    Сообщения:
    3.758
    Симпатии:
    4.167
    Подобное эссе могли забабахать исключительно кнопочники.
    20-30 мм при пересчёте вышеуказанных СК - суперточный результат.
     
    #13
  14. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная
    C этим трудно не согласиться. Причём в большинстве других случаев результат много хуже. А вот фразе: "Среднеквадратическая ошибка преобразования плановых координат между «ПМСК Москвы» и «МСК Москвы» на пунктах каркасной сети составляет 0.006 м." охотно верю. Году в 2005-2006 попал ко мне фрагмент предварительно уравненного каталога координат МГГТ из более чем 1300 пунктов. И была возможность получить 7 параметров.
    В Pinnacle картинка выглядела так
    7ParamMosMGGT001+Resid-datum.png
     
    #14
  15. Дядя Вова

    Форумчанин

    Регистрация:
    29 дек 2009
    Сообщения:
    3.758
    Симпатии:
    4.167
    Очень возможно, что координаты пунктов МСК Москвы переуравнивали, отсюда такой достойный результат при получении параметров. А вот МСК 50 (он же модифицированный СК-42) никто не переуравнивал 100%, соответственно ни о каких 6 мм речи не может быть. (к сожалению)
     
    #15
  16. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная
    Перемножив и сложив матрицы получил
    Table.jpg
    Нашёл ошибки и у себя в утилитке и в тех параметрах, которые ввёл в Pinnacle.::biggrin24.gif::
     
    #16
    ITT, Grandpa, genfild и 2 другим нравится это.
  17. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    830
    Симпатии:
    689
    Адрес:
    Любимый Иркутск - середина земли.
    Извините дядьку, я не понял! О каком перемножении, приводящем к изменению параметров идет речь? У вас для величин из 3-й и 4-й колонок, являющихся производными от исходных величин из1-й и 2-й колонки отличия наблюдаются только у двух первых. У других 5-ти изменений нет. Как то странно. По моим банальным вычислениям, для первых двух параметров третьей колонки они должны быть 434,905м и 140,706м. Ну и для 4-й колонки также, только со знаком минус.

    И еще вопрос. А некоторой эмпирической оценкой параметров связи между ПМСК Москвы и WGS-84 вы обладаете? У меня такие оценки есть и я бы хотел сравнить их с выводами других специалистов без "нарушений режима".
     
    #17
  18. В.Шуфотинский

    В.Шуфотинский Модератор
    Команда форума Форумчанин

    Регистрация:
    10 дек 2008
    Сообщения:
    17.358
    Симпатии:
    4.976
    Без нарушения коммерческой тайны? Или о каком "режиме" Вы говорите?
     
    #18
  19. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.938
    Адрес:
    Златоглавая и Белокаменная
    В векторном виде первое равенство B=D+s×R×A => A=(1/s)×inv(R)[B-D]=-(1/s)×inv(R)×D+(1/s)×inv(R)×B;
    inv(R)=trans(R) — обратная матрица равна транспонированной.
    Если бы элементы вектора столбца D были малы и малы были бы внедиагональные элементы матрицы вращения R, тогда этим произведением ( -(1/s)×[trans(R)-E]×D; Е-единичная матрица ) можно было пренебречь. Но посмотрите, чему равно Δz!
    В конце концов, практика есть критерий истины, как учили классики. Попробуйте перевести туда-сюда с одним набором и с другим.
     

    Вложения:

    #19
  20. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    830
    Симпатии:
    689
    Адрес:
    Любимый Иркутск - середина земли.
    Посмотрел формулы. Увидел, что уважаемый stout вывел корректирующие выражения, снимающие недостатки преобразования подобия, возникающие из-за нарушения условия малости углов и/или переносов начала координат. В принципе, при величинах переносов более 500м(5Е+6мм) и даже при малых углах(2-6 сек.дуги – 1Е-5 - 3Е-5) преобразование подобия должно реализовываться в полной мере, т.е. без замены синусов малых углов на сами углы, а косинусов малых углов на единички. При больших переносах для достижения мм точности преобразования координат необходимо знать диагональные и внедиагональные элементы матрицы преобразования до 7-го -8-го знаков, а это означает, что необходимо строго реализовать матрицу R(Ex,Ey,Ez)=R1(Ex)*R2(Ey)*R3(Ez) и ей обратную для прямого и обратного преобразований координат. Здесь R1, R2 и R3 матрицы вращения вокруг соответствующих осей координат на углы, стоящие в скобках. При этом сами параметры преобразования не изменяются, усложняется вид матрицы вращения – она перестает быть столь утрированной.
    --- Сообщения объединены, 10 ноя 2013, Оригинальное время сообщения: 10 ноя 2013 ---
    Хочу добавить следующее.


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

    Ребра этого графа имеют направления и номера, что позволяет вычислять параметры при нужном направлении перехода. В частности, ребра 1, 2 и 3 отражают информацию из соответствующих колонок таблицы из поста #16. Пользоваться этим графом можно подобно схеме нивелирной сети – вычислять значения параметров при задании фиксированных величин и даже уравнивать значения при сочетании фиксированных и предварительно определенных параметров.
    Например, наличие в графе МСК-50, худо-бедно связанной с ГСК-42, у которой есть гостовские параметры связи с WGS-84, позволяет сделать банальные вычисления, позволяющие найти параметры связи 4,5 и 6.
    По моему графу получается следующее.

    1 2 3 4 5 6
    Dx 118,754 -316,151 -434,905 411,335 -23,57 292,581
    Dy 61,782 -78,924 -140,706 281,656 140,95 219,874
    Dz 93,237 -589,65 -682,887 762,687 79,8 669,45
    Ex 2,40896 1,57273 -0,83623 0,83623 0 -1,57273
    Ey 3,47502 2,69209 -0,78293 1,13293 0,35 -2,34209
    Ez -1,29688 2,34693 3,64381 -2,85381 0,79 -1,55693
    dm -6,5177 -8,4507 -1,933 2,153 0,22 8,6707


    Здесь параметры из колонок 1, 2(из Руководства по ПМСК Москвы) и 5(ГОСТ Р 51794-2008) взяты в качестве исходных, а другие вычислены с помощью логики графа.
    Интересно отметить, что данные верхней части графа являются открытыми, а нижней части ( в частности значения ребра 6)- закрытыми, входят в секретную часть документа, устанавливающего связи ПМСК Москвы. И в этом мой ответ на вопрос В.Шуфотинского.
     
    #20
    ITT, Grandpa и genfild нравится это.

Поделиться этой страницей

  1. Этот сайт использует файлы cookie. Продолжая пользоваться данным сайтом, Вы соглашаетесь на использование нами Ваших файлов cookie.
    Скрыть объявление