<div dir="ltr">Hi,<div><br></div><div>sorry for late answer, I had overlooked this email<br><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Mar 7, 2014 at 7:30 AM, Hamish <span dir="ltr"><<a href="mailto:hamish_b@yahoo.com" target="_blank">hamish_b@yahoo.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div class="">maning wrote:<br>
<br>
> I'm testing r.hazard.flood and noticed that it computes the flood and<br>
> mti layers on the full region of the elevation raster instead of the<br>
> pre-defined region settings.<br>
><br>
> The relvant code I found from the r.hazard.flood is this:<br>
>    # Detect cellsize of the DEM<br>
>     info_region = grass.read_command('g.region', flags = 'p',<br>
> rast =<br>
> '%s' % (r_elevation))<br>
>     dict_region = grass.parse_key_val(info_region, ':')<br>
>     resolution = (float(dict_region['nsres']) +<br>
> float(dict_region['ewres']))/2<br>
>     grass.message("Cellsize : %s " % resolution)<br>
><br>
> Would it be possible to either respect the current region settings or<br>
> add a flag to choose between the dem region settings or current region<br>
> settings?<br>
<br>
<br>
</div>if it needs to detect the original raster cell resolution (usually that's only needed to avoid aliasing artifacts in certain situations) it should use <a href="http://r.info" target="_blank">r.info</a> to get the answer. I'm guessing due to the averaging of the ew and ns resolutions avoiding aliasing isn't the case.<br>


<br>
if a module wants to change the region (and almost none should ever do that except for g.region by itself) it should set up a temporary WIND_OVERRIDE first, in the case of python scripting there is an easy grass python function to make that and clean it up at the end. (otherwise parallel jobs get their regions messed up mid-run, and region changes without you expecting it will)<br>


<br>
I suspect grass.raster's raster_info() is the better way for r.hazard.flood to do what it's trying to do now. But also just querying the current g.region info without changing anything is probably even better, as that is what the other raster commands will expect to use. If the user should align perfectly with the input map first, it should be noted in the help page for them to do it manually before running the module.<br>

</blockquote><div><br></div><div>Thanks for suggestion, applied in r59879 </div><div><br></div><div>Cheers,</div><div>madi</div></div><div><br></div>-- <br><div dir="ltr"><div><font color="#666666">Best regards,</font></div>

<div><font color="#666666"><br></font></div><div><font color="#666666">Dr. Margherita DI LEO    </font></div><div><span style="color:rgb(102,102,102);font-family:arial;font-size:small">Scientific / technical project officer</span><br>

</div><div><font color="#666666"><br></font></div><div><font color="#666666">European Commission - DG JRC </font></div><div><font color="#666666">Institute for Environment and Sustainability (IES)</font></div><div><font color="#666666">Via Fermi, 2749</font></div>

<div><font color="#666666">I-21027 Ispra (VA) - Italy - TP 261</font></div><div><font color="#666666">       </font></div><div><font color="#666666">Tel. +39 0332 78 3600   </font></div><div><font color="#666666"><a href="mailto:margherita.di-leo@jrc.ec.europa.eu" target="_blank">margherita.di-leo@jrc.ec.europa.eu</a></font></div>

<div><font color="#666666"><br></font></div><div><font color="#666666">Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission.</font></div>

</div>
</div></div></div>