Shapely: opérations géométriques

objectif

Expliquer et montrer comment buffer, intersection, union et distance avec Shapely.

code minimal

from shapely.geometry import Point, Polygon
p = Point(0,0).buffer(1.0)
q = Polygon([(0,0),(2,0),(2,2),(0,2)])
inter = p.intersection(q)
float(inter.area), float(p.distance(q))

utilisation

# union et difference
u = p.union(q)
d = q.difference(p)
float(u.area), float(d.area)

variante(s) utile(s)

# simplification géométrique
ps = p.simplify(0.1, preserve_topology=True)
float(ps.area)

notes

• Shapely travaille en coordonnées cartésiennes; reprojetez en mètres si besoin.
• buffer rayon dépend de l’unité du CRS.