Outils pour analyser et visualiser des MNT (à partir de GDAL 1.7.0).
Introduit dans la version 1.7.0.
Usage :
Pour générer une carte de reliefs ombragés à partir de n’importe quel raster d’élévation géré par GDAL :
gdaldem hillshade input_dem output_hillshade
[-z ZFactor (défaut=1)] [-s scale* (défaut=1)]"
[-az Azimuth (défaut=315)] [-alt Altitude (défaut=45)]
[-alg ZevenbergenThorne] [-combined]
[-compute_edges] [-b Band (défaut=1)] [-of format] [-co "NAME=VALUE"]* [-q]
Pour générer une carte des pentes à partir de n’importe quel raster d’élévation géré par GDAL :
gdaldem slope input_dem output_slope_map"
[-p use percent slope (défaut=degrees)] [-s scale* (défaut=1)]
[-alg ZevenbergenThorne]
[-compute_edges] [-b Band (défaut=1)] [-of format] [-co "NAME=VALUE"]* [-q]
Pour générer une carte d’aspect à partir de n’importe quel raster d’élévation géré par GDAL. Créé en sortie un raster en float de 32 bit avec des valeurs de pixels de 0 à 360 indiquant l’azymuth :
gdaldem aspect input_dem output_aspect_map"
[-trigonometric] [-zero_for_flat]
[-alg ZevenbergenThorne]
[-compute_edges] [-b Band (défaut=1)] [-of format] [-co "NAME=VALUE"]* [-q]
Pour générer une carte de relief en couleur à partir de n’importe quel raster d’élévation géré par GDAL où color_text_file contient les lignes au format “elevation_value rouge vert bleu” :
gdaldem color-relief input_dem color_text_file output_color_relief_map
[-alpha] [-exact_color_entry | -nearest_color_entry]
[-b Band (défaut=1)] [-of format] [-co "NAME=VALUE"]* [-q]
Pour générer une carte d’Index de Rugosité du Terrain (Terrain Ruggedness Index (TRI)) à partir de n’importe quel raster d’élévation géré par GDAL :
gdaldem TRI input_dem output_TRI_map
[-compute_edges] [-b Band (défaut=1)] [-of format] [-q]
Pour générer une carte d’Index de Position Topographique (Topographic Position Index (TPI)) à partir de n’importe quel raster d’élévation géré par GDAL :
gdaldem TPI input_dem output_TPI_map
[-compute_edges] [-b Band (défaut=1)] [-of format] [-q]
Pour générer une carte de rugosité à partir de n’importe quel raster d’élévation géré par GDAL :
gdaldem roughness input_dem output_roughness_map
[-compute_edges] [-b Band (défaut=1)] [-of format] [-q]
Note
Scale est le rapport des unités verticales sur celles horizontales pour Feet:Latlong utilisez scale=370400, pour Meters:LatLong utilisez scale=111120
Cette commande possède 7 modes différents :
Les options générales suivantes sont disponibles :
Pour tous les algorithme, sauf color-relief, une valeur nodata dans le jeu de données cible sera émise si au moins un pixel définie dans la valeur nodata est trouvé dans une fenêtre 3x3 centrée autour du pixel source. La conséquence de cela est qu’il y aura un bord d’un pixel autour de chaque image définie à la valeur nodata. À partir de GDAL 1.8.0, si l’option -compute_edges est définie, ‘’gdaldem’’ calculera les valeurs au bord de l’image ou si une valeur nodata est trouvée dans la fenêtre 3x3, en interpolant les valeurs manquantes.
Cette commande renvoie un raster de 8 bit avec un bel effet de relief ombragé. Cela est très utile pour visualiser le terrain. Vous pouvez en option définir l’azymuth et l’altitude la source de lumière, un facteur d’exagération vertical et un facteur d’échelle à tenir compte pour les différences entre les unités verticales et horizontales.
La valeur 0 est utilisée comme valeur nodata en sortie.
Les options spécifiques suivantes sont disponibles :
Cette commande utilise un raster MNT et renvoie un raster en float 32 bit avec des valeurs de pente. Vous avez la possibilité de définir le type de pente que vous voulez : degrés ou pourcentage. Dans le cas où les unités horizontales diffèrent des unités verticales vous pouvez également fournie un facteur d’échelle.
La valeur -9999 est utilisé comme valeur nodata en sortie.
Les options spécifiques suivantes sont disponibles :
Cette commande renvoie un raster float 32 bit avec des valeurs entre 0° et 360° représentant l’azymuth dont les pentes font face. La définition de l’azymuth est : 0° signifie que la pente est face au Nord, 90° face à l’Est, 180° face au sud et 270° face à l’Ouest (en supposant que le haut du raster en entrée est orienté au Nord). La valeur -9999 de l’aspect est utilisé comme valeur nodata pour indiquer un aspect indéfinie dans les zones plates avec la pente = 0.
Les options spécifiques suivantes sont disponibles :
En utilisant ces deux options, l’aspect renvoyé par ‘’gdaldem aspect’’ doit être identique à celui de la commande ‘’r.slope.aspect’’ de GRASS. Autrement il est identique à celui de la commande aspect.cpp de Matthew Perry.
Cette commande renvoie un raster à 3 bande (RVB) ou à 4 bandes (RVBA) avec des valeurs calculées à partir de l’élévation et d’un fichier de configuration de couleur au format texte, contenant l’association entre les différentes valeurs d’élévation et la couleur désirée correspondante. Par défaut, les couleurs entre les valeurs d’élévation données sont mélangées en douceur et le résultat est un beau MNT coloré. Les options -exact_color_entry ou -nearest_color_entry peuvent être utilisées pour éviter cette interpolation linéaire pour les valeurs qui n’ont pas de correspondance avec un index de couleur du fichier de configuration.
Les options spécifiques suivantes sont disponibles :
Le mode color-relief est le seul mode gérant le format VRT en sortie. Dans ce cas, il traduira le fichier de configuration des couleurs en éléments <LUT> appropriés. Notez que les élévations définie en pourcentage seront traduit en valeur absolue, ce qui doit être pris en compte lorsque les statistiques du raster source diffère de celui qui a été utilisé lors de la construction du VRT.
Le fichier de configuration des couleurs au format texte contient généralement 4 colonnes par ligne : la valeur de l’élévation et les composants correspondants de Rouge, Vert, Bleu (entre 0 et 255). La valeur d’élévation peut être une valeur en virgule flottante, ou le mot-clé nv pour la valeur nodata. L’élévation peut aussi être exprimée en pourcentage : 0 % étant la valeur minimale trouvé dans le raster, 100 % la valeur maximale.
Une colonne supplémentaire peut être ajouté optionnellement pour le composant alpha. S’il n’est pas définie, l’opacité complète (255) est supposée.
Différents séparateurs de champs sont acceptés : virgule, tabulation, espaces, ‘:’.
Les couleurs généralement utilisées par GRASS peuvent également être spécifiées en utilisant leur nom, au lieu du triplet RVB. La liste des noms gérés est : white, black, red, green, blue, yellow, magenta, cyan, aqua, grey/gray, orange, brown, purple/violet et indigo.
Depuis GDAL 1.8.0, les fichiers de palette .cpt GMT sont également géré (COLOR_MODEL = RGB suelement).
Note
La syntaxe du fichier de configuration de couleur est dérivé de celui géré par la commande r.colors de GRASS. Les fichiers (.clr) de table de couleur HDR d’ESRI correspondent également à cette syntaxe. Le composent alpha et la gestion des tabulations et virgules comme séparateurs sont des extensions spécifiques à GDAL.
Par exemple :
3500 white
2500 235:220:175
50% 190 185 135
700 240 250 150
0 50 180 50
nv 0 0 0 0
Cette commande renvoie un raster à une seule bande avec des valeurs calculées à partir de l’élévation. TRI signifie Terrain Ruggedness Index, qui est définie comme la différence moyenne entre un pixel central et ses cellules l’entourant (voir Wilson et al 2007, Marine Geodesy 30:3-35).
La valeur -9999 est utilisé comme valeur nodata en sortie.
Il n’y a pas d’options spécifiques.
Cette commande renvoie un raste à une seule bande avec des valeurs calculées à partir de l’élévation. TPI signifie Topographic Position Index, qui est définie comme la différence entre un pixel central et la moyenne des cellules l’entourant (voir Wilson et al 2007, Marine Geodesy 30:3-35).
La valeur -9999 est utilisé comme valeur nodata en sortie.
Il n’y a pas d’options spécifiques.
Cette commande renvoie un raster à une seule bande calculé à partir de l’élévation. La rugosité est la plus grande différence inter-cellule d’un pixel central et ses cellules l’entourant, comme définie dans Wilson et al (2007, Marine Geodesy 30:3-35).
La valeur -9999 est utilisé comme valeur nodata en sortie.
Il n’y a pas d’options spécifiques.
Matthew Perry <perrygeo@gmail.com>, Even Rouault <even.rouault@mines-paris.org>, Howard Butler <hobu.inc@gmail.com>, Chris Yesson <chris.yesson@ioz.ac.uk>
Derived from code by Michael Shapiro, Olga Waupotitsch, Marjorie Larson, Jim Westervelt : U.S. Army CERL, 1993. GRASS 4.1 Reference Manual. U.S. Army Corps of Engineers, Construction Engineering Research Laboratories, Champaign, Illinois, 1-425.
Documentation des commandes GRASS connexes :