Форумы GIS-Lab.info

Здравствуйте!
Помогите, пожалуйста, решить проблемму со сдвигом координат.
У меня есть два TAB файла, с системами координат - одна в метрах, другая в градусах.
GDAL говорит, что они такие:
файл A: файл Б: Онлайн сервис для конвертации говорит, что такие:
файл А: файл Б: В файле A у меня полигон, в файле Б точка с координатами в градусах: я хочу ее сконвертировать в систему координат файла А (в метрах).
При конвертации во всех онлайн сервисах и с помощью GDAL и с помощью openlayers она переводится в точку: Что не правильно относительно файла A (точка не попадает в полигон)
И только в MapInfo она переводится в точку: что правильно, относительно полигона.
Вопрос: Почему возникают такие отличия? Может ли это быть из-за того, что
"degree" в GDAL 0.0174532925199433, а не 0.017453292519943295 ?
Или из-за того что фактически LOCAL_CS это не указанная система координат и между ней и EPSG:3857 может быть какой то сдвиг?

Snef писал(а): ↑

08 апр 2020, 11:46

У меня есть два TAB файла, с системами координат - одна в метрах, другая в градусах.

Выкладывайте оба TAB файла или их фрагменты.

Выкладываю файлы.
Так отображается в MapInfo:

А так после конвертации другими способами (онлайн сервисы, GDAL, openlayers):

Snef писал(а): ↑

08 апр 2020, 11:46

LOCAL_CS["Nonearth",UNIT["Meter",1.0]]

Должна быть описана LOCAL_CS, скорее всего - это МСК. В таком виде как у Вас не пойдет.

Спасибо большое! А почему Map Info может правильно преобразовать WGS 84 (которая в градусах) к этой самой LOCAL_CS, если открыть сначала файл с полигонами, а потом с точками? Можно ли как то узнать что это за LOCAL_CS? В .wor файле есть такая информация:

Snef писал(а): ↑

08 апр 2020, 14:36

А почему Map Info может правильно преобразовать WGS 84

Вы уверены в правильности? Скорее всего, наложение произошло случайным образом.

lam писал(а): ↑

08 апр 2020, 14:44

Вы уверены в правильности?

Да, именно в map info точки там где должны быть. Это очень странно и похоже на какую-то погрешность.

Snef писал(а): ↑

08 апр 2020, 11:46

В файле A у меня полигон, в файле Б точка с координатами в градусах:
КОД: ВЫДЕЛИТЬ ВСЁ

13.685983, 6.276257

Ваши координаты в градусах попадают в Африку, так и должно быть?

lam писал(а): ↑

08 апр 2020, 15:33

Ваши координаты в градусах попадают в Африку, так и должно быть?

Да, с этим все правильно (должно быть где-то в Камеруне ) Важно только взаимное расположение объектов слоев.

Snef писал(а): ↑

08 апр 2020, 16:18

Важно только взаимное расположение объектов слоев.

Необходимо знать или попытаться определить параметры системы координат в LOCAL_CS. Только после этого можно производить трансформацию.
Почему у Вас не совпадают результаты трансформаций между онлайн сервисами и мапинфо?
Да потому, что Вы для онлайн сервисов обозначаете LOCAL_CS как:

Snef писал(а): ↑

08 апр 2020, 11:46

WGS 84 / Pseudo-Mercator (EPSG:3857)

В мапинфо же трансформация вообще не происходит по всем правилам, потому что LOCAL_CS указана как "Non-Earth (meters)", а для этой системы координат трансформация вообще не предусмотрена.
Укажите для мапинфо что LOCAL_CS:

Snef писал(а): ↑

08 апр 2020, 11:46

WGS 84 / Pseudo-Mercator (EPSG:3857)

и получите одинаковые результаты с онлайн сервисами.

Спасибо, стало яснее! Скажите, пожалуйста, каким образом Map Info переводит: в или где об этом можно прочитать? Я думаю это распространенный случай?
По LOCAL_CS - в настройках map info стоит "План-схема (метры)".

Данные на самом деле кривые, и скорее всего, два слоя наложены в предположении, что LOCAL_CS это EPSG:3857 (отсюда попадание в Африку) - просто не стали разбираться с системой координат. Затем уже в системе WGS 84 нарисовали в Африке слой с точками, так, чтобы они визуально совпадали со слоем полигонов.
Т.е например была точка полигона в LOCAL_CS 1518467.7, 696453.07 и просто сказали, что это не LOCAL_CS-план схема, а EPSG:3857 и точка в Африке.

Было бы очень нужно узнать как MapInfo переводит 13.685983, 6.276257 в свои 1521811.2, 697887.6, чтобы просто отобразить ее (на другой карте) так же как в MapInfo.

Запутали Вы нас со своей Африкой
Как я и предположил ранее, у Вас МСК-10 зона 1(https://mapbasic.ru/msk10)
С помощью мапинфо сохраняете обратно точки, которые у Вас в WGS 84, в "План-схема (метры)".
Потом присваиваете и полигонам и точкам параметры МСК-10 зона 1 и пересчитываете в WGS84 во всех онлайн сервисах и с помощью GDAL и с помощью openlayers.

Snef писал(а): ↑

09 апр 2020, 09:06

нужно

Разработчики Мапинфы сделали глупость, и вместо того, чтобы ругаться на попытки трансформировать такие системы координат, их трансформируют, почти наверняка считая, что это градусы, и скорее всего на WGS84. Гадать не нужно, нужно писать нормальную СК.

gamm писал(а): ↑

11 апр 2020, 11:13

Разработчики Мапинфы сделали глупость, и вместо того, чтобы ругаться на попытки трансформировать такие системы координат, их трансформируют, почти наверняка считая, что это градусы, и скорее всего на WGS84. Гадать не нужно, нужно писать нормальную СК.

Полностью с Вами согласен, "План-схему (метры)" надо закрывать для возможности координатных трансформаций, это уже не первые слои TAB которые мне попадались с такими градусами на WGS84.

gamm, lam огромное вам спасибо за помощь!