Las self-intersection en polígonos pueden ser determinadas en PyQGIS con el método "saga:polygonselfintersection" que se encuentra en la Processing Tool Box. Es equivalente a la 'unary_union' del módulo python shapely.ops que se encuentra descrito aquí. Para probar el procedimiento se usó el shapefile de la imagen siguiente:
El código empleado se encuentra a continuación:
1 2 3 4 5 6 7 | import processing polys = QgsVectorLayer("/home/zeito/pyqgis_data/layer.shp","intersection","ogr") output = "/home/zeito/pyqgis_data/intersection.shp" processing.runAndLoadResults("saga:polygonselfintersection", {'POLYGONS':polys, 'INTERSECT':output}) |
Después de ejecutado el código anterior en la Python Console de QGIS, el resultado visualizado en la imagen siguiente es tal como se esperaba.
No hay comentarios:
Publicar un comentario