Obtenir des polygones voronoi qui s'étendent jusqu'à la limite d'un autre polygone


10

J'ai certains hubs (de couleur jaune) dans un quartier (de couleur bleu clair). J'ai besoin de trouver des polygones indiquant la plaque tournante la plus proche de chaque point du district.

Les polygones voronoi que j'ai obtenus en utilisant QGIS sont colorés en marron dans l'image. Cependant, j'ai besoin des polygones voronoi pour remplir la frontière extérieure comme j'ai besoin des polygones pour remplir tout le district.

entrez la description de l'image ici

L'ajout d'une région tampon tout en créant les polygones voronoi, puis en utilisant la fonction d'intersection pour recadrer les polygones à la limite du district me donne des polygones incorrects, comme le montre l'image suivante:

entrez la description de l'image ici

Comment générer les polygones voronoi dont j'ai besoin?


1
Sauf s'ils l'ont corrigé, qgis ne crée pas le bon voronoi s'il est tamponné. voir le lien ici gis.stackexchange.com/questions/11866/…

avez-vous essayé l'outil rdeldir dans GME? spatialecology.com/gme/rdeldir.htm
mwil

2
@Dan a raison: chaque bit du chiffre inférieur est faux (les segments de ligne doivent couper les points sur leurs côtés, mais la plupart ne le font pas). Une solution de contournement consiste à placer quatre points artificiels bien au-delà de l'étendue du polygone (aux coins d'une grande boîte englobante) et à les inclure dans le calcul du polygone de Voronoï, puis à intersecter le résultat avec le polygone.
whuber

Si vous avez ArcMap, l'implémentation que j'ai fournie le fera, sinon, incluez des "points à l'infini" comme le suggère Bill

Réponses:


4

Je pense que j'ai peut-être une solution à votre problème. Grâce au plugin Sextante, vous pouvez accéder aux modules à partir d'un programme appelé SAGA GIS. Dans le module Points se trouve un outil appelé Thiessen Polygons. Ceci est un autre nom pour les polygones de Voronoi. Cet outil devrait fonctionner pour vous. Je viens de l'essayer sur un ensemble de pluviomètres et d'obtenir exactement la même zone de couverture que ArcGIS et l'herbe (v.voronoi) me donnent. Je me sens donc assez confiant dans la sortie au moins. J'espère que ça aide.


Je n'ai pas trouvé l'outil Thiessen Polygons dans le module Points de la boîte à outils de géotraitement QGIS. Cependant, j'ai trouvé dominoc925.blogspot.com.au/2013/07/… décrit comment le faire en utilisant SAGA GIS directement, ce qui fonctionne très bien.
AndrewHarvey

0

En résumé des commentaires précédents ...

Ajoutez des points à l'infini, puis utilisez la fonction de coupe. Cette vidéo le montre assez clairement:

https://www.youtube.com/watch?v=eos1nR-5Tio

J'espère que cela aide quelqu'un!


La procédure expliquée dans la vidéo est incorrecte. Comme le montre la deuxième image de la question, l'ajout d'une zone tampon puis l'utilisation de l'intersection donnent des polygones inexacts.
shreyasgm
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.