Importante
La traducción es un esfuerzo comunitario puede unirse. Esta página está actualmente traducida en |progreso de traducción|.
17.15. Cortando y uniendo capas raster
Nota
En esta lección veremos otro ejemplo de preparación de datos espaciales, para continuar usando geoalgoritmos en escenarios del mundo real.
Para esta lección, vamos a calcular una capa de pendiente para un área que rodea el área de una ciudad, que se da en una capa vectorial con un solo polígono. El MDE base se divide en dos capas ráster que, juntas, cubren un área mucho mayor que la que rodea la ciudad con la que queremos trabajar. Si abre el proyecto correspondiente a esta lección, verá algo como esto.
Estas capas tienen dos problemas:
Cubren un área que es demasiado grande para lo que queremos (estamos interesados en una región más pequeña alrededor del centro de la ciudad)
Están en dos archivos diferentes (los límites de la ciudad caen en una sola capa ráster, pero, como se ha dicho, queremos un área adicional a su alrededor).
Ambos se pueden resolver fácilmente con los geoalgoritmos adecuados.
Primero, creamos un rectángulo que define el área que queremos. Para hacerlo, creamos una capa que contiene el cuadro delimitador de la capa con los límites del área de la ciudad, y luego lo almacenamos en búfer, para tener una capa ráster que cubra un poco más de lo estrictamente necesario.
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.
Advertencia
La sintaxis cambió en versiones recientes; establezca la distancia y vértice del arco a .25
Aquí está el cuadro delimitador resultante obtenido usando los parámetros que se muestran arriba
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.
Nota
Las versiones recientes tienen una interfaz más compleja. Asegúrese de seleccionar al menos un método de compresión.
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.
Nota
Puede ahorrar tiempo fusionando primero y luego recortando, y evitará llamar dos veces al algoritmo de recorte. Sin embargo, si hay varias capas para fusionar y tienen un tamaño bastante grande, terminará con una capa más grande que luego puede ser difícil de procesar. En ese caso, es posible que tenga que llamar al algoritmo de recorte varias veces, lo que puede llevar mucho tiempo, pero no se preocupe, pronto veremos que hay algunas herramientas adicionales para automatizar esa operación. En este ejemplo, solo tenemos dos capas, por lo que no debería preocuparse por eso ahora.
Con eso, obtenemos el MDE final que queremos.
Ahora es el momento de calcular la capa de pendiente.
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.
Con el nuevo MDE, ahora se puede calcular la pendiente.
Y aquí está la capa de pendiente resultante.
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.
Los procesos de reproyección pueden haber causado que la capa final contenga datos fuera del cuadro delimitador que calculamos en uno de los primeros pasos. Esto se puede solucionar recortándolo de nuevo, como hicimos para obtener el MDE base.