[GRASS-user] Building topologically correct stream network

Jarek Jasiewicz jarekj at amu.edu.pl
Thu Jul 23 11:45:42 EDT 2009



Markus Metz pisze:
>
> Margherita Di Leo wrote:
>> Hi all
>> Hi Markus,
>> thanks for your answer. I found that the problem is in the topology 
>> of the river network. I think that v.strahler produces many '0' 
>> because it finds many little 'circles' in the stream vector. I tried 
>> to produce the stream network both with r.watershed -> r.to.vect and 
>> with r.strahler.sh -> r.to.vect: the results obtained in both ways 
>> are very similar. Is there a way to eliminate these circles, not by 
>> hand ?
> r.strahler.sh tries to if I understand the script right, but there 
> must be side effects: what were parallel streams with a ladder may be 
> merged into one stream.
>
> Maybe a new module to extract both raster and vector stream networks 
> from any flow accumulation map would indeed be useful.
Yes, r.to.vect will almost never give correct directional network. It 
require, as Marcus pointed different approach. r.flow ma be a good 
starting point to create new module
Jarek
>
> Markus M
>
>>
>> Markus Metz wrote:
>>>
>>> Margherita Di Leo wrote:
>>>> I'm sorry but the procedure that I proposed does not produce good 
>>>> results. In fact the ascii file produced by v.strahler has very 
>>>> strange orders, and also many 0 (that does not have sense in 
>>>> strahler order!). I think the problem is generated in finding the 
>>>> outlets 
>>> Outlet points are where stream is not null and drainage direction is 
>>> negative, e.g.
>>>
>>> r.mapcalc "outlet_points = if(!isnull(streams) && drainage < 0, 
>>> streams, null())"
>>>
>>> That would give you outlet points whose cell value is the stream 
>>> segment number leading to the outlet point.
>>>
>>>> and/or confluences.
>>> All nodes of the stream vector, excluding outlet points?
>>>
>>> Markus M
>>>
>>>> Waiting for the response of the List (please)..
>>>>
>>>> margherita
>>>>
>>>>
>>>> Jarek Jasiewicz wrote:
>>>>> Hi
>>>>>
>>>>> maybe this will help:
>>>>>
>>>>> http://gis.vsb.cz/GIS_Ostrava/GIS_Ova_2009/sbornik/Lists/Papers/060.pdf 
>>>>>
>>>>>
>>>>> it require another gis program but gives great results!
>>>>>
>>>>> Jarek
>>>>>
>>>>>
>>>>> Margherita Di Leo pisze:
>>>>>> Hi Stephen,
>>>>>> I actually am fighting against almost the same problem, so let's 
>>>>>> do it together :)
>>>>>> I found that r.strahler.sh works well to delineate streams and 
>>>>>> correct topology, (I use a treshold of 1km^2, it works well but 
>>>>>> don't ask me why :) ), but (in my little experience) does not 
>>>>>> work to order streams, so I must run v.strahler after r.strahler.sh.
>>>>>> I know it's a very small hint but hope this helps.
>>>>>>
>>>>>> Best
>>>>>> margherita
>>>>>>
>>>>>>
>>>>>> I have a low relief area in the costal plain of georgia that I am 
>>>>>> trying to produce a topologically correct stream network. There 
>>>>>> are many breaks in the flow accumulation grid and I would like to 
>>>>>> delineate the streams and "correct" the topology after r.thin and 
>>>>>> r.to.vect have been preformed so that I can then use v.strahler 
>>>>>> to order the streams. Thanks for any help and I will provide and 
>>>>>> data that would be helpful. thanks,
>>>>>>
>>>>>>
>>>>>> -- Stephen Sefick Let's not spend our time and resources thinking 
>>>>>> about things that are so little or so large that all they really 
>>>>>> do for us is puff us up and make us feel like gods. We are 
>>>>>> mammals, and have not exhausted the annoying little problems of 
>>>>>> being mammals. -K. Mullis
>>>>>>
>>>>>
>>>>
>>>>
>>>
>>
>>


More information about the grass-user mailing list