Страница 1 из 1

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

Добавлено: 12 ноя 2014, 09:27
MI user
Уважаемые коллеги, нужна ваша подсказка в вопросе перепроецирования из EPSG:3857 в местную систему координат (МСК).
Проблема в том, что один и тот же файл перепроецируется из EPSG:3857 в МСК:
  • в MapInfo корректно
  • с помощью GDAL (или QGIS) со смещением
Вот как это выглядит:
2 layers - Mapinfo.png
Пример из MapInfo
2 layers - Mapinfo.png (31.33 КБ) 11090 просмотров
2 layers - QGIS.png
Пример из QGIS (файл преобразован в GDAL)
2 layers - QGIS.png (113.35 КБ) 11090 просмотров
Исходные данные:
  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)?
Спасибо.

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

Добавлено: 12 ноя 2014, 09:35
ericsson
Точно измерить сдвиг в метрах по X,Y и сравнить с параметрами towgs пробовали?

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

Добавлено: 12 ноя 2014, 09:36
Игорь Белов
Поменяйте знаки параметров разворота:

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

+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22

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

Добавлено: 12 ноя 2014, 14:01
MI user
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 КБ) 11044 просмотра

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

Добавлено: 12 мар 2017, 14:16
t[0]p
Игорь Белов писал(а):Поменяйте знаки параметров разворота:

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

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

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

Добавлено: 12 мар 2017, 14:48
gamm
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