¿Polígonos Voronoi que se agotan hasta el infinito?

17

He estado trabajando en la implementación del script voronoi.py (traducción de Bill Simon del código de Steve Fortune). De acuerdo con el enlace QGIS es básicamente el mismo código con el que estoy trabajando, pero la forma de generar los polígonos puede extenderse hasta el infinito difiere entre los dos. En el caso de QGIS, supongo que algún rectángulo delimitador se genera en cierta medida más allá de la extensión del punto del sitio. En mi caso, genero puntos de sitio falsos como un círculo que rodea los puntos de interés del sitio. Los puntos del círculo tienen un radio el doble del diámetro de la extensión de los puntos del sitio. Los resultados entre los dos se muestran en la figura adjunta, la implementación de QGIS a la izquierda y mi implementación a la derecha. .

Comosepuedever,laimplementacióndeQGISproducepolígonosextrañosquenotienenpuntosdesitio(veaelsitio473comoejemploenambasimplementaciones).Estoesobvioqueestámal,perotampocoestoysegurodesimiimplementaciónescorrecta.ComonotengoaccesoaunalicenciadeArcInfoparaArcMap,nopuedocompararloqueArcMapgeneradeformanativaynopuedoverlamaneradeadjuntarelshapefiledepuntosdelsitioparaquealguienlopruebe.Asíquesialguientienealgúncomentariosobremiusodeun"círculo infinito" para ayudar a forzar polígonos de borde en lugar de un "rectángulo infinito", lo apreciaría. Mejor aún, si alguien está interesado en probar la salida en alguna otra plataforma GIS, puedo proporcionar el shapefile de los puntos del sitio que usé.

Parece que parte del problema es que los polígonos más allá de los puntos del sitio se superponen, por lo tanto, cuando se simbolizan sin un relleno, aparecen como polígonos distintos. La siguiente figura muestra lo que sucede cuando se selecciona el polígono sobre el sitio 473 y 415. Así que las preguntas originales permanecen, pero hay una nueva arruga.

Es un error (o una "falta de definición"), los resultados correctos en QGIS si se utiliza un búfer de 0%. Lo que aún queda por resolver es si un patrón denso de "puntos infinitos" circulares dará el mismo resultado que los "bordes infinitos" rectangulares.

    
pregunta PolyGeo 06.07.2011 - 14:15

1 respuesta

9

Pude confirmar que mi implementación se compara al menos con ArcGIS (pero requiere una licencia de ArcInfo). La implementación de QGIS es errónea si se utiliza una distancia de búfer. Para aquellos que usan ArcGIS con una licencia de ArcView, el conjunto de herramientas para Voronoi y / o Delaunay se puede encontrar como Herramientas de triangulación .

    
respondido por el user681 12.07.2011 - 19:39

Lea otras preguntas en las etiquetas