Outdated version of the documentation. Find the latest one here.

8.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.

8.3.1. basic Follow Along: Calcul d’un Ombrage

Le DEM que vous avez en ce moment sur votre carte doit vous montrer l’élévation du terrain, mais il peut parfois sembler un peu abstrait. Il contient toute l’information 3D sur le terrain dont vous avez besoin, mais il ne ressemble pas à un objet 3D. Pour obtenir un meilleur aspect du terrain, il est possible de calculer un ombrage, qui est un raster qui cartographie le terrain en utilisant la lumière et l’ombre pour créer une image 3D.

Pour travailler avec des DEM, vous devez utiliser l’outil d’analyse QGIS tout-en-un MNE (Modèles de terrain).

  • Cliquez sur l’élément du menu Raster ‣ Analyse ‣ MNE (Modèles de terrain).

  • Dans la boîte de dialogue qui apparaît, assurez-vous que le Fichier source est la couche MNE.

  • Mettez comme Fichier de sortie hillshade.tif dans le répertoire exercise_data/residential_development.

  • Vérifiez également que l’option Mode a sélectionné Hillshade.

  • Cochez la case à côté de Charger dans la carte une fois terminé.

  • Vous pouvez laisser toutes les autres options inchangées.

  • Cliquez sur OK pour générer l’ombrage.

  • Lorsque l’on vous dit que le processus est achevé, cliquer sur OK dans le message pour le fermer.

  • Cliquez sur Fermer sur la boîte de dialogue principale MNE (Modèles de terrain).

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

../../../_images/hillshade_raster.png

That looks nice and 3D, but can we improve on this? On its own, the hillshade looks like a plaster cast. Can’t we use it together with our other, more colorful rasters somehow? Of course we can, by using the hillshade as an overlay.

8.3.2. basic 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.

  • Changez la symbologie du DEM d’origine pour utiliser le schéma Pseudo-couleur comme dans l’exercice précédent.

  • Cachez toutes les couches, exceptées les couches DEM et hillshade.

  • Cliquez et faites glisser le DEM pour être sous la couche hillshade dans la Légende de la carte.

  • Configurez la couche hillshade comme transparente en ouvrant ses Propriétés de la couche et rendez-vous à l’onglet Transparence.

  • Configurez la Transparence globale à 50% :

  • Cliquez sur OK dans la boîte de dialogue Propriétés de la couche. Vous recevrez un résultat comme ça :

    ../../../_images/hillshade_pseudocolor.png
  • Désactivez la couche hillshade puis réactivez-la dans la Légende de la carte pour voir la différence que cela fait.

Using a hillshade in this way, it’s possible to enhance the topography of the landscape. If the effect doesn’t seem strong enough to you, you can change the transparency of the hillshade layer; but of course, the brighter the hillshade becomes, the dimmer the colors behind it will be. You will need to find a balance that works for you.

N’oubliez pas de sauvegarder votre carte quand vous avez terminé.

Note

Pour les deux prochains exercices, utilisez une nouvelle carte. Chargez le jeu de données raster DEM dans la carte (exercise_data/raster/SRTM/srtm_41_19.tif). C’est pour simplifier les choses pendant que vous travaillez avec les outils d’analyse raster. Sauvegardez la carte comme exercise_data/raster_analysis.qgs.

8.3.3. moderate 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 faire cela, vous devez utiliser le mode Pente de l’outil MNE (Modèles de terrain) .

  • Ouvrez l’outil comme avant.

  • Sélectionnez l’option ModePente :

    ../../../_images/dem_slope_dialog.png
  • Définissez l’emplacement de la sauvegarde à to exercise_data/residential_development/slope.tif

  • Cochez la case à coche Charger dans la carte....

  • Cliquez sur OK et fermez les boîtes de dialogue lorsque le processus est terminé, et cliquez sur Fermer pour fermer la boîte de dialogue. Vous verrez un nouveau raster chargé dans votre carte.

  • Avec le nouveau raster sélectionné dans la Légende de la carte, cliquez sur le bouton Étendre l’histogramme à tout le jeu de données. Vous verrez maintenant la pente du terrain, les pixels noirs représentant un terrain plat et les pixels blancs un terrain pentu :

    ../../../_images/slope_raster.png

8.3.4. moderate Try Yourself Calcul de l’aspect

L’aspect du terrain fait référence à la direction à laquelle il fait face. Comme cette étude prend place dans l’hémisphère sud, les propriétés devraient idéalement être construites sur une pente orientée au Nord afin qu’elles puissent rester dans la lumière du soleil.

  • Utilisez le mode Aspect de l’outil MNE (Modèles de terrain) pour calculer l’aspect du terrain.

Vérifiez vos résultats

8.3.5. moderate Follow Along: Utilisation de la Calculatrice Raster

Think back to the estate agent problem, which we last addressed in the Vector Analysis lesson. Let’s imagine that the buyers now wish to purchase a building and build a smaller cottage on the property. In the Southern Hemisphere, we know that an ideal plot for development needs to have areas on it that are north-facing, and with a slope of less than five degrees. But if the slope is less than 2 degrees, then the aspect doesn’t matter.

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.

  • Cliquez sur Raster > Calculatrice raster... pour enclencher cet outil.

  • To make use of the aspect dataset, double-click on the item aspect@1 in the Raster bands list on the left. It will appear in the Raster calculator expression text field below.

Le Nord est à 0 (zéro) degrés, donc pour le terrain à face nord, son aspect doit être plus grand que 270 degrés et plus petit que 90 degrés.

  • Dans le champ Expression de la calculatrice raster, entrez cette expression :

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

  • Mettez le fichier de sortie comme aspect_north.tif dans le répertoire exercise_data/residential_development/.

  • Assurez-vous que la case Ajouter le résultat au projet est cochée.

  • Cliquez sur OK pour démarrer le processus.

../../../_images/raster_calculator.png

Votre résultat doit être cela :

../../../_images/aspect_result.png

8.3.6. moderate Try Yourself

Maintenant que vous avez fait l’aspect, créez deux nouvelles analyses séparées pour la couche MNE.

  • La première identifiera toutes les zones où la pente est plus petite ou égale à 2 degrés.

  • La seconde est similaire, mais la pente devra être plus petite ou égale à 5 degrés.

  • Sauvegardez-les sous exercise_data/residential_development/ comme slope_lte2.tif et slope_lte5.tif.

Vérifiez vos résultats

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

Vous avez maintenant trois nouvelles analyses raster sur la couche MNE :

  • guilabel:aspect_north: le terrain face au nord

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

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

Là où les conditions de ces couches se rencontrent, elles sont égales à 1. Ailleurs, elles sont égales à 0. Cependant, si vous multipliez un de ces rasters par un autre, vous obtiendrez les zones où les deux sont égales à 1.

Les conditions à remplir sont : égal ou inférieur à 5 degrés de pente, le terrain doit être orienté au Nord ; mais égal ou inférieur 2 degrés de pente, la direction vers laquelle le terrain est orienté n’a pas d’importance.

Cependant, vous devez trouver des zones où la pente est égale ou inférieure à 5 degrés ET`le terrain orienté au Nord, :kbd:`OU la pente est égale ou inférieure à 2 degrés. Un tel terrain serait approprié pour le développement.

Pour calculer les zones qui satisfont ces critères :

  • Ouvrez à nouveau votre Calculatrice raster

  • Utilisez la liste Bandes de raster, les boutons Opérateurs, et votre clavier pour construire cette expression dans la zone de texte Expression de la calculatrice raster :

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

  • Sauvegardez la sortie sous exercise_data/residential_development/ comme all_conditions.tif.

  • Cliquez sur OK dans la Calculatrice raster. Vos résultats :

    ../../../_images/development_analysis_results.png

8.3.8. moderate Follow Along: Simplification du Raster

Comme vous pouvez le voir dans l’image du bas, l’analyse combinée nous a laissé avec beaucoup de très petites zones où les conditions sont remplies. Mais elles ne sont pas très utiles pour nos analyses, puisqu’elles sont trop petites pour y construire quelque chose. Éliminons toutes ces toutes petites zones inutiles.

  • Ouvrez l’outil Tamiser (Raster ‣ Analyse ‣ Tamiser).

  • Mettez comme Fichier source all_conditions, et comme Fichier de sortie all_conditions_sieve.tif (sous exercise_data/residential_development/).

  • Mettez pour les deux valeurs Seuil and Connexions de pixel à 8, puis exécuter l’outil.

../../../_images/raster_seive_dialog.png

Une fois le processus exécuté, la nouvelle couche sera chargée dans le canevas. Mais quand vous essayez d’utiliser l’outil histogramme complet pour visualiser les données, voici ce qui se passe :

../../../_images/seive_result_incorrect.png

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

  • Regardez les métadonnées sous l’onglet Métadonnées de la boîte de dialogue Propriétés de la couche. Regardez dans la section en bas Propriétés.

../../../_images/seive_metadata.png

Alors que ce raster, comme celui duquel il est dérivé, devrait seulement présenter les valeurs 1 et 0, il possède comme valeur STATISTICS_MINIMUM un très grand nombre négatif. Une enquête sur les données montre que ce nombre agit comme une valeur nulle. Comme nous recherchons seulement les zones qui n’ont pas été éliminées, mettons ces valeurs nulles à zéro.

  • Ouvrez à nouveau la Calculatrice raster, et construisez cette expression :

    (all_conditions_sieve@1 <= 0) = 0

    Cela maintiendra toutes les valeurs zéro existantes, tout en définissant les nombres négatifs à zéro ; ce qui éliminera toutes les zones avec une valeur 1 intacte.

  • Sauvegardez la sortie sous exercise_data/residential_development/ comme all_conditions_simple.tif.

Votre sortie ressemble à cela :

../../../_images/raster_seive_correct.png

This is what was expected: a simplified version of the earlier results. Remember that if the results you get from a tool aren’t what you expected, viewing the metadata (and vector attributes, if applicable) can prove essential to solving the problem.

8.3.9. In Conclusion

Vous avez vu comment tirer toutes les sortes de produits d’analyse d’un MNE. Cela inclut l’ombrage, la pente et les calculs d’aspect. Vous avez également vu comment utiliser la calculatrice raster pour analyser plus loin et combiner les résultats.

8.3.10. 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.