TP QGIS : Introduction GRASS

Introduction GRASS avec Quantum

GIS
Ce document fait partie dune srie de tutoriels destins au logiciel Quantum GIS dans
sa version 1.8 Lisboa. Il ne remplace pas le manuel utilisateur du logiciel, il est donc
conseill de sy rfrer pour plus de dtails sur les fonctionnalits.
Auteur : Marie SILVESTRE marie.silvestre@upmc.fr
CNRS FR3020 FIRE Fdration Ile-de-France de Recherche sur lEnvironnement
Ce document est sous la licence CC BY-SA 3.0 : Creative Commons Paternit Partage
des conditions initiales lidentique.
Pour suivre ce document dune manire pratique, lorsque des actions sont demandes aux
lecteurs, les paragraphes dinstructions sont signals par une barre verticale.
Notions abordes :
Cration dun jeu de donnes GRASS,
Import de couches raster et vecteur,
Cration dun rseau hydrographique partir dun MNT,
Ralisation de profils,
Visualisation en 3D avec NVIZ.
GRASS est un logiciel SIG libre trs puissant permettant notamment deffectuer des
traitements dimages pousss (classification dimage de tldtection, etc.). Ce logiciel est
indpendant de QGIS mais il a t intgr QGIS sous la forme dune extension principale
(installe par dfaut).
Nous allons travailler sur la zone des sources du Drac (affluent de lIsre) qui se situe
dans le dpartement des Hautes-Alpes. Lobjectif est dapprendre crer un jeu de donnes
GRASS et de sinitier quelques fonctions du logiciel : cration dun rseau hydrographique,
profil en long, visualisation 3D.
Lextension GRASS au sein de QGIS sactive via le Gestionnaire dextensions.

1
1.1

Cration dun jeu de donnes GRASS

Principe de structuration dune base de donnes GRASS

Toutes les donnes GRASS se stockent dans un rpertoire, souvent appel grassdata.
Ce rpertoire est ensuite subdivis en sous-rpertoires, chacun correspondant un projet
(SECTEUR ou LOCATION en anglais). Un secteur se caractrise par un systme de coordonnes et une emprise spatiale. Enfin, un secteur se dcompose en diffrents Jeux de
donnes (ou MAPSET en anglais). Ce systme permet de dfinir des droits de lecture et
dcriture diffrents selon les jeux de donnes et permet plusieurs utilisateurs de travailler sur une mme base de donnes centralise. Pour un utilisateur isol, ce dcoupage
en jeux de donnes permet surtout dorganiser les donnes par sous-rgions et / ou par
thmatiques de travail.
Pour chaque secteur, GRASS cre automatiquement un jeu de donnes appel PERMANENT.
Dans un environnement multi-utilisateurs, en lui dfinissant des droits de lecture seule, il

TP QGIS : Introduction GRASS

permet de stocker les donnes de bases, susceptibles dtre utilises dans tous les jeux
de donnes. Dans le cas dun utilisateur isol, il sagit dun jeu de donnes comme les
autres.
Pour plus dinformations sur la structuration dune base de donnes GRASS, consultez
cet article de Martin Laloux sur le PortailSIG.

1.2

Cration dun jeu de donnes GRASS via QGIS

Lextension GRASS de QGIS permet de faciliter ltape de cration dune base de donnes GRASS en se basant sur les caractristiques du projet courant (systme de coordonnes
et emprise).
Crez un projet QGIS ayant pour SCR le RGF93/Lambert 93 (EPSG 2154).
Chargez le MNT du drac et zoomez sur son emprise.
Lextension GRASS active propose sa propre barre doutils :

Figure 1 Barre doutils GRASS.

Activez lextension GRASS et lancez la cration dun nouveau jeu de donnes.
Voici les tapes de cration dun jeu de donnes GRASS via QGIS :
1. Cration du rpertoire contenant lensemble des donnes GRASS (Godatabase ou
base de donnes).
2. Nom du secteur crer.
3. Dfinition du systme de coordonnes du secteur.
4. Dfinition de lemprise du secteur : ici, pour rentrer automatiquement lemprise
de la fentre QGIS, il faut cliquer sur le bouton "Fixer lemprise courante de
QGIS".
5. Nom du jeu de donnes crer.
Crez une base de donnes GRASS avec :
- le secteur "Drac" ayant pour projection le RGF93/Lambert 93 (EPSG 2154) et
pour emprise la zone du MNT du Drac charge,
- le jeu de donnes "tp_intro".

TP QGIS : Introduction GRASS

1 : cration de la Godatabase

2 : nom du secteur

3 : cration de la Godatabase

4 : nom du secteur

5 : nom du jeu de donnes

6 : rcapitulatif

TP QGIS : Introduction GRASS

Import de donnes dans GRASS

Un jeu de donnes nouvellement cr est vide. Le logiciel GRASS dispose de ses propres
formats de donnes raster et vecteur. Pour utiliser des donnes dans GRASS il faut donc
les importer.
Prsentation de la fentre doutils GRASS
Laccs aux fonctionnalits de GRASS se fait en cliquant sur le bouton "Ouvrir les
outils GRASS". La fentre qui souvre alors prsente trois onglets :
"Arborescence des modules" donne accs un ensemble (non exhaustif) de fonctions GRASS,
"Liste des modules" propose laccs aux mme fonctions mais permet de les filtrer
rapidement,
"Parcourir" permet dexplorer les jeux de donnes du secteur ouvert.

2.1

Import dun raster dans un jeu de donnes GRASS

GRASS permet dimporter tous les formats de donnes vecteur et raster supports par
la librairie GDAL/OGR et propose dimporter des donnes charges dans QGIS comme
par exemple pour un raster via loutil GRASS "Gestion de fichier / Importer dans
GRASS / Importer un raster dans GRASS / Importer un raster dans GRASS depuis
QGIS / Importer un raster dj charg" (r.in.gdal.qgis).
Lorsque lon clique sur un outil, un nouvel onglet apparat automatiquement qui lui-mme
se dcompose en trois sous-onglets :
"Options" : paramtres (principaux) de loutil,
"Rendu" : rcapitule la commande lance et son excution,
"Manuel" : aide sur loutil et liste exhaustive de ses paramtres.
Importez le MNT du drac dans le jeu de donnes "tp_intro".
Visualisation dun raster GRASS
Dans longlet "Parcourir", pour mettre jour la liste des couches des jeux de donnes,
il faut cliquer sur le bouton Rafrachir. Les couches y sont tries par grand type : raster et
vecteur. Pour ajouter une couche dun jeu de donnes GRASS QGIS, il faut slectionner
la couche puis cliquer sur le bouton "Ajoute la carte slectionne la carte".
Supprimez le raster MNT_drac.tif de QGIS.
Visualisez le raster import au jeu de donnes "tp_intro".
Par dfaut, un raster na pas de symbologie et apparat en gris uniforme. Loutil GRASS
"Raster / Gestion des cartes de couleurs / Dfinir une table de couleurs raster
depuis une table" (r.colors.table) permet dappliquer une symbologie prdfinie
un raster.
Appliquez la table de couleur "Maps relatives ranges of rater values to elevation color
ramp" au MNT (testez-en aussi dautres !).

TP QGIS : Introduction GRASS

2.2

Import de donnes vectorielles dans un jeu de donnes GRASS

Limport de couche vecteur se fait via loutil GRASS "Gestion de fichier / Importer
dans GRASS / Importer un vecteur dans GRASS / Importer un vecteur dj charg"
(v.in.ogr.qgis).
Nous allons
Nous allons
table est de
CLC.
Importez la

importer la couche doccupation du sol (CLC_drac) dans GRASS :

avoir besoin du code CLC sous forme numrique or le champ CODE_06 de la
type texte : crez un nouveau champ de type entier pour y stocker les codes
couche dans GRASS et affichez-la. Supprimez la couche .shp du projet.

Cration dun raster dombrage

Pour amliorer le rendu 2D des donnes nous allons crer partir du MNT un raster
dombrage.
Paramtres de la rgion GRASS
Avant de se lancer dans les traitements raster avec GRASS il faut toujours vrifier les
paramtres de rgion : tendue et rsolution. Ils sont stocks dans le fichier WIND qui se
trouve la racine de chaque jeu de donnes. Ils sont utiliss par GRASS pour dfinir la
taille (tendue et rsolution) des rasters crer en sortie de traitement.
Ltendue de la rgion se matrialise dans QGIS par un rectangle rouge et correspond
par dfaut ltendue du jeu de donnes. Loutil "Paramtre de rgion / Rduire la
rgion courante jusqu ce quelle rencontre des donnes non vides dans le raster"
(g.region.zoom) permet de caler ltendue de la rgion sur celle dun raster existant.
Modifiez ltendue de la rgion pour quelle concide avec celle du MNT.
Pour vrifier et modifier la rsolution de la rgion, il faut cliquer sur le bouton "diter
la rgion courante GRASS" de la barre doutils GRASS.
Vrifiez et modifiez la rsolution de la rgion pour quelle concide avec celle du MNT.
Ombrage
Un raster dombrage se gnre avec loutil GRASS "Raster / Analyse spatiale /
Analyse de terrain / Crer une carte dombrage" (r.shaded.relief).
Crez un raster dombrage avec une altitude de 45 degrs et un azimut de 315 degrs.
Affichez lombrage en transparence au dessus de toutes les autres couches.

Cration dun profil en long

Nous allons gnrer un profil daltitude et un profil doccupation du sol le long dune
ligne que nous allons crer dans QGIS.
Crez un shapefile comportant une ligne qui traverse une valle.
Importez cette ligne dans le jeu de donnes GRASS.

TP QGIS : Introduction GRASS

Loutil permettant de calculer un profil demande en entre un fichier de coordonnes et

le raster sur lequel on souhaite faire le profil. Il faut donc convertir notre couche de ligne
en un fichier de coordonnes. Loutil dexport au format ascii se trouve dans "Gestion de
fichier / Exporter depuis GRASS / Exporter un vecteur depuis GRASS / Convertir
les vecteurs binaires GRASS en ASCII GRASS" (v.out.ascii).
Exportez la couche de ligne au format ASCII :
- spcifiez le chemin et le nom du fichier (sans extension),
- prcisez le Format de sortie standard,
- dans les options avances, rentrez un espace comme sparateur de champ.
Le fichier de coordonnes ne doit comporter que les coordonnes X et Y des points
spares par un espace (une paire de coordonnes par ligne), pas de ligne de titre.
Nettoyez le fichier de points de manire adapte.
La console GRASS
Loutil de profil GRASS ne dispose pas dinterface graphique dans QGIS. Pour lutiliser,
il faut passer par la console GRASS. Celle-ci souvre via la liste des modules puis en cliquant
sur shell.
Ouvrez la console GRASS. Dplacez-vous jusqu votre rpertoire de travail, o est
enregistr votre fichier de points (par exemple cd C:\TD_FIRE_SIG\donnees\grass).
Lancez la commande :
cat fichier_points | r.profile -c -g input=nom_mnt output=profil_mnt
Dans cette instruction, -c indique quil faut exporter le code RGB des points et -g
leurs coordonnes.
Ouvrez le fichier profil_mnt dans un tableur et tracez le profil daltitude.

Figure 2 Profil daltitude.

TP QGIS : Introduction GRASS

Profil doccupation du sol

De la mme manire nous allons crer un profil doccupation du sol partir de la
couche CLC importe. Loutil de conversion raster se trouve dans "Gestion de fichier /
Conversion de type de carte / Convertir un vecteur en raster GRASS / Convertir
le vecteur en raster en utilisant les valeurs dattributs" (v.to.rast.attr).
Convertissez la couche vecteur doccupation du sol en raster en prcisant le Champs
dattributs contenant le code CLC.
Appliquez au raster la table de couleur "EU Corine land cover colors".
Ralisez le profil doccupation du sol sur les mmes points que prcdemment.

Extraction dun rseau hydrographique partir dun MNT

Nous allons extraire le rseau hydrographique partir du MNT. Loutil "Raster /

Modles spatiaux / Modlisation hydrologique / Analyse de bassins de retenue
deau" (r.watershed) permet de gnrer :
le raster de direction des coulements,
le raster daccumulation : nombre de cellules qui scoulent dans dautres cellules,
le raster des segments de cours deau : la finesse du rseau dpend de la taille minimale
fixe pour chaque bassin,
le raster des bassins versants (un bassin correspondra un segment de cours deau
et un seul).
Gnrez les rasters de cours deau et de bassins versants. Choisissez tout dabord une
taille minimale de bassin de 5 km2 sachant quun pixel mesure 25 m par 25m.
Testez ensuite dautres valeurs.
Extraction de lamont dun point
Nous souhaitons extraire la surface de bassin situe en amont dun point prcis, dont
nous connaissons les coordonnes. Pour cela il faut utiliser loutil GRASS "Raster /
Modles spatiaux / Modlisation hydrologique / Crer des bassins dinondation"
(r.water.outlet), lui fournir un raster de direction des coulements et les coordonnes
de lexutoire choisit.
En vous basant sur le raster des cours deau gnr par loutil r.watershed et en utilisant
loutil QGIS "Saisie de coordonnes", relevez les coordonnes du point exutoire de
votre choix.
Ralisez ensuite lextraction de lamont de ce point.
Conversion en vecteur
Nous allons convertir le raster des cours deau en couche vecteur. Avant cela il faut
procder une tape de nettoyage sommaire du trac des cours deau. Loutil "Raster /
Analyse spatiale / Transformer des entits / Fines cellules non nulles indiquant
la prsence de lignes" (r.thin) permet daffiner le trac.
Affinez le trac des cours deau.

TP QGIS : Introduction GRASS

Loutil dexport des donnes raster en vecteur se trouve dans "Gestion de fichier /
Conversion de type de carte / Convertir un raster en vecteur avec GRASS / Convertir
le raster en lignes vectorielles" (r.to.vect.line).
Crez la couche vectorielle des cours deau.
Avec loutil de "conversion en zones vecteurs" (r.to.vect.area), convertissez la
couche des bassins versants. Affectez lui une symbologie sans remplissage.

Visualisation en 3D avec NVIZ

Loutil "Visionneur 3D (NVIZ)" permet reprsenter en trois dimensions un jeu de

donnes. Loutil permet de choisir un raster dlvation et un autre raster qui sera drap sur
ce raster dlvation. Une couche vectorielle peut galement tre ajoute la visualisation.
Lancez le visualisateur NVIZ avec le MNT sur lequel sera drap le raster doccupation
du sol. Utilisez galement la couche vecteur des rivires.
Explorez linterface pour vous dplacer dans le bloc 3D, changez lorientation du soleil,
ajoutez une flche du nord, etc.

Figure 3 Visualisation 3D avec NVIZ.