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.
Pour calculer l’emprise rectangulaire, nous pouvons utiliser l’algorithme Créer un polygone à partir de l’emprise d’une couche.
Pour lui ajouter un tampon, nous utilisons l’algorithme Tampon à distance fixe, avec les valeurs de paramètres suivantes.
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.
Il s’agit d’un rectangle arrondi, mais nous pouvons facilement obtenir l’emprise équivalente avec des angles droits, en lançant sur elle l’algorithme Créer un polygone à partir de l’emprise de la couche. Nous aurions pu créer un tampon sur les limites de la ville d’abord et ensuite calculer l’emprise rectangulaire, ce qui nous aurait économisé une étape.
Vous remarquerez que les rasters ont une projection différente de celle du vecteur. Nous devons donc les reprojeter avant de poursuivre, en utilisant l’outil Projection (reprojection).
Note
Les versions récentes ont des interfaces plus complexes. Sélectionnez au moins une méthode de compression.
Avec cette couche qui contient l’emprise rectangulaire de la couche raster que nous voulons obtenir, nous pouvons recadrer les deux couches rasters, en utilisant l’algorithme Découper un raster selon une couche de masque.
Une fois que les couches ont été recadrées, elles peuvent être fusionnées en utilisant l’algorithme Mosaic raster layers de SAGA.
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.
Une couche de pente peut être calculée avec l’algorithme Pente, Aspect, Courbes, mais le MNT obtenu dans la dernière étape n’est pas approprié comme entrée, puisque les valeurs d’élévation sont en mètre mais la taille de cellule n’est pas exprimée en mètre (la couche utilise un SCR avec des coordonnées géographiques). Une reprojection est nécessaire. Pour reprojeter une couche raster, l’algorithme Projection (reprojection) peut encore être utilisé. Nous reprojetons dans un SCR avec des mètres comme unités (par exemple 3857), donc nous pouvons ensuite calculer correctement la pente, avec soit SAGA soit GDAL.
Avec ce nouveau MNT, la pente peut maintenant être calculée.
Et voici la couche de pente de résultat.
La pente produite par l’algorithme Pente, Aspect, Courbes peut être exprimée en radians, mais les degrés sont une unité plus pratique et plus commune. L’algorithme Conversions métriques nous aidera à faire la conversion (mais dans le cas où vous ne sauriez pas que cet algorithme existe, vous pouvez toujours utiliser la calculatrice raster que nous avons déjà manipulée).
En reprojetant la couche de pente reconvertie avec la Reprojection de couche raster, nous obtenons la couche finale que nous souhaitons.
Avertissement
A faire : Ajouter une image
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.