17.15. Découpage et fusion de couches raster¶

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.

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

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 Fusionner de GDAL.

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.

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.