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

Обсуждение материалов сайта: вопросы, замечания, предложения
Ответить
Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2229
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1501
Откуда: Казань

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

Сообщение Игорь Белов » 21 ноя 2014, 20:25

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

Статья опубликована.
The purpose of computing is insight, not numbers

Аватара пользователя
vveco
Активный участник
Сообщения: 195
Зарегистрирован: 20 май 2009, 06:40
Репутация: 3
Откуда: Саратов
Контактная информация:

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

Сообщение vveco » 21 мар 2015, 18:53

Очень актуальная статья. Спасибо!
Не могу понять в чем причина, но не получается воспроизвести тестовый пример..
Файл данных 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

Что не так?

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 21 мар 2015, 19:43

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

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2229
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1501
Откуда: Казань

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

Сообщение Игорь Белов » 22 мар 2015, 15:58

Александр Мурый писал(а):Вот ссылка на бинарники актуальной версии 4.8.0
Это первое.

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

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

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 22 мар 2015, 21:49

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

Charaunica
Интересующийся
Сообщения: 35
Зарегистрирован: 11 мар 2013, 11:03
Репутация: 1

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

Сообщение Charaunica » 08 май 2015, 12:23

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

Пересохранила нужные слои чертежа объекта (в МСК города) из 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
Вложения
dxf.rar
dxf , который перепроецировался в МСК объекта
(586.84 КБ) 786 скачиваний
Проверка_QGIS.jpg
Проверка проекции созданной проекции в QGIS
Проверка_QGIS.jpg (67.18 КБ) 22299 просмотров
Проверка_shell.jpg
Процесс работы утилиты
Проверка_shell.jpg (86.77 КБ) 22299 просмотров

Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

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

Сообщение Донецков » 08 май 2015, 13:16

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

Charaunica
Интересующийся
Сообщения: 35
Зарегистрирован: 11 мар 2013, 11:03
Репутация: 1

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

Сообщение Charaunica » 08 май 2015, 13:34

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

Донецков
Гуру
Сообщения: 3058
Зарегистрирован: 19 май 2010, 19:44
Репутация: 189

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

Сообщение Донецков » 08 май 2015, 14:22

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

Charaunica
Интересующийся
Сообщения: 35
Зарегистрирован: 11 мар 2013, 11:03
Репутация: 1

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

Сообщение Charaunica » 08 май 2015, 14:52

Донецков писал(а): МСК города у Вас "каталожное" или сами вычисляли? , т.к. +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/ по созданию разных МСК, вот и вникаю во всю эту красоту)

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2229
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1501
Откуда: Казань

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

Сообщение Игорь Белов » 08 май 2015, 19:40

Charaunica писал(а):координаты улетели
Ну, «улетели» сильно сказано. Так, отползли :)

Действительно, смещение наблюдается. Проблема, по-видимому, в том, что в СК аэродрома не прописан явно датум. В принципе, задача решается корректно, если разбить её на два этапа:
  • открыть слои в городской системе и пересохранить в шейпы в 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
Вложения
airdrome.zip
(17.9 КБ) 854 скачивания
The purpose of computing is insight, not numbers

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 2229
Зарегистрирован: 04 янв 2011, 22:00
Репутация: 1501
Откуда: Казань

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

Сообщение Игорь Белов » 03 май 2020, 15:25

Статья приведена в соответствие с сегодняшними реалиями. Система координат создаётся в формате WKT в двух вариантах косой проекции Меркатора.
The purpose of computing is insight, not numbers

Ответить

Вернуться в «Материалы сайта»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 13 гостей