Сдвинуть растр на 360 градусов восточнее

[GWolf]

Обрезал ЦМР ETopo1 по границам СССР.

ETopo_USSR.png (56.09 КБ) 8538 просмотров

Теперь надо бы тот кусочек, который расположен восточнее долготы 180 и который одиноко болтается слева перенести на 360 градусов на восток, чтобы он соединился с остальной частью СССР (надеюсь санкции на меня за это не наложат).

Работаю в QGIS и SAGA. В QGIS есть модуль Affine Transformations, но он работает только с векторными объектами.
Может кто знает как это сделать?

[Petruxin]

А может проще ?

[ericsson]

или так http://habrahabr.ru/post/235283/

[GWolf]

А может проще ?

или так http://habrahabr.ru/post/235283/

Я так понимаю, что вы предлагаете перейти к проекции в которой бы эти два кусочка были рядом. Но мне нужно получить цельный растр физически а не визуально. Чтобы потом по нему считать морфометрию. Хотя может я не понял, что вы хотели сказать.

[trir]

а у вас физически получилось два растра?

[ericsson]

А что, по-вашему, получается, если сменить проекцию на такую, где СССР не разбивается 180-м меридианом, а потом экспортировать этот растр?

[GWolf]

а у вас физически получилось два растра?

Да, я сначала вырезал границей СССР основную часть, а потом сдвинув ее на 360 градусов западнее вырезал тот маленький кусок. Но это не принципиально. Уже переделал.

А что, по-вашему, получается, если сменить проекцию на такую, где СССР не разбивается 180-м меридианом, а потом экспортировать этот растр?

Получается вот так:

ETopo_Albers-Siberia.png (189.04 КБ) 8405 просмотров

Это в проекции Albers-Siberia, описание которой есть по ссылке, которую дал Petruxin. Действительно кусочек оказался там, где он должен быть, но теперь белый разрыв все портит.

Причем в оригинальном не спроецированном файле ETopo1 данные на область этого разрыва есть, в чем можно убедиться например взглянув на остров Врангеля. В спроецированном варианте у него отсутствует западная половина, а в не спроецированном она есть.

ETopo_WGS84_arrow.png (106.24 КБ) 8405 просмотров

Как с этим бороться?

[GWolf]

Справиться с этой проблемой невозможно?

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

[ericsson]

Картинка похожа на результат ошибки в софте.
Может попробуете любую другую проекцию, решающую те же задачи? Обычно, глюки подобного рода связаны с какой-то конкретной проекцией.

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

Зачем же так уродовать исходные данные? Растр дан в географических координатах, так что логично использовать цилиндрическую проекцию, лучше равнопромежуточную. Пример:

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

+proj=eqc +lat_0=0 +lon_0=100 +lat_ts=60 +x_0=0 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs