7.3. Lesson: Analyse de terrain

Certains types de raster vous permettent de gagner plus de perspicacité dans le terrain que ce qu’ils représentent. Les Modèles Numériques d’Élévation (MNE) sont particulièrement utiles dans ce sens. Dans cette leçon, vous utiliserez les outils d’analyse de terrain pour en savoir plus sur la zone d’étude pour le développement résidentiel proposé plus tôt.

Objectif de cette leçon : Utiliser les outils d’analyse de terrain pour avoir plus d’information sur le terrain.

7.3.1. Follow Along: Calcul d’un Ombrage

Nous allons utiliser la même couche DEM que dans la leçon précédente. Si vous commencez ce chapitre à partir de zéro, utilisez le panneau navigateur et chargez le fichier raster/SRTM/srtm_41_19.tif.

La couche DEM vous montre l’élévation du terrain, mais elle peut parfois sembler un peu abstraite. Elle contient toutes les informations 3D sur le terrain dont vous avez besoin, mais elle ne ressemble pas à un objet 3D. Pour avoir une meilleure impression du terrain, il est possible de calculer un ombrage, qui est un raster qui cartographie le terrain en utilisant la lumière et les ombres pour créer une image d’aspect 3D.

Nous allons utiliser des algorithmes dans le menu Raster ▶ analyse de terrain.

1. Cliquez sur le menu ombrage

2. L’algorithme vous permet de spécifier la position de la source lumineuse : Azimuth a des valeurs de 0 (Nord) à 90 (Est), 180 (Sud) et 270 (Ouest), tandis que Angle verticale définit la hauteur de la source lumineuse (0 à 90 degrés). Nous utiliserons les valeurs par défaut :

   

3. Enregistrez le fichier dans un nouveau dossier raster_analysis dans le dossier exercise_data avec le nom ombrage

4. Pour finir, cliquez sur Exécuter

Vous aurez maintenant une nouvelle couche appelée hillshade qui ressemble à cela :

Le résultat semble sympathique et en 3D. Mais pouvons-nous l’améliorer ? L’ombrage semble un peut en plâtre. Pouvons-nous le combiner avec nos autres rasters plus colorés ? C’est effectivement possible en utilisant l’ombrage en tant que couverture.

7.3.2. Follow Along: Utilisation d’un ombrage comme une couverture

Un ombrage peut fournir des informations très utiles sur la lumière du soleil à un moment donné de la journée. Mais il peut aussi être utilisé à des fins esthétiques, pour donner à la carte une meilleure apparence. La clé pour cela est de configurer l’ombrage pour qu’il soit presque complètement transparent.

1. Changer la symbologie de la couche srtm_41_19 originale pour utiliser le schéma Pseudocolor comme dans l’exercice précédent

2. Cacher toutes les couches sauf les couches srtm_41_19 et ombrage

3. Cliquez et faites glisser le srtm_41_19 pour qu’il se trouve sous la couche ombrage dans le panneau couches

4. Réglez la couche ombrage pour qu’elle soit transparente en cliquant sur l’onglet Transparence dans les propriétés de la couche

5. Définissez l” Opacité globale à 50%.

   Vous obtiendrez le résultat suivant :

   

6. Désactivez et réactivez la couche ombrage dans le panneau couches pour voir la différence que cela fait.

En utilisant un ombrage de cette façon, il est possible d’améliorer la topographie du paysage. Si l’effet ne vous semble pas suffisamment fort, vous pouvez changer la transparence de la couche hillshade ; mais bien sûr que plus l’ombrage devient brillant, plus les couleurs derrière lui seront sombres. Vous devrez trouver un équilibre qui fonctionne pour vous.

N’oubliez pas de sauvegarder le projet lorsque vous avez terminé.

7.3.3. Follow Along: Calcul d’une Pente

Une autre chose utile à savoir à propos du terrain est comment est sa pente. Si, par exemple, vous voulez construire des maisons ici, alors vous avez besoin d’un terrain relativement plat.

Pour ce faire, vous devez utiliser l’algorithme pente du menu terrain ▶ analyse terrain raster.

1. Ouvrir l’algorithme

2. Choisissez srtm_41_19 comme Couche d’élévation.

3. Enregistrez la sortie sous la forme d’un fichier portant le nom « pente » dans le même dossier que le « ombrage ».

4. Cliquez sur executer.

Vous allez maintenant voir la pente du terrain, les pixels noirs étant des terrains plats et les pixels blancs, des terrains escarpés :

7.3.4. Try Yourself Calculer l’aspect

Aspect est la direction de la boussole à laquelle fait face la pente du terrain. Un aspect de 0 signifie que la pente est orientée vers le nord, 90 vers l’est, 180 vers le sud et 270 vers l’ouest.

Comme cette étude se déroule dans l’hémisphère sud, les propriétés devraient idéalement être construites sur une pente orientée vers le nord afin qu’elles puissent rester à la lumière du soleil.

Utilisez l’algorithme Aspect du menu traitement ▶ analyse terrain raster pour obtenir la couche.

Vérifiez vos résultats

7.3.5. Follow Along: Utilisation de la Calculatrice Raster

Repensez au problème des agents immobiliers, que nous avons abordé pour la dernière fois dans la leçon sur l’analyse vecteur. Imaginons que les acheteurs souhaitent maintenant acheter un immeuble et construire un chalet plus petit sur la propriété. Dans l’hémisphère sud, nous savons qu’un terrain idéal pour le développement doit comporter des zones orientées vers le nord, avec une pente inférieure à cinq degrés. Mais si la pente est inférieure à 2 degrés, alors l’aspect n’a pas d’importance.

Heureusement, vous avez déjà des rasters qui vous montrent la pente aussi bien que l’aspect, mais vous n’avez aucune façon de savoir où les deux conditions sont immédiatement satisfaites. Comment cette analyse peut-elle être faite ?

La réponse réside dans la Calculatrice Raster.

QGIS dispose de différents calculateurs rasters :

• Raster ▶ Calculatrice Raster

• En cours de traitement

  • Analyse raster ▶ Calculatrice raster

  • GDAL ▶ Raster miscellaneous ▶ calculatrice raster

  • SAGA ▶ Raster calculus ▶ calculatrice raster

Chaque outil conduit aux mêmes résultats, mais la syntaxe peut être légèrement différente et la disponibilité des opérateurs peut varier.

Nous utiliserons analyse rasteer ▶ calculatrice raster dans la boîte à outils traitement

1. Ouvrez l’outil en double-cliquant dessus.

   • La partie supérieure gauche de la boîte de dialogue répertorie toutes les couches raster chargées sous la forme name@N, où name est le nom de la couche et N est la bande.

   • Dans la partie supérieure droite, vous verrez un grand nombre d’opérateurs différents. Arrêtez-vous un instant pour vous rappeler qu’un raster est une image. Vous devriez la voir comme une matrice 2D remplie de chiffres.

2. Le nord est à 0 (zéro) degré, donc pour que le terrain soit orienté vers le nord, son aspect doit être supérieur à 270 degrés et inférieur à 90 degrés. La formule est donc la suivante: :

   aspect@1 <= 90 OR aspect@1 >= 270
   

3. Vous devez maintenant configurer les détails du raster, comme la taille de la cellule, l’étendue et le CRS. Cela peut être fait manuellement ou automatiquement en choisissant une couche de référence. Choisissez cette dernière option en cliquant sur le bouton … à côté du paramètre couche de référence.

4. Dans le dialogue, choisissez la couche aspect, car nous voulons obtenir une couche avec la même résolution.

5. Sauvegardez la couche sous aspect_north.

   Le dialogue devrait ressembler à :

   

6. Pour finir, cliquez sur Exécuter

Votre résultat doit être cela :

Les valeurs de sortie sont 0 ou 1. Qu’est-ce que cela signifie ? La formule que nous avons écrite contient l’opérateur conditionnel OR. Par conséquent, le résultat final sera False (0) et True (1).

7.3.6. Try Yourself Plus de critères

Maintenant que vous avez fait l’aspect, créez deux nouvelles couches à partir du DEM.

• La première identifie les zones où la pente est inférieure ou égale à 2 degrés

• La seconde est similaire, mais la pente doit être inférieure ou égale à 5 degrés.

• Enregistrez-les sous exercise_data/raster_analysis sous les noms slope_lte2.tif et slope_lte5.tif.

Vérifiez vos résultats

7.3.7. Follow Along: Combinaison des résultats de l’analyse raster

Vous avez maintenant généré trois couches raster à partir du DEM :

• aspect_north : terrain orienté au nord

• slope_lte2 : pente égale ou inférieure à 2 degrés

• slope_lte5 : pente égale ou inférieure à 5 degrés

Lorsque les conditions sont remplies, la valeur du pixel est 1. Ailleurs, elle est de 0. Par conséquent, si vous multipliez ces raster, les pixels qui ont une valeur de 1 pour toutes les rasters obtiendront une valeur de 1 (les autres obtiendront 0).

Les conditions à remplir sont les suivantes :

• à 5 degrés de pente ou moins, le terrain doit être orienté vers le nord

• à 2 degrés de pente ou en dessous, la direction à laquelle le terrain fait face n’a pas d’importance.

Par conséquent, vous devez trouver des zones où la pente est égale ou inférieure à cinq degrés AND le terrain est orienté vers le nord, OR la pente est égale ou inférieure à 2 degrés. Un tel terrain serait propice au développement.

Pour calculer les zones qui satisfont ces critères :

1. Ouvrir à nouveau la Calculatrice raster.

2. Utilisez cette expression dans Expression: :

   ( aspect_north@1 = 1 AND slope_lte5@1 = 1 ) OR slope_lte2@1 = 1
   

3. Définissez le paramètre couche reference à aspect_north (peu importe si vous en choisissez un autre - ils ont tous été calculés à partir de srtm_41_19)

4. Enregistrez le résultat sous exercise_data/raster_analysis/ sous all_conditions.tif

5. Cliquez sur Exécuter

Le résultat :

7.3.8. Follow Along: Simplification du Raster

Comme vous pouvez le voir sur l’image ci-dessus, l’analyse combinée nous a laissé de nombreuses et très petites zones où les conditions sont remplies. Mais ces zones ne sont pas vraiment utiles pour notre analyse, car elles sont trop petites pour qu’on puisse s’appuyer sur quoi que ce soit. Débarrassons-nous de toutes ces minuscules zones inutilisables.

1. Ouvrez l’outil tamis (GDAL ▶ Analyse raster dans la boite a outil traitement)

2. Mettez le fichier Input à all_conditions, et le fichier tamis à all_conditions_sieve.tif (sous exercise_data/raster_analysis/).

3. Réglez le Threshold à 8 (minimum huit pixels contigus), et vérifiez le Use 8-connectedness.

   

   Une fois le traitement terminé, la nouvelle couche sera chargée.

   

   Que se passe-t-il ? La réponse se trouve dans les métadonnées du nouveau fichier raster.

4. Visualisez les métadonnées sous l’onglet Information du dialogue proprietes couche. Regardez la valeur STATISTICS_MINIMUM :

   

   Ce raster, comme celui dont il est dérivé, ne doit comporter que les valeurs 1 et 0, mais il comporte également un très grand nombre de négatif. L’examen des données montre que ce nombre agit comme une valeur nulle. Puisque nous ne recherchons que les zones qui n’ont pas été filtrées, fixons ces valeurs nulles à zéro.

5. Ouvrez la calculatrice raster, et construire cette expression: :

   (all_conditions_sieve@1 <= 0) = 0
   

   Cela permettra de maintenir toutes les valeurs non négatives, et de mettre les nombres négatifs à zéro, laissant toutes les zones avec la valeur 1 intactes.

6. Enregistrez le résultat sous exercise_data/raster_analysis/ sous le nom de all_conditions_simple.tif.

Votre sortie ressemble à cela :

C’est ce qui était attendu : une version simplifiée des résultats précédents. Souvenez-vous que si les résultats que vous obtenez à partir d’un outil ne sont pas ceux que vous attendiez, visualiser les métadonnées (et les attributs vectoriels, si applicable) peut s’avérer essentiel pour résoudre le problème.

7.3.9. Follow Along: Reclassification raster

Nous avons utilisé la calculatrice raster pour effectuer des calculs sur des couches raster. Il existe un autre outil puissant que nous pouvons utiliser pour extraire des informations des couches existantes.

Retour à la couche « aspect ». Nous savons maintenant qu’elle a des valeurs numériques dans une fourchette de 0 à 360. Ce que nous voulons faire, c’est reclasser cette couche en d’autres valeurs discrètes (de 1 à 4), selon l’aspect :

• 1 = Nord (de 0 à 45 et de 315 à 360)

• 2 = Est (de 45 à 135)

• 3 = Sud (de 135 à 225)

• 4 = Ouest (de 225 à 315)

Cette opération peut être réalisée avec la calculatrice raster, mais la formule deviendrait très très grande.

L’outil alternatif est l’outil reclassifier par table dans le menu analyse raster dans la boite outils.

1. Lancez l’outil

2. Choisissez aspect comme Couche raster d'entrée.

3. Cliquez sur le … de Table de reclassement. Un dialogue semblable à un tableau apparaîtra, dans lequel vous pourrez choisir les valeurs minimales, maximales et les nouvelles valeurs pour chaque classe.

4. Cliquez sur le bouton Ajouter ligne et ajoutez 5 lignes. Remplissez chaque rangée comme suit et cliquez sur OK :

   

   La méthode utilisée par l’algorithme pour traiter les valeurs seuils de chaque classe est définie par la limite de gamme.

5. Enregistrez la couche sous le nom de file:reclassified.tif dans le dossier exercise_data/raster_analysis/

   

6. Cliquez sur executer.

Si vous comparez la couche native aspect avec celle reclassifiée, il n’y a pas de grandes différences. Mais en regardant la légende, vous pouvez voir que les valeurs vont de 1 à 4.

Donnons à cette couche un meilleur style.

1. Ouvrez le panneau style couches

2. Choisissez Paletted/Unique values, au lieu de Singleband gray.

3. Cliquez sur le bouton classer pour récupérer automatiquement les valeurs et leur attribuer des couleurs aléatoires :

   

La sortie doit ressembler à ceci (vous pouvez avoir différentes couleurs étant donné qu’elles ont été générées de manière aléatoire) :

Grâce à cette reclassification et au style palettisé appliqué à la couche, vous pouvez immédiatement différencier les zones d’aspect.

7.3.10. Follow Along: Interroger raster

Contrairement aux couches vecteur, les couches raster n’ont pas de table d’attributs. Chaque pixel contient une ou plusieurs valeurs numériques (rasters à bande unique ou multibandes).

Toutes les couches rasters que nous avons utilisées dans cet exercice ne consistent qu’en une seule bande. Selon la couche, les valeurs des pixels peuvent représenter des valeurs d’élévation, d’aspect ou de pente.

Comment pouvons-nous interroger la couche raster pour obtenir la valeur d’un pixel ? Nous pouvons utiliser le bouton ^{Identifier des entités} !

1. Sélectionnez l’outil dans la barre d’outils Attributs.

2. Cliquez sur un emplacement aléatoire de la couche srtm_41_19. identifier résultat apparaîtra avec la valeur de la bande à l’emplacement cliqué :

   

3. Vous pouvez changer la sortie du panneau Résultat de l’identification du mode arbre actuel à un mode table en sélectionnant Table dans le menu Vue en bas du panneau :

   

Cliquer sur chaque pixel pour obtenir la valeur du raster pourrait devenir ennuyeux après un certain temps. Nous pouvons utiliser le plugin outil valeur pour résoudre ce problème.

1. Allez à Extensions ▶ Gérer/Installer des extensions…

2. Dans l’onglet tous, tapez value t dans la zone de recherche

3. Sélectionnez l’extension Value Tool, appuyez sur Installer l’extension et ensuite Fermer la fenêtre.

   

   Le nouveau panneau Value Tool apparaîtra.

   Astuce

   Si vous fermez le panneau, vous pouvez le rouvrir en l’activant dans la section Vue ▶ Panneaux ▶ Value Tool ou en cliquant sur l’icône dans la barre d’outils.

4. Pour utiliser le plugin, il suffit de cocher la case Activer et de s’assurer que la couche srtm_41_19 est active (cochée) dans le panneau Couches.

5. Déplacez le curseur sur la carte pour voir la valeur des pixels.

   

6. Mais il y a plus. Le plugin Value Tool vous permet d’interroger toutes les couches raster actives dans le panneau couches. Activez à nouveau les couches aspect et pente et passez la souris sur la carte :

   

7.3.11. In Conclusion

Vous avez vu comment dériver toutes sortes de produits d’analyse à partir d’un DEM. Il s’agit notamment des calculs de la hauteur, de la pente et de l’aspect. Vous avez également vu comment utiliser la calculatrice matricielle pour analyser et combiner ces résultats. Enfin, vous avez appris comment reclasser une couche et comment interroger les résultats.

7.3.12. What’s Next?

Maintenant que vous avez deux analyse : l’analyse vectorielle qui vous montre les parcelles susceptibles de convenir, et l’analyse raster qui vont montre le terrain susceptibles de convenir. Comment peut-on les combiner pour arriver à un résultat final pour ce problème ? C’est le sujet de la prochaine leçon, qui commence dans le module suivant.