Hi Rich,<br><br><div class="gmail_quote">On Sat, Apr 7, 2012 at 12:23 AM, Rich Shepard <span dir="ltr"><<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I have several questions about r.basin that I hope can be quickly<br>
answered.<br>
<br>
1.) Is there an example that can be added to the manual page? I would like<br>
to better understand the prefix and threshold parameters.<br></blockquote><div><br></div><div>r.basin produces several maps for every run, which are flow direction, accumulation, Horton ordered streams etc.. Prefix parameter is simply a string given by the user in order to distinguish all the maps produced by a single run of the program, i.e. every set of coords for outlet. For example, </div>
<div><br></div><div><a href="http://r.basin.py">r.basin.py</a> map=elevation@PERMANENT prefix=out easting=636654.791181 northing=218824.126649 threshold=20</div><div><br></div><div>will produce output maps named:</div><div>
<br></div><div># raster</div><div><div>out_elevation_accumulation </div><div>out_elevation_aspect </div><div>out_elevation_dist2out </div><div>out_elevation_distance </div>
<div>out_elevation_drainage </div><div>out_elevation_hack</div><div>out_elevation_hillslope_distance</div><div>out_elevation_horton</div><div>out_elevation_shreve</div><div>out_elevation_slope</div><div>out_elevation_strahler</div>
<div><br></div><div># vector</div><div>out_elevation_basin </div><div>out_elevation_mainchannel </div><div>out_elevation_network </div><div>out_elevation_outlet </div></div><div><br></div><div>Threshold parameter is the same of r.watershed. 'Autothreshold ' flag' simply uses a threshold area of 1 km^2 and is intended as a tentative.</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
2.) There are several module dependencies, but r.stream.basin, <a href="http://r.wf.py" target="_blank">r.wf.py</a>,<br>
and <a href="http://r.ipso.py" target="_blank">r.ipso.py</a> cannot be downloaded and installed using g.extension; they're<br>
not found in the addon repository.<br></blockquote><div><br></div><div>For what concerns r.stream.basin, it is a typo in the documentation, correct is r.stream.basins, thanks for pointing me out, I just updated it in the description file.</div>
<div><a href="http://r.wf.py">r.wf.py</a> and <a href="http://r.ipso.py">r.ipso.py</a> are in the addOns repository:</div><div><br></div><div><a href="http://grass.osgeo.org/wiki/GRASS_AddOns#r.wf">http://grass.osgeo.org/wiki/GRASS_AddOns#r.wf</a></div>
<div><a href="http://grass.osgeo.org/wiki/GRASS_AddOns#r.ipso">http://grass.osgeo.org/wiki/GRASS_AddOns#r.ipso</a></div><div><br></div><div>I can't see any obvious error in the Makefile..... maybe an experienced dev could please check why they are not working with g.extension.</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
3.) The r.stream.extract module is not working for my projects' data using<br>
version 6.5svn on linux. I have detailed vector maps of the stream networks<br>
in each of the subbasins I need to describe and model. If I run v.to.rast on<br>
these maps can I use the outlet cell as the values for the easting and<br>
northing parameters to the command?<br></blockquote><div><br></div><div>Please note that r.basin is intended to work with only one outlet at time. This means you should script it in a for cycle if you run it for several outlets. You can't give to r.basin the coordinates directly from a vector map.</div>
<div>I have indicated to run r.stream.extract before r.basin because most of the times the calculated stream network doesn't match with the natural one, so that your coordinates should be adjusted to match the calculated stream network in order to get a result from r.basin. By the way, in a first run you can give a try with your coordinates.</div>
<div><br></div><div>HTH,</div><div>madi</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Thanks,<br>
<br>
Rich<br>
<br>
______________________________<u></u>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/grass-user</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Ing. Margherita Di Leo, Ph.D.<br>