Obtenir un double graphique de QGIS?

J'ai un fichier .shp, et j'ai besoin d'obtenir le double graphique de ce fichier (c'est-à-dire que chaque polygone centroïde devrait devenir un sommet, et deux sommets sont liés si leurs polygones associés sont adjacents).

J'ai réussi à obtenir les centroïdes des polygones de mon fichier.

Comment puis-je ajouter des lignes entre les sommets qui devraient en effet être connectés?

Vous pouvez utiliser Python uniquement (versions 2.7.x et 3.x) sans QGIS:

1) Avec Fiona , bien fait et 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) avec GeoPandas , galbées et 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) Je suppose que vous pouvez faire la même chose avec PyQgis et itertools (Python 2.7 pour les versions 2.x et 3.x pour les versions 3.x)