By the way. Instead of being setting everytime g.region xxxx -p can I save a region definition and use it everytime? If so, how do I load region definition in a script?<div><br></div><div>Thanks</div><div><br>Best regards</div>
<div>Pedro<br><br><div class="gmail_quote">On Fri, Feb 5, 2010 at 6:50 PM, Glynn Clements <span dir="ltr">&lt;<a href="mailto:glynn@gclements.plus.com">glynn@gclements.plus.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
Pedro Roma wrote:<br>
<br>
&gt; Ok but for instance. If my region is defined with a Spatial Resolution of<br>
&gt; 1kilometer and, if I do some r.mapcalc with two LANDSAT images, the output<br>
&gt; image will have which one of the resolutions? the Region or the 2 landsat<br>
&gt; images?<br>
<br>
</div>The region.<br>
<br>
Most raster modules (including r.mapcalc) generate maps with the<br>
bounds and resolution of the current region. Input maps are cropped,<br>
padded and/or resampled to the current region, so all processing<br>
within the module takes place on the grid defined by the current<br>
region.<br>
<br>
The main exceptions are the r.in.* modules, which import files<br>
cell-for-cell, using the bounds specified in the file&#39;s metadata if<br>
present. The current region is ignored.<br>
<br>
Modules which perform their own resampling (e.g. r.resamp.*, r.proj)<br>
create output maps according to the current region but read their<br>
input maps without cropping or resampling.<br>
<font color="#888888"><br>
--<br>
</font><div><div></div><div class="h5">Glynn Clements &lt;<a href="mailto:glynn@gclements.plus.com">glynn@gclements.plus.com</a>&gt;<br>
</div></div></blockquote></div><br></div>