Réponse courte

Une carte de pente correcte dans QGIS dépend surtout de ce que vous faites avant d'appuyer sur exécuter : confirmer que le MNT est dans un SCR projeté en mètres, régler le facteur z pour que les unités verticales et horizontales concordent, et décider si vous voulez des degrés ou des pourcents. Lancez ensuite Traitement → Analyse de terrain raster → Pente (ou l'algorithme GDAL Pente), classez la sortie avec un rendu Pseudo-couleur à bande unique, et vérifiez quelques valeurs par rapport à un terrain que vous connaissez. La palette de couleurs représente les 5 derniers pour cent du travail.

La pente est une couche de mesure, pas une image. Si les unités sous-jacentes sont fausses, la plus jolie palette du monde rapporte tout de même des absurdités.

Étape 1 — Vérifier le MNT avant toute chose

Chargez le MNT et ouvrez Propriétés de la couche → Information. Vous recherchez trois choses :

  • SCR. Il devrait s'agir d'un système projeté en mètres — par exemple EPSG:32630 (zone UTM 30N) ou une grille nationale. S'il indique EPSG:4326 (WGS84) ou tout SCR géographique, les unités sont des degrés et la pente sortira erronée. Reprojetez d'abord (voir étape 2).
  • Taille de pixel. Devrait indiquer quelque chose comme 30, -30 (mètres). Si elle indique 0.000277..., -0.000277..., ce sont des degrés — un autre signe que vous êtes dans un SCR géographique.
  • Unités verticales. Confirmez que l'altitude est en mètres. Les données américaines sont souvent en pieds, ce qui change le facteur z.

Confirmez aussi que le fichier est bien une altitude, et non une image d'ombrage que quelqu'un aurait enregistrée en GeoTIFF. L'outil Identifier devrait renvoyer des valeurs d'altitude (par ex. 412,7), et non des triplets RVB (par ex. 188, 188, 188).

Étape 2 — Reprojeter si nécessaire

Si le MNT est géographique, reprojetez-le. Raster → Projections → Warp (Reprojeter), ou dans Traitement lancez GDAL → Warp :

  • SCR cible : votre zone UTM locale, par ex. EPSG:32630.
  • Méthode de rééchantillonnage : Bilinéaire pour l'altitude (jamais Plus proche voisin — elle crée des artefacts en marches d'escalier qui deviennent de fausses pointes de pente ; Cubique convient aussi).
  • Résolution de sortie : définissez une taille de pixel nette, par ex. 30 m.

L'équivalent en ligne de commande :

gdalwarp -t_srs EPSG:32630 -r bilinear -tr 30 30 dem.tif dem_utm.tif

Étape 3 — Lancer l'algorithme Pente

Deux bonnes options figurent dans la boîte à outils de Traitement :

Natif QGIS : Analyse de terrain raster → Pente. Il produit la pente en degrés, applique un facteur z et utilise la méthode de Horn sur une fenêtre 3×3. Réglez :

  • Facteur Z : 1 si les unités verticales et horizontales sont toutes deux en mètres ; 0.3048 si l'altitude est en pieds et la grille en mètres.

GDAL Pente (GDAL → Analyse raster → Pente) offre plus de contrôle :

  • « Pente exprimée en pourcent au lieu de degrés » — cochez ceci pour le pourcentage de dénivelé.
  • « Échelle (rapport des unités vert. aux horiz.) » — c'est le facteur z déguisé : 1 pour mètre/mètre, 0.3048 pour pied/mètre.
  • « Calculer les bords » — activez pour éviter une bordure NoData d'un pixel de large.

Équivalent en ligne de commande :

gdaldem slope dem_utm.tif slope_deg.tif        # degrés
gdaldem slope -p dem_utm.tif slope_pct.tif     # pourcent

Degrés vs pourcent — choisir délibérément

Ils mesurent la même chose sur des échelles différentes et sont faciles à confondre :

  • Les degrés vont de 0 à 90. Intuitifs pour les seuils d'aléa (par ex. les pentes sensibles aux glissements de terrain souvent signalées au-delà de 25–35 degrés) et pour s'adresser à des interlocuteurs non-SIG.
  • Le pourcentage de dénivelé est (dénivelé / distance) × 100. Une pente de 45 degrés est de 100 pour cent, pas de 100 degrés ; une pente de 100 pour cent est raide, pas verticale. Courant en ingénierie routière et de drainage.

Conversion : degrés = atan(pourcent / 100) × 180 / π. Quel que soit votre choix, inscrivez l'unité dans le titre de la légende — « Pente (degrés) » — afin que le prochain lecteur ne puisse pas la mal interpréter.

Étape 4 — Classer et styliser

Clic droit sur la couche de pente → Propriétés → Symbologie → Pseudo-couleur à bande unique. Plutôt qu'une rampe continue, utilisez des classes discrètes ou classifiées liées à la décision que la carte soutient. Un schéma courant de stabilité de terrain en degrés :

ClassePlage (deg)Signification
Plat0–5Constructible, faible érosion
Doux5–15Gérable
Modéré15–25Prudence
Raide25–35Contraint
Très raide>35Aléa élevé

Les classes discrètes communiquent les seuils bien mieux qu'une rampe lisse bleu-rouge, qui masque l'emplacement de la rupture significative.

Passe d'AQ pratique

Avant que la carte ne quitte votre machine :

  • Vérification de l'histogramme. Dans la Symbologie, chargez le min/max et affichez l'histogramme. Pour la plupart des paysages, l'essentiel se situe sous 40 degrés avec une longue queue. Un pic marqué à 89–90 degrés signale une erreur d'unités ou de SCR.
  • Contrôle ponctuel d'une pente connue. Utilisez l'outil Identifier sur une entité que vous pouvez vérifier — un front de carrière, un talus de route nivelé, une berge de lac plate (devrait indiquer près de 0).
  • Regardez les bords, pas seulement le centre. Les jointures de dalles, les bordures NoData et les zones comblées dans le SRTM/lidar produisent des artefacts en anneau ou en bandes qui n'apparaissent qu'aux marges.
  • Méfiez-vous du bruit sur le lidar haute résolution. Un MNT à 1 m amplifie le micro-relief (retours de végétation, ornières de véhicules). Si vous voulez une pente à l'échelle des formes du relief, rééchantillonnez vers une grille plus grossière ou appliquez un léger filtre passe-bas avant de calculer la pente, et indiquez-le dans les métadonnées.

Pièges courants et leurs causes

  • Pente bloquée près de 90 partout. Le MNT est dans un SCR géographique — les degrés servent de distance horizontale. Reprojetez vers l'UTM. C'est de loin l'échec le plus fréquent.
  • Pente trop faible d'un facteur constant. Les unités verticales sont en pieds mais le facteur z est resté à 1. Réglez le facteur z à 0,3048.
  • Pente mouchetée et bruitée. Lidar haute résolution avec de vraies micro-entités, ou rééchantillonnage au plus proche voisin lors de la reprojection. Utilisez le rééchantillonnage bilinéaire et envisagez un lissage.
  • La légende indique « pourcent » mais les valeurs ressemblent à des degrés (ou inversement). L'option -p de GDAL a été activée de manière incohérente entre les exécutions. Relancez et confirmez l'unité.
  • Traiter un ombrage comme une pente. L'ombrage est une visualisation en relief ombré éclairé sous un angle ; ce n'est pas une mesure et il ne peut pas être classé en seuils de pente.

Le point de vue Bathyl

Nous construisons les rasters de pente comme des couches de preuve, ce qui signifie que le SCR, le facteur z, la résolution source et le choix d'unité sont consignés à côté du fichier, et non gardés dans la tête de quelqu'un. Une carte de pente n'est fiable que lorsqu'un autre analyste peut voir exactement comment elle a été produite et reproduire les mêmes chiffres à partir du même MNT.

Lectures Bathyl complémentaires

Sources