[GRASS-user] intersection of network-lines and areas
    achim 
    ak7 at jupiter.uni-freiburg.de
       
    Mon Mar  2 11:02:05 EST 2009
    
    
  
Thanks a lot!!
I will try this.
achim
Markus Metz schrieb:
> 
> achim wrote:
>> I attach an example-jpg to illustrate my point of view...
>>   
> OK, I think I understand.
> You could get the start points with r.mapcalc
> r.mapcalc "start_points = if(!isnull(stream_segments) && drainage < 0,
> 1, null())"
> 
> convert the start points to vector
> r.to.vect in=start_points out=start_points feature=point
> 
> get coordinates for each point with v.to.db map=start_points -p option=coor
> ideally pipe the output of v.to.db to something that generates a script
> to call r.water.outlet with each set of coords, otherwise it's manual
> editing
> 
    
    
More information about the grass-user
mailing list