Определение параметров СК

[DDim]

Есть чертеж на съемке в автокаде, вот с такой вот сеткой:


Сетка с кружочками, это МСК61 с праметрами для Mapinfo:

Код: Выделить всё

"МСК-61 зона 2", 8, 1001, 7, 40.98333333333, 0, 1, 2300000, -4811057.628

Вторая сетка это вроде как МЖСК - то есть это ITRF2008 близкий к WGS84.
Сам чертеж сделан в координатах как раз МЖСК, чтобы нормально подгрузить в Mapinfo я копирую с базовой точкой объекты, потом вокруг неё кручу чтобы сетка стала вертикальной, и импортирую это всё в Mapinfo через dxf.

Вот возник вопрос: можно как то имея на руках такой чертеж вычислить параметры проекции МЖСК для MapInfo?

Нашел здесь на форуме ссылку на программу от Ракурса для вычисления 7 параметров, она мне вычислила вот такие параметры, но мне кажется я что то не так делаю, по недопониманию (я эколог, а не геодезист), потому что с этими параметрами не ложится ничего никуда.


На чертеже у одной координаты лишняя 7 дописана вначале, так как по таблицам с координатами точек которые у меня есть в этой системе - там её не должно быть.

Так же прикладываю два файла с координатами точек в разных проекциях.

[trir]

Таки это банльный Гаусс-Крюгер зона 7

"sk42 z7", 8, 1001, 7, 39, 0, 1, 500000, 0

который, кстати, есть в MapInfo

P.S. А что такое МЖСК?

[DDim]

Местная железнодорожная система координат в своем составе содержит пространственную местную железнодорожную систему координат (далее - ПМСК) и местную плоскую железнодорожную систему координат (далее - МЖСК).
При создании МЖСК за основу принимают международную геоцентрическую систему координат ITRF2008.

Это отсюда.

Что касается Гаусса-Крюгера, то при предложенных вами параметрах в Мапинфо я получаю вот такое:


При этом, вроде как из цитаты распоряжения РЖД следует, что там вряд ли элипсойд Крассовского, который в датуме 1001 в MapInfo. Насколько я понял, элипсойд там скорее WGS84.

[trir]

ну так замени 1001 на 104

[Донецков]

По виду координаты МЖСК похожи на UTM84 (для Ростовской обл. соответственно зона) 37N,
проверьте с ней, она по умолчанию есть во всех ГИСах ...

[DDim]

Да я уже играл в эту "угадайку", и 104, и 1001, и 1012-1014(какая то из них уже не помню, ближе всего попадает, ну относительно). Только вот не ложится ровно и всё.

Поэтому и спросил, можно ли как то имея на руках координаты в двух разных проекциях и зная что одна точно МСК61, а вторая на элипсойде вроде как WGS84 узнать параметры проекции. Я так понимаю мне нужно узнать для датума МЖСК параметры: dX, dY, dZ, EX, EY, EZ, m, НулевойМеридиан - в терминологии MapInfo.

Вот результат если 1001, на 104 заменить, как вы предложили:

[DDim]

По виду координаты МЖСК похожи на UTM84 (для Ростовской обл. соответственно зона) 37N,
проверьте с ней, она по умолчанию есть во всех ГИСах ...

похоже, согласен, но вот что то не ложится у меня, может что не так делаю

[Донецков]

Попробуйте

Код: Выделить всё

"МЖСК Зона 37 (WGS 84)", 8, 104, 7, 39, 0, 1, 500000, 0

[Boris]

Да я уже играл в эту "угадайку", и 104, и 1001, и 1012-1014(какая то из них уже не помню, ближе всего попадает, ну относительно). Только вот не ложится ровно и всё.

Поэтому и спросил, можно ли как то имея на руках координаты в двух разных проекциях и зная что одна точно МСК61, а вторая на элипсойде вроде как WGS84 узнать параметры проекции.

Если все как вы предполагаете, то можно. Заходите на сайт компании Ракурс, скачиваете программу "Вычисление 7 параметров" вводите известные координаты, в предположенных вами СК, получаете параметры трансформации между датумами. Заносите их в мапинфо, возможно углы поворота с противоположным знаком (не помню уже про знаки), и получите свою локальную СК.
Для надежности пересчитайте из МСК в любую глобальную СК на то же эллипсоиде, т.к. по моему, МСК программа не поддерживает, а вот Гаусса-Крюгера - поддерживает. И еще там надо с координатами разобраться, помниться там сперва серев шел, а потом восток.

[DDim]

Попробуйте

Код: Выделить всё

"МЖСК Зона 37 (WGS 84)", 8, 104, 7, 39, 0, 1, 500000, 0

я уже выше выложил результат этого - не ложится

Заходите на сайт компании Ракурс, скачиваете программу "Вычисление 7 параметров" вводите известные координаты, в предположенных вами СК, получаете параметры трансформации между датумами.

Я в первом сообщении выложил то что вышло в этой программе, видимо, не то что то делал, потому как и параметры вышли страшные, да и не ложится ничего - попробую повнимательнее всё проделать, может ошибся где, спасибо.

[Игорь Белов]

IMHO в соответствии с документом Росжелдорпроекта Методические указания по построению сетей для ВСМ задача ставится так.

Создаём ПМСК, определённую в ITRF2008, т. е. практически WGS 84.
На её оси вешаем эллипсоид Красовского и создаём МЖСК как проекцию Гаусса-Крюгера. В данном случае получается такой CoordSys:

Код: Выделить всё

"МЖСК (ПМСК)", 8, 999, 3, 0, 0, 0, 7, 39, 0, 1, 500000, 0

По итогам создания ПМСК можно (и нужно) определить параметры трансформирования из местной реализации Пулково 1942 в ITRF2008. Для этого используют координаты/высоты пунктов ГС в обеих СК. Решим эту задачу, взяв вместо пунктов координаты крестов с нулевыми высотоми. Но для утилиты Ракурса придётся подготовить данные. Координаты МСК-61 пересчитаем в СК-42; координаты МЖСК дополним номером зоны 7.

В настройках утилиты для обеих СК укажем тип координат "Проекция Гаусса-Крюгера (X-Y-H)", а в параметрах выберем эллипсоид Красовского и зададим осевой меридиан 39.

Настройки

img_20151225_084308.png (8.47 КБ) 12191 просмотр

Загрузим координаты и вычислим параметры.

Окно программы

img_20151225_084347.png (29.72 КБ) 12191 просмотр

Вычислим масштабный параметр как (Масштаб,единица − 1)×10⁶. Строка CoordSys для МСК-61 готова:

Код: Выделить всё

"МСК-61 зона 2 (Пулково 1942 -> ПМСК)", 8, 9999, 3, 348.5415, 297.4061, 531.6406, 1.75, 2.012, 7.66, -110.255, 0, 7, 40.9833333333333, 0.1, 1, 2300000, -4800000

Все данные МСК должны быть определены в нём. И всем местным данным в СК-42 и СК-63 также должны быть приписаны эти параметры датума. НО!

ВНИМАНИЕ! Несмотря на красивый результат, параметры трансформирования далеки от совершенства. Во-первых, у точек нет высот. Во-вторых, точки занимают крошечную область 50×50 м, что приводит к неустойчивому решению (т. е. если добавить новую точку, параметры существенно изменятся).

[Ariki]

На geodesist.ru есть некоторая информация о том, как эти железнодорожные системы строятся. Там эллипсоид Красовского, система привязана к ITRF2008 трёхпараметрическим преобразованием. Я попробовал подобрать (при условии использования параметров ГОСТ для МСК-61), получилось вот что:

Код: Выделить всё

+proj=tmerc +ellps=krass +lon_0=39 +lat_0=0 +x_0=500000 +y_0=0 +k_0=1 +towgs84=-27.71,-77.36,-83.48 

Это параметры для PROJ.4, вариант для MapInfo приводить не рискну, так как не уверен в знаках параметров, давно с ним не работал.

P. S. Ой, ErnieBoyd меня опередил

[DDim]

ErnieBoyd, Ariki, спасибо огромное за пояснения, попробую поиграться с количеством точек и взять более удалённые, сюда выложил просто маленький кусок, меня точность +-1 метр, да даже и 5 метров вполне устроит

[ Сообщение с мобильного устройства ]

[DDim]

На geodesist.ru есть некоторая информация о том, как эти железнодорожные системы строятся. Там эллипсоид Красовского, система привязана к ITRF2008 трёхпараметрическим преобразованием. Я попробовал подобрать (при условии использования параметров ГОСТ для МСК-61), получилось вот что:

а не подскажете, как вы подбирали параметры (какая то программа или что то ещё), я вот пробовал поискать как вычислять эти параметры для трехпараметрического преобразования и не нашел ничего, только вот программу которая вычисляет 7.

[Ariki]

У меня есть более-менее универсальный скрипт, который это делает:
https://github.com/Ariki/guessproj
Хотя три параметра можно легко подобрать вручную, если пересчитать исходные точки в геоцентрические декартовы координаты.

Поскольку в мое шпаргалке с параметрами проекций была большая путаница со знаками гостовских параметров, попробовал подобрать параметры МЖСК заново. Получилось вот что:

Код: Выделить всё

+proj=tmerc +ellps=krass +lat_0=0 +lon_0=39 +x_0=500000 +y_0=0 +k_0=1 +towgs84=-31.994,-53.130,-94.217 +no_defs

Если есть возможность, проверьте: без обратной связи я легко могу допустить ошибку и не заметить.