Répondre au commentaire

Bonjour,

Je débute aussi et me sert de ces bouts de codes pour pour faire un script dans lequel je croise les profils en long des rivières avec la géologie.

J'ai moi aussi bloqué sur le passage des limites géologiques, ne réussissant pas à extraire les coordonnées x et y des points d'intersection.

Pour cela, j'ai transformé la couche ptintersect au format shapely en utilisant la commande suivante

ptintersec_shp = loads(ptintersec.ExportToWkb())

Il suffit ensuite de récupérer les coordonnées en faisant comme dans la partie précédente, c'est à dire :

xptint,yptint=ptintersec_shp.x,ptintersec_shp.y

Peut être que l'auteur a une solution plus simple.

J'y retourne, il faut maintenant que je trouve comment plotter ces intersections dans mon graph en convertissant les coordonnées xy vers la distance.