"Загрубить" растр
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
"Загрубить" растр
Здравствуйте. Надобно (как всегда "быстренько-до-конца-дня") сделать из исходного растра (ИР) новый растр (НР) с размером_ячейки_НР = 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 - это то, что надо в нашем случае - иль нет?
В современном 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: "Загрубить" растр
gdal_warp и параметры xRes, yRes
Например, в Python, загрубляем в два раза:
Например, в 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,)
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: "Загрубить" растр
Эдуард, спасибо большое!Эдуард Казаков писал(а): ↑04 июн 2018, 14:56gdal_warp и параметры xRes, yRes
Например, в Python, загрубляем в два раза:
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: "Загрубить" растр
Ну да. Не работает?
Код: Выделить всё
r.new <- raster::aggregate(r.old, fact=N, fun=mean)
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
-
- Гуру
- Сообщения: 964
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 34 гостя