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

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

Войти

Программа для пересчета координат из МСК города

Тема в разделе "Выбор и сравнение геодезического ПО", создана пользователем Usopp, 13 фев 2017.

  1. Usopp

    Регистрация:
    19 янв 2017
    Сообщения:
    3
    Симпатии:
    1
    Добрый день!
    Прошу посоветовать программы для пересчета из МСК-города в СК-42, МСК-субъекта. В наличии имеются параметры МСК, установленные на населенные пункты и представленные в виде: координат Х, Y в СК-42, значение осевого меридиана, значения высоты, значение начала отсчета местной СК. Эти параметры необходимо внести в программу и затем по ним пересчитать координаты из МСК-города в СК-42, МСК-субъекта. Желательно чтобы программа позволяла производить импорт/экспорт MIF/MID файлов (MapInfo) с координатами земельных участков. Присматриваюсь к CREDO Транскор, но есть ли в данной программе возможность внесения параметров местных СК и дальнейшего пересчета из местных СК в СК-42, МСК-субъекта? Буду благодарен за помощь!
     
    #1
    sDrawMan нравится это.
  2. trir

    Форумчанин

    Регистрация:
    25 ноя 2014
    Сообщения:
    3.252
    Симпатии:
    931
    Адрес:
    gnomtrir@mail.ru
    #2
    sDrawMan нравится это.
  3. Geo_major

    Форумчанин

    Регистрация:
    28 мар 2010
    Сообщения:
    4.404
    Симпатии:
    5.457
    Адрес:
    Военный городок Москва-74
    http://geodesist.ru/forum/resources/archaeosys.8/

    Перед началом работ внимательно прочитайте мануал к этой программе, иначе программа будет неактивна.
     
    #3
    sDrawMan и Oigi нравится это.
  4. ЮС

    Форумчанин

    Регистрация:
    28 фев 2010
    Сообщения:
    4.563
    Симпатии:
    5.052
    Есть такие возможности в Транскоре (их даже намного больше).
    Можно, например, несколькими методами определять неизвестные параметры МСК (ключи) по общим пунктам в обеих системах.
     
    #4
    sDrawMan нравится это.
  5. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    821
    Симпатии:
    684
    Адрес:
    Любимый Иркутск - середина земли.
    Возможности есть, да вот не выходит "каменный цветок". Я когда расковыривал вопрос по связи городских систем с ГСК (см.http://geodesist.ru/forum/threads/mestnye-sistemy-koordinat-gorodov.12174/) так и не смог заставить адекватно работать Транскор. Причем "неадеквашки" наблюдались как в случае определения параметров, так и в случае использования известных ключей. Поэтому буду признателен, если увижу примеры, или подсказки по этому поводу.
     
    #5
    Иоан4 и sDrawMan нравится это.
  6. ЮС

    Форумчанин

    Регистрация:
    28 фев 2010
    Сообщения:
    4.563
    Симпатии:
    5.052
    А можно глянуть на конкретный пример, где "не выходит цветок"? Интересно попробовать свежую версию Транскора.
    Я не утверждаю, что Транскор всесилен, однако зачастую его способности недооцениваются или используются не в полной мере.
     
    #6
    sDrawMan нравится это.
  7. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.932
    Адрес:
    Златоглавая и Белокаменная
    Уважаемый Вячеслав Трофимович!
    Действительно было бы интересно взглянуть на примеры.
     
    #7
    sDrawMan нравится это.
  8. прохожий

    Форумчанин

    Регистрация:
    15 сен 2008
    Сообщения:
    820
    Симпатии:
    154
    Адрес:
    CREDO
    Usopp, если параметры СК действительно известны, то пересчет можно выполнить в любом из продуктов на платформе CREDO III (кроме Вьювера).
    Импорт MIF/MID в них тоже есть, как и прямой импорт кадастровых XML.
     
    #8
    sDrawMan нравится это.
  9. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    821
    Симпатии:
    684
    Адрес:
    Любимый Иркутск - середина земли.
    Да вот, собственно, свеженький пример решения задачи определения параметров связи ГорСК и СК-42.
    Здесь не удается определить долготу осевого меридиана ГорСК. Это Сахалин, п.Макаров. Анализ показал, что причиной неуспеха является косячок от Кредо-Транскор(у меня версия 2.30.437). Он глючит с вычислениями координат в 24-й зоне.

    Посмотрите, может у вас, что получится.
     

    Вложения:

    #9
    sDrawMan нравится это.
  10. Slinger

    Форумчанин

    Регистрация:
    13 окт 2010
    Сообщения:
    1.270
    Симпатии:
    623
    ВЯЗ,

    Долгота 138˚ 40' .У меня такого плана параметры вышли. Не через транксор, правда.

    Раскрыть Спойлер
    "МСК г. Макаров", 8, 1001, 7, 138.666666667, 0, 1, -293527.679, -5385367.23


    Но ошибка будь здоров, 1.7 метра. На такой небольшой город это слегка многовато. Хотя, для некоторых целей может быть и терпимо.
     
    #10
    sDrawMan нравится это.
  11. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.932
    Адрес:
    Златоглавая и Белокаменная
    С такими параметрами проекции ТМ
    SemiMajorAxis=6378245
    RecipFlattening=298.3
    CMeridian=138.664550467965°
    scale=0.99888379195035
    FalseNorthing=-5379355.380
    FalseEasting=-293345.848
    Получились такие невязки
    Pic001.png
    Хотя сама "согласованность" наборов координат намного выше
    Pic002.png
    Видно, что координаты в МСК получены пересчётом из ГГС.
    Т.о. это указывает на то, что МСК не является проекцией ТМ в прямом смысле.
     
    #11
    sDrawMan нравится это.
  12. ЮС

    Форумчанин

    Регистрация:
    28 фев 2010
    Сообщения:
    4.563
    Симпатии:
    5.052
    Да, не надо искать чёрную кошку (ТМ) там, где её нет.
    Первое, что я всегда делаю при поиске ключа, - тупо делаю прикидку по Хельмерту:
    QIP Shot - Screen 1269.jpg
    По всей видимости так же "тупо" и была образована данная МСК.

    В Транскоре есть ещё несколько методов поиска ключа.
    Аффинное преобразование:
    QIP Shot - Screen 1270.jpg
    МНК в ортографической проекции.
    QIP Shot - Screen 1271.jpg
    В наборе есть ещё несколько методов, однако они дают худшие результаты.
     
    #12
    Иоан4 и sDrawMan нравится это.
  13. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    821
    Симпатии:
    684
    Адрес:
    Любимый Иркутск - середина земли.
    ЮС, спасибо! Я все это, разумеется, опробовал. Но у меня в принципе Транскор не давал таких результатов. Там были десятки метров в невязках. К сожалению мои картинки будут доступны только в понедельник. У вас какая версия Транскора?

    Slinger. Да, я также вне Транскора определил долготу осевого.

    Stout. В ваших вычислениях настораживает огромный масштабный коэффициент. Это можно как-то интерпретировать?
     
    #13
    sDrawMan нравится это.
  14. ЮС

    Форумчанин

    Регистрация:
    28 фев 2010
    Сообщения:
    4.563
    Симпатии:
    5.052
    Транскор 2.31.
    Можно и так, если "ключ" представить с учётом высоты поверхности относимости (расчётная Нпк=58 м).
    Отчёт в процессе расчёта ключа:
    QIP Shot - Screen 1274.jpg
    И так выглядят устанавливаемые параметры:
    QIP Shot - Screen 1275.jpg
     
    #14
    sDrawMan нравится это.
  15. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.932
    Адрес:
    Златоглавая и Белокаменная
    Так не от хорошей жизни это. Если бы весь софт кроме стандартных проекций умел ещё делать комбинации, в частности, на завершающем этапе осуществлять преобразование по Гельмерту, ещё лучше, имел бы какой-нибудь встроенный скриптовый язык или возможность вызова внешних программ, то и надобности в этом не возникло бы. Именно по этой причине я всегда был противником использования проекции Коугия.
    Так и я практически всегда с этого начинаю. Только в моём случае преобразование по Гельмерту тождественно использованию гармонического полинома первой степени, я же начинаю со второй степени.

    Если честно, то пока затрудняюсь в осмыслении результата. Собственно сам масштабный коэффициент в формулах присутствует только в произведении на большую полуось. Но это и понятно. Это как измерять длину не в метрах а в футах. Между коэффициентами гармонического полинома и параметрами проекции существует строгая взаимно однозначная связь. Я использовал формулы из Морозова, там где он определяет коэффициенты перехода из зоны в зону. Только слегка уточнил их вид в MAXIMA. Потом нашёл их же в диссертации Подшивалова.
    В программе можно принудительно изменять масштаб. Если его принять равным 1, то невязки становятся как у Slinger.
    Интересно также, что произведение большой полуоси на масштаб даёт величину ≈6371125. То есть почти средний радиус Земли.
     
    #15
  16. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    821
    Симпатии:
    684
    Адрес:
    Любимый Иркутск - середина земли.
    Это решение, по моему мнению, лучшим образом соответствует логике установления городских систем координат. В нем используется условная поверхность относимости, которая, вероятнее всего, взята равной 60м, а появление масштабного коэффициента в 9-ть единиц шестого знака обусловлено вероятностным характером определения этой величины. Долгота осевого меридиана, в принципе может быть "корявым"(не округленным числом), если начало ГорСК совмещено с пунктом ГГС, который брался в качестве физического закрепления этого начала.
    Заметил только одну странность - в автонастройках Транскора на ГорСК номер зоны выбран равным "-1".
    Оригинальное решение Stout'a в моей головушке не укладывается.
    Продолжаем разговор на затронутую тему.
    Вот вам еще один пример. Он, правда синтетический, т.е. сочетает в себе моделирование с "приглядкой" на реально существующую ГорСК. На нем я также обломал зубки Транскора.
     

    Вложения:

    #16
  17. stout

    Форумчанин

    Регистрация:
    5 янв 2008
    Сообщения:
    4.172
    Симпатии:
    11.932
    Адрес:
    Златоглавая и Белокаменная
    Интересно. Именно пункт 7с выбивается из стройного ансамбля остальных пунктов на пару см.
     
    #17
  18. ЮС

    Форумчанин

    Регистрация:
    28 фев 2010
    Сообщения:
    4.563
    Симпатии:
    5.052
    Пробую варианты расчёта ключа для нового набора пунктов.

    Есть и такой Хельмерт.
    QIP Shot - Screen 1276.jpg
    Да, пункт 7с явно выбивается.
    QIP Shot - Screen 1277.jpg
    Исключаю из набора пункт 7с и несогласованность не превышает 1 мм.
    QIP Shot - Screen 1278.jpg

    Полиноминальный тип преобразования, уклонения в пределах 0,6 мм (7с исключён).
    QIP Shot - Screen 1279.jpg

    МНК в проекции ТМ, уклонения в пределах 3,5 мм (7с исключён).
    QIP Shot - Screen 1280.jpg

    МНК в ортографической проекции, уклонения в пределах 1,3 мм (7с исключён).
    QIP Shot - Screen 1281.jpg

    Просмотр вариантов ключа в ТМ с учётом поверхности относимости.
    QIP Shot - Screen 1282.jpg
    Выбираю начальным пункт 2 и получаю ведомость установления ключа.
    QIP Shot - Screen 1283.jpg
    Так параметры этого ключа выглядят в Транскоре.
    QIP Shot - Screen 1284.jpg

    Вариант определения ключа в ТМ с учётом поверхности относимости + дополнительный доворот (часто применяется для стройсеток или небольших нас. пунктов). Даю сразу ведомость.
    QIP Shot - Screen 1285.jpg
    И так эти параметры выглядит в Транскоре.
    QIP Shot - Screen 1286.jpg , в доп. параметрах МСК добавлен уголок доворота ∝.

    В принципе, любой из полученных ключей можно применять для пересчёта координат.
     
    #18
    Андрей Мороз и stout нравится это.
  19. ВЯЗ

    Форумчанин

    Регистрация:
    17 май 2012
    Сообщения:
    821
    Симпатии:
    684
    Адрес:
    Любимый Иркутск - середина земли.
    Мужики, с Днем Защитника Отечества, Вас! Всего наилучшего!

    Появилось, наконец-то, время, чтобы закончить свою "предъяву" в адрес Кредо-Транскор.
    Весьма признателен ЮСу за вычисления, которые замечательно дополняют мои заключения.
    Суть моей претензии, высказанной в сообщении №5, заключается в том, что Транскор не позволяет строго реализовать вычисления координат в рамках перехода между ГСК и МСКгорода.
    Действительно, если вспомнить вопрос "топикстартера" и мой второй пример, то получается, что имеющихся у меня ключей перехода(Lo,Xo,Yo,xo,yo,Ho) не достаточно для решения задачи средствами Транскора. Их просто не удастся внедрить в качестве исходных данных, чтобы начать ими пользоваться. Сравните мои ключи, приведенные во втором примере и заключительную табличку уважаемого ЮС в сообщении №18. Разве это одно и тоже? Получается, что для введения параметров МСКгорода в Транскор нужно сначала наиграться с опциями этого ПО по вероятностному решению задачи определения ключей перехода, найти какой-то подходящий вариант и уже потом заняться решением основной задачи. Для описания этой интеллектуальной процедуры мне нравится определение - доставать правое ухо левой рукой из под правой коленки. Ну не очень удобно.
    Далее. Мой второй пример содержит данные (за исключением пункта 7с, в котором я ошибся при наборе числа) по формулам, которые систематизировал из имеющихся источников, анонсировал в том же сообщении №5, и протестировал их как мог. В моих тестах три алгоритма вычислений давали мм сходимость между собой, а также при прямых и обратных вычислениях. Но повторить эти тесты с Кредо-Транскор я не смог. Подобно вычислениям ЮС я получал разные параметры перехода, не совпадающие с известными ключами. В итоге я пришел к заключению, что в Транскоре используются свои алгоритмы. В техническом описании я увидел просто ошибочные формулы...Видимо, описание и реальные алгоритмы не одно и то же.
    Теперь последнее. Оно касается вероятностных(определенных по группе пунктов с координатами в паре СК) и детерминированных(заложенных при образовании МСК) ключей перехода. Вероятностные ключи можно использовать только для приближенной настройки СК. Например, при создании опорной геодезической сети, или ПВО, которые будут уравниваться с фиксацией координат пунктов в МСК. Их также можно использовать для настройки работы в режиме РТК, но в ограниченной области и с обязательными проверками на известных пунктах.
    Для решения ответственных задач массового преобразования координат пунктов из ГСК в МСКгорода и обратно, подобно тому, как они поставлены топикстартером, необходимо пользоваться исключительно детерминированными ключами. И Кредо-Транскор здесь пока не помощник.

    Всем успехов!
     
    #19
    voron73 нравится это.
  20. ЮС

    Форумчанин

    Регистрация:
    28 фев 2010
    Сообщения:
    4.563
    Симпатии:
    5.052
    Всё позволяет, надо лишь задавать корректные исходные данные.
    Что мы имеем?
    QIP Shot - Screen 1302.jpg Предположим, что ...920,230 - грубая ошибка, но и ...920,023 вылетает из общего списка около 20 мм. Ладно бы только это, однако и другие "Известные параметры ключа" не совсем соответствуют тому набору координат, что представлены в двух СК.
    1. Например, если пересчитать из ГСК в ТМ с осевым меридианом 107°33'42.75" (известные данные), то наблюдается общий разворот полученной системы относительно МСКГор на 1.96". Вроде мелочь, однако при расстояниях между пунктами до 15 км (5-6) возникает ошибка смещения (разворота) около 14 мм. Для исключения разворота у МСКГор должен быть осевой меридиан 107°33'40.256264".
    1. Средний радиус кривизны элл. Красовского для пункта 7с = 6383267.503 м. Для пересчёта на поверхность относимости масштабный коэффициент должен быть = 1.00008146. Однако по факту (по исходным данным) средний масштаб между ТМ и МСКГор = 1.00008118, то есть разница в 0.3 мм/км, что может показаться несущественным, но между пунктами 5 и 6 это даёт ошибку взаимного положения порядка 4-5 мм.
    Так что представленные "известные параметры"
    в принципе не могли давать миллиметровую сходимость.

    Теперь о Транскоре.
    Да, когда параметры МСК неизвестны, приходится немного "поиграть" для выбора наилучшего варианта. Кстати, это занимает времени меньше, чем, например, вручную вбить координаты 13 пунктов в двух системах.
    Если же параметры известны, тогда ещё проще. Только не надо путать Xo и Yo, которые даны как исходные в 18 зоне, с теми Xo(ТМ) и Yo(ТМ), что в карточке настройки МСК.
    Xo(ТМ) и Yo(ТМ) в карточке даны для зоны МСК (ос.мер.107°33'42.75"). Получить их просто, надо лишь пересчитать Xo и Yo из 18-й зоны в зону МСК.
    Далее заполняется карточка МСК:
    QIP Shot - Screen 1300.jpg и всё, можно работать. Всё строго по "Известным параметрам ключа".

    Однако, как уже говорилось, из-за не совсем корректных исходных данных, пересчёт будет выполняется с определёнными ошибками. Сравните, что было задано в условии задачи (слева) и что получено в результате строгого пересчёта, используя так называемые "известные параметры".
    QIP Shot - Screen 1301.jpg Координаты начала (7с) пересчитаны верно, а всё остальное имеет существенные отклонения из-за ошибки разворота (1.96"), неточного масштаба и "отскока" пункта 7с от общей массы.

    Транскор же, опираясь на всё 13 пунктов и вычисляя свой ключ, даёт значительно лучшее решение.
     
    #20

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

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