[GRASS-user] voronoi for polygons
Martin Rutzinger
Martin.Rutzinger at uibk.ac.at
Wed Aug 9 11:22:35 EDT 2006
hi martin!
i'm not sure what is wrong with v.select, but i can show a second approach how
to derive voronoi polygons (restriction: doesn't work for adjacent polygons):
- convert the selected polygons to raster (v.to.rast)
- switch null() values to 1 and pixels with values to null() (r.mapcalc)
- thin raster (r.thin) several iterations
- convert thinned raster to vector and make polygons (r.to.vect,...)
this works for me
/martin
********************
Hello,
I want to calculate voronoi areas for polygons, but I am stuck in
the "query-part".
using spearfish data I did:
v.extract input=fields at PERMANENT output=fields_extract list=41,36,46,21
v.to.points input=fields_extract output=fields_extract_points type=boundary
v.voronoi input=fields_extract_points output=fields_extract_points_voron
now I wanted to
- group all voronoi polygons which belong to points of one polygon (GIS query
overlap)
- calculate area of merged voronoi polygons and assign them to each original
polygon
but:
v.select operator=overlap gives:
v.select binput=fields_extract_points_voron_clean ainput=fields_extract_clean
atype=area btype=area output=test operator=overlap
Processing ainput areas ...
100%
Writing attributes ...
Layer 1
v.select complete.
Building topology ...
0 primitives registered
0 areas built
0 isles built
Attaching islands:
Attaching centroids: Topology was built.
Number of nodes : 0
Number of primitives: 0
Number of points : 0
Number of lines : 0
Number of boundaries: 0
Number of centroids : 0
Number of areas : 0
Number of isles : 0
anything wrong with the command?
I did it as well without b/atype=area but same result.
The problem how to calculate areas is not yet tackled ,-)
I appreciate any hint, Martin
________________________________________________________________________________
Mag. Martin Rutzinger
rutzinger at alps-gmbh.com
alpS - Natural Hazard Management
http://www.alps-gmbh.com
martin.rutzinger at uibk.ac.at
Institute for Geography, University of Innsbruck
http://geowww.uibk.ac.at
More information about the grass-user
mailing list