[GRASS-user] intersection of network-lines and areas

Markus Metz markus.metz.giswork at googlemail.com
Mon Mar 2 10:59:47 EST 2009

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 

More information about the grass-user mailing list