Oперации с набором растров одновременно

[rhot]

Есть ли возможность в GRASS проводить операции с целым набором растров одновременно? Типа функции raster:brick в R.

[Александр Мурый]

А как действует функция raster:brick в R? И какие именно операции интересуют?

[Александр Мурый]

Для работы с пространственно-временными наборами данных (если это вдруг оно) есть целый фреймворк: https://grass.osgeo.org/grass70/manuals ... intro.html
Как раз там реализована работа с набором растров.

[rhot]

raster:brick действует так же, как и i.group

Набор растров - это данные Sentinel 2.
Интересует oберезание по границам другого растра и сохранение слоя по маске.

[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')
}