J'ai un fichier raster sur l'utilisation des terres dans le nord de l'Italie (cultures / zones urbaines, etc.)
J'ai défini des zones tampons autour de certains points donnés d'une couche vectorielle et je voudrais calculer le% de ces zones occupées par chaque valeur du fichier raster afin de faire une estimation sur les utilisations moyennes des zones sélectionnées (comme 70% cultures, 30% urbain, etc.).
Je sais que dans ArcGIS il y a une fonction Spatial Analyst pour le faire appelée " Zone de tabulation ", mais je ne trouve pas une telle fonctionnalité dans QGIS.
Si quelqu'un peut m'aider, je lui en serais très reconnaissant.
Réponses:
essayez d'utiliser "Cross-Classification and Tabulation" de SAGA GIS à partir du plugin Sextante.
J'espère que cette aide
la source
J'ai codé un plugin pour QGIS appelé LecoS (regardez dans le Plugin Downloader, plugins expérimentaux, installez Scipy d'abord!). Il peut atteindre ce que vous voulez (pourcentages de couverture terrestre pour un raster donné), mais jusqu'à présent, je n'ai pas trouvé le temps de coder une superposition de polygones en python et vous devez donc découper votre raster individuellement par tampon et entité vectorielle au préalable. Voir un cas d'utilisation ici .
la source
Dans ce cas, vous pouvez utiliser la classification croisée et la tabulation comme mentionné ci-dessus. Pour l'utiliser, vous devez d'abord avoir une couche raster avec des valeurs int. Selon votre couche d'utilisation, vous devrez effectuer un reclassement ou un recodage pour que votre raster soit prêt.
L'autre prémisse est que vos polygones (vos zones tampons) doivent avoir une valeur int. Sinon, vous devez le créer (avec la formule de calcul $ id + 1).
Pixellisez votre couche vectorielle en utilisant le code int comme valeur raster. Ici, vous devez mettre certaines conditions: utiliser la même valeur de pixel et la même extension que le raster dont vous souhaitez extraire les informations.
Utilisez maintenant la classification croisée et la tabulation à partir des outils SAGA. Vous devez considérer que cela va créer une table xy avec le nombre de pixels de chaque catégorie d'utilisation dans chaque polygone. Ici, vous devez insérer le "nombre maximum de classes", ce qui signifie que, par exemple, si vous avez un raster avec 10 catégories et 25 classes de polygones, votre maximum sera de 25.
Voici le résultat final: un tableau où vous avez le nombre de pixels de chaque catégorie dans chaque polygone. Si vous devez le convertir en surface, multipliez simplement le nombre de cellules par la résolution de votre raster.
la source
L'outil Zones tabulées se trouve dans la boîte à outils Spatial Analyst de la section Outils zonaux pour ArcGIS 10. Vous ne savez pas quelle version vous exécutez. Notez que vous devez disposer d'une licence distincte pour l'extension Spatial Analyst pour utiliser ces outils. On dirait que cet outil fera ce que vous cherchez.
J'espère que cela t'aides
la source