16.2. Analyse Raster
16.2.1. Calculatrice Raster
La du menu 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.
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:
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 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
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 et cliquez sur 
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
Fig. 16.21 Sélection des options de ré-échantillonnage raster
Dans la fenêtre Aligner les raster, vous pouvez 
Modifier les paramètres de traitement d’un fichier ou 
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.
Fig. 16.22 Alignement de rasters