Перепроецирование из ГСК в проекцию Альберса

Ответить
bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Перепроецирование из ГСК в проекцию Альберса

Сообщение bolotoved » 11 окт 2012, 16:58

При попытке перепроецировать из географической системы координат в проекцию Альберса:

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

gdalwarp -s_srs EPSG:4326 -t_srs "+proj=aea +lat_1=52 +lat_2=64 +lat_0=0 +lon_0=105 +x_0=18500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs" average_temp_gcs.tif average_albers.tif
Конкретно этот файл глючит (с другими данная операция проходила нормально). В чем может быть дело?


Версия GDAL/OGR при сборке 1.9.1
Вложения
average_temp_albers.png
average_temp_albers.png (119.93 КБ) 7640 просмотров
average_temp_gcs.7z
(148.83 КБ) 399 скачиваний
Последний раз редактировалось bolotoved 12 окт 2012, 04:16, всего редактировалось 1 раз.

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Перепроецирование из ГСК в проекцию Альберса

Сообщение Александр Мурый » 11 окт 2012, 17:29

При перепроецировании <gdalwarp> (GDAL 1.9.1) выдаёт ошибку "ERROR 1: tolerance condition error". Быстрый гуглинг показал, что проблема давно известная; только я не понял, решённая ли.

При открытии в QGIS, например, растр отображается нормально:
rast.png
rast.png (57.93 КБ) 7634 просмотра
Редактор материалов, модератор форума

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Re: Перепроецирование из ГСК в проекцию Альберса

Сообщение bolotoved » 11 окт 2012, 17:35

amuriy писал(а):При открытии в QGIS, например, растр отображается нормально:
Картинка выше - тоже в QGIS, проекция проекта = проекции полученного растра

Версия QGIS
1.9.0-Master
Ревизия
f262caa

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9128
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 747
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: Перепроецирование из ГСК в проекцию Альберса

Сообщение Максим Дубинин » 11 окт 2012, 21:45

приложите куда-нибудь растр, пожалуйста.
пристегивайтесь, турбулентность прямо по курсу

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Перепроецирование из ГСК в проекцию Альберса

Сообщение Александр Мурый » 11 окт 2012, 21:49

Максим Дубинин писал(а):приложите куда-нибудь растр, пожалуйста.
Уже приложен в первом посте темы.
Редактор материалов, модератор форума

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9128
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 747
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Re: Перепроецирование из ГСК в проекцию Альберса

Сообщение Максим Дубинин » 11 окт 2012, 22:20

Ой, не заметил, пардон, а пример который не глючит? :) У неглючащего такой же охват?
пристегивайтесь, турбулентность прямо по курсу

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Перепроецирование из ГСК в проекцию Альберса

Сообщение Александр Мурый » 11 окт 2012, 23:58

Максим Дубинин писал(а):а пример который не глючит?
average_albers.tif.7z
(271.58 КБ) 450 скачиваний
Редактор материалов, модератор форума

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Понял, в чем дело

Сообщение bolotoved » 12 окт 2012, 04:25

Растры, которые я перепроецировал представляют собою вырезки по границам россиии до 180 долготы. Перепроецирование происходит нормально, когда разрешение растра достаточно высокое, например 3672 х 976 пикселей. В случае, если растр имеет более крупные пиксели, например, 918 х 244, то они сильно выходят за 180 долготу, что воспринимается GDAL как выход за допустимые границы (долгота > 180 градуса), что приводит к глюкам перепроецирования.
В приложении два файла - растеризованный слой границ РФ с разрешением 918 х 244 и 3672 х 976. Первый при перепоецировании ломается, второй, хотя при его перепроецировании и выдается ошибка выхода за пределы границ, перепроецируется нормально.
Вложения
rus.7z
(40.46 КБ) 434 скачивания

Ответить

Вернуться в «GDAL/OGR»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 16 гостей