Important

La traduction est le fruit d’un effort communautaire auquel vous pouvez vous joindre. Cette page est actuellement traduite à 66.67%.

17.15. Découpage et fusion de couches raster

Note

Dans cette leçon, nous verrons un autre exemple de préparation de données spatiales, pour continuer à utiliser des géoalgorithmes dans des scénarios du monde réel.

Pour cette leçon, nous allons calculer une couche de pente pour une zone entourant une zone de ville, qui est donnée dans une couche vecteur avec un seul polygone. Le MNT de base est divisé en deux couches raster qui, ensemble, couvrent une zone beaucoup plus grande que celle autour de la ville avec laquelle nous voulons travailler. Si vous ouvrez le projet correspondant à cette leçon, vous verrez quelque chose comme ça.

Ces couches ont deux problèmes :

• Elles couvrent une aire qui est trop grande pour ce que l’on veut (nous nous intéressons à une plus petite région autour du centre ville)

• Elles sont dans deux fichiers différents (les limites de la ville tombent dans une seule couche raster, mais, comme cela a été dit, nous voulons des zones supplémentaires autour d’elle).

Ces deux problèmes sont facilement résolvables avec les géoalgorithmes appropriés.

Premièrement, nous créons un rectangle définissant la zone que nous voulons. Pour faire cela, nous créons une couche contenant l’emprise rectangulaire de la couche avec les limites de la zone de la ville, et ensuite nous lui ajoutons un tampon, afin d’obtenir une couche raster qui couvre un peu plus que le strict nécessaire.

To calculate the bounding box , we can use the Extract layer extent algorithm

To buffer it, we use the Buffer algorithm, with the following parameter values.

Avertissement

La syntaxe a changé dans les récentes versions; définissez la Distance et l’Arc de vertex à .25

Voici l’emprise résultante obtenue en utilisant les paramètres montrés au-dessus.

It is a rounded box, but we can easily get the equivalent box with square angles, by running the Extract layer extent algorithm on it. We could have buffered the city limits first, and then calculate the extent rectangle, saving one step.

You will notice that the rasters has a different projection from the vector. We should therefore reproject them before proceeding further, using the Warp (reproject) tool.

Note

Les versions récentes ont des interfaces plus complexes. Sélectionnez au moins une méthode de compression.

With this layer that contains the bounding box of the raster layer that we want to obtain, we can crop both of the raster layers, using the Clip raster by mask layer algorithm.

Once the layers have been cropped, they can be merged using the SAGA Mosaic raster layers algorithm.

Note

Vous pouvez gagner du temps en fusionnant d’abord et en recadrant ensuite, et vous éviterez d’appeler deux fois l’algorithme de découpage. Cependant, s’il y a plusieurs couches à fusionner et qu’elles ont une assez grande taille, vous finirez avec une grande couche qui pourra être difficile à traiter plus tard. Dans ce cas, vous devriez appeler l’algorithme de découpage plusieurs fois, ce qui peut prendre beaucoup de temps, mais ne vous inquiétez pas, nous verrons bientôt qu’il existe des outils supplémentaires pour automatiser cette opération. Dans cet exemple, nous n’avons que deux couches, donc vous n’aurez pas besoin de vous soucier de cela maintenant.

Avec cela, nous obtenons le MNT final que nous voulons.

Il est maintenant temps de calculer la couche de pente.

A slope layer can be computed with the Slope algorithm, but the DEM obtained in the last step is not suitable as input, since elevation values are in meters but cellsize is not expressed in meters (the layer uses a CRS with geographic coordinates). A reprojection is needed. To reproject a raster layer, the Warp (reproject) algorithm can be used again. We reproject into a CRS with meters as units (e.g. 3857), so we can then correctly calculate the slope, with either SAGA or GDAL.

Avec ce nouveau MNT, la pente peut maintenant être calculée.

Et voici la couche de pente de résultat.

The slope produced by the Slope algorithm can be expressed in degrees or radians; degrees are a more practical and common unit. In case you calculated it in radians, the Metric conversions algorithm will help us to do the conversion (but in case you didn’t know that algorithm existed, you could use the raster calculator that we have already used).

Reprojecting the converted slope layer back with the Reproject raster layer, we get the final layer we wanted.

Les processus de reprojection ont peut-être généré une couche finale qui contient des données en dehors de l’emprise que nous avons calculée dans une des premières étapes. Cela peut être résolu en la découpant à nouveau, comme nous l’avons fait pour obtenir le MNT de base.