<div dir="ltr">Hi Kurt!<br>I work/study Hydrology and I use GRASS in my work... so, for my short experinece, I can tell you:<br>the &quot;stream map&quot;, generated by r.watershed (or r.watershed.fast) is the map of the reaches, so the colors you can see and the numbers (categories) associated with each reach (you can see it if you query the map) is simply an &quot;id&quot; of the reach (it is useful to distinguish the reaches...). <br>

r.watershed extracts these stream segments from the DEM basing his analysis on the elevation data... so, depending on the &quot;threshold&quot; value you give him in input, the module will return a different result (so it happens that if you give a threshold value greater you will see in your map only principal reaches, or if you give a lesser threshold value you will obtain more and more reaches - more dense network -). <br>

<br>A way to isolate some reaches -remainig in raster format- could be using r.reclass:<br><br>type &quot;r.reclass&quot; in the shell: it requires you the name of the map to reclassify, the name of a new map (reclassified) and the values of the categories to change (e.g. I want category 30 to become category 1... so I type: 30 = 1). <br>

<br>To do this you can query the map for the reaches you interest and read the actual category (e.g. 30,32,34) then <span style="background-color: rgb(255, 255, 51);">type</span>:<br><span style="background-color: rgb(255, 255, 51);">&gt; r.reclass</span><br>

name of the new map: OPTION:&nbsp;&nbsp; Raster map to be reclassified<br>&nbsp;&nbsp;&nbsp;&nbsp; key: input<br>&nbsp; format: name<br>required: YES<br><br>Enter the name of an existing raster file<br>Enter &#39;list&#39; for a list of existing raster files<br>

Hit RETURN to cancel request<br><span style="background-color: rgb(255, 255, 51);">&gt; streams</span>&nbsp;&nbsp;&nbsp; #from r.watershed<br>&lt;streams&gt;<br><br>OPTION:&nbsp;&nbsp; Name for output raster map<br>&nbsp;&nbsp;&nbsp;&nbsp; key: output<br>&nbsp; format: name<br>

required: YES<br><br>Enter a new raster file name<br>Enter &#39;list&#39; for a list of existing raster files<br>Hit RETURN to cancel request<br><span style="background-color: rgb(255, 255, 51);">&gt; WhatYouWant</span>&nbsp;&nbsp; #the name of the new map I&#39;m going to create<br>

&lt;WhatYouWant&gt;<br><br>OPTION:&nbsp;&nbsp; File containing reclass rules<br>&nbsp;&nbsp;&nbsp;&nbsp; key: rules<br>&nbsp; format: name<br>required: NO<br><br>Enter the name of an existing input file<br>Hit RETURN to cancel request<br><span style="background-color: rgb(255, 255, 51);">&gt;&nbsp;<span style="background-color: rgb(255, 255, 255);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #I have not a file containig the rules</span></span>...<br>

&lt;&gt;<br><br>OPTION:&nbsp;&nbsp; Title for the resulting raster map<br>&nbsp;&nbsp;&nbsp;&nbsp; key: title<br>required: NO<br><span style="background-color: rgb(255, 255, 51);">enter option &gt;&nbsp; </span><span style="background-color: rgb(255, 255, 51);"><span style="background-color: rgb(255, 255, 255);">#I choose not to give a title to the map</span></span> I&#39;m going to create<br>

Enter rule(s), &quot;end&quot; when done, &quot;help&quot; if you need it<br>Data range is 2 to 66<br><span style="background-color: rgb(255, 255, 51);">&gt; 30 = 1<span style="background-color: rgb(255, 255, 255);">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span></span>#...I prefer entering the rules by myself<br style="background-color: rgb(255, 255, 51);">

<span style="background-color: rgb(255, 255, 51);">&gt; 32 = 2</span><br style="background-color: rgb(255, 255, 51);"><span style="background-color: rgb(255, 255, 51);">&gt; 34 = 3</span><br style="background-color: rgb(255, 255, 51);">

<span style="background-color: rgb(255, 255, 51);">&gt; end</span><br><br>so you can open your &quot;WhatYouWant&quot; raster map where you will see only the reaches you selected by category.<br><br>Hoping I&#39;ve been helpful...<br>

<br>Annalisa<br><br><br><br><br><div class="gmail_quote">2008/10/14 Kurt Springs <span dir="ltr">&lt;<a href="mailto:ferret_bard@mac.com" target="_blank">ferret_bard@mac.com</a>&gt;</span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

Thanks to Hamish and Markus Metz,<br>
<br>
What I am doing is comparing two or three (possibly four) types of megaliths in my study area and trying to understand the relation ship between both contemporary monuments and later monuments. &nbsp;While I was looking at old GIS maps, I was looking at the stream and basin maps and had overlaid the monuments on these maps.<br>


<br>
First off, the vector river and stream maps, representing the current water ways are a bit different then the stream map. &nbsp;The physical streams do seem to follow the stream maps, but there are many more branches on the stream map than there are actual streams.<br>


<br>
Some of the megalithic tombs types seem to fall into the web of the streams, while others don&#39;t. &nbsp;I think this is important. &nbsp;I am, unfortunately, an archaeologist, not a hydrologist, so I am not up on watershed lingo.<br>


<br>
What I need to know is what information on these maps may tell me something important. &nbsp;For example, are the stream systems shown on the watershed stream maps just important of them selves, or might the values the colors represent tell me something? &nbsp;Is there a way to isolate them?<br>


<br>
I guess what I am hoping for is a hydrologist who uses grass and can explain what I am looking at in the real world.<br><font color="#888888">
<br>
Kurt Springs_______________________________________________<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/mailman/listinfo/grass-user</a><br>
</font></blockquote></div><br></div>