Erhalten Sie Dual Graph von QGIS?

Ich habe eine .shp-Datei, und ich muss das duale Diagramm aus dieser Datei abrufen (dh jeder Polygonschwerpunkt sollte ein Scheitelpunkt werden, und zwei Scheitelpunkte werden verknüpft, wenn die zugehörigen Polygone benachbart sind).

Ich habe es geschafft, die Zentroide aus den Polygonen meiner Akte zu erhalten.

Wie kann ich Linien zwischen den Eckpunkten hinzufügen, die tatsächlich verbunden werden sollen?

Sie können Python nur (Versionen 2.7.x und 3.x) ohne QGIS verwenden:

1) Mit Fiona , formschönen und itertools

import fiona
from shapely.geometry import MuliPolygon, LineString, mapping
Multi = MultiPolygon([shape(poly['geometry']) for poly in fiona.open("polygons.shp")])
# creation of the dual graph shapefile
import itertools
# schema of the dual graph shapefile
schema = {'geometry': 'LineString','properties': {'test': 'int'}}
with fiona.open('dual_graph.shp','w','ESRI Shapefile', schema) as e:
   for poly1,poly2 in  itertools.combinations(Multi, 2):
       if poly1.touches(poly2):
           e.write({'geometry':mapping(LineString([poly1.centroid, poly2.centroid])), 'properties':{'test':1}})

2) mit GeoPandas , formschönen und itertools

import geopandas as gpd 
polygon = gpd.GeoDataFrame.from_file("polygons.shp")
geoms = polygon['geometry'].tolist()
dual = gpd.GeoDataFrame(gpd.GeoSeries([LineString([poly1.centroid, poly2.centroid]) for poly1,poly2 in  itertools.combinations(geoms, 2) if poly1.touches(poly2)]),columns=['geometry'])
dual.to_file("dual_graph2.shp")

3) Ich nehme an, Sie können dasselbe mit PyQgis und itertools tun (Python 2.7 für 2.x-Versionen und 3.x für 3.x-Versionen).