Форумы GIS-Lab.info

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

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

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

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

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

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

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

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