21 Aug
2020
21 Aug
'20
noon
Message posté par : Alexpiration (alexisqgis(a)gmail.com)
----------------------------------------
Bonjour à tous !
Plusieurs jours que je planche sur un script Python me permettant de réaliser une
opération qui, avec mes connaissance actuelles du logiciel QGIS, me parait compliqué !
Explication de mon but :
J'ai 2 sources de données :
- La première est une couche de centroïdes, dont chaque centroïdes est associés à une
commune de l'Ile de France. Dans la table attributaire de cette couche se trouve le
code INSEE associé à la commune ainsi que son nom
- La deuxieme est une base de donnée au format SQLITE contenant des informations
concernant divers projets en cours sur des communes. Dans cette base de donnée, on y
retrouve le nom du projet, la commune dans laquelle se trouve le projet, le code insee
associé à la commune, et divers autres informations.
J'ai d'abord pensé à créer une simple jointure pour faire en sorte de générer un
centroïde correspondant à chaque projet associé à sa commune. Cependant, ma table
attributaire des centroïdes ne gerent pas les projets qui ont une commune en commun,
uniquement un seul projet est rattaché et par conséquent, des projets de ma base de donnée
ne sont pas pris en compte.
J'ai donc voulu débloquer cette situation avec mes connaissances en Python et avec
PyQGIS.
La démarche de mon script est la suivante :
Faire une boucle sur le champ Code Insee de ma table des centroïdes et dans cette même
boucle, faire une autre boucle sur le champ Code Insee de ma base de donnée.
Avec une structure conditionnelle, s'il y a correspondance de Code Insee, je récupère
la géométrie du centroïdes correspondant , je récupère également les différents attributs
correspondant dans ma base de donnée, et j'ajoute le tous dans une nouvelle couche
Shapefile de points.
Cependant je n'arrive à comprendre ou debugger mon script sur la partie structure
conditionnelle. A la fin de mon script, j'ai bien ma nouvelle couche Shape avec les
champs de ma base de donnée, mais aucune entrée n'est enregistrée dedans, je suppose
donc que mon programme ne trouve aucune correspondance de code insee, ce qui n'est pas
le cas.
Voici mon script :
-----------------
Code :
# Create new temporary layer
new_layer = QgsVectorLayer("Point", "Suivi_Topo", "memory")
# Get 2 layers that want
center_layer = QgsProject.instance().mapLayersByName('centroide_commune_idf')
center_layer = center_layer[0]
suivi_projet_db = QgsProject.instance().mapLayersByName('suivi_projet')
suivi_projet_db = suivi_topo_db[0]
# Get attributes of each layers
center_layer_attr = center_layer.getFeatures()
suivi_topo_db_attr = suivi_topo_db.getFeatures()
# Add Fields of db in the new layer
suivi_topo_field = suivi_topo_db.fields()
new_layer_co = new_layer.dataProvider()
new_layer_co.addAttributes(suivi_topo_field)
new_layer.updateFields()
# Loop in each features, if same attrs, add features
for center_feature in center_layer_attr:
# get geometry of the feature
geom_center = center_feature.geometry()
for db_feature in suivi_topo_db_attr:
if db_feature.attribute('insee') ==
center_feature.attribute('code_insee'):
db_feature_lst = db_feature.attributes()
new_feature = QgsFeature()
new_feature.setGeometry(geom_center)
new_feature.setAttributes(db_feature_lst)
new_layer_co.addFeature(new_feature)
new_layer.updateExtents()
QgsProject.instance().addMapLayer(new_layer)
-----------------
En vous remerciant tous d'avance pour votre contribution :)
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=334932#p334932
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