Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
introsig [2018/03/16 15:38]
qcbs [Exercice 6 - Utiliser GRASS dans QGIS]
introsig [2019/02/28 17:34] (current)
qcbs [Exercice 2 - Tampons et analyses simples]
Line 6: Line 6:
 Guillaume Larocque, Professionnel de recherche. (guillaume.larocque@mcgill.ca) Guillaume Larocque, Professionnel de recherche. (guillaume.larocque@mcgill.ca)
  
-et mars 2017, Université ​de Montréal+29 février ​et 1er mars 2019, Université ​McGill
  
 ===== Présentation ===== ===== Présentation =====
Line 15: Line 15:
 ===== 1 - Installation de QGIS et GRASS ===== ===== 1 - Installation de QGIS et GRASS =====
  
-QGIS 3 et GRASS 7 seront utilisés pour cet atelier. Veuillez vous assurer d'​avoir les logiciels proprement installés sur votre ordinateur avant l'​atelier. ​+QGIS 3.4 et GRASS 7 seront utilisés pour cet atelier. Veuillez vous assurer d'​avoir les logiciels proprement installés sur votre ordinateur avant l'​atelier. ​
  
 Installez simplement QGIS et GRASS avec le programme d'​installation [[http://​qgis.org/​en/​site/​forusers/​download.html|disponible ici]]. ​ Installez simplement QGIS et GRASS avec le programme d'​installation [[http://​qgis.org/​en/​site/​forusers/​download.html|disponible ici]]. ​
Line 61: Line 61:
   * [[http://​www.saga-gis.org/​|SAGA GIS]] - (System for Automated Geoscientific Analyses) Un SIG qui a de nombreuses fonctions pour l'​analyse de données spatiales avancée. ​   * [[http://​www.saga-gis.org/​|SAGA GIS]] - (System for Automated Geoscientific Analyses) Un SIG qui a de nombreuses fonctions pour l'​analyse de données spatiales avancée. ​
   * [[http://​52north.org/​communities/​ilwis|ILWIS Open]] - Un SIG et logiciel pour la télémétrie qui intègre des données matricielles,​ vectorielles et thématiques. ​   * [[http://​52north.org/​communities/​ilwis|ILWIS Open]] - Un SIG et logiciel pour la télémétrie qui intègre des données matricielles,​ vectorielles et thématiques. ​
-  * [[http://​www.gvsig.org/|gvSIG]] - Un logiciel qui est présentement en développement qui a un interface convivialqui accepte les formats de fichiers les plus communs et qui offre de nombreux outils ​pour effectuer des requêtes, créer des modèles, faire des analyses, traiter des réseaux, etc+  * [[http://​www.uoguelph.ca/​~hydrogeo/​Whitebox/|Whitebox GAT]] - Powerful software for geo-spatial analysismodeling and remote sensing.  
 +  * [[https://​www.orfeo-toolbox.org/​|ORFEO toolbox]] - Plate-forme ouverte ​pour le traitement d'​images issues de la télédétection
   * [[http://​www.openjump.org/​|OpenJUMP GIS]] - La version actuelle de OpenJUMP peut lire et écrire des shapefiles et autres fichiers vectoriels, et offre peu de possibilités pour le traitement d'​images. Dans son état actuel, ce logiciel peut être utilisé pour visualiser des données SIG ou pour l'​édition simple de la géométrie et des attributs de données vectorielles. ​   * [[http://​www.openjump.org/​|OpenJUMP GIS]] - La version actuelle de OpenJUMP peut lire et écrire des shapefiles et autres fichiers vectoriels, et offre peu de possibilités pour le traitement d'​images. Dans son état actuel, ce logiciel peut être utilisé pour visualiser des données SIG ou pour l'​édition simple de la géométrie et des attributs de données vectorielles. ​
  
Line 94: Line 95:
 Sauvez-le ensuite avec un extension .csv. Dans QGIS, ajoutez ce fichier comme couche vectorielle en sélectionnant l'​option '​Ajoutez une couche de texte délimitée'​. Utilisez Longitude pour la colonne des X et Latitude pour la colonne des Y et indiquez que le fichier est séparé par des virgules. Dans la fenêtre suivante, spécifiez Système de Coordonnée geographique>​WGS 84 comme SCR. Vous devriez maintenant voir apparaître deux points indiquant la localisation de la Réserve Molson et de l'​Arboretum Morgan. ​ Sauvez-le ensuite avec un extension .csv. Dans QGIS, ajoutez ce fichier comme couche vectorielle en sélectionnant l'​option '​Ajoutez une couche de texte délimitée'​. Utilisez Longitude pour la colonne des X et Latitude pour la colonne des Y et indiquez que le fichier est séparé par des virgules. Dans la fenêtre suivante, spécifiez Système de Coordonnée geographique>​WGS 84 comme SCR. Vous devriez maintenant voir apparaître deux points indiquant la localisation de la Réserve Molson et de l'​Arboretum Morgan. ​
  
-**Étape 2**: Ajoutez une couche '​Google Satellite'​ en Cliquant sur Explorateur (en bas à gauche) et sur XYZ Tiles. Cliquez sur le bouton de droite et sur Nouvelle connexion. Pour Nom, mettre '​Google satellite'​. Pour URL mettrehttp://​mt0.google.com/​vt/​lyrs=y&​hl=en&​x={x}&​y={y}&​z={z}&​s=Ga ​Notez que le SCR du canevas est maintenant '​Google Mercator'​ (WGS 84 / Pseudo mercator). Déplacez cette couche en dessous du canevas en cliquant-glissant le nom de la couche vers le bas dans la table des couches à gauche. Autre source de couches (https://​leaflet-extras.github.io/​leaflet-providers/​preview/​.+**Étape ​2**: Cliquez sur Extensions et ajoutez l'​extension QuickMapServices. Activez cette extension en la cochant dans la liste des extensions disponibles. Ensuite, cliquez sur Internet>​QuickMapServices>​Google>​Google Hybrid pour ajouter un fond de carte Google.  
 + 
 +**Alternative à l'​étape ​2**: Ajoutez une couche '​Google Satellite'​ en Cliquant sur Explorateur (en bas à gauche) et sur XYZ Tiles. Cliquez sur le bouton de droite et sur Nouvelle connexion. Pour Nom, mettre '​Google satellite'​. Pour URL mettre ​ 
 + 
 +<​file>​ 
 +http://​mt0.google.com/​vt/​lyrs=y&​hl=en&​x={x}&​y={y}&​z={z}&​s=Ga 
 +</​file>​ 
 + 
 +Notez que le SCR du canevas est maintenant '​Google Mercator'​ (WGS 84 / Pseudo mercator). Déplacez cette couche en dessous du canevas en cliquant-glissant le nom de la couche vers le bas dans la table des couches à gauche. Autre source de couches (https://​leaflet-extras.github.io/​leaflet-providers/​preview/​).
  
 **Étape 3**: Cliquez sur le menu de changement de SCR en bas à droite. Changez le SCR pour NAD83/UTM Zone 18N.  **Étape 3**: Cliquez sur le menu de changement de SCR en bas à droite. Changez le SCR pour NAD83/UTM Zone 18N. 
  
-**Étape 4**: Numérisez la Réserve Molson et l'​Arboretum Morgan. Vous devez d'​abord ajouter une nouvelle couche ​shapefile ​de type polygone (SCR NAD 83/UTM 18N) et ajouter ​un nouveau champ pour le nom (Donnée Texte, ​largeur ​80). Donnez à ce fichier ​shapefile ​un nom approprié (e.g. Reserves_McGill). Activer le mode édition (­>​Couche­­­­>​­Basculer en mode édition)­­­­­­­­­­ et numérisez les deux réserves en cliquant sur l'​icône {{::​capture_du_2013-09-25_10_58_40.png|}}. Pour vous aider à numériser les deux réserves, utilisez les points que vous avez importés à partir du fichier CSV. Pour la réserve Molson, numérisez la forêt située entre le boulevard résidentiel (Boul Perrot) et l'​autoroute 20. Pour l'​Arboretum,​ numérisez la région boisée située autour du point central. L'​Arboretum fait approximativement 250 hectares. Ajoutez maintenant des polygones adjacents pour indiquer de quelle façon vous prévoyez étendre les réserves. Pour ce faire, il est important d'​activer l'​option d'​accrochage ​(toggle snapping, ​icône d'​aimant rougeet de spécifier un tolérance de 20 m. +**Étape 4**: Numérisez la Réserve Molson et l'​Arboretum Morgan. Vous devez d'​abord ajouter une nouvelle couche ​Geopackage ​de type polygone (Couche>​Créer une couche>​Créer une couche Geopackage) avec un SCR NAD 83/UTM 18N. Ajouter ​un nouveau champ pour le nom (Donnée Texte, ​longueur ​80). Donnez à ce fichier ​Geopackage ​un nom approprié (e.g. Reserves_McGill). Activer le mode édition (­>​Couche­­­­>​­Basculer en mode édition)­­­­­­­­­­ et numérisez les deux réserves en cliquant sur l'​icône {{::​capture_du_2013-09-25_10_58_40.png|}}. Pour vous aider à numériser les deux réserves, utilisez les points que vous avez importés à partir du fichier CSV. Pour la réserve Molson, numérisez la forêt située entre le boulevard résidentiel (Boul Perrot) et l'​autoroute 20. Pour l'​Arboretum,​ numérisez la région boisée située autour du point central. L'​Arboretum fait approximativement 250 hectares. Ajoutez maintenant des polygones adjacents pour indiquer de quelle façon vous prévoyez étendre les réserves. Pour ce faire, il est important d'​activer l'​option d'​accrochage ​en cliquant sur l'icône d'​aimant rouge et de spécifier un tolérance de 20 m. Si vous ne voyez pas l'​aimant rouge, activez la barre d'​outils de numérisation sous Vue>​Barre d'​outils
  
 **Étape 5**: Sortez du mode édition en cliquant sur l'​icône de crayon et sauvez cette couche. ​ **Étape 5**: Sortez du mode édition en cliquant sur l'​icône de crayon et sauvez cette couche. ​
Line 104: Line 113:
 **Étape 6**: Désactivez la couche Google Satellite. ​ **Étape 6**: Désactivez la couche Google Satellite. ​
  
-**Étape 7**: À partir du menu '​Propriétés de la couche',​ changez les couleurs des deux réserves à partir de l'​onglet ​Style. Choisissez '​Catégorisé' ​et spécifiez la colonne contenant les noms des réserves avant de cliquer sur Classer. Depuis l'​onglet étiquettes,​ ajoutez ensuite des étiquettes spécifiant le nom des réserves et définissez des étiquettes pour les extensions des réserves que vous proposez. ​+**Étape 7**: À partir du menu '​Propriétés de la couche',​ changez les couleurs des deux réserves à partir de l'​onglet ​Symbologie. Choisissez '​Catégorisé'​spécifiez la colonne contenant les noms des réserves, et choisissez une Palette "​Random colors", ​avant de cliquer sur Classer. Cliquez sur la case de couleur à côté de chacune des réserves pour choisir une couleur à votre goût. Depuis l'​onglet étiquettes,​ ajoutez ensuite des étiquettes ​(simples) ​spécifiant le nom des réserves et définissez des étiquettes pour les extensions des réserves que vous proposez.  
 + 
 +**Étape 8**: Ajoutez-y les routes (Routes.shp),​ les autres régions boisées (region_boise.shp) et les étendues d'eau (Region_hydrique.shp) à partir du fichier .zip. Définissez des couleurs et des styles de lignes appropriées pour chaque couche en utilisant l'​onglet Symbologie dans le menu >​Propriétés de la couche.  
 + 
 +**Étape 9**: Cliquez sur Projet>​Nouvelle mise en page. Créez une carte qui montre les réserves que vous avez numérisés,​ les routes et les cours d'eau. Note: pour ajouter la carte sur la page blanche, vous devez cliquer sur l'​icône '​Ajouter une nouvelle carte' et cliquer-glisser sur la feuille blanche pour définir l'​étendue de la carte sur la feuille. 
 + 
 +=-=-=-=-=-=-=-=-=-=-=-=- 
 + 
 +**DÉFI 1**: Ajoutez les pipelines et emprises de lignes hydro-électriques (powerlines) en utilisant le plugin Quick OSM
  
-**Étape 8**: Ajoutez-y les routes ​(Routes.shp), les autres régions boisées (region_boise.shp) et les étendues d'eau (Region_hydrique.shp). Définissez des couleurs ​appropriées pour chaque couche en utilisant l'​onglet Style dans le menu >​Propriétés ​de la couche+**DÉFI 2**: Ajoutez ​à votre carte les milieux humides en format Shapefile téléchargés depuis le portail de données ouvertes de la ville Montréal ​(http://​donnees.ville.montreal.qc.ca/). Mettre une palette de couleurs ​qui est fonction de la grandeur ​de la zone humide
  
-**Étape 9**:  ​Créez une jolie carte avec le Composeur d'​impression,​ en ajoutant une Flèche indiquant le nord, des étiquettes et une légende. Sauvegardez cette carte en fichier PNG+**DÉFI 3**: Mettre ​le fond de carte suivant: https://​leaflet-extras.github.io/​leaflet-providers/​preview/#​filter=Stamen.Watercolor.
  
-**Étape 9**: Cliquez sur Projet>​Nouvelle mise en page. Créez une carte qui montre les réserves que vous avez numérisés. Note: pour ajouter la carte sur la page blanche, vous devez cliquer sur l'icône '​Ajouter une nouvelle carte' et cliquer-glisser sur la feuille blanche pour définir l'​étendue de la carte sur la feuille.+**DÉFI 4**: Visualisez votre carte dans un fureteur ​en utilisant ​l'application gis2web
  
-**DÉFI**: Ajoutez les pipelines et emprises de lignes hydro-électriques (powerlines),​ que vous aurez d'​abord téléchargé depuis Open Street Maps en format .osm.  
 ====== Exercice 2 - Tampons et analyses simples ====== ====== Exercice 2 - Tampons et analyses simples ======
  
 **Objectif: Trouver la proportion de voteurs vivant à moins de 500 m d'une région boisée plus grande que 10 hectares dans le district électoral de Jacques-Cartier dans l'​Ouest de l'île de Montréal. ** **Objectif: Trouver la proportion de voteurs vivant à moins de 500 m d'une région boisée plus grande que 10 hectares dans le district électoral de Jacques-Cartier dans l'​Ouest de l'île de Montréal. **
  
-**Étape 1**: Démarrez un nouveau projet et ajoutez ​les couches nommées ​'​sections_vote_31h5.shp' et '​regions_boise.shp'. Convertissez ​ces fichiers ​en NAD 83/UTM 18N en les sauvant sous de nouveaux fichiers (­Bouton de droite de la souris sur le nom de la couche>​Sauvegarder sous... et sélectionnez ​le SCR ). +**Étape 1**: Démarrez un nouveau projet et ajoutez ​la couche nommée ​'​sections_vote_31h5.shp'​. Convertissez ​ce fichier ​en NAD 83/UTM 18N en les sauvant sous de nouveaux fichiers (­Bouton de droite de la souris sur le nom de la couche>​Exporter>​Sauvegarder ​les entités ​sous... ​Sélectionnez le format Geopackage ​et appelez ​le fichier '​Exercice 1' et cliquez sur les (...pour choisir un endroit approprié pour le mettre sur votre ordinateur. Mettre un nom de couche approprié et spécifier le bon SCR. Répétez ces étapes pour ajouter la couche '​region_boisés.shp'​ au même fichier Geopackage dans le même SCR.
  
 **Étape 2**: Spécifiez NAD83/UTM 18N pour le SCR du canevas. Ajoutez-y les deux nouveaux fichiers que vous avez sauvés. Enlevez le fichiers en WGS 84.  **Étape 2**: Spécifiez NAD83/UTM 18N pour le SCR du canevas. Ajoutez-y les deux nouveaux fichiers que vous avez sauvés. Enlevez le fichiers en WGS 84. 
  
-**Étape 3**: Ouvrez la table d'​attribut de la couche régions boisées (en cliquant sur le nom de la couche... table d'​attributs) et ajoutez-y une colonne spécifiant la superficie de chaque région en hectares. Pour ce faire, passez en mode édition, utilisez la calculatrice de champs (icône de calculatrice en bas) et utilisez l'​opérateur Géométrie->​$area. Note: le système de coordonnées UTM est en mètres et 1 hectare= 10 000 mètres carrés. Important: spécifiez Nombre décimal (Réel) et augmentez la précision à 2. Sortez du mode édition. ​+**Étape 3**: Ouvrez la table d'​attribut de la couche régions boisées (en cliquant sur le nom de la couche... table d'​attributs) et ajoutez-y une colonne spécifiant la superficie de chaque région en hectares. Pour ce faire, passez en mode édition ​(icône crayon), utilisez la calculatrice de champs (icône de calculatrice en haut) et utilisez l'​opérateur Géométrie->​$area. Note: le système de coordonnées UTM est en mètres et 1 hectare= 10 000 mètres carrés. Important: spécifiez Nombre décimal (Réel) et augmentez la précision à 2. Sortez du mode édition. ​
  
 **Étape 4**: Effectuez une requête (>​Couche>​Filtrer) pour isoler les surfaces boisées de plus de 10 hectares. Notez que vous ne pouvez pas effectuer de requêtes en mode édition. ​ **Étape 4**: Effectuez une requête (>​Couche>​Filtrer) pour isoler les surfaces boisées de plus de 10 hectares. Notez que vous ne pouvez pas effectuer de requêtes en mode édition. ​
  
-**Étape 5**: Créez un nouveau fichier ​shapefile ​avec un tampon de 500 m entourant ces surfaces boisées (>​Vecteurs>​Outils ​de géotraitement>​Tampons). +**Étape 5**: Créez un nouveau fichier ​Geopackage ​avec un tampon de 500 m entourant ces surfaces boisées (Boîte à outils ​de géotraitement>​Géométrie vectorielle>​Tampon). Spécifiez le champs distance et laissez les autres champs aux valeurs par défaut. Dans le champs Mise en Tampon, cliquez sur [...] et Enregistrer dans le fichier Geopackage créé plus tôt (Exercice1).  
  
-**Étape 6**: Ajouter une colonne définissant l'aire de chaque polygone dans la table d'​attributs de sections_vote_31h5. Ceci nous donne la superficie totale de chaque zone de vote. +**Étape 6**: Ajouter une colonne définissant l'aire de chaque polygone dans la table d'​attributs de sections_vote_31h5. Ceci nous donne la superficie totale de chaque zone de vote en mètres
  
 **Étape 7**: Effectuez une requête (Filtre) pour ne sélectionner que le district électoral Jacques-Cartier (colonne TRI_CEP). ​ **Étape 7**: Effectuez une requête (Filtre) pour ne sélectionner que le district électoral Jacques-Cartier (colonne TRI_CEP). ​
Line 133: Line 149:
 **Étape 8**: Utilisez la fonction Vecteur>​Outils d'​analyse>​Statistiques basiques pour les champs, pour trouver la somme de la colonne 2008_12_08.  ​ **Étape 8**: Utilisez la fonction Vecteur>​Outils d'​analyse>​Statistiques basiques pour les champs, pour trouver la somme de la colonne 2008_12_08.  ​
  
-**Étape 9**: Utilisez la fonction '​Vecteur...Outils de géotraitement...Coupez (Clip)'​ pour générer une couche vecteur ne contenant que la partie des sections de vote du district électoral Jacques-Cartier qui se situent à l'​intérieur du tampon de 500 mètres. Notez que la couche source doit contenir les sections de votez, et le district ​est la couche de découpage+**Étape 9**: Utilisez la fonction '​Vecteur...Outils de géotraitement...Coupez (Clip)'​ pour générer une couche vecteur ne contenant que la partie des sections de vote du district électoral Jacques-Cartier qui se situent à l'​intérieur du tampon de 500 mètres. Notez que la couche source doit contenir les sections de votes, et le tampon ​est la couche de superposition
  
 **Étape 10**: Ajoutez une nouvelle colonne dans la table d'​attributs contenant l'aire de chaque zone de vote après le découpage, donc qui est située à moins de 500 mètres d'un forêt. ($area sur le résultat de l'​étape 9).  **Étape 10**: Ajoutez une nouvelle colonne dans la table d'​attributs contenant l'aire de chaque zone de vote après le découpage, donc qui est située à moins de 500 mètres d'un forêt. ($area sur le résultat de l'​étape 9). 
  
-**Étape 11**: Afin d'​estimer le nombre de voteurs à l'​intérieur du tampon de 500 m, nous allons supposer que le nombre de voteurs est proportionnel à la superficie de chaque zone et calculer la proportion de chaque zone qui est à l'​intérieur du tampon. Ouvrez la table d'​attributs de la couche shapefile créée à l'​étape 9. Ajoutez une nouvelle colonne contenant le nombre de voteurs (Colonne 2008_12_08) multiplié par la proportion de l'​aire ​ de chaque zone  située à l'​intérieur du tampon. ​+**Étape 11**: Afin d'​estimer le nombre de voteurs à l'​intérieur du tampon de 500 m, nous allons supposer que le nombre de voteurs est proportionnel à la superficie de chaque zone et calculer la proportion de chaque zone qui est à l'​intérieur du tampon. Ouvrez la table d'​attributs de la couche shapefile créée à l'​étape 9. Ajoutez une nouvelle colonne contenant le nombre de voteurs (Colonne 2008_12_08) multiplié par la proportion de l'​aire ​ de chaque zone située à l'​intérieur du tampon. ​
  
 <​code>​ <​code>​
Line 149: Line 165:
 ++++ ++++
  
-**DÉFIS** : Trouvez quel district électoral (Colonne TRI_CEP) sur cette carte contient la plus grande longueur cumulée de routes (fichier '​routes.shp'​). Utilisez les fonctions '​Regrouper'​ et 'Somme des longueurs de lignes'​ (dans le menu >​Vecteurs). ​+=-=-=-=-=-=-=-=-=-=- 
 + 
 +**DÉFI 1** : Trouvez quel district électoral (Colonne TRI_CEP) sur cette carte contient la plus grande longueur cumulée de routes (fichier '​routes.shp'​). Utilisez les fonctions '​Regrouper'​ et 'Somme des longueurs de lignes'​ (dans le menu >​Vecteurs). ​
  
 ++++ Réponse | ++++ Réponse |
 Chateauguay,​ avec 587.294km de routes. ​ Chateauguay,​ avec 587.294km de routes. ​
 +++++
 +
 +
 +**DÉFI 2** : Trouvez la superficie totale de milieux humides de type marécage dans la section de votre no.8124. ​
 +
 +++++ Réponse |
 +667446 m2
 ++++ ++++
  
Line 169: Line 194:
 **Étape 2**: Ajoutez le fichier CSV nommé '​BBS_Routes_QC.csv'​ au canevas en utilisant la fonction '​Ajouter une couche de texte délimitée'​ (Choisissez '​Délimiteurs sélectionnés-­>​virgule) et spécifiez le CRS NAD83 (Latitude, Longitude- Géographique). ​ **Étape 2**: Ajoutez le fichier CSV nommé '​BBS_Routes_QC.csv'​ au canevas en utilisant la fonction '​Ajouter une couche de texte délimitée'​ (Choisissez '​Délimiteurs sélectionnés-­>​virgule) et spécifiez le CRS NAD83 (Latitude, Longitude- Géographique). ​
  
-**Étape 3**: Sauvegardez cette couche comme fichier ​shapefile '​BBS_Routes_QC.shp' avec le SCR: NAD83 Quebec Lambert. ​+**Étape 3**: Sauvegardez cette couche comme couches dans un fichier ​Geopackage créé pour l'exercice 3 avec le SCR: NAD83 Quebec Lambert. ​
  
-**Étape 4**: Ajoutez le fichier '​province.shp'​ au canevas. Sélectionnez la province de Québec et sauvegardez le polygone sous un nouveau shapefile ​avec le SCR Quebec Lambert (NAD83). ​+**Étape 4**: Ajoutez le fichier '​province.shp'​ au canevas. Sélectionnez la province de Québec et sauvegardez le polygone sous un nouvelle couche Geopackage ​avec le SCR Quebec Lambert (NAD83). ​
  
-**Étape 5**: Démarrez un nouveau projet et spécifiez ​le SCR Quebec Lambert (NAD83). ​Ajoutez-y ​les fichiers des étapes 3 et 4+**Étape 5**: Changez ​le SCR du canevas pour Quebec Lambert (NAD83). ​Enlevez ​les couches en format Shapefile WGS 84 
  
-**Étape 6**: Ajoutez la table '​BBS_Ovenbird_QC.csv'​ au canevas en utilisant '>​Couche>​Ajoutez une couche vecteur'​ et en listant tous les types de fichier (oui, cette étape est contre-intuitive!). ​+**Étape 6**: Ajoutez la table non-spatiale qui contient les données d'​abondance de la paruline ​'​BBS_Ovenbird_QC.csv'​ au canevas en utilisant '>​Couche>​Ajoutez une couche vecteur'​ et en listant tous les fichier (oui, cette étape est contre-intuitive!). ​
  
-**Étape 7**: À partir du menu '​Propriétés de la couche'​ de la table du fichier ​'​BBS_Routes_QC.shp', trouvez l'​onglet Jointure et joignez la table BBS_Routes et BBS_Ovenbird_QC en utilisant la colonne partagée '​Route'​ (choisissez '​Route'​ pour les deux tables). Désélectionnez la case qui indique "​Mettre la couche en cache dans la mémoire virtuelle"​. Vous devriez maintenant voir le nombre d'​observation de parulines par année apparaître dans la table des attributs du fichier '​BBS_Routes_QC'​. ​+**Étape 7**: À partir du menu '​Propriétés de la couche'​ de '​BBS_Routes_QC',​ trouvez l'​onglet Jointure et joignez la table BBS_Routes et BBS_Ovenbird_QC en utilisant la colonne partagée '​Route'​ (choisissez '​Route'​ pour les deux tables). Désélectionnez la case qui indique "​Mettre la couche en cache dans la mémoire virtuelle"​. Vous devriez maintenant voir le nombre d'​observation de parulines par année apparaître dans la table des attributs du fichier '​BBS_Routes_QC'​. ​
  
 **Étape 8**: Utilisez la '​Calculatrice de Champs'​ pour ajouter une nouvelle colonne (Nombre entier) à la table BBS_Routes qui contient la somme du nombre d'​oiseaux aperçus entre 1980-1994 et 1995-2010. Chaque champ doit être sélectionné un à un depuis "​Champs et valeurs"​ (désolé!). Le type de colonne doit être '​Nombre entier'​ (longueur 10 ou plus). Les noms des colonnes doivent contenir moins de 10 caractères,​ doivent commencer par une lettre et ne doivent pas contenir d'​espaces ou de caractères spéciaux. ​ **Étape 8**: Utilisez la '​Calculatrice de Champs'​ pour ajouter une nouvelle colonne (Nombre entier) à la table BBS_Routes qui contient la somme du nombre d'​oiseaux aperçus entre 1980-1994 et 1995-2010. Chaque champ doit être sélectionné un à un depuis "​Champs et valeurs"​ (désolé!). Le type de colonne doit être '​Nombre entier'​ (longueur 10 ou plus). Les noms des colonnes doivent contenir moins de 10 caractères,​ doivent commencer par une lettre et ne doivent pas contenir d'​espaces ou de caractères spéciaux. ​
  
-**Étape 9**: Sauvegardez les changements et sortez du mode d'​édition. Les observations par année peuvent être enlevées en désactivant la jointure. Notez que seuls les nouveaux champs créés font vraiment partie ​du fichier ​BBS_Routes_QC. ​+**Étape 9**: Sauvegardez les changements et sortez du mode d'​édition. Les observations par année peuvent être enlevées en désactivant la jointure. Notez que seuls les nouveaux champs créés font vraiment partie ​de la couche ​BBS_Routes_QC. ​
  
 **Étape 10**: Vous allez maintenant créer une surface interpolée (continue) afin de visualiser la répartition des parulines couronnées pour chaque période. Assurez-vous que le canevas utilise le même SCR que celui du fichier BBS_Routes_QC. Trouvez Interpolation IDW dans la boîte à outil de géotraitement. Spécifier BBS_Routes_QC avec la colonne correspondant à la période 1980-1994 comme Attribut d'​interpolation et ajoutez la couche en entrée. Vous devez définir le nombre de colonnes, le nombre de lignes et l'​emprise,​ de façon à obtenir un raster en sortie avec un résolution d'​exactement 2km x 2km, de façon à inclure la plupart des Routes (vous pouvez exclure les points au Nord), mais en restant dans les limites de la couverture des points. ​ Utilisez une calculatrice pour vérifier vos calculs. Cette image peut également vous aider: **Étape 10**: Vous allez maintenant créer une surface interpolée (continue) afin de visualiser la répartition des parulines couronnées pour chaque période. Assurez-vous que le canevas utilise le même SCR que celui du fichier BBS_Routes_QC. Trouvez Interpolation IDW dans la boîte à outil de géotraitement. Spécifier BBS_Routes_QC avec la colonne correspondant à la période 1980-1994 comme Attribut d'​interpolation et ajoutez la couche en entrée. Vous devez définir le nombre de colonnes, le nombre de lignes et l'​emprise,​ de façon à obtenir un raster en sortie avec un résolution d'​exactement 2km x 2km, de façon à inclure la plupart des Routes (vous pouvez exclure les points au Nord), mais en restant dans les limites de la couverture des points. ​ Utilisez une calculatrice pour vérifier vos calculs. Cette image peut également vous aider:
Line 187: Line 212:
 {{:​raster_res.png?​200x140|}} {{:​raster_res.png?​200x140|}}
  
-**Étape 11**: Utilisez la fonction '​Découper selon une couche de masque'​ sous Raster>​Extraction pour enlever de la carte interpolée les sections qui se trouvent à l'​extérieur du Québec. Utilisez ​le fichier ​province.shp en Quebec Lambert comme couche de masquage. ​+**Étape 11**: Utilisez la fonction '​Découper selon une couche de masque'​ sous Raster>​Extraction pour enlever de la carte interpolée les sections qui se trouvent à l'​extérieur du Québec. Utilisez ​la couche ​province en Quebec Lambert comme couche de masquage. ​
  
 Répétez les étapes 10-11 pour la période 1995-2010. Répétez les étapes 10-11 pour la période 1995-2010.
  
-**Étape 12**: Cliquez sur la carte avec l'​outil '​Identifier les entités'​ pour étudier les valeurs de la carte et prenez note des valeurs extrêmes. Dans la fenêtre des propriétés de la première couche interpolée,​ développez une palette en associant des couleurs aux valeurs de la couche (Style>Type de rendu>​pseudo-couleurs à bande unique). Choisissez une interpolation linéaire, cliquez sur le symbole '​+'​ et définissez trois nombres couvrant l'​étendu des valeurs sur la carte. Choisissez trois couleurs correspondantes qui offrent un bon contraste. Sauvegardez ce style et donnez-lui un nom. Pour le raster couvrant la période de 1995-2010, chargez le style que vous venez de sauvegarder. La même palette de couleurs devrait apparaître pour les deux périodes. ​+**Étape 12**: Cliquez sur la carte avec l'​outil '​Identifier les entités'​ pour étudier les valeurs de la carte et prenez note des valeurs extrêmes. Dans la fenêtre des propriétés de la première couche interpolée,​ développez une palette en associant des couleurs aux valeurs de la couche (Symbologie>Type de rendu>​pseudo-couleurs à bande unique). Choisissez une interpolation linéaire, cliquez sur le symbole '​+'​ et définissez trois nombres couvrant l'​étendu des valeurs sur la carte. Choisissez trois couleurs correspondantes qui offrent un bon contraste. Sauvegardez ce style en cliquant sur Style>​Sauvegarder le style, ​et donnez-lui un nom. Pour le raster couvrant la période de 1995-2010, chargez le style que vous venez de sauvegarder. La même palette de couleurs devrait apparaître pour les deux périodes. ​
  
 Notez-vous une différence dans l'aire de répartition des parulines couronnées entre les deux périodes? Notez-vous une différence dans l'aire de répartition des parulines couronnées entre les deux périodes?
  
-**DÉFI 1**: Ajoutez une colonne au fichier BBS Routes contenant la température annuelle moyenne (de la couche Quebec_mat_tenths.txt en degrés C x 10) dans un rayon de 20 km autour de chaque route. Vous aurez besoin ​d'outil de Statistiques de zones. ​+**DÉFI 1**: Ajoutez une colonne au fichier BBS Routes contenant la température annuelle moyenne (de la couche ​raster ​Quebec_mat_tenths.txt en degrés C x 10) dans un rayon de 20 km autour de chaque route. Vous aurez besoin ​de l'outil de Statistiques de zones. ​
  
-**DÉFI 2**: Recréez vos raster interpolés avec la distribution de parulines couronnées en utilisant l'​outil "​Multi-level B-Spline Interpolation"​ dans Saga. Utilisez ​le "​Processing toolbox"​ (Traitement). Comparez ces cartes avec celles que vous avez obtenues avec la pondération en fonction de l'​inverse de la distance (Inverse Distance Weighting).+**DÉFI 2**: Recréez vos raster interpolés avec la distribution de parulines couronnées en utilisant l'​outil "​Multi-level B-Spline Interpolation"​ dans Saga. Utilisez ​la boite à outils de géotraitement. Comparez ces cartes avec celles que vous avez obtenues avec la pondération en fonction de l'​inverse de la distance (Inverse Distance Weighting).
  
 ====== Exercice 4 - Téléchargement de fichiers, récapitulation et défi! ====== ====== Exercice 4 - Téléchargement de fichiers, récapitulation et défi! ======
Line 203: Line 228:
 Pour cet exercice, vous devez extraire l'​élévation moyenne ainsi que la couverture du sol à des sites d’occurrence de [[http://​www.gbif.org/​|Global Biodiversity Information Facility (GBIF)]] se trouvant dans le territoire autochtone de Wemindji, dans la région de la Baie-James au Québec. Vous allez vouloir travailler avec le **système de référence UTM 17N NAD83**. Pour se faire, vous devez suivre les étapes suivantes: ​ Pour cet exercice, vous devez extraire l'​élévation moyenne ainsi que la couverture du sol à des sites d’occurrence de [[http://​www.gbif.org/​|Global Biodiversity Information Facility (GBIF)]] se trouvant dans le territoire autochtone de Wemindji, dans la région de la Baie-James au Québec. Vous allez vouloir travailler avec le **système de référence UTM 17N NAD83**. Pour se faire, vous devez suivre les étapes suivantes: ​
  
-  * Sur le portail [[https://​open.canada.ca/​fr|Gouvernement ouvert Canada]], les données raster [[http://​ftp.geogratis.gc.ca/​pub/​nrcan_rncan/​elevation/​cdem_mnec/​|d'​élévation]] ​ (Modèle numérique d'​élévation du Canada) pour les zones 33D et 33E à l'​échelle 1:250,000.+  * Sur le portail [[https://​open.canada.ca/​fr|Gouvernement ouvert Canada]], ​téléchargez ​les données raster [[http://​ftp.geogratis.gc.ca/​pub/​nrcan_rncan/​elevation/​cdem_mnec/​|d'​élévation]] ​ (Modèle numérique d'​élévation du Canada) pour les zones 33D et 33E à l'​échelle 1:250,000.
   * Sur ce même site, téléchargez les Limites législatives des terres autochtones du Canada en format Shapefile (al_ta_ca_shp_gdm_fra.zip). Notez que la région de Wemindji se trouve à la frontière entre les zones 33D et 33E, et le fichier qui vous intéresse a la mention MODIFIE dans le nom.    * Sur ce même site, téléchargez les Limites législatives des terres autochtones du Canada en format Shapefile (al_ta_ca_shp_gdm_fra.zip). Notez que la région de Wemindji se trouve à la frontière entre les zones 33D et 33E, et le fichier qui vous intéresse a la mention MODIFIE dans le nom. 
   * Téléchargez ce {{::​qc_land_use_33de.zip|fichier ZIP}} qui contient un fichier raster GeoTiff de couverture du sol pour la région, ainsi qu'un fichier .qml de style/​couleurs associé. Ouvrez le fichier tif dans QGIS et chargez la palette de couleurs dans l'​onglet de Style des propriétés de la couche. ​   * Téléchargez ce {{::​qc_land_use_33de.zip|fichier ZIP}} qui contient un fichier raster GeoTiff de couverture du sol pour la région, ainsi qu'un fichier .qml de style/​couleurs associé. Ouvrez le fichier tif dans QGIS et chargez la palette de couleurs dans l'​onglet de Style des propriétés de la couche. ​
   * Télécharger {{::​occurrence.txt|ce fichier}} contenant la liste des occurrences présentes dans [[http://​www.gbif.org/​|GBIF]] dans la région de Wemindji avec des localisations en format latitude, longitude WGS84. Notez que ce fichier est délimité par des "​TAB"​. ​   * Télécharger {{::​occurrence.txt|ce fichier}} contenant la liste des occurrences présentes dans [[http://​www.gbif.org/​|GBIF]] dans la région de Wemindji avec des localisations en format latitude, longitude WGS84. Notez que ce fichier est délimité par des "​TAB"​. ​
-  * Vous pouvez ensuite fusionner les fichiers d'​élévation raster (.dem) en un seul fichier raster. ​+  * Vous pouvez ensuite fusionner les fichiers d'​élévation raster (.tifdans zones D et E en un seul fichier raster. ​
   * En utilisant les coordonnées latitude, longitude, ajouter les occurrences de GBIF au canevas. ​   * En utilisant les coordonnées latitude, longitude, ajouter les occurrences de GBIF au canevas. ​
   * "​Découpez"​ les occurrences pour n'​obtenir que celles présentes dans le territoire de Wemindji. ​   * "​Découpez"​ les occurrences pour n'​obtenir que celles présentes dans le territoire de Wemindji. ​
-  * Utiliser l'​extension 'point sampling tool' pour extraire le nom de chaque espèce en latin, la couverture du sol du fichier raster et l'​altitude à chaque endroit où il y a une occurrence GBIF. Notez que certains points contiennent un grand nombre d'​occurrences. ​+  * Utiliser l'​extension 'point sampling tool' pour extraire le nom de chaque espèce en latin, la couverture du sol du fichier raster et l'​altitude à chaque endroit où il y a une occurrence GBIF. Pour ce faire, cliquez sur le nom des colonnes désirés dans la longue liste de couches et de colonnes. Notez que certains points contiennent un grand nombre d'​occurrences. ​
  
-**DÉFI**: Utilisez l'​outil "​Kernel Density Estimation"​ dans SAGA (Processing toolbox) pour obtenir un raster affichant la densité des occurrences d'​oiseaux dans la région de Wemindji.+=-=-=-=-=-=-=-=-=-=- 
 + 
 +**DÉFI ​1**: Utilisez l'​outil "​Kernel Density Estimation"​ dans SAGA (Outils de géotraitement) pour obtenir un raster affichant la densité des occurrences d'​oiseaux dans la région de Wemindji. 
 + 
 +**DÉFI 2**: Utilisez l'​extension Data plotly pour générer un diagramme à barre affichant l'​altitude en fonction du règne (Kingdom) des occurrences observées à Wemindji. ​
 ====== Exercice 5 - Manipulation de rasters et introduction à l'​imagerie satellitaire ====== ====== Exercice 5 - Manipulation de rasters et introduction à l'​imagerie satellitaire ======
  
Line 245: Line 274:
 ++++ ++++
  
-**DÉFI**: Aller sur le site web du USGS EarthExplorer http://​earthexplorer.usgs.gov/​ et créez un compte en cliquant sur '​register'​. Ensuite, choisissez une région n'​importe où dans le monde de la taille de l'île de Montréal, où vous croyez qu'il y a eu beaucoup de changements au cours des 30 dernières années. Recherchez une image Collection 1 - Landsat 4-5TM image datant de quelque part entre 1980-1984 pour cette région et prise au mois de juillet ou août, qui contient moins de 10% de nuages. Téléchargez ce fichier en format "Level 1 product"​. Notez le path/row de cette image et téléchargez une image équivalente pour 2010-2014 (même mois, row et path). Ouvrez-les dans QGIS et calculez le NDVI pour chaque période. Soustrayez le NDVI de la carte ultérieure moins la première carte NDVI. Visualisez la carte résultante avec cette {{:​ndvidiff.qml|palette}} pour voir les zones où les changements dans la végétation on été les plus importants.  +=-=-=-=-=-=-=-=-=-=- 
-** + 
-DÉFI 2**: Utilisez la fonction "​Cluster Analysis for Grids" dans Outil de géotraitement (Processing Toolbox)>​SAGA pour effectuer une classification non-supervisée de vos images en utilisant les bandes 3,4,5 et 7.+**DÉFI ​1**: Aller sur le site web du USGS EarthExplorer http://​earthexplorer.usgs.gov/​ et créez un compte en cliquant sur '​register'​. Ensuite, choisissez une région n'​importe où dans le monde de la taille de l'île de Montréal, où vous croyez qu'il y a eu beaucoup de changements au cours des 30 dernières années. Recherchez une image Collection 1 - Landsat 4-5TM image datant de quelque part entre 1980-1984 pour cette région et prise au mois de juillet ou août, qui contient moins de 10% de nuages. Téléchargez ce fichier en format "Level 1 product"​. Notez le path/row de cette image et téléchargez une image équivalente pour 2010-2014 (même mois, row et path). Ouvrez-les dans QGIS et calculez le NDVI pour chaque période. Soustrayez le NDVI de la carte ultérieure moins la première carte NDVI. Visualisez la carte résultante avec cette {{:​ndvidiff.qml|palette}} pour voir les zones où les changements dans la végétation on été les plus importants. ​ 
 + 
 +**DÉFI 2**: Utilisez la fonction "​Cluster Analysis for Grids" dans Outil de géotraitement (Processing Toolbox)>​SAGA pour effectuer une classification non-supervisée de vos images en utilisant les bandes 3,4,5 et 7.
 ====== Exercice 6 - Utiliser GRASS dans QGIS ====== ====== Exercice 6 - Utiliser GRASS dans QGIS ======
 +
 +Sous Windows, vous devez démarrer QGIS avec l'​option "QGIS with GRASS 7 support"​ dans le menu Démarrer. ​
  
 **Étape 1**: Assurez-vous que l'​extension GRASS est activée et que les icônes GRASS sont bien visibles dans votre menu. Cliquez sur l'​icône '​Nouveau jeu de données'​ et créez un nouvelle géodatabase GRASS, un nouveau secteur, choisissez le SCR NAD83/UTM 18N, spécifiez la région GRASS par défaut en utilisant l'​étendue QGIS courante et choisissez un nom pour le jeu de données. ​ **Étape 1**: Assurez-vous que l'​extension GRASS est activée et que les icônes GRASS sont bien visibles dans votre menu. Cliquez sur l'​icône '​Nouveau jeu de données'​ et créez un nouvelle géodatabase GRASS, un nouveau secteur, choisissez le SCR NAD83/UTM 18N, spécifiez la région GRASS par défaut en utilisant l'​étendue QGIS courante et choisissez un nom pour le jeu de données. ​
  
-**Étape 2**: Sauvez le fichier routes.shp comme une nouvelle couche avec le SCR NAD83 / UTM 18N. Faites de même avec le shapefile Region_Hydrique. ​+**Étape 2**: Sauvez le fichier routes.shp comme une nouvelle couche ​Geopackage ​avec le SCR NAD83 / UTM 18N. Faites de même avec le shapefile Region_Hydrique. ​
  
-**Étape 3**: Vous devez maintenant importer les fichiers dans GRASS. ​Cliquez ​sur '​Ouvrir les outils GRASS', ​cliques ​sur 'arborescence des modules'>'​Gestion de fichiers'>'​Importer dans GRASS'>'​Importer un vecteur dans GRASS'>'​v.in.ogr.qgis'​. Sélectionnez le fichier avec les routes créé à l'​étape 2 et spécifiez un nom pour le fichier en sortie. Répétez cette étape pour le fichier Région_Hydrique. Ajouter les fichiers nouvellement crées au canevas en cliquant sur l'​icône '​Ajouter une couche vectorielle GRASS'​.+**Étape 3**: Vous devez maintenant importer les fichiers dans GRASS. ​Dans le menu Extensions>​GRASS,​ cliquez ​sur '​Ouvrir les outils GRASS', ​cliquez ​sur 'Modules'>'​Gestion de fichier'>'​Importer dans GRASS'>'​Importer un vecteur dans GRASS'>'​v.in.ogr.qgis'​. Sélectionnez le fichier avec les routes créé à l'​étape 2 et spécifiez un nom pour le fichier en sortie. Répétez cette étape pour le fichier Région_Hydrique. Ajouter les fichiers nouvellement crées au canevas en cliquant sur l'​icône '​Ajouter une couche vectorielle GRASS'​.
  
 **Étape 4**: Trouvez la fonction nommée r.in.gdal et importez le fichier 31h05dem.tif dans GRASS. Ajoutez cette nouvelle couche au canevas QGIS.  **Étape 4**: Trouvez la fonction nommée r.in.gdal et importez le fichier 31h05dem.tif dans GRASS. Ajoutez cette nouvelle couche au canevas QGIS. 
Line 260: Line 293:
 **Étape 5**: Trouvez la fonction g.region.multiple.raster et utiliser le fichier importé à l'​étape 4 pour définir la région courante (entrez le nom exact de la couche dans la boîte). ​ **Étape 5**: Trouvez la fonction g.region.multiple.raster et utiliser le fichier importé à l'​étape 4 pour définir la région courante (entrez le nom exact de la couche dans la boîte). ​
  
-**Étape 6**: Convertissez la couche des routes en format raster en utilisant la fonction v.to.rast.constant. Faites la même chose avec les régions hydriques. ​+**Étape 6**: Convertissez la couche des routes en format raster en utilisant la fonction v.to.rast.constant. Acceptez la valeur de 1 par défaut pour le Raster. Faites la même chose avec les régions hydriques. ​
  
-**Étape 7**: Utilisez la fonction r.grow.distance pour créer une surface raster continue dans laquelle la valeur de la distance à la route la plus proche est assignée à chaque pixel. Utilisez le fichier raster des routes (Étape 6) comme fichier en entrée et spécifiez un nom pour la carte raster de distance en sortie (laissez la couche de valeurs vide). ​+**Étape 7**: Utilisez la fonction r.grow.distance pour créer une surface raster continue dans laquelle la valeur de la distance à la route la plus proche est assignée à chaque pixel. Utilisez le fichier raster des routes (Étape 6) comme fichier en entrée et spécifiez un nom pour la carte raster de distance en sortie (laissez la couche de valeurs vide). Je vous encourage à visualiser cette couche dans QGIS et comprendre ce qu'​elle représente
  
 **Étape 8**: Utilisez la fonction r.null.to pour remplacer la valeur NULL par 0 dans le fichier raster Région hydrique. ​ **Étape 8**: Utilisez la fonction r.null.to pour remplacer la valeur NULL par 0 dans le fichier raster Région hydrique. ​
  
-**Étape 9**: Utilisez r.**mapcalculator** **N'​EXISTE PLUS** ​pour enlever les zones couvertes d'eau du fichier de l'​étape 7. Spécifiez le fichier raster Region Hydrique comme couche A et le fichier créé à l'​étape 7 comme couche B. Cliquez sur "​Utiliser la région de cette couche"​ à côté de la couche A. Dans le boîte de formule, entrez: ​+**Étape 9**: Trouvez la console GRASS "GRASS Shell" dans la liste des fonctions et écrivez la commande suivante ​pour enlever les zones couvertes d'eau du fichier de l'​étape 7.
  
 <​file>​ <​file>​
-((A-1)*-1)*B+r.mapcalc expression="​distance_sans_eau=((region_hydrique-1)*-1)*distance_routes" ​
 </​file>​ </​file>​
  
-et donnez un nom au fichier ​en sortie+en spécifiant les bons noms pour les couches
  
 **Étape 10**: Pour exclure les zones situées à moins de 1km des routes, créez un fichier texte (dans Notepad our Text Edit) avec ce contenu: ​ **Étape 10**: Pour exclure les zones situées à moins de 1km des routes, créez un fichier texte (dans Notepad our Text Edit) avec ce contenu: ​
Line 281: Line 314:
 </​file>​ </​file>​
  
-**Étape 11**: Utilisez la fonction r.reclass pour re-classifier le raster créé à l'​étape 9, et en utilisant les fichier ​de l'​étape 10 pour définir les règles de re-classification. ​+**Étape 11**: Utilisez la fonction r.reclass pour re-classifier le raster créé à l'​étape 9, et en utilisant les fichiers ​de l'​étape 10 pour définir les règles de re-classification. ​
  
 **Étape 12**: Utilisez la fonction r.clump pour donner un identifiant unique à toutes les zones contiguës. ​ **Étape 12**: Utilisez la fonction r.clump pour donner un identifiant unique à toutes les zones contiguës. ​
Line 293: Line 326:
 remplacez ID avec l'​identifiant de la zone la plus grande, et Etape_12 avec le nom que vous avez donné au raster de l'​Étape 12. Cette commande va créer une fichier nommé Largest_patch dans lequel la zone avec la superficie la plus grande sera isolée avec une valeur de 1. Vous pouvez ensuite afficher ce raster et changer la palette pour voir cette zone.  remplacez ID avec l'​identifiant de la zone la plus grande, et Etape_12 avec le nom que vous avez donné au raster de l'​Étape 12. Cette commande va créer une fichier nommé Largest_patch dans lequel la zone avec la superficie la plus grande sera isolée avec une valeur de 1. Vous pouvez ensuite afficher ce raster et changer la palette pour voir cette zone. 
  
-**Étape 15**: Convertissez ensuite ce fichier en vecteur avec la fonction r.to.vect.area. ​+**Étape 15**: Convertissez ensuite ce fichier en vecteur avec la fonction r.to.vect.area ​et visualisez le dans QGIS
  
 ====== Terminologie GRASS====== ====== Terminologie GRASS======