Перепроектирование ASTER/Terralook средствами GDAL

Ответить
Михаил Парилов
Интересующийся
Сообщения: 21
Зарегистрирован: 10 ноя 2007, 16:55
Репутация: 0
Откуда: Архара, Амурская область
Контактная информация:

Перепроектирование ASTER/Terralook средствами GDAL

Сообщение Михаил Парилов » 23 ноя 2007, 08:49

Пытаюсь использовать GDAL для перепроектировки Терралук Астер.
Хочу на выходе получить файл спроектированный в Albers Equal Area conic, Russia- Far East. Прописываю следующую команду:

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

gdalwarp -s_srs "+proj=longlat +datum=WGS84 +no_defs" -t_srs "+proj=longlat +lat_0=52 +lon_0=135 +ellps=krass +datum=Pulkovo_1942 +x_0=23500000 +y_0=0 +units=m +no_defs" -tr 15 15 -rn  -srcnodata 0 -dstnodata 0 E:\images\AST_L1B_00309262002022630_20071109174321_13397.jpg A09262002.tif
Выдается сообщение:

"Creating output file that is 0P x 0L.
ERROR1: Attempt to create 0x0x3 TIFF file, but width, height and bands must be positive"

Понятно, что речь идет о параметрах растра, но непонятно где они должны быть указаны...
Михаил Парилов

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

Re: Перепроектирование ASTER/Terralook средствами GDAL

Сообщение Максим Дубинин » 23 ноя 2007, 20:21

Миша, в команде довольно много глюков, начиная с
-t_srs "+proj=longlat
выходная проекция - Альберс! +proj=aea

не вдаваясь в подробности, вот правильная команда

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

gdalwarp -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84" -t_srs "+proj=aea +lat_1=52 +lat_2=64 lat_0=0 +lon_0=135 +ellps=krass +datum=Pulkovo_1942 +x_0=23500000 +y_0=0" -tr 15 15 -rn -srcnodata 0 -dstnodata 0 AST_L1B_00309262002022630_20071109174321_13397.jpg test.tif
есть тонкости с датумом, но это на потом
пристегивайтесь, турбулентность прямо по курсу

Михаил Парилов
Интересующийся
Сообщения: 21
Зарегистрирован: 10 ноя 2007, 16:55
Репутация: 0
Откуда: Архара, Амурская область
Контактная информация:

Re: Перепроектирование ASTER/Terralook средствами GDAL

Сообщение Михаил Парилов » 24 ноя 2007, 17:34

Максим, еще не всё :)
GDAL выдает
creating output file...
а потом
Error 1: Invalid nBandCount or bMaskIsFloat argument in SourceNoDataMask
Честно говоря, все слова понимаю, а смысл- нет:)
Остается надеяться, что этот пример будет использоваться другими пользователями программы GDAL, потому что, насколько я понимаю, хорошей инструкции по GDAL на русском языке еще не существует, или я не прав?
Даже тот Tutorial который есть на английском на сайте gdal отнюдь не исчерпывающий
Михаил Парилов

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

Re: Перепроектирование ASTER/Terralook средствами GDAL

Сообщение Максим Дубинин » 24 ноя 2007, 20:39

прикольная проблемка, вот что я выяснил

выражение абсолютно правильное и работает у меня без проблем (версия GDAL 1.4.2), бьюсь об заклад, что у тебя версия 1.4.3, в которой создатели GDAL что-то нахимичили и стала появляться эта ошибка (вот тут подробнее про это).

Советую "обновиться" до 1.4.2 и попробовать еще раз. Судя по трекеру, они эту ошибку в 1.4.4 уже исправили, но ее пока не выпустили.

По поводу русских материалов, есть только 2 статьи у нас на сайте:
  • Использование GDAL для перепроектировки данных на базе эллипсоида на сферу
    Использование GDAL для привязки растровых материалов
+ статья про rgdal
пристегивайтесь, турбулентность прямо по курсу

Михаил Парилов
Интересующийся
Сообщения: 21
Зарегистрирован: 10 ноя 2007, 16:55
Репутация: 0
Откуда: Архара, Амурская область
Контактная информация:

Re: Перепроектирование ASTER/Terralook средствами GDAL

Сообщение Михаил Парилов » 25 ноя 2007, 16:58

Максим, ты как всегда прав!
Скачал предпоследнюю версию GDAL по ссылке http://download.osgeo.org/gdal/win32/1. ... exe142.zip
и все пошло как по маслу!
Единственное неудобство, что в результирующем файле при попытке посмотреть "image information" ERDAS выдает "The spheroid "<NULL>" is not recognized" ,т.е. Эрдас не видит датума и проекции в результирующем файле.
Зато как быстро GDAL перепроектирует растры, Эрдасу такая скорость и не снилась!
Михаил Парилов

Ответить

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

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

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