Le portail francophone de la géomatique

Bonjour à tous,

Avant de vous expliquez se que je veut faire et où je rencontre des problèmes. Je travaille sur QGIS 2.16.3, GRASS GIS 7 et SAGA 2.1.2. Je suis assez habitué à travailler sur ArcGis mais c'est la première fois que je me sers de QGIS.

Le besoin : créer une carte par commune sur l'ensemble d'une département avec l'IGN scan 25 en fond et un shapefile vecteur par dessus

Je dispose donc d'une carte IGN au format GPKG, un shapefile polygone des limites des communes du département et le shapefile vecteur à traiter

Mon idée : Vu qu'il y a plus d'une centaine de commune sur le département et que le shapefile vecteur est évolutif dans le temps, je me suis dit qu'il serait plus malin de faire un modèle. Je ne maitrise pas le langage python, du coup, je suis parti sur le modeleur graphique de QGIS (je me servais de son équivalent sur ArcGis).

J'avais préparer un image du brouillon de mon modèle pour faciliter la compréhension mais je n'arrive pas l'ajouter... Du coup, je vais essayer d'être clair dans mon explication d'où j'en suis dans le modèle.

Je n'ai aucun problème pour réaliser le traitement des vecteurs que j'ai effectué de la manière suivante :

1 - Une extraction par attribut dans la couche Commune pour extraire les limites de la commune appelé "Tartempion"
2 - Utilisation de l'outil Coupe pour découper le shapefile vecteur par rapport à "Tartempion" et appelé "Vector_Clip"
3 - Utilisation de l'outil "set style for vector layer" pour appliquer le style que je voulais à "Vector clip"

Jusque là ça tourne sans problème.

Au départ, j'avais pensé à seulement définir un masque sur l'ensemble de la commune "Tartempion" et ça marche très bien quand je le fais avec l'extension Mask de QGis mais cette option n'est pas disponible dans le modeleur... Du coup, j'ai essayé avec l'option l'outil mask de GRASS mais ça me dit que la couche n'a pas été générée correctement. Du coup, à paritr de là j'ai changé de solution.

Du moment, que j'ai commencé à essayer de découper la couche raster par rapport à "Tartempion", les seuls résultats que j'ai eu c'est plantage sur plantage de QGIS avec comme erreur : GPKG : bad application_id 0x47504B47.

Aprés avoir fait plein de tests qui ne pouvait pas fonctionner (et j'ai compris pourquoi). J'ai trouvé une nouvelle approche que je pensait fonctionnelle mais qui ne marche pas :

1 - utiliser "r.rgb" de GRASS pour découper mon raster en 3 rasters "red","blue","green"
2 - utiliser "clip raster with polygon" de SAGA pour découper les 3 rasters avec "Tartempion"
3 - utiliser "r.composite" de GRASS pour reformer mon raster RVB

Mais ça ne marche pas et toujours avec l'erreur ci-dessus. Je commence à me poser des questions sur le GPKG qui est un format que je ne connais pas trop et dont j'ai pas l'habitude de travailler avec (je travaille plutôt en TIFF). J'ai fait des recherches mais ce format est normalement bien supporté par les outils dont je me sers.

Voilà mes questions :
- Pensez vous que ma méthode d'approche puisse fonctionner pour avoir le résultat que j'escompte ?
- Est-il possible que le problème vient du format ? j'ai essayé de changer le format en TIFF mais j'ai eu la même que ci-dessus avec plantage d'ArcGis.
- Avez-vous déjà ce genre de problème et comment avez-vous fait pour le régler ?

D'habitude je préfère trouver la réponse par moi-même mais là j'avoue que je reste un peu dubitatif sur ce problème.

Merci d'avance pour votre aide et les pistes que vous me donnerez

Bonne journée.

Dernière modification par Sankreet (Mon 05 October 2020 13:23)