I'm reading the wiki [1] on the r.stream.basins module, and a question has come to mind.
I want to create the drainage areas for many specific pour points. The sample in the wiki page, using r.reclass on the streams raster, reads:
echo '42 = 42
252 = 252
188 = 188
* = NULL' >tmp
r.reclass input=streams output=sel_streams_1 <tmp1
r.stream.basins dir=dirs stream=streams basins=bas_basin_3
Should the last line not be:
r.stream.basins dir=dirs stream=sel_streams_1 basins=bas_basins_3
^^^^^^^^^^^^^^^
??
Also, for clarity, perhaps "tmp" should be "tmp3" in both cases above
Thanks,
Micha
[1]`` [http://grasswiki.osgeo.org/wiki/R.stream.*#r.stream.basins](http://grasswiki.osgeo.org/wiki/R.stream.*#r.stream.basins)
ยทยทยท
--
Micha Silver
GIS Consultant, Arava Development Co.
[http://www.surfaces.co.il](http://www.surfaces.co.il)
On Fri, Dec 27, 2013 at 10:47 AM, Micha Silver <micha@arava.co.il> wrote:
I'm reading the wiki [1] on the r.stream.basins module, and a question has
come to mind.
I want to create the drainage areas for many specific pour points. The
sample in the wiki page, using r.reclass on the streams raster, reads:
echo '42 = 42
252 = 252
188 = 188
* = NULL' >tmp
r.reclass input=streams output=sel_streams_1 <tmp1
r.stream.basins dir=dirs stream=streams basins=bas_basin_3
Should the last line not be:
r.stream.basins dir=dirs stream=sel_streams_1 basins=bas_basins_3
^^^^^^^^^^^^^^^
??
Also, for clarity, perhaps "tmp" should be "tmp3" in both cases above
I have attempted to fix that, please review
http://grasswiki.osgeo.org/grass-wiki/index.php?title=R.stream.*&diff=20069&oldid=17583
For further corrections you may update directly the Wiki since it is a Wiki 
Best
Markus