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

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

Добавлено: 25 янв 2015, 19:30
GWolf
Обрезал ЦМР ETopo1 по границам СССР.
ETopo_USSR.png
ETopo_USSR.png (56.09 КБ) 8373 просмотра
Теперь надо бы тот кусочек, который расположен восточнее долготы 180 и который одиноко болтается слева перенести на 360 градусов на восток, чтобы он соединился с остальной частью СССР (надеюсь санкции на меня за это не наложат).

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

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

Добавлено: 25 янв 2015, 19:47
Petruxin
А может проще ?

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

Добавлено: 25 янв 2015, 20:44
ericsson

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

Добавлено: 25 янв 2015, 21:52
GWolf
Petruxin писал(а):А может проще ?
ericsson писал(а):или так http://habrahabr.ru/post/235283/
Я так понимаю, что вы предлагаете перейти к проекции в которой бы эти два кусочка были рядом. Но мне нужно получить цельный растр физически а не визуально. Чтобы потом по нему считать морфометрию. Хотя может я не понял, что вы хотели сказать.

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

Добавлено: 25 янв 2015, 22:33
trir
а у вас физически получилось два растра?

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

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

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

Добавлено: 26 янв 2015, 13:06
GWolf
trir писал(а):а у вас физически получилось два растра?
Да, я сначала вырезал границей СССР основную часть, а потом сдвинув ее на 360 градусов западнее вырезал тот маленький кусок. Но это не принципиально. Уже переделал.
ericsson писал(а):А что, по-вашему, получается, если сменить проекцию на такую, где СССР не разбивается 180-м меридианом, а потом экспортировать этот растр?
Получается вот так:
ETopo_Albers-Siberia.png
ETopo_Albers-Siberia.png (189.04 КБ) 8240 просмотров
Это в проекции Albers-Siberia, описание которой есть по ссылке, которую дал Petruxin. Действительно кусочек оказался там, где он должен быть, но теперь белый разрыв все портит.

Причем в оригинальном не спроецированном файле ETopo1 данные на область этого разрыва есть, в чем можно убедиться например взглянув на остров Врангеля. В спроецированном варианте у него отсутствует западная половина, а в не спроецированном она есть.
ETopo_WGS84_arrow.png
ETopo_WGS84_arrow.png (106.24 КБ) 8240 просмотров
Как с этим бороться?

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

Добавлено: 27 янв 2015, 19:51
GWolf
Справиться с этой проблемой невозможно?

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

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

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

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

Добавлено: 28 янв 2015, 08:57
Игорь Белов
Зачем же так уродовать исходные данные? Растр дан в географических координатах, так что логично использовать цилиндрическую проекцию, лучше равнопромежуточную. Пример:

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

+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