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

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

Добавлено: 11 окт 2012, 16:58
bolotoved
При попытке перепроецировать из географической системы координат в проекцию Альберса:

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

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

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

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

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

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

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

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

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

Добавлено: 11 окт 2012, 21:45
Максим Дубинин
приложите куда-нибудь растр, пожалуйста.

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

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

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

Добавлено: 11 окт 2012, 22:20
Максим Дубинин
Ой, не заметил, пардон, а пример который не глючит? :) У неглючащего такой же охват?

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

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

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

Добавлено: 12 окт 2012, 04:25
bolotoved
Растры, которые я перепроецировал представляют собою вырезки по границам россиии до 180 долготы. Перепроецирование происходит нормально, когда разрешение растра достаточно высокое, например 3672 х 976 пикселей. В случае, если растр имеет более крупные пиксели, например, 918 х 244, то они сильно выходят за 180 долготу, что воспринимается GDAL как выход за допустимые границы (долгота > 180 градуса), что приводит к глюкам перепроецирования.
В приложении два файла - растеризованный слой границ РФ с разрешением 918 х 244 и 3672 х 976. Первый при перепоецировании ломается, второй, хотя при его перепроецировании и выдается ошибка выхода за пределы границ, перепроецируется нормально.