But with this I can subset one band data at a time. Is it somehow possible
to use i.group and region definition to obtain subsets of all bands?
Scripting might be an answer but I have no much idea about it.
what about working with different computational regions?
http://grass.osgeo.org/grass64/manuals/html64_user/rasterintro.html :
"[...]
As a general rule in GRASS:
* Raster output maps have their bounds and resolution equal to those
of the current computational region.
* Raster input maps are automatically cropped/padded and rescaled
(using nearest-neighbour resampling) to match the current region.
* Raster input maps are automatically masked if a raster map named
MASK exists. The MASK is only applied when reading maps
from the disk.
[...]"
so you can do calculations/run modules in different computational region without
extracting of many subsets of your data before.
best regards
Helmut
___________________________________________________________
WEB.DE DSL Doppel-Flat ab 19,99 €/mtl.! Jetzt mit
gratis Handy-Flat! http://produkte.web.de/go/DSL_Doppel_Flatrate/2