Страница 1 из 1
Oперации с набором растров одновременно
Добавлено: 14 ноя 2017, 13:05
rhot
Есть ли возможность в GRASS проводить операции с целым набором растров одновременно? Типа функции raster:brick в R.
Re: Oперации с набором растров одновременно
Добавлено: 14 ноя 2017, 15:24
Александр Мурый
А как действует функция raster:brick в R? И какие именно операции интересуют?
Re: Oперации с набором растров одновременно
Добавлено: 14 ноя 2017, 15:38
Александр Мурый
Для работы с пространственно-временными наборами данных (если это вдруг оно) есть целый фреймворк:
https://grass.osgeo.org/grass70/manuals ... intro.html
Как раз там реализована работа с набором растров.
Re: Oперации с набором растров одновременно
Добавлено: 14 ноя 2017, 15:44
rhot
raster:brick действует так же, как и i.group
Набор растров - это данные Sentinel 2.
Интересует oберезание по границам другого растра и сохранение слоя по маске.
Re: Oперации с набором растров одновременно
Добавлено: 16 ноя 2017, 18:00
rhot
В общем, не нашел ничего.
Проблему решать циклами (R, Python). Например, сохранить все 10 каналов Sentinel по маске леса в R:
Код: Выделить всё
for (i in 1:10){
print(paste('Processing',i))
exp <- paste0('sent.',i,'=if(sentinel.',i,' && forest, sentinel.',i,', null())')
execGRASS('r.mapcalc', expression=exp, flags='overwrite')
}