Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I have certain hubs (coloured yellow) in a district (coloured light blue). I need to find polygons indicating the closest hub to every point in the district.

The voronoi polygons that I obtained using QGIS are coloured brown in the image. However, I need the voronoi polygons to fill the outer boundary as I need the polygons to fill the entire district.

enter image description here

Adding a buffer region while making the voronoi polygons and then using the intersect function to crop the polygons to the district boundary gives me wrong polygons as shown by the following image:

enter image description here

How do I generate the voronoi polygons that I require?

share|improve this question
1  
Unless they have fixed it, qgis doesn't create the proper voronoi if buffered. see the link here gis.stackexchange.com/questions/11866/… –  Dan Patterson Sep 4 '13 at 16:43
    
have you tried the rdeldir tool in GME? spatialecology.com/gme/rdeldir.htm –  mwil Sep 4 '13 at 19:33
2  
@Dan is right: every bit of the lower figure is wrong (the line segments should bisect the points at their sides but most of them don't). A workaround is to place four artificial points far beyond the extent of the polygon (at corners of a large bounding box) and include them in the Voronoi polygon calculation, then intersect the result with the polygon. –  whuber Sep 4 '13 at 20:13
    
If you have ArcMap, the implementation I provided will do this, otherwise, include some "points at infinity" as Bill suggests –  Dan Patterson Sep 4 '13 at 21:20
add comment

1 Answer

I think I may have a solution to your problem. Through the Sextante plugin you can access modules from a program called SAGA GIS. In the Points module is a tool called Thiessen Polygons. This is another name for Voronoi polygons. This tool should work for you. I just tried it out on a set of rain gauges and get the exact same coverage area that ArcGIS and grass (v.voronoi) give me. So I feel pretty confident in the output at least. I hope this helps.

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.