Система координат города Краснодар

lil_saviour, В свое время с помощью скрипта удалось сконвертировать файл mapinfow.prj (на тот момент рабочей ГИС системой была именно mapinfo) в систему proj4 и занести все СК в кугис. А поскольку в mapinfow.prj скопилось проекций весьма много, штук 300, то впоследствии выяснилось, что кугис нереально тормозит при открытии меню "Пользовательские проекции". То бишь секунд 15-20 уходило на открытие окна. Поэтому было принято решение добавлять проекции в формат proj4 по мере необходимости

Если же говорить об отдельных файлах проекций .prj, которые чаще всего идут вместе с шейпами, дабы проекция была подвачена кугисом/аркгисом, то вопроса массовой конвертации из таких prj в proj4 не возникало. Не было необходимости

 

Принял. Я разобрался, где была ошибка. Пришлось вникнуть в синтаксис proj4.
В параметрах prj, что мне предложили (8, 1001, 7, 39, 0, 1, 21800.013, -4974499.820) указан стандартный датум с кодом 1001 - Pulkovo 1942. По умолчанию, этому датуму соответсвует эллипсоид Красовского и следующие семь переходных параметров от WGS84: 24, -123, -94, -0,02, 0.25, 0.13, 1,1 (dX, dY, dZ, EX, EY, EZ, m). При таком переводе точки уходили мимо. Но, когда Вы мне прислали ее же, но в proj4 , то в ее синтаксисе был фактически описан другой датум, тоже на основе Красовского, но семь параметров другие: 23.57, -140.95, -79.8, 0, 0.35, 0.79, -0.22. В этих параметрах, все работает прекрасно. Единственное, что для получения такого proj4 из prj, согласно синтаксисам из англоязычных мануалов, prj должен был выглядеть так: 8, 9999, 3, 23.57, -140.95, -79.8, 0, 0.35, 0.79, -0.22, 7, 39, 0, 1, 21800.013, -4974499.820. То есть не стандартный датум, а пользовательский, с параметрами, которые были отображены в proj4. Как то так, в общем.
Проблема, в общем, решена. Мне только теперь интересно понять, как будучи не геодезистом в будущем адекватно вычислияь эти 7 параметров (dX, dY, dZ, EX, EY, EZ, m) для неизвестных мне проекций...желательно из массива точек с известными коориднатами в двух системах. Но это, видиом, как-то потом...

 

* helmert3d
* GeoFindKey
* Таблицы MS Excel для преобразования координат

 

lil_saviour,

Да, это параметры датума 1001 из Mapinfo. Разумеется он некорректный и при пересчёте в проекции, использующие другой датум, либо в системы координат, которые не являются проекциями (тот же EPSG:4326, к примеру), естественно не подходит

Параметры ГОСТа. Сначала ГОСТа 51794-2008, затем ГОСТ 32453-2017. Свежее, вроде бы, нету

Правильнее будет 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, 7, 39, 0, 1, 21800.013, -4974499.820

Обратите внимание на знаки перед значениями поворота в датуме.

Mapinfo, Arcgis, Autocad используют преобразование Бурша-Вольфа и знаки перед значением поворота там такие же, как и в ГОСТе.
QGIS, Global Mapper используют преобразование Гельмерта, поэтому там обратные знаки.

В 97% Вам не придётся их вычислять. Параметры ГОСТа и были разработаны, чтобы в камеральных условиях производить вычисления с помощью калькуляторов. Всё, что Вам необходимо будет в случае чего вычислить - это осевой меридиан и оба смещения
2% я бы оставил на проекции, где используется датум СК-95. В основном это объекты Транснефти, например.
1% - это кастомные датумы. Которые используются, например, на Сахалине

Ну и да, геодезисты на местности занимаются локализацией и устранением погрешностей, поэтому на местности тоже могут быть вычислены свои местные, локальные датумы

 

Ну, prj в qgis и не используется, поэтому вроде как нормально все. А при переводе proj4 все съедается корректно с теми же параметрами.
А вот геодезиста для локализации тут у меня не завелось, к сожалению. Может, сами освоим чуть позже. Так еще и геодезисты, сидящие на заводах, как правило, ни шиша не дают тех параметров, в которых держат свои генпланы (у кого секретность, а кто просто не знает). И приходится все переводить с помощью excel и такой-то матери.

Спасибо за помощь!

 

Вставлю 5 копеек про ArcGIS. Там можно создать преобразование и по Гельмерту (Coordinate Frame) и Бурша-Вольфа (Position Vector). Там гибко сделано. Не говоря про преобразования Молоденского, сетки NTv2 и т. д. Метод пересчёта описан в отдельном файле, не помню расширение - gfc или gcf. А вот proj4, при 7-параметрическом использует именно Бурша-Вольфа.

 

А вот здесь секундочку. Коли твой ArcGIS может "и то и это", то не проверишь ли, те ли знаки я выбрал в последней версии helmert3d? Я так внятного ответа и не добился на форуме.

 

Проверить могу, не раньше следующей недели. Стоит учесть, что ArcGIS параметры не вычисляет, пользователь их туда сам вводит и выбирает метод. Это на всякий случай.

 

Ну так оно и надо: Взять набор координат, взять строку proj4 и получить трансформированные координаты. А в helmert3d по двум наборам координат надо получить исходную строку proj4.