GDAL MapInfo > proj4, подскажите новичку

Ответить
gadkiy_uut
Активный участник
Сообщения: 139
Зарегистрирован: 04 апр 2011, 18:32
Репутация: 0
Откуда: Москва

GDAL MapInfo > proj4, подскажите новичку

Сообщение gadkiy_uut » 20 мар 2018, 19:19

Может быть недостаточно читал документацию, подскажите, есть такой пример:

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

$ gdalsrsinfo -o proj4 osr/data/lcc_esri.prj
Можно ли подготовить на основе этих данных:
http://www.mapbasic.ru/msk50
Файл формата *.prj Что бы получить из описания СК в формате MapInfo описание в proj4? Пока что делаю это вручную, в текстовом редакторе, но сейчас утыкаюсь в непонимание как прописывать towgs84, хотел проверить это через GDAL. У меня получается вот так:

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

+proj=tmerc +lat_0=0 +lon_0=38.48333333333 +k=1 +x_0=2250000 +y_0=-5712900.566 +ellps=krass +towgs84=23.57, -140.95, -79.8, 0, 0.35, 0.79, -0.22 +units=m +no_defs
Правила преобразования есть тут:
http://wiki.gis-lab.info/w/%D0%A2%D1%80 ... %B8_PROJ.4
Но у меня пока не получается до конца разобраться с towgs. Приходится дополнительно спрашивать подсказки на форуме.

Аватара пользователя
Denis Rykov
Гуру
Сообщения: 3376
Зарегистрирован: 11 апр 2008, 21:09
Репутация: 529
Ваше звание: Author
Контактная информация:

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение Denis Rykov » 20 мар 2018, 19:48

Насколько я помню, тут есть уже готовый список МСК: http://nextgis.ru/nextgis-qgis/
Spatial is now, more than ever, just another column- The Geometry Column.

trir
Гуру
Сообщения: 5271
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1013
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение trir » 20 мар 2018, 20:08

Правила преобразования есть тут:
плохо читаете - Ссылки п.9

trir
Гуру
Сообщения: 5271
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1013
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение trir » 20 мар 2018, 20:15

Но у меня пока не получается до конца разобраться с towgs.
а в чём затык? их обсуждали в обсуждении

gadkiy_uut
Активный участник
Сообщения: 139
Зарегистрирован: 04 апр 2011, 18:32
Репутация: 0
Откуда: Москва

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение gadkiy_uut » 21 мар 2018, 10:02

trir писал(а):
20 мар 2018, 20:08
плохо читаете - Ссылки п.9
Попробовал вариант из п. 9, ссылок, файл res4.txt. На мой взгляд данные для МСК-50 там представлены некорректно и вызывает сомнение вся информация в res4.txt. Кто составитель данного файла?

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

+proj=tmerc +ellps=krass +towgs84=-118.754,-61.782,-93.237,2.40896,3.47502,-1.29688,6.5177 +units=m +lon_0=38.48333333333 +lat_0=0 +k_0=1 +x_0=2250000 +y_0=-5712900.566
Выглядит вот так (на глаз вижу что не корректно):
Снимок1.PNG
Снимок1.PNG (138.19 КБ) 13716 просмотров
Старый вариант выглядит более корректным (опять же на глаз) и совпадает с данными http://www.mapbasic.ru/msk50:

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

+proj=tmerc +lat_0=0 +lon_0=38.48333333333 +k=1 +x_0=2250000 +y_0=-5712900.566 +ellps=krass +towgs84=23.57, -140.95, -79.8, 0, 0.35, 0.79, -0.22 +units=m +no_defs
Снимок2.PNG
Снимок2.PNG (106.73 КБ) 13716 просмотров
И всё же вопрос, с помощью gdalsrsinfo можно преобразовать данные строк http://www.mapbasic.ru в proj4 автоматом? Каждый раз методом тыка подбирать towgs сложно, хотелось бы наработанный алгоритм.

trir
Гуру
Сообщения: 5271
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1013
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение trir » 21 мар 2018, 10:14

Попробовал вариант из п. 9, ссылок, файл res4.txt
там параметры отсюда, но брал я их очень давно и эти параметры явно были не правильные - раз их заменили в исходнике
Каждый раз методом тыка подбирать towgs сложно, хотелось бы наработанный алгоритм.
идёте в поле и снимаете пункты, а потом их уравниваете и расчитываете правильные параметры. А без поля можно только ГОСТ'ом пользоватся
Последний раз редактировалось trir 21 мар 2018, 13:23, всего редактировалось 1 раз.

trir
Гуру
Сообщения: 5271
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1013
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение trir » 21 мар 2018, 10:16

И всё же вопрос, с помощью gdalinfo можно преобразовать данные строк http://www.mapbasic.ru в proj4 автоматом?
можно этим - Ссылки п.8

gadkiy_uut
Активный участник
Сообщения: 139
Зарегистрирован: 04 апр 2011, 18:32
Репутация: 0
Откуда: Москва

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение gadkiy_uut » 21 мар 2018, 10:27

trir писал(а):
21 мар 2018, 10:16
можно этим - Ссылки п.8
Скачал. Как этим пользоваться?

trir
Гуру
Сообщения: 5271
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1013
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение trir » 21 мар 2018, 10:32

компелируешь и нажимаешь кнопку Button1

gadkiy_uut
Активный участник
Сообщения: 139
Зарегистрирован: 04 апр 2011, 18:32
Репутация: 0
Откуда: Москва

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение gadkiy_uut » 21 мар 2018, 10:57

В какой среде?

trir
Гуру
Сообщения: 5271
Зарегистрирован: 09 апр 2010, 19:30
Репутация: 1013
Ваше звание: просто мимо прохожу
Откуда: Ё-бург

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение trir » 21 мар 2018, 11:03


Oval
Участник
Сообщения: 61
Зарегистрирован: 30 дек 2014, 14:40
Репутация: 14

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение Oval » 21 мар 2018, 13:21

gadkiy_uut писал(а):
21 мар 2018, 10:02
Выглядит вот так (на глаз вижу что не корректно):
есть уверенность что картографическая основа точна?

gadkiy_uut
Активный участник
Сообщения: 139
Зарегистрирован: 04 апр 2011, 18:32
Репутация: 0
Откуда: Москва

Re: GDAL MapInfo > proj4, подскажите новичку

Сообщение gadkiy_uut » 21 мар 2018, 15:02

Oval писал(а):
21 мар 2018, 13:21
есть уверенность что картографическая основа точна?
Проверял с OSM подосновой. Понимание есть что точна, это видно на глаз, в первом случае точность туда-сюда 2-3 метра, во втором улетает метров на 20. Можно наложить и на другие подосновы (openlayers плагин к сожалению уже не работает в версии 3.0 но я успел проверить в версии 2.8, так же бъётся на слоях фотоснимков)
Единственное на вторм скриншоте случайно отключил один слой с дорогами.
И вот на мой взгляд корректные параметры:

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

+proj=tmerc +lat_0=0 +lon_0=38.48333333333 +k=1 +x_0=2250000 +y_0=-5712900.566 +ellps=krass +towgs84=23.57, -140.95, -79.8, 0, 0.35, 0.79, -0.22 +units=m +no_defs
Уточнял эти параметры тут:
viewtopic.php?f=27&t=21881&p=149332#p149328
Помогли участники форума. Спасибо Giser, помог уточнить.

Ответить

Вернуться в «GDAL/OGR»

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

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