Импорт и обработка спроецированной растровой карты

[Странник]

Нужно загрузить растровое изображение вот такой карты http://www.zin.ru/animalia/coleoptera/i ... arctic.gif
в ГИС (ArcView, ArcInfo, QGis)
Эти области нужно оцифровать по линиям и наложить на нормальную карту мира.
Проблема в том, что данная карта растровая и не понятно в какой проекции. Я боюсь, что при деформации карты вручную по опорным точкам она может неправильно исказиться.

Каков порядок действий, что нужно сделать и как?

[Petruxin]

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

[ Сообщение с мобильного устройства ]

[Странник]

Тут нелинейное преобразование требуется, точками я запарюсь его добиваться

[gamm]

Каков порядок действий, что нужно сделать и как?

план А) узнать у тех, кто делал карту, проекцию (например, найти печатный официальный вариант этой или аналогичной карты из той же конторы, на которой есть проекция)

план В) поскольку есть линии градусной сетки, оцифровать их, и автоматически наставить вдоль них много точек, получим ячейки градусной сетки с детальными сторонами. Дальше интерполировать линии сторон отдельных ячеек, например построив на них полиномы низкой степени, и линейно интерполировать коэффициенты между линиями, чтобы сгустить сетку (получится плавный переход верхней линии ячейки в нижнюю, и аналогично для сторон). Глобальные полиномы лучше не использовать.

план С) поскольку проекций полушарий (типов проекций) не так много, то взять из плана (В) точки линий, и подобрать проекцию и ее параметры скриптов в R или Pyton, займет минут двадцать счета. Если точно не подберется, то невязки убрать по плану (В)

[Странник]

Спасибо! Не ожидал, что так сложно. В Фотошопе предварительно попробовал трансформировать, нелинейной деформацией, но там только вручную растянут можно и центр изображения указать нельзя. Точность конечно никакая.
Я подозреваю, что эти линии были нарисованы экспертно на глазок, в Фотошопе поверх готовой проецированной карты, которую взяли откуда-то в электронном виде.

[Sergey Astakhov]

На глаз проекция на вот эту похожа - http://en.wikipedia.org/wiki/Mollweide_projection

[Странник]

Скачал и установил FlexProjector. Куча возможностей, но ничего не могу при помощи него сделать. Ни файл загрузить, чтобы подобрать проекцию и визуально его перепроецировать, ни в слепую это выполнить не получается. Читаю документацию, но ничего из того, что написано в документации с растровыми файлами сделать не получается. Пишет либо "не могу загрузить", либо "ошибка", либо что-то ему не хватает, либо пропорции файла не нравятся. Брр... Зачем так пишут, проще надо. Как им пользоваться?

[Странник]

Похоже в нем только из географических координат в разные проекции можно карты преобразовать, но не наоборот.

[Странник]

На глаз проекция на вот эту похожа - http://en.wikipedia.org/wiki/Mollweide_projection

Я не нашел эту проекцию в QGIS, при создании слоя. Поликоническую тоже не нашел.
Пока, в надежде на светлое будущее, векторизую карту как есть, указав систему координат WGS84 равнопромежуточную, элипсоид . Потом надо будет переназначить правильную.
В меню QGIS обозначены какие то преобразования проекций, но я ещё не разбирался.

Буду благодарен за полезные советы.

ps...

Вот, кстати, только что получил от авторов уточненные данные по проекции:
World_Robinson
Authority: Custom

Projection: Robinson
False_Easting: 0,0
False_Northing: 0,0
Central_Meridian: 100,0
Linear Unit: Meter (1,0)

Geographic Coordinate System: GCS_WGS_1984
Angular Unit: Degree (0,0174532925199433)
Prime Meridian: Greenwich (0,0)
Datum: D_WGS_1984
Spheroid: WGS_1984
Semimajor Axis: 6378137,0
Semiminor Axis: 6356752,314245179
Inverse Flattening: 298,257223563

В связи с этим вопросы:
1) Как эти данные добавить в проект, и не нарушится ли при этом геометрия уже сделанных шейпов?
2) Приведёт ли добавление этих данных к правильному масштабированию и позиционированию? Я что-то не понимаю принципа, карта-то обрезанная.
3) Как эти данные использовать для трансформации?

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

Правильная последовательность действий такова.

Создаёте пользовательскую проекцию с такими параметрами:

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

+proj=robin +lon_0=100 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs

Назовите эту проекцию, скажем, "WGS 84 / Robinson 100E".

Трансформировать исходный растр глупо, поскольку это только испортит совершенный оригинал. Нужно просто для palearctic.gif сформировать правильный файл привязки. У меня получился такой файл palearctic.wld:

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

3847.9
0
0
-3847.9
-11914228
8820535

Помещаете этот файл в папку, где лежит растр. Открываете растр; по запросу QGIS указываете для него созданную Вами координатную систему. Ту же самую систему назначаете для проекта.

не нарушится ли при этом геометрия уже сделанных шейпов?

Если геометрия нарушится, значит, она неправильная.

[Странник]

Вот это то, что надо, спасибо!
Шейпов пока не много, так что рискну.
Циферки привязки, правда, непонятны для меня.

[Донецков]

Циферки привязки, правда, непонятны для меня

найдите и прочитайте ппро спецификацию world-файла, на этом сайте такое было...
gis-lab.info/qa/tfw.html

[udav]

На этой схеме (это не карта, а схема) всего-то линий столько, что их можно пересчитать по пальцам рук.
Поэтому еще более правильная последовательность действий такая: возьмите какую-нибудь удобную привязанную подложку и нарисуйте на ней заново. Вы не только получите простейшим и быстрейшим путем неискаженную схему, но и сделаете ее лучше оригинала Ведь легко заметить, что авторы схемы проводили рубежи по отлично видным физико-географическим рубежам. А вы не просто проведете по ним же, а проведете с еще большей точностью.

Кстати, с зоогеографической точки зрения продукт так себе. Они даже не удосужились придумать, куда относить переходные зоны. Автор явно не в ладах со своей же темой.