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

Вопросы по нескольким пакетам сразу, или вопросы, которые непонятно к какой ГИС отнести
Ответить
nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

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

Сообщение nickleb » 04 июн 2018, 14:27

Здравствуйте. Надобно (как всегда "быстренько-до-конца-дня") сделать из исходного растра (ИР) новый растр (НР) с размером_ячейки_НР = 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 - это то, что надо в нашем случае - иль нет?

Аватара пользователя
Эдуард Казаков
Гуру
Сообщения: 546
Зарегистрирован: 23 апр 2014, 17:11
Репутация: 532
Откуда: Planet Earth
Контактная информация:

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,)

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

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

Сообщение nickleb » 04 июн 2018, 14:58

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

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

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

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

Сообщение rhot » 04 июн 2018, 15:17

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

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

r.new <- raster::aggregate(r.old, fact=N, fun=mean)
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

nickleb
Гуру
Сообщения: 964
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

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

Сообщение nickleb » 04 июн 2018, 15:52

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

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

Ответить

Вернуться в «Общий - ПО»

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

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