<div dir="ltr">Dear Moritz,<div><br></div><div style>Thanks a lot for your quick answer and sorry for my own "latency". Yes, indeed, the Medial Axis Transform (or Skeleton) of Polygons dataset can be computed via the Voronoi diagram of the corresponding set of points (discretizing the polygons boundaries). However, I assume that this is not a scalable solution with a huge set of polygons (insofar as the v.to.points resolution has to be of enough fine grain). I was thinking to a morphological thinning based solution... Maybe I'm wrong.</div>
<div style><br></div><div style>With GRASS GIS 6.4.0+42329, and the enclosed (small) dataset, the tessel.shp output file does not match expectations. Could you please tell me where I'm wrong?</div><div style><br></div>
<div style><div><span class="" style="white-space:pre">   </span>g.remove vect=buildings,newbuildings,vor,tessel</div><div><span class="" style="white-space:pre">    </span></div><div><span class="" style="white-space:pre">   </span>v.in.ogr dsn=/tmp/src/src.shp output=buildings</div>
<div><span class="" style="white-space:pre">    </span>g.region vect=buildings</div><div><br></div><div><span class="" style="white-space:pre">   </span>v.category input=buildings output=newbuildings option=add type=boundary</div>
<div><span class="" style="white-space:pre">    </span>v.out.ogr input=newbuildings type=boundary dsn=/tmp/src olayer=build layer=1 format=ESRI_Shapefile</div><div><br></div><div><span class="" style="white-space:pre">        </span>v.to.points -i --overwrite input=newbuildings output=newbuildings_pts dmax=0.1</div>
<div><span class="" style="white-space:pre">    </span>v.out.ogr input=newbuildings_pts type=point dsn=/tmp/src olayer=pts layer=1 format=ESRI_Shapefile</div><div><span class="" style="white-space:pre">  </span></div><div><span class="" style="white-space:pre">   </span>v.voronoi -l --overwrite input=newbuildings_pts output=vor</div>
<div><span class="" style="white-space:pre">    </span>v.out.ogr input=vor type=line dsn=/tmp/src olayer=vor layer=1 format=ESRI_Shapefile</div><div><br></div><div><span class="" style="white-space:pre">       </span>v.dissolve input=vor output=tessel</div>
<div><br></div><div><span class="" style="white-space:pre">   </span>v.out.ogr input=tessel type=line dsn=/tmp/src olayer=tessel layer=1 format=ESRI_Shapefile</div><div><br></div><div style>Thanks a lot once again for your help.</div>
<div style><br></div><div style>Thomas Leduc </div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Sep 2, 2013 at 3:13 PM, Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On 29/08/13 17:38, Thomas Leduc wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Let's consider a set of input polygons (represented by dark gray<br>
polygonal footprints in the enclosed screenshot [1]). I wonder whether<br>
there already exists a (set of) GRASS instructions to produce the<br>
"surrounding tessellation". I mean, could you please tell me how to<br>
compute the set of multicolor polygons in the enclosed screenshot?<br>
<br>
As an example, the blue polygon is the set of points that are:<br>
   - overlayed by the dark gray polygon labeled #1<br>
   - or closer to the dark gray polygon labeled #1 than to any other<br>
input polygon.<br>
<br>
Thanks a lot for any help/advice.<br>
</blockquote>
<br></div>
You could try the following workflow:<br>
<br>
- add category values to your boundaries (if that is important make them equal to the cat values of your centroids, otherwise just using v.category op=add type=boundary)<br>
- use v.to.points to transform the boundaries to points (each point will have the cat value of the corresponding boundary)<br>
- v.voronoi to create voronoir diagrams for these points<br>
- v.dissolve to disolve boundaries between polyons with the same category numbers<br>
<br>
Moritz<br>
______________________________<u></u>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/grass-user</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Thomas LEDUC
</div>