Hi Jens,
Sorry, my fault. I should have asked you about your platform first.
Assuming you are on Windows I guess python would be most suitable language for
scripting. Unfortunately I am not that good at python, but I am sure others
will fill in.
You might need something like:
points = grass.read_command("v.to.db", map='YOUR_VECTOR_POINTS", columns='',
option='coor', separator=',', flags='p')
for p in points.splitlines():
x=p.split(',')[2]
y= p.split(',')[3]
grass.run_comand('r.stream.basin', dirs=..., coors="%d,%d" % (x, y),
...)
Good luck!
Stefan
-----Original Message-----
From: Jens Hegg [mailto:jensenh...@gmail.com]
Sent: 10. februar 2015 21:22
To: Blumentrath, Stefan
Cc: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] r.stream.basins vector outlets help (Jens Hegg)
Sefan,
Thanks for the help. I'm still not getting it to run, but I am not at all
familiar with anything beyond running basic commands on the command console,
I've mostly just used the GUI. Should these be input into the command console
or the Python shell? I've tried both, the console seems to not recognize most
of it, the Python shell gives me lots of syntax errors (which seems like a good
start, I just don't know Python at all).
Here is what I have:
points=$(v.db.select -c map=Fish_and_Water_PointsEPSG4140 column=cat)) for p in
$points do
v.extract --overwrite input=Fish_and_Water_PointsEPSG4140
output=temp_point_map cat=(p) type=point
r.stream.basins direction=HenrysDrainage@HenrysFork
points=temp_point_map basins=Basin_(p) done
The first error is a syntax error at the dollar sign in points=$(v.db.select...
Am I getting close? The idea is that temp_point_map gets overwritten at each
loop so all that remains is the last one. I'm wanting to name each new basin by
the category so I can go back and figure out it's attribute data.
Jens
On Feb 9, 2015, at 11:04 PM, Blumentrath, Stefan wrote:
> Hei Jens,
>
> If you want to have watersheds which overlap, you have to have them in
> different maps as neither raster nor vector data with topology (as in GRASS
> GIS) allow for overlapping areas.
> That again means you would have to run r.stream.basins for your vector points
> individually, e.g. in a loop like this:
>
> points=$(v.db.select -c map=mypoints column=cat) for p in $points do
> v.extract ...
> r.stream.basisns ...
> done
>
> Cheers
> Stefan
>
> -----Original Message-----
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Jens Hegg
> Sent: 10. februar 2015 05:32
> To: grass-user@lists.osgeo.org
> Subject: Re: [GRASS-user] r.stream.basins vector outlets help (Jens
> Hegg)
>
> Ok, I figured out how to make it work for the vector file. I was including
> both the flow direction and the stream network along with the vector file. It
> turns out that only the direction and vector outlets files are needed and
> then it runs perfectly. The manual examples really should be updated to
> include a vector points example.
>
> I do still have one problem that someone might be able to help with. Many of
> the points I use for outlets are just downstream of one another so it creates
> watersheds that are just small slivers between points. Is it possible for it
> to create overlapping watersheds as you move down the stream, or do I need to
> go back and do this manually?
>
> Jens
>
>
>
> On Feb 9, 2015, at 6:08 PM, grass-user-requ...@lists.osgeo.org wrote:
>
>> Send grass-user mailing list submissions to
>> grass-user@lists.osgeo.org
>>
>> To subscribe or unsubscribe via the World Wide Web, visit
>> http://lists.osgeo.org/mailman/listinfo/grass-user
>> or, via email, send a message with subject or body 'help' to
>> grass-user-requ...@lists.osgeo.org
>>
>> You can reach the person managing the list at
>> grass-user-ow...@lists.osgeo.org
>>
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of grass-user digest..."
>>
>>
>> Today's Topics:
>>
>> 1. Import/Export of ESRI Arc ascii grids (Thomas Adams) 2. GRASS
>> GIS 7.0.0 RC2 released (Markus Neteler) 3. Re: Import/Export of ESRI
>> Arc ascii grids (Markus Neteler) 4. Re: Import/Export of ESRI Arc
>> ascii grids
>> (C?sar Augusto Ram?rez Franco)
>> 5. Re: v.generalize: does it take forever? (F?bio Dias) 6.
>> r.stream.basins vector outlets help (Jens Hegg)
>>
>>
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
