Message posté par : Olivier Pompier (pompierolivier(a)hotmail.com) ---------------------------------------- Merci JD pour votre réponse. Effectivement, cela a diminué le temps de traitements significativement. Le code fonctionnel ci-après: ----------------- Code : index = QgsSpatialIndex(features) feats = [feat for feat in layer.getFeatures()] for feat in feats: for id in index.intersects(feat.geometry().boundingBox()): if feat.id() > id: intersection = feat.geometry().intersection(feats[id].geometry()) elt = QgsFeature(fields) elt.setAttribute('nb_sup',compteur) elt.setGeometry(intersection) ----------------- ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=330175#p330175 Pour y répondre : qgis_fr(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net