<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
On 07/29/2010 06:38 PM, David Townshend wrote:
<blockquote
cite="mid:AANLkTinGJaGddVuwHwLvGvb-7Vr62hjK6qyKuR72=jc2@mail.gmail.com"
type="cite">
<pre wrap="">Hi
I am trying the generate watersheds for several outlets by looping
over r.water.outlet in a python script (attached). The basin procedure
is that for each outlet point (read using v.info -t), my script runs
r.water.outlet, r.to.vect and v.patch, to combine the watersheds into
a simgle vector map. I am connecting the watersheds to the outlet
points by setting the category value in the output from r.to.vect,
i.e. before patching. The trouble is that some of the watersheds
overlap, e.g. where an outlet lies inside another watershed. v.patch
doesn't seem to handle this too well, and the result is that I have a
vector map with a few areas, a few boundaries, a few centroids and a
lot of errors.
Is there a better way of approaching this? And how do I allow overlapping areas?
</pre>
</blockquote>
No, since GRASS enforces topology, there's no such thing as overlapping
polygons. If you run r.water.outlet twice in the same stream network,
you will get one basin that's a sub-catchment of a larger basin. If you
try to patch these two vector maps, assuming the patch succeeds, and
there are no topology errors, you'll get the smaller basin as one
polygon, and the larger basin with the smaller one "cut out" as the
second polygon.<br>
<br>
I've encountered a similar problem to yours: many small overlaps or
gaps along the boundary between two adjacent basins created by two runs
of r.water.outlet. <br>
I don't know of any easy solution, other than using the options
included with "v.clean" to remove small areas and snap boundaries that
are within a given threshold. THen after v.clean, there will often
still be a few errors that need to be ironed out manually with v.digit.<br>
<br>
HTH,<br>
-- <br>
Micha<br>
<blockquote
cite="mid:AANLkTinGJaGddVuwHwLvGvb-7Vr62hjK6qyKuR72=jc2@mail.gmail.com"
type="cite">
<pre wrap="">Regards
David
This mail was received via Mail-SeCure System.
</pre>
<pre wrap=""><fieldset class="mimeAttachmentHeader"></fieldset>
_______________________________________________
grass-user mailing list
<a class="moz-txt-link-abbreviated"
href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>
<a class="moz-txt-link-freetext"
href="http://lists.osgeo.org/mailman/listinfo/grass-user">http://lists.osgeo.org/mailman/listinfo/grass-user</a>
This mail was received via Mail-SeCure System.
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Micha Silver
Arava Development Co. +972-52-3665918
<a class="moz-txt-link-freetext" href="http://surfaces.co.il">http://surfaces.co.il</a>
</pre>
</body>
</html>