Outdated version of the documentation. Find the latest one here.

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 MNE 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.

../../../_images/medfordarea.png

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 la boîte englobante de la couche avec les limites de la zone de la ville, et ensuite nous lui mettons un tampon, afin d’obtenir une couche raster qui couvre un peu plus que le strict nécessaire.

Pour calculer la boîte englobante, nous pouvons utiliser l’algorithme Créer un polygone à partir de l’emprise de la couche.

../../../_images/bbox.png

Pour lui mettre un tampon, nous utilisons l’algorithme Tampon à distance fixe, avec les valeurs de paramètres suivantes.

../../../_images/buffer_dialog.png

Voici la boîte englobante de résultat obtenue en utilisant les paramètres montrés au-dessus.

../../../_images/buffer1.png

C’est une boîte ronde, mais nous pouvons facilement obtenir la boîte équivalente avec des angles droits, en lançant pour elle l’algorithme Créer un polygone à partir de l’emprise de la couche. Nous pourrions avoir fait un tampon sur les limites de la ville d’abord, et ensuite calculer l’emprise du rectangle, en sauvegardant à chaque étape.

../../../_images/buffer_squared.png

Vous remarquerez que les rasters ont une projection différente que celle du vecteur. Nous devons donc les reprojeter avant de poursuivre, en utilisant l’outil Projection (reprojection).

../../../_images/warp1.png

Avec cette couche qui contient la boîte englobante de la couche raster que nous voulons obtenir, nous pouvons recadrer les deux couches raster, en utilisant l’algorithme Découper une grille avec des polygones.

../../../_images/clip1.png

Once the layers have been cropped, they can be merged using the Merge raster layers algorithm

../../../_images/merge1.png

A cellsize is needed for the merged layer. We will use the same one of the input ones. You do not need to know it in advance before calling the algorithm. Just click on the button in the right–hand size of the text field and you will have a dialog to enter small mathematical formulas, and a list of frequently used values, among them the cellsizes and bounding coordinates of all available layers.

Note: You can save time merging first and then cropping, and you will avoid calling the clipping algorith twice. However, if there are several layers to merge and they have a rather big size, you will end up with a large layer than it can later be difficult to process. In that case, you might have to call the clipping algorithm several times, which might be time consuming, but don’t worry, we will soon see that there are some additional tools to automate that operation. In this example, we just have two layers, so you shouldn’t worry about that now.

With that, we get the final DEM we want.

../../../_images/finaldem.png

Now it is time to compute the slope layer.

A slope layer can be computed with the Slope, Aspect, Curvature 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, so we can then correctly calculate the slope.

Here is the reprojected DEM.

../../../_images/warpeddem.png

With the new DEM, slope can now be computed.

../../../_images/slope.png

And here is the resulting slope layer.

../../../_images/slopereproj.png

The slope produced by the Slope, Aspect, Curvature algorithm is expressed in radians, but degrees are a more practical and common unit. 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).

../../../_images/metricconversions.png

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

docs/training_manual/processing/img/cutting_merging/reproject_back.png

The reprojection processes have caused the final layer to contain data outside the bounding box that we calculated in one of the first steps. This can be solved by clipping it again, as we did to obtain the base DEM.