Растр в Растр с другой grid-сеткой средствами gdal/ogr
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Растр в Растр с другой grid-сеткой средствами gdal/ogr
Как перевести растр одной grid-сетки в растр с другой grid-сеткой средствами gdal/ogr? Помнится, что в GoldenSoftWare Surfer для таких целей использовал slice-процедуру или interpolate (получение значения Z в каждой точке X, Y новой сетки по существующему grid' у). А в gdal/ogr как это отскриптовать?
- Дмитрий Барышников
- Гуру
- Сообщения: 2572
- Зарегистрирован: 17 ноя 2009, 19:17
- Репутация: 261
- Откуда: Москва
Re: Растр в Растр с другой grid-сеткой средствами gdal/ogr
Пояcните, что вы имеете в виду. У растра вроде нет координаты z. Или речь о dem?
[ Сообщение с мобильного устройства ]
[ Сообщение с мобильного устройства ]
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Растр в Растр с другой grid-сеткой средствами gdal/ogr
Дмитрий, к примеру: у меня файл глубин в формате ASC (ArcGIS grid-файл с 1-минутным разрешением по широте и долготе). У меня есть регулярная сетка в декартовых неизвестно какой картографической проекции, но географические координаты от автора этой проекции я все-же заполучил. Так, что даже "не растр-в-растр" - как я объявил, а в цикле, точка-за-точкой с географ. координатами - это будет универсальнее, получать значение глубины для каждой точки из описанного выше ASC'а . Какой-нибудь процедурой bilinear'ной интерполяции, как делал в свое время в Surfer'е.
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Растр в Растр с другой grid-сеткой средствами gdal/ogr
...наверное, для этого надо использовать gdalwarp...
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Растр в Растр с другой grid-сеткой средствами gdal/ogr
...или, что-то другое?..
- Дмитрий Барышников
- Гуру
- Сообщения: 2572
- Зарегистрирован: 17 ноя 2009, 19:17
- Репутация: 261
- Откуда: Москва
Re: Растр в Растр с другой grid-сеткой средствами gdal/ogr
Алгоритмы перепроецирования работают несколько иначе - в результирующем растре берется точка и ищется какое значение соответствует ей в исходном. Разумеется, обычно точка в исходном точно не попадает в центр пикселя, поэтому используется интерполяция.
Да, если я правильно вас понимаю, вам должен помочь gdalwarp.
При этом ключ r указывает метод интерполяции. Как бы ASCII grid тоже поддерживается: http://gdal.org/frmt_various.html#AAIGrid
Да, если я правильно вас понимаю, вам должен помочь gdalwarp.
При этом ключ r указывает метод интерполяции. Как бы ASCII grid тоже поддерживается: http://gdal.org/frmt_various.html#AAIGrid
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: Растр в Растр с другой grid-сеткой средствами gdal/ogr
Дмитрий, спасибо.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 14 гостей