Python >> Tutoriel Python >  >> Python Tag >> NumPy

convertir une liste de géométrie galbée en tableau numpy

J'ai créé un fichier de formes avec vos coordonnées de points et le code suivant produit ces tableaux numpy :

import fiona
from shapely.geometry import shape
import numpy as np

path = '/home/zeito/pyqgis_data/points_list.shp' #your points

points = fiona.open(path)

geoms = [ shape(feat["geometry"]) for feat in points ]

list_arrays = [ np.array((geom.xy[0][0], geom.xy[1][0])) for geom in geoms ]

for array in list_arrays:
    print array

Le résultat après l'avoir exécuté est :

[ 656822.07964268 -185003.72434373]
[ 656917.75456579 -184985.66727045]
[ 656997.7888964  -185001.57868412]