Местная система координат линейного объекта

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

Написал статью, приглашаю к конструктивному обсуждению. Спасибо Чаровнице за идею.

Статья опубликована.

[vveco]

Очень актуальная статья. Спасибо!
Не могу понять в чем причина, но не получается воспроизвести тестовый пример..
Файл данных inv.txt:
51N 22E 50N 20E
команда
geod -I -f "%.10f" -F "%f" +ellps=WGS84 +units=m inv.txt
результат:
-127.3190807751 51.1375476134 180292.395160

Ну и с остальными командами также.

Использую этот пакет:
proj446_win32_bin.zip: Prebuilt Win32 executables, DLL including NAD27 grid shift files

Что не так?

[Александр Мурый]

Дело в версии PROJ.4 - у вас очень старая версия ("proj446"), там действительно почему-то считается не как в статье. Какие-то другие параметры внутри, наверное.
Вот ссылка на бинарники актуальной версии 8.0. Проверил - тамошний geod.exe считает как в статье.

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

Вот ссылка на бинарники актуальной версии 4.8.0

Это первое.

Во-вторых, если установлен QGIS, утилиты PROJ.4 находятся вместе с прочими в директории bin. Можно добавить путь к ним в переменную PATH, но гораздо авантажнее воспользоваться ярлычком "MSYS Shell" на рабочем столе. Тогда Вы получаете ещё и отличный эмулятор UNIX-терминала со всеми вкусными инструментами.

В-третьих, разработчики PROJ.4 выпустили версию 9, в которую, как и обещалось, интегрированы библиотеки GeographicLib. Теперь утилита geod вычисляет геодезические задачи на эллипсоиде. По теме статьи это делает достаточно точным первоначальное определение параметра α, и итерации в большинстве случаев становятся ненужными. В соответствии с этим я подправлю статью, когда сборки от OSGeo и NextGIS будут включать новую версию. А пока для пользователей Windows я скомпилировал из исходников версии 4.9.1 бинарные файлы.

[Александр Мурый]

Я на всякий случай залил архив с бинарниками на сервер, чтобы не потерялось.

[Charaunica]

Столкнулась со следующими проблемами при использовании созданной системы линейного объекта. В чем суть:

Пересохранила нужные слои чертежа объекта (в МСК города) из dwg в dxf, импортировала в QGIS (в созданную по методике статьи СК обьекта). Чертеж-10 развернуло как нужно, но координаты улетели. Самое забавное, что при тестировании проекции в QGIS она не села буквально на 15-20 см по сравнению с исходными данными.

Опять же, не могу утверждать, что все сделала правильно, поэтому готова исходниками поделиться.

МСК города
+proj=tmerc +lat_0=0 +lon_0=27.5705643917 +k_0=1 +x_0=911.2154 +y_0=-5975255.9586 +towgs84=-8.6433,-136.8101,-59.7188,3.353,3.251,5.904,0.629 +ellps=krass

МСК объекта
+proj=omerc +lat_0=53.9200334055 +lonc=27.712314778 +alpha=17.8625874234 +gamma=90 +k_0=0.9999963+x_0=0 +y_0=0 +ellps=WGS84

Протяженность объекта от 0 координат - 2168,22

Приблизительные координаты торцов, полученные после перепроецирования: -108,3 -2116,2 ; -104,4 53,2

[Донецков]

Если не секрет, то в чем цель была установления МСК линейного объекта? не проще ли было совместить исходные данные с МСК города (т.к. это очень локальный участок)...

[Charaunica]

Если не секрет, то в чем цель была установления МСК линейного объекта? не проще ли было совместить исходные данные с МСК города (т.к. это очень локальный участок)...

Абсолютно не секрет: мне нужна МСК аэродрома, она активно используется в Беларуси, особенно в отчетности по аэронавигационным препятствиям. Пишу магистерскую по ГИС-обеспечению аэродрома, как одна из задач получить данную МСК и анализировать препятствия уже в ГИС.

[Донецков]

О, перепутал, что длина объекта не 2168.22 км, а 2168.22 м ...
МСК города у Вас "каталожное" или сами вычисляли? , т.к. +towgs84=-8.6433,-136.8101,-59.7188,3.353,3.251,5.904,0.629 - явно "доморощенное"
в МСК объекта +ellps=WGS84 - это необходимость?

[Charaunica]

МСК города у Вас "каталожное" или сами вычисляли? , т.к. +towgs84=-8.6433,-136.8101,-59.7188,3.353,3.251,5.904,0.629 - явно "доморощенное"
в МСК объекта +ellps=WGS84 - это необходимость?

Находила через ПО Ракурс "Вычисление 7 параметров" (по известным координатам 5 пунктов в WGS84, СК-42 и МСК. Насчет +ellps=WGS84 - это же инструмент proj4, без него никак,

А вообще все делалось по рецептам с http://gis-lab.info/ по созданию разных МСК, вот и вникаю во всю эту красоту)

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

координаты улетели

Ну, «улетели» сильно сказано. Так, отползли

Действительно, смещение наблюдается. Проблема, по-видимому, в том, что в СК аэродрома не прописан явно датум. В принципе, задача решается корректно, если разбить её на два этапа:

• открыть слои в городской системе и пересохранить в шейпы в WGS84;
• переоткрыть эти шейпы и сохранить их копии в СК аэродрома.

Однако правильное решение — заменить +ellps=WGS84 на +datum=WGS84.

Кроме того, "+no deft" следует исправить на "+no_defs". И ещё я бы подправил масштабный коэффициент под расстояние 2168.22 м.

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

+proj=omerc +lat_0=53.9200334055 +lonc=27.712314778 +alpha=17.8625874234 +k=0.9999222461 +x_0=0 +y_0=0 +gamma=90 +datum=WGS84 +units=m +no_defs

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

Статья приведена в соответствие с сегодняшними реалиями. Система координат создаётся в формате WKT в двух вариантах косой проекции Меркатора.