Растр в Растр с другой grid-сеткой средствами gdal/ogr

[nickleb]

Как перевести растр одной grid-сетки в растр с другой grid-сеткой средствами gdal/ogr? Помнится, что в GoldenSoftWare Surfer для таких целей использовал slice-процедуру или interpolate (получение значения Z в каждой точке X, Y новой сетки по существующему grid' у). А в gdal/ogr как это отскриптовать?

[Дмитрий Барышников]

Пояcните, что вы имеете в виду. У растра вроде нет координаты z. Или речь о dem?

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

[nickleb]

Дмитрий, к примеру: у меня файл глубин в формате ASC (ArcGIS grid-файл с 1-минутным разрешением по широте и долготе). У меня есть регулярная сетка в декартовых неизвестно какой картографической проекции, но географические координаты от автора этой проекции я все-же заполучил. Так, что даже "не растр-в-растр" - как я объявил, а в цикле, точка-за-точкой с географ. координатами - это будет универсальнее, получать значение глубины для каждой точки из описанного выше ASC'а . Какой-нибудь процедурой bilinear'ной интерполяции, как делал в свое время в Surfer'е.

[nickleb]

...наверное, для этого надо использовать gdalwarp...

[nickleb]

...или, что-то другое?..

[Дмитрий Барышников]

Алгоритмы перепроецирования работают несколько иначе - в результирующем растре берется точка и ищется какое значение соответствует ей в исходном. Разумеется, обычно точка в исходном точно не попадает в центр пикселя, поэтому используется интерполяция.
Да, если я правильно вас понимаю, вам должен помочь gdalwarp.
При этом ключ r указывает метод интерполяции. Как бы ASCII grid тоже поддерживается: http://gdal.org/frmt_various.html#AAIGrid

[nickleb]

Дмитрий, спасибо.