Le portail francophone de la géomatique

Bonjour,
L'erreur indique un mauvais paramètre et comme les paramètres attendus sont essentiellement des couches, c'est que l'une d'elles ne l'est pas.
Perso, ce que je fais souvent, c'est d'exécuter mes codes dans la console python de QGIS pas à pas pour voir ce qui est renvoyé.
Et dans le code proposé, layer1 ne renvoie pas une couche à la sortie mais une entité (la dernière parcourue dans la boucle). print ne fait qu'imprimer/afficher les évènements, il ne sélectionne pas les entités.
Essaie en remplaçant

#on demande de selectionner  dans la table que les adresses ayant dans leur type d'occupation le mot "Logement"
for f in layer1.getFeatures():
        print f['TYPE_D_OCC' is '%Logement%']

# on recupere la selection (variable f) dans une autre couche layer1
layer1=f

par

layer1.setSelectedFeatures([]) # pour s'assurer d'avoir une sélection vide au départ
for feature in layer1.getFeatures(): #  chaque entité de la couche 
    if 'Logement' in feature['TYPE_D_OCC'] :  # dont le champ occupation contient 'Logement'
        layer1.select(feature.id())   # est ajoutée à la sélection

Ce code  sélectionne directement dans la couche les entités dont le type d'occupation contient 'Logement'. Ça devrait marcher.
plus d'infos par ici

Dernière modification par SANTANNA (Thu 11 June 2015 14:23)

Bonjour,
As-tu essayé d'exécuter le nouveau code?
Par défaut, les algorithmes de Traitement, sauf erreur de ma part, fonctionnent sur les entités sélectionnées dans la couche et s'il n'y a pas de sélection, sur toute la couche. il suffit donc de sélectionner les entités qui t'intéressent (ce que fait le bout de code proposé) et appliquer le même code d'intersection entre sur les couches layer1 et layer2. Pas besoin de récupérer ailleurs rien que les entités sélectionnées.

Bonjour,

en intégrant juste les modifs proposées, il a marché ou pas? Et si ça n'a pas marché, qu'est-ce qui est renvoyé comme erreur ou résultat?

Ensuite, pourquoi avoir changé aussi "radicalement" le code? L'objectif a atteindre a-t-il changé? Parce qu'avec les propositions de modification apportées à ton premier code, j'obtiens de mon coté le résultat attendu. Pourquoi donc ce choix de repartir sur d'autres fonctions?
Enfin, quelques règles sous Python :
- les "import ..." et "from ... import ..." doivent préférentiellement être au début de ton code; et il n'est pas forcément nécessaire d'importer toutes les bibliothèques possibles (je pense que dans le lot, il y en a qqs uns qui ne sont pas nécessaires pour ton analyse)
- si ce bout de code ne fait pas partie d'un ensemble de traitement avec d'autres fonctions, il n'est pas vraiment nécessaire de définir une fonction. Mais ta fonction définie ne sélectionne pas non plus toutes les entités "Logement"...
- Si tu veux récupérer le résultat d'une fonction, il te faut utiliser un "return" dans ta fonction

ok, c'est parce que lors du parcours de la table, certaines entités n'ont pas de type d'occupation de renseigné (valeur NULL) donc comparaison impossible. La boucle "bug" donc s'arrête. On va donc ignorer les cas de valeur nulle. En décomposé, ça donnerait quelque chose comme ceci...

layer1.setSelectedFeatures([]) # pour s'assurer d'avoir une sélection vide au départ
for feature in layer1.getFeatures(): #  chaque entité de la couche 
    if feature['TYPE_D_OCC'] <> NULL: dont le champ occupation est non nulle 
        if 'Logement' in feature['TYPE_D_OCC'] :  # et contient le texte 'Logement'
            layer1.select(feature.id())   # est ajoutée à la sélection

Bonjour, dans la continuité de mon sujet je reviens vers vous avec un autre problème (encore plus difficile).
Vu que les géométries des mes deux couches ne s'intersecte pas; l'idée est de faire un BUFFER à partir de la même sélection (type d'occupation=logement) et ensuite:

-Récupérer dans un fichier 1 (csv ou excel si possible) les entités intersectées par le buffer et dont ID_RIL=ID_buffer

-Récupérer dans un fichier 2 (csv ou excel si possible) les entités intersectées par le buffer et dont ID_RIL DIFFERENT ID_buffer

-Récupérer dans un fichier 3 (csv ou excel si possible) les entités qui ne sont PAS intersectées et dont ID_RIL EXISTE dans ID_buffer

-Récupérer dans un fichier 4 (csv ou excel si possible) les entités qui ne sont NI intersectées et dont dont ID_RIL N'EXISTE PAS dans ID_buffer.

J'arrive pas à récupérer ma selection (type_occupation_logement) dans une couche temporaire.
et lorsque j’exécute le script dans la console python, la couche intersection en sortie est vide.

Voici mon début de code (j'avoue qu'il n'est pas vraiment terrible):

# -*-coding: utf8 -*-
import os  #importer la librairie os
import sys #importer la librairie du systeme (sys)
from osgeo import ogr # importer la librairie ogr
import processing #importer l'ensemble des algorithmes
from qgis.core import *
from qgis.gui import * 
#from qgis.analysis import * #importe les deux libraires  d'analyse de couche et de geometry
#on charge la couche adresse à jour
layer1 = iface.addVectorLayer("H:\SIG\Bases adresses\Shapefiles/adresse20150201.shp", "adresse2015", "ogr")
layer1.setSelectedFeatures([]) # pour s'assurer d'avoir une sélection vide au départ
for feature in layer1.getFeatures(): #  chaque entité de la couche 
    if feature ['TYPE_D_OCC']<> NULL:
        if 'Logement' in feature['TYPE_D_OCC'] :  # dont le champ occupation contient 'Logement'
            layer1.select(feature.id())   # est ajoutée à la sélection
from qgis.analysis import QgsGeometryAnalyzer
QgsGeometryAnalyzer().buffer(layer1, "H:\SIG\Bases adresses\Shapefiles/output.shp", 0.08, False, False, -1)
layer2 = iface.addVectorLayer("H:\SIG\Bases adresses\Shapefiles/output.shp", "buffer", "ogr")# on ajoute le buffer
#on charge la couche RIL avec laquelle on va faire l'intersection
layer3= iface.addVectorLayer("H:\SIG\Bases adresses\Shapefiles/ADRESSE_RIL_2015.shp", "ADRESSE_RIL_2015", "ogr")
#from qgis.analysis import * #importe les deux libraires  d'analyse de couche et de geometry
from qgis.analysis import QgsOverlayAnalyzer       
overlayAnalyzer = QgsOverlayAnalyzer()
overlayAnalyzer.intersection(layer2, layer3, r"H:\SIG\Bases adresses\Shapefiles\intersection.shp") # false pour qu'il nus retourne que les entités selectionnées
layer3= iface.addVectorLayer("H:\SIG\Bases adresses\Shapefiles/intersection.shp", "intersection", "ogr")