[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