Конвертация план-схемы в WGS-84 в Mapinfo

[Institor]

Такой вот экстент карты

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

Bounds (6931.72539153, -1099.33453179) (14823.5243692, 25228.1022486)

не очень удобен, разная точность по осям получается. Я бы добавил в конец строки еще

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

, -100000, -100000, 100000, 100000

Тогда погрешность будет одинаковая и равная 0,1 мм.

[delman]

Попробуйте CoordSys без Bounds:

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

CoordSys Earth Projection 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, "m", 60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481

Получилось наконец, но карта со смещением на 10 метров и поворотом против часовой стрелки. Как это можно поправить?

[Донецков]

Значит

60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481

не совсем, то что нужно

[delman]

Такой вот экстент карты

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

Bounds (6931.72539153, -1099.33453179) (14823.5243692, 25228.1022486)

не очень удобен, разная точность по осям получается. Я бы добавил в конец строки еще

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

, -100000, -100000, 100000, 100000

Тогда погрешность будет одинаковая и равная 0,1 мм.

Тоесть добавить подряд
Bounds (6931.72539153, -1099.33453179) (14823.5243692, 25228.1022486), -100000, -100000, 100000, 100000

?

И почему в скобках?

[delman]

Значит

60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481

не совсем, то что нужно

Теперь я это понял, для части города подходит, а для всего города по краям большая погрешность, нужно искать другую строку переноса.

[Institor]

Тоесть добавить подряд
Bounds (6931.72539153, -1099.33453179) (14823.5243692, 25228.1022486), -100000, -100000, 100000, 100000

?

И почему в скобках?

Нет. слова "bounds ( ) ( )" заменить на "-100000, -100000, 100000, 100000"

[delman]

Нет. слова "bounds ( ) ( )" заменить на "-100000, -100000, 100000, 100000"

Тогда со строкой

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

CoordSys Earth Projection 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, "m", 60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481, -100000, -100000, 100000, 100000

выдает

2015-10-02_144053.jpg (20.85 КБ) 8715 просмотров

Подскажите автоматические способы смещения групп объектов относительно карты? По хорошему создать бы более точную строку в mif или в файл mapinfo.prj, но чем точнее это сделать не нашел пока.

[Игорь Белов]

со строкой

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

CoordSys Earth Projection 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, "m", 60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481, -100000, -100000, 100000, 100000

выдает ...

Правильная строка в MIF:

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

CoordSys Earth Projection 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, "m", 60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481 Bounds (-100000, -100000) (100000, 100000)

Подскажите автоматические способы смещения групп объектов относительно карты?

Проще всего слой сдвигать изменением параметров false_easting и false_northing (37637.627999999997 и 10157.481 в данном случае). Можно также менять central_meridian и latitude_of_origin (60.535844865282 и 55.049659241555 соответственно).

[Institor]

А, точняк. Забыл сказать, что Projection 8 надо заменить на Projection 2008. Из справки:

Вы можете указать предельные границы для координатной системы в файле MAPINFOW.PRJ. Чтобы сделать это, добавьте 2000 к номеру проекции и укажите границы после списка параметров проекции. Общая форма следующая:

name, projectionnum + 2000, projection parameters, x1, y1, x2, y2

У меня они все в mapinfow.prj забиты с границами, поэтому и забыл про данный нюанс.
Мне так удобнее - можно просто копировать нужную строчку из mapinfow.prj

[delman]

А, точняк. Забыл сказать, что Projection 8 надо заменить на Projection 2008. Из справки:

Вы можете указать предельные границы для координатной системы в файле MAPINFOW.PRJ. Чтобы сделать это, добавьте 2000 к номеру проекции и укажите границы после списка параметров проекции. Общая форма следующая:

name, projectionnum + 2000, projection parameters, x1, y1, x2, y2

У меня они все в mapinfow.prj забиты с границами, поэтому и забыл про данный нюанс.
Мне так удобнее - можно просто копировать нужную строчку из mapinfow.prj

Можно увидеть пример вашего mapinfow.prj с предельными границами?

[delman]

Правильная строка в MIF:

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

CoordSys Earth Projection 8, 9999, 3, 23.57, -140.95, -79.8, 0, -0.35, -0.79, -0.22, 0, "m", 60.535844865282, 55.049659241555, 0.9959922445, 37637.627999999997, 10157.481 Bounds (-100000, -100000) (100000, 100000)

Ничего не изменилось при добавлении после Bounds

Проще всего слой сдвигать изменением параметров false_easting и false_northing (37637.627999999997 и 10157.481 в данном случае). Можно также менять central_meridian и latitude_of_origin (60.535844865282 и 55.049659241555 соответственно).

А есть параметры, отвечающие за кручение (против часовой повернуть все обьекты)? и параметр немного увеличить (растянуть) ее?

[Игорь Белов]

А есть параметры, отвечающие за кручение (против часовой повернуть все обьекты)? и параметр немного увеличить (растянуть) ее?

Во-первых, перенос осевого меридиана (central_meridian 60.535844865282) приведёт к развороту, примерно равному ∆λ sin φ. Сжимайте/растягивайте масштабным коэффициентом (scale_factor 0.9959922445).

Во-вторых, полную свободу в этом вопросе даёт аффинное преобразование (добавляется строка Affine с параметрами).

Внимание! Важно помнить, что предложение CoordSys имеет ограничение на длину, и одновременно впихнуть явное определение датума + Affine + Bounds не удастся.

[Донецков]

Интересно

Внимание! Важно помнить, что предложение CoordSys имеет ограничение на длину, и одновременно впихнуть явное определение датума + Affine + Bounds не удастся.

не встречал об этом ссылок...

[Institor]

Есть ограничение, между 200 и 300 символов, подозреваю, что 256
Но впихнуть можно, просто не нужно гнаться за огромным количеством значащих цифр в параметрах.
Аффинное это отдельная тема. Тоже можно добавить в строку проекции, просто прибавив к номеру еще 1000.

Вы также можете указать координатную систему с границами и аффинным преобразованием. В этом случае добавьте 3000 к номеру проекции, и перечислите границы после параметров аффинного преобразования. Общая форма следующая:

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

Название, номер проекции + 3000, параметры проекции, название единиц измерения, A, B, C, D, E , F, x1, y1, x2, y2

Взять три точки с известными координатами, взять их координаты по "неправильной" строчке проекции, вычислить коэффициенты преобразования между ними в экселе(или где-нибудь еще) по формулам из статьи http://gis-lab.info/qa/affine-math.html, и дело в шляпе. Только это надо на небольшую площадь делать.

Можно увидеть пример вашего mapinfow.prj с предельными границами?

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

"++СК42 Зона 23 GGC", 2008, 9999, 3, 23.9, -141.3, -80.9, 0, -0.37, -0.85, -0.12, 0, 7, 135, 0, 1, 23500000, 0, 23000000, 5000000, 24000000, 6000000

например. Можно еще сузить. Но я считаю, что точности в 0,5 мм хватит всем. Дело вкуса.

[Донецков]

Есть ограничение, между 200 и 300 символов, подозреваю, что 256

- не сталкивался, но нужно проверить...