<figure><img src="image_1.png" /></figure>

cirad
AGRICULTURAL RESEARCH
FOR DEVELOPMENT

<figure><img src="image_2.png" /></figure>

tetis
TERRITOIRE ENVIRONNEMENT TÉLÉDÉTECTION
INFORMATION SPATIALE

# Prise en main du logiciel QGIS pour la mise en place et l'utilisation d'un Système d'Information Géographique
## QGIS V3.4.14

Annelise TRAN, Cirad UMR TETIS, Saint-Denis, Réunion, France

Stéphane DUPUY, Cirad UMR TETIS, Montpellier, France

<figure><img src="image_3.png" /></figure>

Ce(tte) œuvre est mise à disposition selon les termes de la [Licence Creative Commons Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/).
Ce document a été élaboré pour la formation d'initiation aux Systèmes d'Information Géographique (QGIS 3.4.14) des agents du Conseil Départemental de Mayotte.

© Cirad, Département Environnements et Sociétés, Unité Mixte de Recherche Territoires Environnement Télédétection et Information Spatiale, Direction régionale Réunion Mayotte, janvier 2020
# SOMMAIRE

## TABLE DES MATIERES

<table><tr><td>1</td><td>Prise en main de QGIS</td><td>5</td></tr><tr><td>1.1</td><td>Installation et démarrage</td><td>5</td></tr><tr><td>1.1.1</td><td>Téléchargement</td><td>5</td></tr><tr><td>1.1.2</td><td>Installation</td><td>5</td></tr><tr><td>1.1.3</td><td>Démarrage</td><td>6</td></tr><tr><td>1.2</td><td>Découverte de différents outils, création d'une première carte</td><td>6</td></tr><tr><td>1.2.1</td><td>Interface graphique</td><td>6</td></tr><tr><td>1.2.2</td><td>Création d'un projet</td><td>7</td></tr><tr><td>1.2.3</td><td>Ajouter une couche au format vecteur</td><td>7</td></tr><tr><td>1.2.4</td><td>Afficher / désafficher une couche</td><td>7</td></tr><tr><td>1.2.5</td><td>Naviguer dans la carte</td><td>7</td></tr><tr><td>1.2.6</td><td>Sélectionner et lire les informations d'une couche</td><td>8</td></tr><tr><td>2</td><td>Symbologie des couches géographiques</td><td>9</td></tr><tr><td>2.1</td><td>Modifier la symbologie d'une couche</td><td>9</td></tr><tr><td>2.1.1</td><td>Symbole unique</td><td>9</td></tr><tr><td>2.1.2</td><td>Catégories</td><td>9</td></tr><tr><td>2.1.3</td><td>Carte choroplèthe</td><td>10</td></tr><tr><td>2.1.4</td><td>Ajouter des étiquettes</td><td>10</td></tr><tr><td>2.2</td><td>Ajouts de fonds de cartes</td><td>11</td></tr><tr><td>2.2.1</td><td>Open Street Map</td><td>11</td></tr><tr><td>2.2.2</td><td>Google Earth</td><td>11</td></tr><tr><td>3</td><td>Création de nouvelles couches géographiques</td><td>12</td></tr><tr><td>3.1</td><td>Intégration de points GPS dans QGIS</td><td>12</td></tr><tr><td>3.1.1</td><td>A partir d'un fichier GPS (.gpx)</td><td>12</td></tr><tr><td>3.1.2</td><td>A partir d'un fichier texte (csv, txt, ...)</td><td>12</td></tr><tr><td>3.2</td><td>Création de nouvelles couches</td><td>13</td></tr><tr><td>3.2.1</td><td>Création de points</td><td>13</td></tr><tr><td>3.2.2</td><td>Création de lignes</td><td>14</td></tr><tr><td>3.2.3</td><td>Création de polygones</td><td>14</td></tr><tr><td>4</td><td>Mise en page</td><td>15</td></tr><tr><td>4.1</td><td>Ouvrir le composeur d'impression</td><td>15</td></tr><tr><td>4.2</td><td>Habillage de la carte</td><td>15</td></tr><tr><td>4.3</td><td>Exporter la carte</td><td>16</td></tr><tr><td>5</td><td>Tables attributaires</td><td>17</td></tr><tr><td>5.1</td><td>Requêtes attributaires</td><td>17</td></tr></table>
<table><tr><td>5.1.1</td><td>Sélection selon les attributs</td><td>17</td></tr><tr><td>5.1.2</td><td>Exporter les résultats d'une sélection</td><td>17</td></tr><tr><td>5.2</td><td>Calculer les valeurs d'un champ</td><td>18</td></tr><tr><td>5.3</td><td>Supprimer ou ajouter un champ</td><td>19</td></tr><tr><td>5.3.1</td><td>Ajouter un nouveau champ</td><td>19</td></tr><tr><td>5.3.2</td><td>Supprimer un champ</td><td>19</td></tr><tr><td>5.4</td><td>Joindre des tables attributaires</td><td>19</td></tr><tr><td>6</td><td>Analyse spatiale – introduction</td><td>22</td></tr><tr><td>6.1</td><td>Outils d'analyse spatiale format vecteur (exemples)</td><td>22</td></tr><tr><td>6.1.1</td><td>Outil Clip</td><td>22</td></tr><tr><td>6.1.2</td><td>Zone tampon (buffer)</td><td>23</td></tr><tr><td>6.2</td><td>Travailler avec le format raster</td><td>24</td></tr><tr><td>6.2.1</td><td>Ajouter une couche raster</td><td>24</td></tr><tr><td>6.2.2</td><td>Projeter un raster</td><td>24</td></tr><tr><td>6.2.3</td><td>Découper un raster</td><td>25</td></tr><tr><td>6.2.4</td><td>Calculer un raster de pente à partir d'un MNT</td><td>26</td></tr><tr><td>6.3</td><td>Outils d'analyse spatiale vecteur/raster (exemple)</td><td>27</td></tr></table>
# 1 PRISE EN MAIN DE QGIS

## 1.1 Installation et démarrage

### 1.1.1 Téléchargement

* QGIS website : https://qgis.org

* Page de téléchargement : https://qgis.org/fr/site/forusers/download.html

* La version courante est la version 3.10 'Coruna'

* La version stable (LTR : Long Term Release) est la version 3.4.14

<figure><img src="image_4.png" /></figure>

### 1.1.2 Installation

* Cliquer sur le fichier **QGIS-OSGeo4W-3.4.14-1-Setup-x86_64.exe** (version 64 bits) ou **QGIS-OSGeo4W-3.4.14-1-Setup-x86.exe** (version 32 bits) pour installer QGIS 3.4

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>
### 1.1.3 Démarrage

* Utiliser le raccourci bureau

<figure><img src="image_9.png" /></figure>

* Ou Démarrer / Programmes / QGIS 3.4 / QGIS Desktop 3.4.14

## 1.2 Découverte de différents outils, création d'une première carte

### 1.2.1 Interface graphique

<figure><img src="image_10.png" /></figure>

Personnaliser l'interface : clic droit dans la barre d'outils pour ajouter/supprimer de l'interface des outils ou des panneaux.
### 1.2.2 Création d'un projet

* Sauvegardez votre projet (par exemple sous le nom **exercice1.qgis**) à partir du Menu Projet / Sauvegarder ou de la barre d'outils « Projet »

<figure><img src="image_11.png" /></figure>

* **Recommandation : sauvegarder régulièrement votre projet**

* A tout moment vous pouvez fermer QGIS, le redémarrer et retrouver votre projet :

* à partir du Menu Projet / Ouvrir

* en cliquant sur l'icône ouvrir de la barre d'outils « Projets »

* depuis l'explorateur Windows en double-cliquant sur le nom du projet

### 1.2.3 Ajouter une couche au format vecteur

* Cliquer sur l'icône

<figure><img src="image_12.png" /></figure>

* Choisir le fichier **EXERCICE1 / COMMUNE.shp**

(source des données : Geofla®Communes édition 2016 Mayotte
https://geoservices.ign.fr/documentation/diffusion/telechargement-donnees-libres.html#geofla)

* Cliquer sur Ouvrir

* Ajouter d'autres couches vecteurs :

* **OSM_ROADS.shp**

* **OSM_WATERWAYS.shp**

(source des données Open Street Map : https://www.openstreetmap.org)

### 1.2.4 Afficher / désafficher une couche

Dans le panneau « Couches », cocher / décocher la case située devant le nom de la couche.

Ordre des couches : dans le panneau « Couches », vous pouvez modifier l'ordre d'apparence des différentes couches.

### 1.2.5 Naviguer dans la carte

Utiliser la Palette d'outils suivante :

<figure><img src="image_13.png" /></figure>

* Tester les différents outils :

* Déplacement libre

* Zoom avant, zoom arrière

* Zoom sur l'étendue de toutes les couches

* Zoom précédent / suivant
* Contrôler l'échelle et les coordonnées géographiques de la carte.

<figure><img src="image_14.png" /></figure>

### 1.2.6 Sélectionner et lire les informations d'une couche

Utiliser la barre d'outils « Attributs » :

<figure><img src="image_15.png" /></figure>

* **Identifier** : Sélectionner une couche dans le panneau « Couches », puis cliquer sur un objet pour faire apparaître les informations associées.

* Sélectionner graphiquement des objets géographiques

* Tester différentes méthodes de sélection

<figure><img src="image_16.png" /></figure>

* Déselectionner

<figure><img src="image_17.png" /></figure>

* Ouvrir la table attributaire d'une couche

<figure><img src="image_18.png" /></figure>

* Mesurer une distance / surface / angle

<figure><img src="image_19.png" /></figure>

* Créer un signet : les « signets » (书签) permettent de garder en mémoire une localisation pour pouvoir y retourner ensuite

<figure><img src="image_20.png" /></figure>

* Choisir une zone d'intérêt

* Cliquer sur

<figure><img src="image_21.png" /></figure>

* Entrer un nom

* Utiliser les signets : permet de retrouver les signets déjà créés (cliquer sur le signet pour retourner à la vue correspondante) ou de les supprimer

<figure><img src="image_22.png" /></figure>
# 2 SYMBOLOGIE DES COUCHES GEOGRAPHIQUES

Dans cet exercice on utilisera les mêmes couches de données géographiques que dans l'exercice 1.

## 2.1 Modifier la symbologie d'une couche

Depuis le panneau "Couches", faire un clic droit sur la couche pour laquelle on souhaite modifier son apparence dans la carte. Sélectionner en bas du menu contextuel "Propriétés" puis la rubrique "Style" dans la fenêtre "Propriétés" qui s'est ouverte.

Les choix de rendu s'adaptent dans la fenêtre en fonction du type d'entité de la couche : point, ligne ou polygone.

<figure><img src="image_23.png" /></figure>

### 2.1.1 Symbole unique

Tous les objets (=entités) sont représentés de la même manière : choisir dans la liste en haut de la fenêtre "symbole unique", puis paramétrer la couleur (éventuellement l'épaisseur du trait, la couleur du contour, le type de remplissage ...).

### 2.1.2 Catégories

Les entités ont un affichage différent en fonction des valeurs qualitatives dans la table attributaire associée : choisir "Catégorisé".

Sélectionner la colonne contenant les valeurs du champ, puis cliquer sur le bouton "Classer" et le bouton "Appliquer".

Toutes les valeurs d'un champ peuvent être différenciées. Ce mode peut être utilisé aussi bien pour des valeurs de type texte ou numérique.
* Choisissez « Catégorisé » pour représenter différemment chaque commune de Mayotte.

<figure><img src="image_24.png" /></figure>

### 2.1.3 Carte choroplèthe

Les entités ont un affichage différent en fonction des valeurs quantitatives dans la table attributaire associée : choisir "Gradué". "Gradué" ne s'emploie qu'avec des champs aux valeurs de type numérique.

Sélectionner d'abord la colonne, puis cliquer sur le bouton "Classer" et le bouton "Appliquer".

Vous pouvez aussi modifier les paramètres de classification et le rendu de chaque classe.

* Choisissez « Gradué » pour représenter la population pour chaque commune de Mayotte.

### 2.1.4 Ajouter des étiquettes

Depuis le panneau "Couches", faire un clic droit sur la couche / "Propriétés" / "Etiquettes".

Choisir « Etiquettes simples », puis pour « Etiqueter avec » indiquer le champ contenant les valeurs à afficher. [par exemple le nom des communes NOM_COM pour la couche COMMUNE.shp]

Modifier le style de la police et cliquer sur OK.
## 2.2 Ajouts de fonds de cartes

### 2.2.1 Open Street Map

Depuis le Menu Web / Quick Map Services / OSM standard

### 2.2.2 Google Earth

* Depuis le panneau Explorateur / XYZ Tuiles / Right-Click / Nouvelle connexion

* Ajouter le nom de la connexion : Google Satellite

* Ajouter l'URL (faire un copié/collé pour éviter les erreurs !) :
http://www.google.cn/maps/vt?lyrs=s@189&gl=cn&x={x}&y={y}&z={z}

* Click OK : une nouvelle connexion est disponible

* Double-cliquer ou faire un glissé/déposé pour ajouter la couche au projet

<figure><img src="image_25.png" /></figure>
# 3 CREATION DE NOUVELLES COUCHES GEOGRAPHIQUES

## 3.1 Intégration de points GPS dans QGIS

### 3.1.1 A partir d'un fichier GPS (.gpx)

Glisser-déposer le fichier depuis le panneau « explorateur » à la carte.

* Utiliser les fichiers .gpx créés lors de la formation à l'utilisation du GPS

### 3.1.2 A partir d'un fichier texte (csv, txt, ...)

Dans cet exercice nous utilisons la localisation des stations météorologiques de Mayotte (disponibles sur le site https://donneespubliques.meteofrance.fr).

* Cliquer sur l'icône 🔧 (ajouter une couche de texte délimité)

<figure><img src="image_26.png" /></figure>

* Choisir le fichier **EXERCICE3 / FICHE_STATIONS.csv**

* Indiquer le format du fichier : CSV

* Donner la définition de la géométrie : champs X et Y, format, système de coordonnées de référence

* Cliquer sur Ajouter puis Fermer

<figure><img src="image_27.png" /></figure>
# 3.2 Création de nouvelles couches

## 3.2.1 Création de points

### Création d'un nouveau fichier Shapefile

* Couche / Nouveau / Nouvelle couche vectorielle

* Ou cliquer sur l'icône Nouvelle couche vectorielle

* La fenêtre suivante s'ouvre

<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

* Nom de fichier : choisir un répertoire (/EXERCICE3) et le nom du fichier qui sera créé (par exemple : exercice_points.shp)

* Type de géométrie : choisir Points

* Indiquer le système de coordonnées de référence (au choix – par exemple UTM zone 38 S)

* Créer les champs de la table attributaire souhaités (par exemple NOM, format texte)

* Cliquer sur OK -> la nouvelle couche est créée et ajoutée à la carte

### Numérisation

* Utiliser un fonds de carte Google satellite (§2.2.2) et zoomer à l'échelle 1 :2500

* Utiliser la boîte à outil Numérisation

<figure><img src="image_30.png" /></figure>
* Cliquer sur l'outil 📉 et numériser la position de quelques maisons ou bâtiments

* Renseigner les champs de la table attributaire au fur et à mesure de la saisie, ou bien à partir de la table attributaire.

<figure><img src="image_31.png" /></figure>

### 3.2.2 Création de lignes

* Procéder de la même manière que précédemment (§3.2.1) avec le type de géométrie « Ligne »

<figure><img src="image_32.png" /></figure>

* Utiliser l'outil 📉 pour numériser quelques routes

<figure><img src="image_33.png" /></figure>

### 3.2.3 Création de polygones

* Procéder de la même manière que précédemment (§3.2.1) avec le type de géométrie « Polygone »

<figure><img src="image_34.png" /></figure>

* Utiliser l'outil 📉 pour numériser quelques grands bâtiments

<figure><img src="image_35.png" /></figure>
# 4 MISE EN PAGE

## 4.1 Ouvrir le composeur d'impression

* Projet / Nouveau composeur d'impression / Indiquer un Titre pour la nouvelle composition

* Icône

<figure><img src="image_36.png" /></figure>

* La fenêtre du composeur d'impression s'ouvre :

<figure><img src="image_37.png" /></figure>

## 4.2 Habillage de la carte

* Tester les différents outils

* Ajouter une nouvelle carte

<figure><img src="image_38.png" /></figure>

* Ajouter une image

<figure><img src="image_39.png" /></figure>

* Ajouter du texte (par exemple : titre, copyright, etc).

<figure><img src="image_40.png" /></figure>

* Ajouter une légende

<figure><img src="image_41.png" /></figure>
* Ajouter l'échelle

<figure><img src="image_42.png" /></figure>

* Ajouter une forme

<figure><img src="image_43.png" /></figure>

* Ajouter une flèche

<figure><img src="image_44.png" /></figure>

* Ajouter une table d'attribut

<figure><img src="image_45.png" /></figure>

* Ajouter du HTML

<figure><img src="image_46.png" /></figure>

* Outils de zoom :

<figure><img src="image_47.png" /></figure>

* Outils pour verrouiller, grouper / dégrouper, aligner/répartir les éléments de la carte :

<figure><img src="image_48.png" /></figure>

### 4.3 Exporter la carte

* Exporter votre carte dans différents formats : PDF, image, SVG

<figure><img src="image_49.png" /></figure>
# 5 TABLES ATTRIBUTAIRES

## 5.1 Requêtes attributaires

### 5.1.1 Sélection selon les attributs

Dans cet exercice il s'agit de sélectionner les communes de plus de 10000 habitants :

* Ajouter la couche COMMUNES.shp

* Ouvrir la table attributaire de la couche commune

* Cliquer sur l'icône

<figure><img src="image_50.png" /></figure>

<figure><img src="image_51.png" /></figure>

* Dans l'espace « Expression », entrer la requête souhaitée (utiliser les champs et opérateurs proposés dans le bloc central)

* Cliquer sur Sélection

<figure><img src="image_52.png" /></figure>

### 5.1.2 Exporter les résultats d'une sélection

Dans cet exercice il s'agit de sélectionner les routes de type « primaire » et « secondaire » et d'exporter les résultats de la sélection dans un nouveau fichier shapefile : routes_principales.shp

* Ajouter la couche OSM_ROADS.shp

* Ouvrir la table attributaire de la couche Routes

* Cliquer sur l'icône

<figure><img src="image_53.png" /></figure>

<figure><img src="image_54.png" /></figure>

* Dans l'espace « Expression », entrer la requête souhaitée (utiliser les champs et opérateurs proposés dans le bloc central)
<figure><img src="image_55.png" /></figure>

* Cliquer sur Sélection

* Exporter les résultats : Click-droit sur la couche OSM_ROADS / Exporter / Enregistrer les entités sélectionnées sous...

## 5.2 Calculer les valeurs d'un champ

Dans cet exercice il s'agit de calculer la densité de population (nb d'habitants / km²) en fonction des valeurs contenues dans les champs « POPULATION » (nombre d'habitants par commune) et « SUPERFICIE » (superficie de la commune, en hectares) de la table attributaire de la couche COMMUNES.shp

* Ouvrir la table attributaire de la couche COMMUNES.shp

* Cliquer sur l'icône (calculatrice de champ)

* Cocher la case *Créer un nouveau champ*

* Définir le nom et le type du nouveau champ (par exemple, DENSITEPOP, Integer)

* Entrer l'expression du calcul dans la zone d'expression

* Rappel $1 \text{ km}^2 = 100 \text{ ha}$

* Cliquer sur OK

* Visualiser la densité de population par commune à Mayotte en créant une carte choroplète (§2.1.3)
<figure><img src="image_56.png" /></figure>

## 5.3 Supprimer ou ajouter un champ

### 5.3.1 Ajouter un nouveau champ

Avec l'outil

<figure><img src="image_57.png" /></figure>

Attention : le type de champ (nombre entier, nombre décimal, texte, date) est défini à la création du champ et ne peut pas être modifié ensuite.

### 5.3.2 Supprimer un champ

Avec l'outil

<figure><img src="image_58.png" /></figure>

## 5.4 Joindre des tables attributaires

Des données de recensement de population plus récentes (2017) sont disponibles sur le site de l'INSEE (https://www.insee.fr/fr/statistiques/3291775?sommaire=2120838).

Dans cet exercice nous allons effectuer une jointure entre la table attributaire de la couche COMMUNES.shp, et la table POPULATION2017.csv (table seule) qui contient les nouvelles données de recensement.

* Cliquer sur l'icône 📄 (ajouter une couche de texte délimité)

<figure><img src="image_59.png" /></figure>

* Choisir le fichier EXERCICE5/POPULATION2017.csv

* Indiquer le format du fichier : CSV

* Choisir : pas de géométrie (juste la table)
<figure><img src="image_60.png" /></figure>

* Résultat : la table POPULATION2017 est importée dans le projet comme une table seule (couche non géographique)

<figure><img src="image_61.png" /></figure>

POPULATION2017

* La table attributaire cible est la table attributaire de la couche COMMUNES.shp (couche géographique)

<figure><img src="image_62.png" /></figure>

<figure><img src="image_63.png" /></figure>

COMMUNE

* Ouvrir les deux tables attributaires et identifier le nom du champ-clef commun aux deux tables :

* 0 : 1

* 0 : 1

* Ouvrir la fenêtre des Propriétés de la couche Communes (click-droit / Propriétés) et sélectionner le Menu Jointure

* Cliquer sur (+) (ajouter une jointure)

* Définir les propriétés de la jointure

* 0 : La table attributaire à joindre

* 0 : Le champ clef dans la table attributaire à joindre

* 0 : Le champ clef dans la table attributaire cible

Important : les champs clefs doivent être de même type.
<figure><img src="image_64.png" /></figure>

* Cliquer sur OK / Appliquer

* Visualiser le résultat de la jointure :

* Dans la table attributaire des communes

* En réalisant une carte choroplète de la population par commune en 2017
# 6 ANALYSE SPATIALE – INTRODUCTION

## 6.1 Outils d'analyse spatiale format vecteur (exemples)

### 6.1.1 Outil Clip

* Dans cet exercice nous allons « découper » (« clipper ») la couche du réseau hydrographique (OSM WATERWAY) selon les contours de la commune de Mamoudzou

* Créer un nouveau projet et ajouter les couches EXERCICE6/MAMOUDZOU.shp et EXERCICE6/OSM WATERWAY.shp

* Ouvrir l'outil du Menu Vecteur / Outils de géotraitement / Découper

* Renseigner :

* la couche en entrée

* la couche de découpe

* (option) le nom du fichier en sortie

<figure><img src="image_65.png" /></figure>

* Visualiser le résultat

<figure><img src="image_66.png" /></figure>
### 6.1.2 Zone tampon (buffer)

* Dans cet exercice nous allons créer une zone tampon de 100 mètres de large autour des rivières

* Ouvrir l'outil du Menu Vecteur / Outils de géotraitement / Zone tampon

* Renseigner :

* La couche en entrée (la couche géographique avec les rivières de la commune de Mamoudzou créée précédemment)

* La distance tampon (50 m pour avoir une zone de 100 m de large)

<figure><img src="image_67.png" /></figure>

* Visualiser le résultat

<figure><img src="image_68.png" /></figure>
# 6.2 Travailler avec le format raster

## 6.2.1 Ajouter une couche raster

* Cliquer sur l'icône

<figure><img src="image_69.png" /></figure>

* Choisir le fichier **EXERCICE6/ srtm_46_15.tif**

Il s'agit d'un modèle numérique de terrain (MNT) au format raster : la valeur d'un pixel correspond à son altitude, en mètres (source des données : http://srtm.csi.cgiar.org/srtmdata/)

* Ouvrir les propriétés du raster (click-droit/Propriétés) et contrôler le système de coordonnées de la couche (Onglet Information) : ....................................................

* Modifier la symbologie du raster

## 6.2.2 Projeter un raster

Le MNT a pour système de coordonnées le Système de coordonnées géographiques WGS84. Afin de pouvoir découper le MNT selon les contours de Mayotte (à partir de la couche des communes, qui est en coordonnées projetées UTM38S), il faut d'abord créer un nouveau raster dans la même projection (UTM 38 S).

* Ouvrir le Menu Raster / Projections / Projection

* Sélectionner le raster en entrée : srtm_46_15.tif

* Indiquer les systèmes de coordonnées source (géographique WGS84) et cible (UTM 38 S)

* Indiquer la résolution spatiale du raster de sortie (90 m)

* Indiquer le nom du fichier raster en sortie (par ex. DEM_UTM38S)

<figure><img src="image_70.png" /></figure>

* Visualiser le résultat
### 6.2.3 Découper un raster

Dans cet exercice il s'agit de découper le MNT selon les contours de Mayotte.

* Ajouter la couche COMMUNES.shp

* Menu Raster / Extraction / Découper Raster / Découper Raster selon couche de masque

* Raster en entrée : le DEM en coordonnées UTM 38 S créé précédemment

* Couche de masque : la couche COMMUNES

* Indiquer le nom de la couche raster en sortie (par ex. DEM_MAYOTTE)

<figure><img src="image_71.png" /></figure>

* Visualiser le résultat

<figure><img src="image_72.png" /></figure>
### 6.2.4 Calculer un raster de pente à partir d'un MNT

* Menu Raster / Analyse / Pente

<figure><img src="image_73.png" /></figure>
## 6.3 Outils d'analyse spatiale vecteur/raster (exemple)

### Statistiques zonales

Dans cet exercice il s'agit de calculer l'altitude moyenne pour chaque commune de Mayotte

* Ouvrir la boîte à outils de traitement > Menu Analyse raster > Statistiques zonales

* Raster en entrée : le MNT de Mayotte

* La couche vecteur avec le zonage : les communes

* Préfixe de la colonne en sortie : DEM_

* Statistique à calculer : moyenne, minimum, maximum, écart-type... au choix

<figure><img src="image_74.png" /></figure>

* Visualiser les résultats : dans la table attributaire de la couche « COMMUNE », des champs ont été ajoutés.
<figure><img src="image_75.png" /></figure>

cirad
LA RECHERCHE AGRONOMIQUE
POUR LE DÉVELOPPEMENT

Département Environnements et Sociétés, Unité Mixte de Recherche Territoires Environnement
Télédétection et Information Spatiale, Direction régionale Réunion Mayotte

Maison de la Télédétection, 500 rue Jean-François Breton, 34093 Montpellier cedex 5

www.cirad.fr