[GRASS-user] Python-Loop over points

Micha Silver micha at arava.co.il
Tue Feb 22 14:59:56 EST 2011

Hello Johannes:

On 02/22/2011 04:57 PM, stephen sefick wrote:

>> Am 22.02.2011 um 02:07 schrieb Daniel Victoria:
>>> I'm not sure but from what I read in the r.stream.basins description,
>>> you can have more than one point in your vector file in roder to
>>> delineate your basins. Now, in case you have one basin inside another,
>>> what I did once was:
>>> 1) get x,y coordinates of each point using v.out.ascii
>>> 2) for each point run r.watershed...
>> Yes, thats what I want as well, but I want to automatize it and not to run for each
>> point manually. I also want to run
>> for each point a grass-function (not necessarily r.stream etc.)
>> I just want to know how it basically works to loop over points.
>> to things I thought:
>> 1) loop over the ascii if that is possible

Here's an outline of a method I once used to get all points where a 
stream network crosses a road, then loop thru those points to create 
drainage basins for each of those points. Maybe you can distill from 
this what you need.

Micha Silver
Arava Development Co. +972-52-3665918

More information about the grass-user mailing list