Dapatkan grafik ganda dari QGIS?

Saya memiliki file .shp, dan saya perlu mendapatkan grafik ganda dari file ini (artinya, setiap poligon centroid harus menjadi simpul, dan dua simpul dihubungkan jika poligon yang terkait berdekatan).

Saya berhasil mendapatkan centroid dari poligon file saya.

Bagaimana saya bisa menambahkan garis antara simpul yang memang harus dihubungkan?

Anda dapat menggunakan hanya Python (versi 2.7.x dan 3.x) tanpa QGIS:

1) Dengan Fiona , rupawan dan 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) dengan GeoPandas , rupawan dan 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) Saya kira Anda dapat melakukan hal yang sama dengan PyQgis dan itertools (Python 2.7 untuk versi 2.x dan 3.x untuk versi 3.x)