Страница 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](./download/file.php?id=5158&sid=716183bac6e0ac2af85d6e39a28486a2)
- 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
Максим Дубинин
Ой, не заметил, пардон, а пример который не глючит?
![Smile :)](./images/smilies/icon_smile.gif)
У неглючащего такой же охват?
Re: Перепроецирование из ГСК в проекцию Альберса
Добавлено: 11 окт 2012, 23:58
Александр Мурый
Максим Дубинин писал(а):а пример который не глючит?
Понял, в чем дело
Добавлено: 12 окт 2012, 04:25
bolotoved
Растры, которые я перепроецировал представляют собою вырезки по границам россиии до 180 долготы. Перепроецирование происходит нормально, когда разрешение растра достаточно высокое, например 3672 х 976 пикселей. В случае, если растр имеет более крупные пиксели, например, 918 х 244, то они сильно выходят за 180 долготу, что воспринимается GDAL как выход за допустимые границы (долгота > 180 градуса), что приводит к глюкам перепроецирования.
В приложении два файла - растеризованный слой границ РФ с разрешением 918 х 244 и 3672 х 976. Первый при перепоецировании ломается, второй, хотя при его перепроецировании и выдается ошибка выхода за пределы границ, перепроецируется нормально.