Hi all!<br><br>As you know landsat data has separate files for different bands. Creating a subset is a wise option to work on region of our interest. Therefore, I used this method to subset my image - <br><div style="margin-left: 40px;">

<span>d.zoom || g.region
</span><br><span>g.region -p
</span><br><span>r.mapcalc <b class="highlight">subset</b>=original
</span><br></div>as suggested in <a href="http://osgeo-org.1803224.n2.nabble.com/Raster-Map-subset-td4523291.html#a4526819">http://osgeo-org.1803224.n2.nabble.com/Raster-Map-subset-td4523291.html#a4526819</a>.<br><br>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.<br>

<br>Thanks and regards,<br><br>Chethan S.<br>