Fusion de couches

Un cas d'école : un collaborateur récupère un grand nombre de fichiers de données et souhaite les exploiter proprement. Il s'agit de plusieurs centaines de fichiers cadastraux allemands.

Téléchargez, dézippez et ouvrez la couche Koln gml.zip.

Le but est d'automatiser :

1. L'ouverture d'un grand nombre de fichiers GML
2. Leur conversion en shapes
3. Leur édition de structure
4. Enfin, leur fusion en un seul shape.

import os
import shutil
from os import listdir
from os.path import isfile, join
 
# On supprime les futures variables utilisées, car certaines semblent gêner les exécutions multiples du script
my_var = ['ShapesOuvertsNoms', 'LayersShapesOuverts', 's', 'n', 'mon_shape', 'repertoireShapes', 'pathShape', 'nom', 'writer']
for v in my_var:
if v in globals():
del v
 
monRepertoire = 'C:/Users/Georges/Downloads/Koln/lines/'
 
# On liste les GML en triant d'éventuels fichiers non GML (GFS par exemple)
mesGML_ = [f for f in listdir(monRepertoire) if isfile(join(monRepertoire, f))]
mesGML = [val for val in mesGML_ if val.endswith('.gml')]
 
QgsProject.instance().removeAllMapLayers()
iface.mapCanvas().refresh()
 
NbreGmlValid = 0
NbreGmlNoValid = 0
 
# On va stocker les noms des GML correctement ouverts mais aussi leurs appels en layers
GmlOuvertsNoms = []
LayersGmlOuverts = []
 
# Boucle pour ouvrir tous les GML sur QGIS
for gml in mesGML :
    mon_gml = QgsVectorLayer(monRepertoire+gml, gml, 'ogr')
 
    # On précise la projection
    crs = mon_gml.crs()
    crs.createFromId(25832)
    mon_gml.setCrs(crs)
    QgsProject.instance().addMapLayer(mon_gml)
 
    # On compte et liste les GML correctement ouverts mais aussi leurs appels en layers
    if mon_gml.isValid():
        print(gml, 'est valide ! OKLM.')
        NbreGmlValid += 1
        GmlOuvertsNoms.append(gml)
        LayersGmlOuverts.append(mon_gml)
 
        # On compte les GML en erreur
        if not mon_gml.isValid():
        print(gml, 'est invalide ! Nous le supprimons.')
        QgsProject.instance().removeMapLayer(mon_gml.id())
        NbreGmlNoValid += 1
 
print('Nombre de GML correctement ouverts :',NbreGmlValid)
print('Nombre de GML invalides et supprimés :',NbreGmlNoValid)
 
# EXPORT EN SHAPE
repertoireShapes = 'shapes'
pathShape = os.path.join(monRepertoire, repertoireShapes)
 
# On va stocker les noms des shapes correctement ouverts mais aussi leurs appels en layers
ShapesOuvertsNoms = []
LayersShapesOuverts = []
 
# Création d'un répertoire pour accueillir nos shapes
if not os.path.exists(pathShape):
    os.makedirs(pathShape)
else:
    shutil.rmtree(pathShape)
    os.makedirs(pathShape)
 
# Boucle pour exporter nos GML en shape
for s,n in zip(LayersGmlOuverts, GmlOuvertsNoms):
    nom = n.replace('.gml', '')
    writer = QgsVectorFileWriter.writeAsVectorFormat(s, pathShape+'/' + nom + '.shp', 'utf-8', driverName='ESRI Shapefile', onlySelected=False)
 
    mon_shape = QgsVectorLayer(pathShape+'/' + nom + '.shp', nom, 'ogr')
 
    ShapesOuvertsNoms.append(nom)
    LayersShapesOuverts.append(mon_shape)
 
# OUVRIR LES SHAPES
QgsProject.instance().removeAllMapLayers()
for x in LayersShapesOuverts:
    QgsProject.instance().addMapLayer(x)
 
# On va supprimer tous les champs de tous nos shapes sauf 2 champs
# Lister les champs d'un layer si besoin et tester avec my_vectorlayer = iface.activeLayer()
# D'abord on crée la liste de champs à supprimer
for my_shape in LayersShapesOuverts:
    prov = my_shape.dataProvider()
    field_names = [field.name() for field in prov.fields()]
    listeChampSuppression = []
 
    for a in field_names:
        listeChampSuppression.append(a)
 
    listeChampSuppression.remove('gml_id')
    listeChampSuppression.remove('measuredHe')
 
    # Puis on supprime la liste de champs à supprimer
    for c in listeChampSuppression:
        with edit(my_shape):
            my_field = my_shape.fields().indexFromName(c)
            my_shape.dataProvider().deleteAttributes([my_field])
        my_shape.updateFields()
 
# Fusion des shapes
parameters = {'LAYERS': ShapesOuvertsNoms,'CRS': 'EPSG:25832','OUTPUT': monRepertoire+'results.shp'}
processing.run("qgis:mergevectorlayers", parameters)
maFusion = QgsVectorLayer(monRepertoire+'results.shp', 'Ma jolie fusion', 'ogr')
QgsProject.instance().addMapLayer(maFusion)

Liens ou pièces jointes

chemin_de_fer.zip	[ ]	0 Ko
data_BDTOPO_V3_Dep05_adresse.zip	[ ]	3889 Ko
data_IRIS_2019.zip	[ ]	45905 Ko
decathlon_france.zip	[308 magasins Décathlon français depuis OSM le 27 décembre 2020]	11 Ko
eau.zip	[ ]	0 Ko
glaciers.zip	[ ]	231 Ko
iso_iris.zip	[Des zones isochrones à 15 minutes autour de 308 POIs.]	12125 Ko
Koln gml.zip	[ ]	2818 Ko
peaks.zip	[ ]	14 Ko
peaks_selection.zip	[ ]	1 Ko
simple_countries.zip	[ ]	1880 Ko
sol.zip	[ ]	0 Ko
troncons_routes.zip	[ ]	0 Ko