[GRASS-user] Region definition
luisapena1979 at gmail.com
Fri Mar 25 10:57:52 EDT 2011
I have a raster map with a size like 25000 x 30000 with valid values on all
its extent. With another raster, I defined a sub-extent of this raster, much
smaller where only a small area is valid and the rest it NULL() like this:
r.mapcalc output = if(regional_area at PERMANENT>0,2004 at Reg,null())
In this case, regional_area is a smaller area and I want to select only a
portion of 2004 at Reg that falls inside regional_area.
So I obtained output as a raster map with only some values (4000x3000)
If I define g.region rast=output I still obtain a 25000 x30000 region which
does not make sense since I only have a small valid area.
How can I define a region, using g.region for only output valid pixels?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the grass-user