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...
Cheers daniel On Fri, Feb 18, 2011 at 11:49 AM, Johannes Radinger <[email protected]> wrote: > Hello, > > > -------- Original-Nachricht -------- >> Datum: Mon, 7 Feb 2011 22:25:31 -0200 >> Von: Daniel Victoria <[email protected]> >> An: Johannes Radinger <[email protected]> >> CC: [email protected] >> Betreff: Re: [GRASS-user] Python-Loop over points > >> Hi Johannes, >> >> You will need to use for loop. Something like this (UNTESTED!) >> >> for p in range(5): >> points_map = "point"+str(p) >> out_map = "output"+str(p) >> grass.run_command("r.streams.basins", dir = "flow_direction, >> points=points_map, basins = out_map) >> >> That should work, or at least something similar >> >> Daniel > > Thank you for your suggestions Daniel, I also thought about > a for loop. But in your case the loop is using for each iteration > a input-pointmap called "point"+str(p). > > The difference in my > case is that I've got only one vector map but with several points > in it and I want to iterate over the single points. > How is it possible to iterate over the single points? The points are > stored in a kind of list format, so it should somehow be possible to > loop over this file, but how? > > any suggestions? > > /johannes > >> >> On Wed, Feb 2, 2011 at 9:39 AM, Johannes Radinger <[email protected]> >> wrote: >> > Hello, >> > >> > I am performing e.g. the r.stream.basins calcutlation in python: >> > >> > grass.run_command("r.stream.basins", >> > dir = "flow_direction", >> > points = "point1", >> > basins = "output") >> > >> > that's working perfectly in the case of one single point in "point1". >> > How can I do that multiple times when the point layer consist of e.g. 5 >> > points and I want to create an output for each point (create 5 output >> rasters). >> > >> > How can that be done in a python script? >> > >> > >> > thank you >> > >> > Johannes >> > >> > -- >> > NEU: FreePhone - kostenlos mobil telefonieren und surfen! >> > Jetzt informieren: http://www.gmx.net/de/go/freephone >> > _______________________________________________ >> > grass-user mailing list >> > [email protected] >> > http://lists.osgeo.org/mailman/listinfo/grass-user >> > > > -- > NEU: FreePhone - kostenlos mobil telefonieren und surfen! > Jetzt informieren: http://www.gmx.net/de/go/freephone > _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
