Le portail francophone de la géomatique

Bonjour à la communauté,

Je fais appel à vous car je suis bloquée (malgré mes recherches sur les forum en anglais) sur une action python qui semble pourtant toute bête... :

Je suis dans une première boucle (un shapefile de tous les Etats) et je fais une sélection dans ce shapefile des Etats pour récupérer l'identifiant correspondant à la 1ère entité de ce shapefile :

layer.startEditing()

for feat in layerEtats.getFeatures():
  id_etat = feat.attribute("ID")

  # une fois que j'ai récupéré cet état, je veux faire une sélection dans mon shapefile des unités administratives plus fines, ce que j'appelle "layer"
  expr1 = (u' "id_state" = \'')
  expr2 = id_etat
  expr3 = (u'\'')
  expr = str(expr1) + str(expr2) + str(expr3)
  selection = layer.getFeatures(QgsFeatureRequest().setFilterExpression(expr))
  ids = [i.id() for i in selection]
  layer.setSelectedFeatures(ids)

A ce stade, mes éléments (entités administratives plus fines dont l'Etat correspond au 1er dans la table attributaire du shapefile des Etats) se sélectionnent dans ma couche.

La deuxième étape est que je veux modifier les attributs de chacune de ces entités sélectionnées en itérant sur chacune d'entre elles....et c'est là que je bloque.

Je fais :

  for elem in layer.selectedFeatures():
    # calculs
    # blabla
    layer.changeAttributeValue(feat.id(),writefield,result)

Je dois donc avoir un problème sur la seconde itération (qui est pourtant à l'intérieur de ma 1ère boucle).

Le comportement de son script est étrange car il me modifie les valeurs des 70 premières unités administratives plus fines (soit exactement le nombre d'Etats...) et quand mon script bloque à la fin, j'ouvre ma table attributaire et ma sélection est pourtant bien faite.

--> le problème serait donc que le script arrive à sélectionner les entités mais il écrit sur tout le layer au lieu d'écrire uniquement sur les entités sélectionnées !


Si quelqu'un pouvait me débloquer...