Greetings<div><br></div><div>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:</div>
<div>r.mapcalc output = if(regional_area@PERMANENT&gt;0,2004@Reg,null())</div><div>In this case, regional_area is a smaller area and I want to select only a portion of 2004@Reg that falls inside regional_area.</div><div>So I obtained output as a raster map with only some values (4000x3000)</div>
<div>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.</div><div>How can I define a region, using g.region for only output valid pixels?</div>
<div><br></div><div>Thanks</div><div>Luisa</div>