possibilités supplémentaires

#Création de géométries:
wkt = "POLYGON((....))"
essai = Testpoly(nom="test3", geom=WKTSpatialElement(wkt))
wkt2 = "POLYGON((....))"
essai2 = Testpoly(nom="test4", geom=WKTSpatialElement(wkt))
session.add_all([essai1,essai2])
session.commit()
# et les analyses spatiales:
session.query(Testpoly).filter(Testpoly.geom.touches(uneautregéométrie)).count()
# + buffer, intersection, convex hull..;