Le portail francophone de la géomatique

Bonjour,

Il y a trois étapes :
- La première est de retruscturer toutes tes couches, afin de permettre leur fusion. Le plus simple est de fusionner une copie de ces centaines de couches auxquelles auraient été soustraites les colonnes attributaires
- Une fois cela fait, on procède à la fusion
- Puis enveloppe convexe sur la fusion

Concernant l'étape de restructuration, on aurait pu espérer que le modeleur graphique de QGIS permette de traiter un dossier entier de shapefiles auxquels appliquer l'algo "refactoriser les champs" mais si le modeleur comprend l'option "exécuter par lot" pour toutes les couches ouvertes dans un projet, son mode de fonctionnement ne permet pas d'opérer facilement l'export automatique de couches retructurées aisément car il faut donner, pour chacune des couches, le chemin et le nom du fichier d'export.

Il faut donc programmer.

La ligne suivante, à exécuter dans une fenêtre de commande ou via un batch, si tu es sous windows, te permettra de ne retenir que la colonne géométique des couches. Un champ FID sera automatiquement attribué aux couches résultantes (avec le suffixe _geom)

for /R %F IN (*.shp) DO ogr2ogr -f "ESRI Shapefile" "%~nF_geom.shp" "%~nxF" -dialect sqlite -sql "select geometry from '%~nF'"

Avant de lancer cette commande, il faut te positionner dans le dossier comprenant les shapefiles

Super ! Je n'avais pas remarqué :-)

Par contre, je ne sais pas si le sextante modeler permet d'enlever automatiquement tous les champs des couches afin de ne conserver que la géométrie si celles-ci n'ont pas la même structure ?

Bonjour,

merci pour vos conseils.

L'outil "refactoriser les champs" pose effectivement problème car la structure de mes shapefiles ne sont pas les mêmes. Je n'arrive donc pas à obtenir via un batch des couches ayant une structure identique. Peut-être existe-t-il toutefois une "formule magique" à indiquer dans la colonne "correspondance de champs" de cet outil...je ne l'ai pas trouvé...

Même problématique pour l'algo "supprimer la colonne". Comme indiqué par Santanna, je n'arrive pas à effacer toutes les colonnes via un batch, ni à effacer toutes les colonnes sauf "ID"...

MathieuR, pardonne mon ignorance, mais je ne comprends pas où rentrer ta ligne de code : dans la commande QGIS ? dans la commande Windows ? j'ai tenté sans succès, peux-tu m'indiquer la marche à suivre stp ?

Merci beaucoup pour votre aide

ps : je confirme, en mode batch, l'auto remplissage des noms de fichiers en sorite fonctionne très bien

Merci MathieuR.

Je n'arrive pas à faire fonctionner le code. J'ai le message d'erreur "'ogr2ogr' n'est pas reconnu en tant que commande interne ou externe, un programme exécutable ou un fichier de commabdes".

J'ai réinstallé Python (3.6 32bits) ainsi que GDAL (2.1.4 cp32 win32). Ces installations ont bien fonctionné, mais le même message d'erreur s'affiche toujours. ogr2ogr is not defined...

C'est fait, plus de problème lié à ogr2ogr, mais nouveau message d'erreur pour chaque shp se trouvant dans mon dossier : "ERROR 1 : sql : unrecognised field name geometry"

Merci beaucoup MathieuR, c'est parfait !