Intégration du SIG GRASS¶

Lancer l’extension GRASS¶

Pour pouvoir utiliser les fonctionnalités de GRASS et/ou visionner des données vecteurs ou raster dans QGIS, vous devez sélectionner et charger l’extension GRASS à l’aide du gestionnaire d’extensions. Cliquez sur le menu Extensions ‣ Installer/gérer les extensions, sélectionnez GRASS et cliquez sur [OK].

Vous pouvez maintenant charger des données raster et vecteur depuis un SECTEUR GRASS existant (voir section sec_load_grassdata). Ou alors vous pouvez créer un nouveau SECTEUR GRASS à l’aide de QGIS (voir section Créer un nouveau SECTEUR GRASS) et y importer des données raster et vecteur (voir section Importer des données dans un SECTEUR GRASS) pour réaliser des traitements à l’aide de la boîte à outils GRASS (voir section La Boîte à outils GRASS).

Charger des données GRASS raster et vecteur¶

Avec l’extension GRASS, vous pouvez charger des données raster ou vecteur à l’aide du bouton approprié dans la barre de menu. Ici nous utiliserons comme exemple le jeu de données QGIS sur l’Alaska (voir section Échantillon de données). Il contient un SECTEUR GRASS avec trois couches vectorielles et un raster d’élévation.

1. Créez un nouveau répertoire grassdata, téléchargez le jeu de données QGIS ‘Alaska’ qgis_sample_data.zip depuis http://download.osgeo.org/qgis/data/ et décompressez le dans le répertoire grassdata.

2. Lancez QGIS.

3. Si cela n’a pas déjà été fait dans une précédente session QGIS, chargez l’extension GRASS en cliquant sur Extensions ‣ Gestionnaire d’extensions et sélectionnez GRASS. La barre d’outils GRASS apparaît dans la barre de menu.

4. Dans la barre d’outils GRASS, cliquez sur le bouton ^{Ouvrir le jeu de données} pour ouvrir le gestionnaire de Base de données.

5. Pour Base de données GIS, parcourez puis sélectionnez ou entrez le chemin vers le répertoire nouvellement créé : grassdata.

6. Vous devriez maintenant être capable de sélectionner le SECTEUR alaska et le jeu de données demo.

7. Cliquez sur [OK]. Notez que les outils GRASS sont maintenant accessibles dans la barre d’outils.

8. Cliquez sur ^{Ajouter une couche raster GRASS}, choisissez le fichier gtopo30 et cliquez sur [OK]. Vous isionnez alors la couche d’élévation.

9. Cliquez sur ^{Ajouter une couche vectorielle GRASS}, choisissez la couche alaska et cliquez sur [OK]. La couche vectorielle Alaska s’affiche au-dessus du raster gtopo30. Vous pouvez modifier les propriétés de la couche d’information comme décrit dans le chapitre Fenêtre Propriétés d’une couche vecteur (par exemple modifier la transparence, changer la couleur du contour ou celle du remplissage).

10. Chargez également les deux autres couches vecteur, rivers et airports, et modifiez leurs propriétés.

Comme vous le voyez, il est très simple de charger des couches vecteur et raster de GRASS dans QGIS. Voir les sections suivantes pour l’édition des données GRASS et la création d’un nouveau SECTEUR. D’autres exemples de SECTEURs sont disponibles sur le site web de GRASS à http://grass.osgeo.org/download/sample-data/.

Secteur et Jeu de données GRASS¶

Les données GRASS sont stockées dans un répertoire référencé sous le nom GISDBASE. Ce répertoire, souvent appelé grassdata, doit être créé avant que vous commenciez à travailler avec l’extension GRASS dans QGIS. Dans ce répertoire, les données GRASS sont organisées par projets et stockées dans des sous-répertoires appelés SECTEUR ( LOCATION en Anglais). Chaque SECTEUR est défini par son système de coordonnées, sa projection et son étendue géographique. Chaque SECTEUR peut contenir plusieurs Jeux de données (MAPSETs en Anglais) (sous-répertoires du SECTEUR) qui sont utilisés pour subdiviser le projet en différents thèmes, sous régions ou espaces de travail pour chaque membre d’une équipe (Neteler & Mitasova 2008 Bibliographie). Pour pouvoir analyser des couches raster ou vecteur à l’aide des modules GRASS, vous devez les importer dans un SECTEUR. (Ce n’est pas complètement vrai, car avec les modules GRASS r.external et v.external, vous pouvez lier (en lecture seule) des données externes gérées par GDAL/OGR sans les importer. Mais comme il ne s’agit pas d’une fonctionnalité courante pour les débutants sur GRASS, elle ne sera pas décrite ici).

Figure GRASS location 1:

Données GRASS du SECTEUR Alaska

Créer un nouveau SECTEUR GRASS¶

À titre d’exemple, voici la manière dont le SECTEUR alaska, projeté en Albers Equal Area et ayant pour unité le pied, a été créé pour l’échantillon de données QGIS. Ce SECTEUR alaska sera utilisé pour tous les exemples et exercices GRASS qui suivent. Il est utile de le télécharger et de l’installer sur votre ordinateur (voir Échantillon de données).

1. Démarrez QGIS et assurez vous que l’extension GRASS est chargée.

2. Affichez le shapefile alaska.shp (voir section vector_load_shapefile) du jeu de données QGIS Alaska (voir Échantillon de données).

3. Dans la barre d’outils GRASS, cliquez sur ^{Nouveau jeu de données} pour ouvrir l’assistant de création de Jeux de données.

4. Sélectionnez un répertoire existant de base de données GRASS (GISDBASE) grassdata ou créez en un pour le nouveau SECTEUR avec le gestionnaire de fichiers de votre ordinateur. Cliquez sur le bouton [Suivant].

5. Nous pouvons utiliser cet assistant à la fois pour créer un nouveau Jeu de données dans un SECTEUR existant (voir section Ajouter un nouveau Jeu de données) et pour créer également un nouveau SECTEUR. Cliquez sur le bouton radio Créez un nouveau secteur (voir figure_grass_location_2).

6. Entrez un nom pour le SECTEUR – nous avons utilisé ‘alaska’ – et cliquez sur le bouton [Suivant].

7. Définissez la projection en cliquant sur le bouton radio Projection pour activer la liste des projections.

8. Nous utilisons la projection Albers Equal Area Alaska (pieds). Étant donné que nous savons qu’elle correspond au code EPSG 2964, nous le saisissons dans le champ de recherche. (Note : Si vous souhaitez reproduire la manipulation pour un autre SECTEUR et une autre projection dont vous ne connaissez pas le code EPSG, cliquez sur ^{Statut de la projection} dans le coin inférieur droit de la barre d’état de QGIS (voir section Utiliser les projections)).

9. Saisissez 2964 dans le Filtre pour sélectionner la projection.

10. Cliquez sur [Suivant].

11. Pour définir la région par défaut, nous devons saisir les limites Nord, Sud , Est et Ouest du SECTEUR. Ici il suffit de cliquer sur le bouton [Fixer l’emprise courante de |qg|], pour appliquer l’emprise du shapefile alaska.shp déjà chargé comme emprise par défaut.

12. Cliquez sur [Suivant].

13. Nous avons aussi besoin de définir un Jeu de données dans notre nouveau SECTEUR (étape indispensable lors de la création d’un nouveau SECTEUR). Vous pouvez l’appeler comme vous le souhaitez - nous utiliserons ‘demo’. GRASS crée automatiquement un Jeu de données spécial appelé PERMANENT, conçu pour stocker les données essentielles du projet, son emprise spatiale par défaut et la définition du système de coordonnées (voir Neteler & Mitasova 2008 Bibliographie).

14. Vérifiez le résumé pour vous assurez que tout est correct et cliquez sur [Terminer].

15. Le nouveau SECTEUR ‘alaska’ et les deux Jeux de données ‘démo’ et ‘PERMANENT’ sont créés. Le jeu de données ouvert à ce moment est ‘démo’, tel que vous l’avez défini.

16. Notez que certains outils de la barre d’outils GRASS qui n’étaient pas accessibles le sont maintenant.

Figure GRASS location 2:

Création d’un nouveau SECTEUR ou Jeu de données GRASS dans QGIS

Si ce processus semble long, il s’agit en fait d’un moyen simple et rapide de créer un SECTEUR. Le SECTEUR ‘alaska’ est maintenant prêt pour l’importation de données (voir section Importer des données dans un SECTEUR GRASS). Vous pouvez également utiliser des données raster ou vecteur existantes dans le SECTEUR ‘alaska’ inclues dans le jeu de données QGIS ‘Alaska’ Échantillon de données et continuez avec la section Le modèle vecteur de GRASS.

Importer des données dans un SECTEUR GRASS¶

Cette section donne un exemple d’importation de données raster et vecteur dans le SECTEUR GRASS ‘alaska’ fournit dans le jeu de données QGIS ‘Alaska’. Nous utiliserons la couche raster d’occupation du sol landcover.img et la couche vectorielle au format GML lakes.gml, toutes deux présentes dans le jeu de données ‘Alaska’ (voir Échantillon de données).

1. Démarrez QGIS et assurez vous que l’extension GRASS est chargée.

2. Dans la barre d’outils GRASS, cliquez sur ^{Ouvrir un jeu de données} pour ouvrir l’assistant Jeu de données.

3. Sélectionnez comme base de données GRASS, le répertoire grassdata dans le jeu de données QGIS Alaska, puis le SECTEUR ‘alaska’, le Jeu de donnée ‘demo’ et cliquez sur [OK].

4. Maintenant cliquez sur ^{Ouvrir les outils GRASS}. La boîte à outils GRASS s’ouvre (voir section La Boîte à outils GRASS).

5. Pour importer la couche raster landcover.img, cliquez sur le module r.in.gdal dans l’onglet Arborescence des modules. Ce module GRASS vous permet d’importer les fichiers raster gérés par la librairie GDAL dans un SECTEUR GRASS. La fenêtre r.in.gdal apparaît.

6. Naviguer jusqu’au répertoire raster dans le jeu de données QGIS ‘Alaska’ et sélectionnez le fichier landcover.img.

7. Définissez landcover_grass comme nom de sortie pour le raster et cliquez sur [Lancer]. Dans l’onglet Rendu, vous voyez la commande GRASS en cours r.in.gdal -o input=/path/to/landcover.img output=landcover_grass.

8. Lorsque Terminé avec succès s’affiche, cliquez sur [Vue]. La couche raster landcover_grass est maintenant importée dans GRASS et pourra être affichée dans QGIS.

9. Pour importer le fichier GML lakes.gml, cliquez sur le module v.in.ogr dans l’onglet Arborescence des modules. Ce module vous permet d’importer des données vectorielles gérées par OGR dans un SECTEUR GRASS. La fenêtre v.in.ogr apparaît.

10. Naviguer jusqu’au répertoire gml dans le jeu de données QGIS ‘Alaska’ et sélectionnez le fichier lakes.gml.

11. Définissez lakes_grass comme nom de sortie et cliquez sur [Lancer]. Vous n’avez pas besoin des autres options dans cet exemple. Dans l’onglet Rendu, vous voyez la commande GRASS en cours v.in.ogr -o dsn=/path/to/lakes.gml output=lakes\_grass.

12. Lorsque Terminé avec succès s’affiche, cliquez sur [Vue]. La couche raster :file:`lakes_grass`est maintenant importée dans GRASS et pourra être affichée dans QGIS.

Le modèle vecteur de GRASS¶

Il est important de comprendre le modèle vectoriel de GRASS avant de faire de la numérisation.

En général, GRASS utilise un modèle topologique pour les couches vecteur.

Cela signifie que les surfaces ne sont pas représentées par des polygones fermés et distincts, mais par une ou plusieurs limites. Une limite entre des polygones adjacents n’est numérisée qu’une seule fois et est partagée par les deux surfaces. Les limites doivent être connectées sans trous. Une surface est identifiée (et libellée) via le centroïde de la surface.

Outre les limites et centroïdes, une couche vectorielle peut également contenir des points et des lignes. Tous ces éléments de géométrie peuvent être mélangés dans une couche vectorielle et seront représentés dans différentes ‘sous-couches’ dans une carte vectorielle GRASS. Ainsi, une couche GRASS n’est pas un vecteur ou un raster, mais un niveau à l’intérieur d’une couche vectorielle. Il est important de bien distinguer ceci (même s’il est possible de mélanger des éléments de géométries différentes, c’est inhabituel et même dans GRASS, on l’utilise dans des cas particuliers tel que l’analyse de réseau. Normalement, vous devriez stocker des éléments de géométries différentes dans des couches différentes).

Il est possible de stocker plusieurs ‘sous-couches’ dans une couche vectorielle. Par exemple, des champs, de la forêt et des lacs peuvent être stockés dans une couche vectorielle. Des forêts et des lacs adjacents partagent les mêmes limites, mais ils auront des tables attributaires différentes. Il est aussi possible de faire correspondre une table attributaire aux limites. Par exemple, la limite entre un lac et une forêt peut être une route qui peut avoir une table attributaire différente.

La ‘sous-couche’ est définie dans GRASS par un chiffre. Ce chiffre définit s’il y a plusieurs sous-couches à l’intérieur d’une couche vectorielle (par exemple, il définit s’il s’agit de lac ou de forêt). Pour l’instant, il s’agit d’un nombre, mais dans des versions futures GRASS pourra utiliser des noms pour les sous-couches dans l’interface utilisateur.

Les données attributaires peuvent être stockées dans le SECTEUR au format dBase, SQLite3 ou dans des tables de bases de données externes comme par exemple : PostgreSQL, MySQL, Oracle, etc.

Les données attributaires sont liées à la géométrie par le biais d’un champ ‘category’.

‘Category’ (clé, ID) est un entier attaché à la géométrie, et il est utilisé comme lien vers une colonne de clé dans la table de base de données.

Création d’une nouvelle couche vectorielle GRASS¶

Pour créer une nouvelle couche vectorielle GRASS à l’aide de l’extension GRASS, cliquez sur ^{Créez une nouvelle couche vectorielle GRASS} dans la barre d’outils. Entrez le nom de la couche dans la fenêtre et vous pouvez commencer à digitaliser un point, une ligne ou un polygone en suivant les instructions de la section Numérisation et édition de couche vectorielle GRASS.

Dans GRASS, il est possible de gérer plusieurs types de géométrie (point, ligne et surface) dans une seule couche d’information, car GRASS utilise un modèle vecteur topologique. Vous n’avez donc pas besoin de sélectionner un type de géométrie quand vous créez une couche vectorielle GRASS. Ce comportement est différent de celui de la création de shapefile avec QGIS, car les shapefiles utilisent un modèle vecteur d’entité simple (voir section Créer de nouvelles couches vecteur).

Numérisation et édition de couche vectorielle GRASS¶

Les outils de numérisation pour les couches vectorielles de GRASS sont accessibles via ^{Éditer une couche vectorielle GRASS} dans la barre d’outils. Assurez-vous d’avoir ouvert une couche vectorielle GRASS et sélectionné la sous-couche dans la légende avant d’utiliser l’outil d’édition. La figure figure_grass_digitizing_2 montre la fenêtre GRASS qui s’affiche quand vous cliquez sur l’outil d’édition. Les outils et les paramètres d’édition sont détaillés dans les paragraphes suivants.

Barre d’outils

Sur la figure figure_grass_digitizing_1 vous pouvez voir la barre d’outils d’édition GRASS de l’extension GRASS. Le tableau table_grass_digitizing_1 récapitule les fonctions disponibles.

Figure GRASS digitizing 1:

Barre d’outils d’édition GRASS

Icône	Outil	Fonction
	Nouveau Point	Numérise un nouveau point
	Nouvelle Ligne	Numérise une nouvelle ligne
	Nouveau Contour	Numérise un nouveau contour (terminer la numérisation en sélectionnant un nouvel outil)
	Nouveau Centroïde	Numérise un nouveau centroïde (permet d’étiqueter un polygone existant)
	Déplacer un sommet	Déplace un sommet d’une ligne ou d’un polygone existant et indique sa nouvelle position
	Ajouter un sommet	Ajoute un nouveau sommet à une ligne existante
	Effacer un sommet	Efface un sommet d’une ligne existante (confirmez le sommet sélectionné par un autre clic)
	Déplacer l’élément	Déplacez la limite, la ligne, le point ou le centroïde sélectionné puis cliquez sur la nouvelle position
	Couper la ligne	Coupe une ligne existante en deux parties
	Effacer l’élément	Efface une limite, une ligne, un point ou un centroïde existant (confirmez l’élément sélectionné avec un autre clic)
	Éditer les attributs	Édite les attributs de l’élément sélectionné (notez qu’un seul élément peut représenter plusieurs géométries, voir ci-dessus)
	Fermer	Ferme la session et sauvegarde l’état actuel (reconstruit la topologie après)

Tableau Numérisation avec GRASS 1 : outils d’édition GRASS

Onglet Categorie

L’onglet Catégorie vous permet de définir la manière dont les valeurs du champ category sont assignées au nouvel élément géométrique.

Figure GRASS digitizing 2:

Onglet Catégorie de la fenêtre d’édition GRASS

• Mode : La catégorie qui sera appliquée aux nouveaux éléments.

  • Prochain non utilisé - Applique la valeur suivante non utilisée du champ category à l’élément géométrique.

  • Saisie manuelle - Saisir manuellement la valeur du champ ‘category’ pour l’élément géométrique.

  • Pas de catégorie - Ne pas remplir le champ category. C’est par exemple utilisé pour les surfaces, car les valeurs de catégorie sont stockées via le centroïde.

• Categorie - Un identifiant (ID) associé à chaque objet numérisé. Il est utilisé pour connecter les objets géométriques avec ces attributs.

• Couche - Chaque objet peut être connecté à différentes tables attributaires au travers des différentes sous-couches. Le numéro de sous-couche par défaut est 1.

Onglet Paramètres

L’onglet Paramètres vous permet de définir la tolérance d’accrochage en pixels-écrans. Le seuil définit à partir de quelle distance les nouveaux points ou les nouvelles lignes sont accrochées automatiquement à des noeuds existants. Cela aide à éviter de créer des trous ou des superpositions entre les contours. La valeur par défaut est fixée à 10 pixels.

Figure GRASS digitizing 3:

Onglet Paramètres de la fenêtre d’édition de GRASS

Onglet Convention des signes

L’onglet Convention des signes vous permet d’afficher et modifier la symbologie, la couleur des différentes formes géométriques ainsi que leur statut topologique (par exemple : contour ouvert / fermé).

Figure GRASS digitizing 4:

Onglet Convention des signes de la fenêtre d’édition de GRASS

Onglet Table

L’onglet Table donne des informations sur la table attributaire d’une ‘sous-couche’ donnée. C’est ici que vous pouvez ajouter des colonnes à une table attributaire existante ou créer une nouvelle table attributaire pour une nouvelle couche vectorielle GRASS (voir section Création d’une nouvelle couche vectorielle GRASS).

Figure GRASS digitizing 5:

Onglet Table de la fenêtre d’édition de GRASS

L’outil région GRASS¶

La définition d’une région (définir une emprise spatiale de travail) dans GRASS est très importante pour travailler avec des couches rasters. Le travail d’analyse vecteur n’est, par défaut, pas limitée à une région définie. Mais, tous les rasters nouvellement créés auront l’emprise spatiale et la résolution de la région GRASS en cours d’utilisation, indépendamment de leur emprise et résolution d’origine. La région courante GRASS est stockée dans le fichier $LOCATION/$MAPSET/WIND, et celui-ci définit les limites Nord, Sud, Est et Ouest, le nombre de lignes et de colonnes ainsi que la résolution spatiale horizontale et verticale.

Il est possible d’afficher ou de masquer l’emprise de la région GRASS dans QGIS à l’aide du bouton ^{Afficher la région courante GRASS}.

A l’aide du bouton ^{Éditer la région courante GRASS} vous avez accès à une fenêtre qui vous permet de modifier la région courante ainsi que sa symbologie. Entrez les nouvelles limites et résolution et cliquez sur [OK]. Cette fenêtre vous permet aussi de définir une nouvelle région interactivement à l’aide de la souris. Pour définir ce rectangle d’emprise, cliquez avec le bouton gauche de la souris et définissez un rectangle que vous terminerez en cliquant de nouveau sur le bouton gauche de la souris et cliquez sur [OK].

Le module GRASS g.region propose un grand nombre de paramètres pour définir de façon appropriée les limites et la résolution d’une région pour faire de l’analyse raster. Vous pouvez vous servir de ces paramètres dans la boîte à outils GRASS décrite dans la section La Boîte à outils GRASS.

La Boîte à outils GRASS¶

La fenêtre ^{Ouvrir les outils GRASS} donne accès aux fonctionnalités GRASS qui permettent de travailler sur les données d’un SECTEUR et d’un Jeu de données. Pour utiliser la Boîte à outils GRASS, vous devez ouvrir un SECTEUR et un Jeu de données sur lequel vous avez des droits d’écriture (que vous avez normalement si vous avez créé le Jeu de données). Cela est nécessaire car les rasters et les vecteurs nouvellement créés lors des analyses doivent être écrits dans le SECTEUR et Jeu de données courant.

Figure GRASS Toolbox 1:

La Boîte à outils GRASS et la l’Arborescence des modules

Travailler avec les modules GRASS¶

La console de la Boîte à outils GRASS vous donne accès à pratiquement tous les modules GRASS (plus de 300) en ligne de commande. Afin d’offrir un environnement de travail plus agréable, environ 200 d’entre eux sont disponibles via l’interface graphique de la Boîte à outils GRASS.

La liste des modules GRASS disponibles via la boîte à outils de QGIS 2.2 est détaillée sur le wiki de GRASS : http://grass.osgeo.org/wiki/GRASS-QGIS_relevant_module_list.

Il est aussi possible de personnaliser le contenu de la boîte à outils GRASS. Ceci est décrit dans la section Paramètrer la boîte à outils GRASS.

Comme indiqué sur la figure figure_grass_toolbox_1, vous pouvez chercher le module GRASS approprié en utilisant l’onglet Arborescence des modules ou en utilisant l’onglet Liste des Modules pour faire une recherche.

Lorsque vous cliquez sur un module, un nouvel onglet apparaît proposant trois sous-onglets : Options, Rendu et Manuel.

Options

L’onglet Options propose une interface simplifiée où vous pouvez sélectionner un raster ou un vecteur en cours de visualisation dans QGIS et saisir les paramètres spécifiques au module avant de le lancer.

Figure GRASS module 1:

Boîte à outils GRASS, onglet Options d’un module

Tous les paramètres du module ne sont généralement pas fournis afin de simplifier les fenêtres. Pour utiliser des paramètres qui ne se trouvent pas dans la fenêtre, vous devez utiliser la console GRASS et lancer les modules en lignes de commande.

Une nouvelle fonctionnalité depuis QGIS 1.8.0 est l’ajout d’un bouton Afficher les options avancées >> en-dessous de la fenêtre simplifiée de l’onglet Options. Pour l’instant seul le module v.in.ascii a été adapté afin de servir d’exemple d’utilisation mais d’autres le seront dans les prochaines versions de QGIS. La finalité est de pouvoir recourir à toutes les options de GRASS sans devoir ouvrir la console GRASS.

Rendu

Figure GRASS module 2:

Boîte à outils GRASS, onglet Sortie d’un module

L’onglet Rendu fournit des informations sur l’état de sortie du module. Quand vous cliquez sur le bouton [Lancer], le module passe sur l’onglet Rendu et vous voyez les informations sur le processus en cours. Si tout se passe bien, vous verrez finalement le message Terminé avec succès.

Manuel

Figure GRASS module 3:

Boîte à outils GRASS, onglet Manuel d’un module

L’onglet Manuel montre la page HTML d’aide du module GRASS. Vous pouvez vous en servir pour voir les autres paramètres du module et pour avoir une connaissance plus approfondie de l’objet du module. À la fin de chaque page d’aide d’un module, vous avez des liens vers Main Help index (index principal), Thematic.index (index par thème) et Full.index (index complet). Ces liens vous donnent les mêmes informations que si vous utilisiez directement g.manual.

Astuce

Afficher les résultats immédiatement

Si vous voulez voir immédiatement dans votre fenêtre carte le résultat des calculs du module, vous pouvez utiliser le bouton ‘Vue’ au bas de l’onglet du module.

Exemples de modules GRASS¶

Les exemples suivants décrivent les possibilités de certains modules GRASS.

Création de courbes de niveau¶

Le premier exemple permet de créer une couche vectorielle de courbes de niveau à partir d’un modèle numérique de terrain (MNT). Ici, nous considérerons que le SECTEUR Alaska a été installé comme décrit dans la section Importer des données dans un SECTEUR GRASS.

• Premièrement, ouvrez le secteur en cliquant sur le bouton ^{Ouvrir le jeu de données} et choisissez le secteur Alaska.

• Maintenant chargez le raster gtopo30 en cliquant sur le bouton ^{Ajouter une couche raster GRASS} puis en sélectionnant le raster gtopo30 dans le secteur demo.

• Ouvrez la boite à outils à l’aide du bouton ^{Ouvrir les outils GRASS}.

• Dans la liste des outils double-cliquez sur Raster -> Gestion de surface -> Générer des lignes vectorielles de contours.

• Maintenant, cliquez sur l’outil r.contour, cela ouvrira une fenêtre comme expliqué ci-dessus (voir Travailler avec les modules GRASS). Le raster gtopo30 devrait apparaitre dans le champ Nom du raster en entrée.

• Dans le champ Increment between Contour levels , saisissez la valeur 100. (Cela va créer des courbes de niveau tous les 100 mètres)

• Saisisez dans le champ Nom de la couche vectorielle en sortie, le nom ctour_100.

• Cliquer sur [Lancer] pour lancer le traitement. Attendez quelques instants que le message Terminé avec succès apparaisse à l’écran. Cliquez enfin sur [Vue] puis [Fermer].

Comme il s’agit d’une grande région, cela prendra un certain temps à s’afficher. Une fois l’affichage terminé, vous pouvez ouvrir la fenêtre de propriétés de la couche pour changer la couleur des courbes de niveau afin qu’elles apparaissent clairement au dessus de la couche raster d’élévation comme décrit dans Fenêtre Propriétés d’une couche vecteur.

Zoomez sur une petite région montagneuse du centre de l’Alaska. Avec un zoom important, vous constaterez que les courbes de niveau sont constituées de lignes brisées avec des angles vifs. GRASS offre la possibilité de généraliser les cartes vecteurs à l’aide de l’outil v.generalize, tout en conservant leur forme générale. L’outil utilise différents algorithmes ayant différents objectifs. Certains de ces algorithmes (par exemple : Douglas Peucker et Réduction de Vertex) simplifient les lignes en supprimant des sommets. La couche simplifiée se chargera plus rapidement. Cette commande est utile lorsque vous avez une couche vectorielle très détaillée et que vous créez une carte à petite échelle où les détails ne sont donc pas nécessaires.

Astuce

L’outil de simplification

Vous remarquerez que fTools dispose aussi d’un outil de simplification Simplifier la géometrie qui fonctionne comme l’algorithme Douglas-Peucker de GRASS, v.generalize.

Cependant, le but de cet exemple est différent. Les courbes de niveau créées avec r.contour ont des angles vifs qui doivent être lissés. Parmi les algorithmes de v.generalize, il y a l’algorithme de Chaiken qui fait justement ça (comme Hermite splines). Gardez à l’esprit que ces algorithmes peuvent ajouter des sommets supplémentaires au vecteur, l’amenant à se charger encore plus lentement.

• Ouvrez la Boîte à outils GRASS et double cliquez sur Vecteur -> Développer la carte -> Généralisation. Cliquez alors sur le module v.generalize pour ouvrir sa fenêtre d’options.

• rifier que la couche vectorielle ‘ctour_100’ apparait dans le champ Nom de la couche vectorielle en entrée.

• Dans la liste des algorithmes choisissez Chaiken. Laisser les autres options par défaut et descendez à la dernière ligne pour donner le nom de la couche d’information à créer : Nom de la couche vectorielle en sortie ‘ctour_100_smooth’, et cliquez sur [Lancer].

• Cela peut prendre plusieurs minutes. Lorsque le texte Terminé avec succès apparait, cliquez sur le bouton [Vue] puis sur [Fermer].

• Vous pouvez changer la couleur de cette couche vectorielle pour qu’elle apparaisse clairement sur le raster et qu’elle contraste aussi avec la couche de départ.Vous remarquerez que les nouvelles courbes de niveau ont des angles plus arrondis que l’original tout en restant fidèle à la forme globale d’origine.

Figure GRASS module 4:

Module GRASS v.generalize utilisé pour simplifier une couche vectorielle

Astuce

Autres utilisations de r.contour

La procédure décrite ci-dessus peut être utilisée dans d’autres cas similaires. Si vous disposez d’une couche d’informations raster représentant des précipitations, par exemple, vous pouvez utiliser la même méthode pour créer des isohyètes (lignes reliant des points d’égales quantités de précipitations).

Créer un ombrage avec effet 3D¶

Différentes méthodes sont utilisées pour afficher les modèles numérique de terrain et donner un effet 3D au carte. L’utilisation de courbes de niveau comme décrit ci-dessus est un des moyens souvent utilisés pour produire des cartes topographiques. Un autre moyen de rendre cet effet 3D est d’utiliser l’ombrage. L’ombrage est créé à partir du modèle numérique de terrain (MNT) en calculant d’abord les pentes et les expositions puis en simulant la position du soleil dans le ciel ce qui donne à chaque cellule une valeur de réflectance. Les pentes éclairées par le soleil sont plus claires et les pentes à l’abri du soleil sont plus sombres.

• Commencez par ouvrir la couche raster gtopo30. Ouvrez la Boîte à outils GRASS et dans la catégorie Raster double cliquez sur Analyse spatiale ‣ Analyse de terrain.

• Cliquez ensuite sur r.shaded.relief pour lancer le module.

• Changer l’ Azimuth du soleil par rapport au nord, en degrés et mettez 315.

• Saisissez gtopo30_shade comme nom pour la nouvelle couche d’ombrage et cliquez sur le bouton [Lancer].

• Quand le calcul est terminé, ajoutez le raster d’ombrage à la fenêtre carte. Normalement, il devrait s’afficher en niveau de gris.

• Pour voir les deux couches d’informations ombrage et gtopo30 en même temps, placez la couche ombrage sous la couche gtopo30 dans le gestionnaire de couches et ouvrez la fenêtre Propriétés de la couche gtopo30, allez sur l’onglet Transparence et fixez la transparence à environ 25%.

Vous devriez maintenant avoir la couche gtopo30 en couleur et en transparence, affiché au dessus de la couche d’ombrage en niveau de gris. Pour bien visualiser l’effet d’ombrage, décochez puis recochez la couche gtopo30_shade dans la légende.

Utiliser la console GRASS

L’extension GRASS de QGIS est faite pour les utilisateurs ne connaissant pas GRASS et qui ne sont pas familiers avec les modules et les options. Ainsi, certains modules dans la Boîte à outils n’apparaissent pas avec toutes les options possibles et certains n’apparaissent pas du tout. La console GRASS donne accès à ces modules additionnels qui n’apparaissent pas dans la Boîte à outils et aussi aux options des modules qui n’apparaissent que de façon simplifiés dans la Boîte à outils. Cet exemple montre l’utilisation des options supplémentaires du module r.shaded.relief utilisé ci-dessus.

Figure GRASS module 5:

La console GRASS utilisation du module r.shaded.relief

Le module r.shaded.relief possède un paramètre zmult qui multiplie la valeur de l’altitude (exprimé dans la même unité que les coordonnées X - Y) ce qui a pour effet d’accentuer le relief.

• Ouvrez le raster gtopo30 comme ci-dessus, lancez la Boîte à outils GRASS et ouvrez la console GRASS. Dans la console, entrez la ligne suivante r.shaded.relief map=gtopo30 shade=gtopo30_shade2 azimuth=315 zmult=3 et pressez [Entrée].

• Une fois le calcul terminé, allez sur l’onglet Parcourir et double-cliquez sur le nouveau raster gtopo30_shade2 pour l’afficher dans QGIS.

• Comme expliqué ci-dessus, placez le raster d’ombrage sous le raster gtopo30 puis vérifiez la transparence du raster gtopo30. Vous devriez constater que le relief apparaît plus marqué qu’avec le premier raster d’ombrage.

Figure GRASS module 6:

Affichage du relief ombré créé avec le module GRASS r.shaded.relief

Statistiques raster avec des couches vectorielles¶

L’exemple suivant comment un module GRASS peut aggréger des données raster et ajouter des colonnes de statistiques pour chaque polygone d’une couche vectorielle.

• Encore une fois, nous allons utiliser le jeu de données Alaska. Référez vous à Importer des données dans un SECTEUR GRASS pour importer les shapefiles contenus dans le répetoire shapefiles dans GRASS.

• Un étape intermédiaire est nécessaire : des centroïdes doivent importés afin d’avoir une couche GRASS vecteur complète (qui inclue les contours et les centroïdes).

• Dans la Boîte à outils choisissez Vecteur -> Gestion des entités et ouvrez le module v.centroids.

• Entrez ‘forest_areas’ comme nom de couche en sortie et lancez le module.

• Now load the forest_areas vector and display the types of forests - deciduous, evergreen, mixed - in different colors: In the layer Properties window, Symbology tab, choose from Legend type ‘Unique value’ and set the Classification field to ‘VEGDESC’. (Refer to the explanation of the symbology tab in sec_symbology of the vector section.)

• Réouvrez la Boîte à outils GRASS et ouvrez Vecteur -> Mise à jour vectorielle via d’autres cartes.

• Cliquez sur le module v.rast.stats. Saisissez gtopo30 et forest_areas.

• Un seul paramètre additionnel est requis : Entrez elev pour le column prefix, et cliquez sur le bouton [Lancer]. C’est un opération lourde qui peut durer longtemps (jusqu’à deux heures).

• Pour finir, ouvrez la table attributaire de forest_areas, et vérifiez que plusieurs nouvelles colonnes ont étés ajoutées dont elev_min, elev_max, elev_mean, etc., pour chaque polygone de forêt.

Travailler avec le navigateur GRASS¶

Une autre fonctionnalité utile de la boîte à outils GRASS est le navigateur de SECTEUR GRASS. Sur la figure figure_grass_module_7 vous pouvez voir le SECTEUR en cours avec ses Jeux de données.

Dans la partie gauche de la fenêtre vous pouvez naviguer dans tous les Jeux de données du SECTEUR courant. La partie droite de la fenêtre affiche des informations sur le raster ou le vecteur sélectionné (par exemple la résolution, l’emprise, la source des données, les tables attributaires pour les vecteurs et un historique des commandes).

Figure GRASS module 7:

Navigateur de SECTEUR GRASS

La barre d’outils de l’onglet Parcourir donne accès aux outils de gestion du SECTEUR sélectionné :

• Ajoute la carte sélectionnée à la carte QGIS

• Copie la carte sélectionnée

• Renomme la carte sélectionnée

• Effacer la carte sélectionnée

• Région courante réglée sur la carte choisie

• Rafraîchir

Les commandes Renommer la carte sélectionnée et Effacer la carte sélectionnée ne fonctionnent qu’avec les cartes présentes dans votre Jeu de données sélectionné. Tous les autres outils fonctionnent aussi avec les autres Jeux de données.

<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE qgisgrassmodule SYSTEM "http://mrcc.com/qgisgrassmodule.dtd">

<qgisgrassmodule label="Vector buffer" module="v.buffer">
        <option key="input" typeoption="type" layeroption="layer" />
        <option key="buffer"/>
        <option key="output" />
</qgisgrassmodule>