GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли


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

Обсуждение материалов сайта: вопросы, замечания, предложения

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

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

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

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

Игорь Белов
Гуру
 
Зарегистрирован: 04 янв 2011
Откуда: Казань (на карте)
 
Сообщения: 1241
Репутация: 743
Статьи: 12
Проекты: 1

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

Сообщение vveco » 21 мар 2015, 19: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

Что не так?
  • 0

Аватара пользователя
vveco
Активный участник
 
Зарегистрирован: 20 май 2009
Откуда: Саратов (на карте)
 
Сообщения: 189
Репутация: 3
Проекты: 2

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

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

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

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4818
Репутация: 648
Статьи: 3
Проекты: 5/1

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

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

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

Это первое.

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

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

Игорь Белов
Гуру
 
Зарегистрирован: 04 янв 2011
Откуда: Казань (на карте)
 
Сообщения: 1241
Репутация: 743
Статьи: 12
Проекты: 1

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

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

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

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4818
Репутация: 648
Статьи: 3
Проекты: 5/1

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

Сообщение Charaunica » 08 май 2015, 13: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
  • 0

Вложения
dxf.rar
dxf , который перепроецировался в МСК объекта
(586.84 КБ) Скачиваний: 234
Проверка_QGIS.jpg
Проверка проекции созданной проекции в QGIS
Проверка_QGIS.jpg (67.18 КБ) Просмотров: 8390
Проверка_shell.jpg
Процесс работы утилиты
Проверка_shell.jpg (86.77 КБ) Просмотров: 8390
Charaunica
Интересующийся
 
Зарегистрирован: 11 мар 2013
 
Сообщения: 33
Репутация: 1

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

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

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

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

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

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

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


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

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

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

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

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

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

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

Сообщение Charaunica » 08 май 2015, 15: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/ по созданию разных МСК, вот и вникаю во всю эту красоту)
  • 0

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

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

Сообщение Игорь Белов » 08 май 2015, 20: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
  • 1

Вложения
airdrome.zip
(17.9 КБ) Скачиваний: 251
Игорь Белов
Гуру
 
Зарегистрирован: 04 янв 2011
Откуда: Казань (на карте)
 
Сообщения: 1241
Репутация: 743
Статьи: 12
Проекты: 1


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

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

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


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).