[GRASS-user] Python-Loop over points

[Micha Silver]

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.
http://www.surfaces.co.il/?p=241



-- 
Micha Silver
Arava Development Co. +972-52-3665918
http://surfaces.co.il