Страница 1 из 1

"Загрубить" растр

Добавлено: 04 июн 2018, 14:27
nickleb
Здравствуйте. Надобно (как всегда "быстренько-до-конца-дня") сделать из исходного растра (ИР) новый растр (НР) с размером_ячейки_НР = N * размер _ячейки_ИР
В современном Python, в R через какие библиотеки (и их функции) это сейчас грамотно делают
уважающие науки о Земле и свой авторитет настоящие ГИС-специалисты?
ИР с ячеёй 0.5 градусов lon/lat
(с привязкой
+proj=longlat +a=6371229 +b=6371229 +no_defs)
в НР c ячеёй 2.5 градуса lon/lat

... В R (library ("raster")) - функции aggregatе/disaggregate - это то, что надо в нашем случае - иль нет?

Re: "Загрубить" растр

Добавлено: 04 июн 2018, 14:56
Эдуард Казаков
gdal_warp и параметры xRes, yRes

Например, в Python, загрубляем в два раза:

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

raster = gdal.Open('input_dataset.tif')
input_xRes = raster.GetGeoTransform[1]
input_yRes = raster.GetGeoTransform[4]
N = 2
new_raster = gdal.Warp('new_dataset.tif',raster, format = 'GeoTIFF', xRes = input_xRes*N, yRes = input_yRes*N,)

Re: "Загрубить" растр

Добавлено: 04 июн 2018, 14:58
nickleb
Эдуард Казаков писал(а):
04 июн 2018, 14:56
gdal_warp и параметры xRes, yRes

Например, в Python, загрубляем в два раза:
Эдуард, спасибо большое!

Re: "Загрубить" растр

Добавлено: 04 июн 2018, 15:17
rhot
nickleb писал(а):
04 июн 2018, 14:27
В R (library ("raster")) - функции aggregatе/disaggregate - это то, что надо в нашем случае - иль нет?
Ну да. Не работает?

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

r.new <- raster::aggregate(r.old, fact=N, fun=mean)

Re: "Загрубить" растр

Добавлено: 04 июн 2018, 15:52
nickleb
rhot писал(а):
04 июн 2018, 15:17
Ну да. Не работает?
КОД: ВЫДЕЛИТЬ ВСЁ

r.new <- raster::aggregate(r.old, fact=N, fun=mean)
rhot, то, что доктоR :) прописал! ехал-google'ил-вспоминал... сейчас всё в "реалити" прошло... спасибо!