16.2. Analyse Raster

16.2.1. Calculatrice Raster

La Calculatrice Raster du menu Raster vous permet d’effectuer des calculs sur la base des valeurs des pixels d’un raster existant (voir Fig. 16.20). Le résultat est écrit dans un nouveau raster dans un des formats gérés par GDAL.

../../../_images/raster_calculator1.png

Fig. 16.20 Calculatrice Raster

The Raster bands list contains all loaded raster layers that can be used. To add a raster to the raster calculator expression field, double click its name in the Fields list. You can then use the operators to construct calculation expressions, or you can just type them into the box.

In the Result layer section, you will need to define an output layer. You can:

  • checkbox Create on-the-fly raster instead of writing layer to disk:

    • If unchecked, the output is stored on the disk as a new plain file. An Output layer path and an Output format are required.

    • If checked, a virtual raster layer, i.e. a raster layer defined by its URI and whose pixels are calculated on-the-fly, is created. It’s not a new file on disk; the virtual layer is still connected to the rasters used in the calculation meaning that deleting or moving these rasters would break it. A Layer name can be provided, otherwise the calculation expression is used as such. Removing the virtual layer from the project deletes it, and it can be made persistent in file using the layer Export ► Save as… contextual menu.

  • Define the Spatial extent of the calculation based on an input raster layer extent, or on custom X,Y coordinates

  • Set the Resolution of the layer using columns and rows number. If the input layer has a different resolution, the values will be resampled with the nearest neighbor algorithm.

  • With the checkbox Add result to project checkbox, the result layer will automatically be added to the legend area and can be visualized. Checked by default for virtual rasters.

The Operators section contains all available operators. To add an operator to the raster calculator expression box, click the appropriate button. Mathematical calculations (+, -, *, … ) and trigonometric functions (sin, cos, tan, … ) are available. Conditional expressions (=, !=, <, >=, … ) return either 0 for false or 1 for true, and therefore can be used with other operators and functions.

Indication

Voir aussi l’algorithme Calculatrice raster.

16.2.1.1. Exemples

Conversion en pieds de données d’élévation en mètres

Pour créer un raster d’élévation en pieds à partir d’un raster en mètres, vous devez utiliser le facteur de conversion 3,28. La formule est la suivante :

"elevation@1" * 3.28

Utilisation d’un masque

Si vous souhaitez masquer certaines zones d’un raster, par exemple si vous vous intéressez uniquement aux altitudes supérieures à 0, vous pouvez utiliser la formule suivante qui crée un masque et l’applique au raster en une seule étape.

("elevation@1" >= 0) * "elevation@1"

En d’autres termes, pour tous les pixels supérieurs ou égal à 0, l’expression conditionnelle renvoie 1, ce qui conserve la valeur originelle en la multipliant par 1. Sinon, l’expression conditionnelle renvoie 0, ce qui remplace la valeur originelle par 0. Ceci permet de créer le masque à la volée.

Classify a Raster

Si vous souhaitez classer un raster, par exemple en deux classes d’altitude, vous pouvez utiliser la formule suivante pour créer un raster contenant deux valeurs, 1 et 2, en une seule étape :

("elevation@1" < 50) * 1 + ("elevation@1" >= 50) * 2

En d’autres termes, les pixels de valeur inférieure à 50, ils prennent la valeur 1. Pour tous les pixels supérieurs ou égal à 50, ils prennent la valeur 2.

Or you can use the IF operator.

if ( elevation@1 < 50 , 1 , 2 )

16.2.2. Alignement de rasters

Cet outil permet de prendre plusieurs rasters en entrée et de les aligner exactement :

  • reprojète dans le même SCR

  • échantillonne à la même taille de cellule et décale sur la même grille,

  • découpe une zone d’intérêt,

  • rééchantillonne les valeurs lorsque cela est nécessaire.

Les rasters seront sauvegardés dans de nouveaux fichiers.

Tout d’abord, ouvrez les outils à partir de Raster ► Aligner les Rasters… et cliquez sur symbologyAdd Ajouter un nouveau raster pour choisir un raster existant. Sélectionnez un fichier de sortie pour enregistrer le raster après l’alignement, la méthode de ré-échantillonnage et si nécessaire Rééchelonner les valeurs en fonction de la taille de la cellule. La méthode de ré-échantillonnage est à choisir parmi (voir Fig. 16.21) :

  • Au plus proche voisin

  • Bilinéaire (noyau de 2x2)

  • Cubique (noyau de 4x4) : approximation par convolution cubique

  • B-Spline cubique (noyau de 4x4) : approximation par B-Spline cubique

  • Lanczos (noyau de 6x6) : interpolation Lanczos avec fenêtrage avec la fonction sinus cardinal

  • Moyenne : calcule la moyenne de tous les pixels qui ne sont pas des NODATA

  • Mode : sélectionne la valeur la plus fréquente parmi les pixels

  • Maximum, Minimum, Médiane, Premier Quartile (Q1) or Troisième Quartile (Q3) de tous les pixels qui ne sont pas des NODATA

../../../_images/raster_align_edit.png

Fig. 16.21 Sélection des options de ré-échantillonnage raster

Dans la fenêtre Aligner les raster, vous pouvez symbologyEdit Modifier les paramètres de traitement d’un fichier ou symbologyRemove Supprimer un fichier existant de la liste des couches raster.Vous pouvez également choisir parmi les options suivantes (voir Fig. 16.22) :

  • sélectionner une Couche de référence,

  • re-projeter dans un nouveau SCR,

  • choisir une Taille de cellule différente,

  • choisir un Décalage de la grille,

  • Découper selon l’emprise : elle peut être définie par l’utilisateur, basée sur une des couches ou sur l’emprise de la carte.

  • Taille de sortie,

  • Ajouter les rasters alignés au canevas de la carte.

../../../_images/raster_align.png

Fig. 16.22 Alignement de rasters