Разница в перепроецировании GDAL / MapInfo

Ответить
Аватара пользователя
MI user
Участник
Сообщения: 63
Зарегистрирован: 14 янв 2013, 17:08
Репутация: 33
Откуда: Оренбург

Разница в перепроецировании GDAL / MapInfo

Сообщение MI user » 12 ноя 2014, 09:27

Уважаемые коллеги, нужна ваша подсказка в вопросе перепроецирования из EPSG:3857 в местную систему координат (МСК).
Проблема в том, что один и тот же файл перепроецируется из EPSG:3857 в МСК:
  • в MapInfo корректно
  • с помощью GDAL (или QGIS) со смещением
Вот как это выглядит:
2 layers - Mapinfo.png
Пример из MapInfo
2 layers - Mapinfo.png (31.33 КБ) 10995 просмотров
2 layers - QGIS.png
Пример из QGIS (файл преобразован в GDAL)
2 layers - QGIS.png (113.35 КБ) 10995 просмотров
Исходные данные:
  1. Эталонный слой в формате TAB MapInfo в проекции МСК (прилагается -
    Layers.zip
    Слои
    (377.88 КБ) 390 скачиваний
    ).
  2. Файл в формате TAB MapInfo в проекции EPSG:3857, который надо преобразовать (прилагается
    Layers.zip
    Слои
    (377.88 КБ) 390 скачиваний
    ).
  3. Параметры МСК в формате MapInfo:

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

    "МСК зона 1", 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, 7, 52.03333333333, 0, 1, 1300000, -5309414.70
  4. Параметры той же самой МСК, адаптированные для GDAL:

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

    +proj=tmerc +lat_0=0 +lon_0=52.03333333333 +k=1 +x_0=1300000 +y_0=-5309414.7 +ellps=krass +towgs84=23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22 +units=m +no_defs
    
Вопросы:
  1. Правильно ли я составил строку с проекцией для GDAL исходя из параметров MapInfo?
  2. Если строка составлена правильно, то в чём может быть причина смещения координат при преобразовании в GDAL (QGIS)?
Спасибо.

ericsson
Гуру
Сообщения: 3321
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 748
Ваше звание: Вредитель полей

Re: Разница в перепроецировании GDAL / MapInfo

Сообщение ericsson » 12 ноя 2014, 09:35

Точно измерить сдвиг в метрах по X,Y и сравнить с параметрами towgs пробовали?

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

Re: Разница в перепроецировании GDAL / MapInfo

Сообщение Игорь Белов » 12 ноя 2014, 09:36

Поменяйте знаки параметров разворота:

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

+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22
The purpose of computing is insight, not numbers

Аватара пользователя
MI user
Участник
Сообщения: 63
Зарегистрирован: 14 янв 2013, 17:08
Репутация: 33
Откуда: Оренбург

Re: Разница в перепроецировании GDAL / MapInfo

Сообщение MI user » 12 ноя 2014, 14:01

ErnieBoyd писал(а):Поменяйте знаки параметров разворота:

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

+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22
ОГРОМНОЕ СПАСИБО!
Поменял параметры и всё получилось:
2 layers -OK!.png
2 layers -OK!.png (21.15 КБ) 10949 просмотров

Аватара пользователя
t[0]p
Участник
Сообщения: 93
Зарегистрирован: 15 ноя 2014, 14:35
Репутация: 7
Откуда: Томск

Re: Разница в перепроецировании GDAL / MapInfo

Сообщение t[0]p » 12 мар 2017, 14:16

Игорь Белов писал(а):Поменяйте знаки параметров разворота:

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

+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22
Я что то не совсем понимаю... какие знаки ? Имеется в виду точки на запятые ? Или убрать пробелы?

gamm
Гуру
Сообщения: 4044
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1049
Ваше звание: программист
Откуда: Казань

Re: Разница в перепроецировании GDAL / MapInfo

Сообщение gamm » 12 мар 2017, 14:48

t[0]p писал(а):Я что то не совсем понимаю... какие знаки ? Имеется в виду точки на запятые ? Или убрать пробелы?
знаки - они у чисел, плюсы и минусы. Только разберитесь где смещение, где углы (менять нужно знаки углов).
towgs84 - Datum transformation to WGS84
The seven parameter case uses delta_x, delta_y, delta_z, Rx - rotation X, Ry - rotation Y, Rz - rotation Z, M_BF - Scaling.
PROJ4

Ответить

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

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

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