Tutorial: Delineación de la Cuenca y las corrientes
Sitio: | OpenCourseWare for GIS |
Curso: | QGIS para Aplicaciones Hidrológicas |
Libro: | Tutorial: Delineación de la Cuenca y las corrientes |
Imprimido por: | Guest user |
Día: | lunes, 30 de diciembre de 2024, 19:06 |
Descripción
Tabla de contenidos
- 1. Introducción
- 2. Teoría
- 3. Descarga de mosaicos (Tiles) DEM
- 4. Unir los mosaicos DEM
- 5. Reproyecta el DEM
- 6. Creación de un recorte del DEM
- 7. Rellenado de huecos (depresiones)
- 8. Delineación de corrientes
- 9. Definir un punto de flujo de salida
- 10. Delineación de la Cuenca
- 11. Almacenamiento de datos en un GeoPackage
- 12. Estilizado de la Cuenca resultante
- 13. Conclusiones
1. Introducción
Uno de los usos más importantes del SIG en Hidrología es la delineación de Cuencas y corrientes. esta lección presenta un flujo de trabajo genérico para la delineación de corrientes y cuencas para áreas donde existen datos abiertos disponibles: El modelo de elevación SRTM (Shuttle Radar Topography Mission) 1 Arc Second DEM. El flujo de trabajo será aplicado en el caso de la cuenca Rur en Alemania.
Despúes de esta lección, serás capaz de:
- Encontrar y descargar tiles SRTM desde el sitio web EarthExplorer
- Hacer una sola capa de los mosaicos de capas ráster en un ráster virtual
- Reproyectar capas ráster
- Crear subconjuntos de ráster
- Llenar hundimientos (sinks) en el DEM
- Calcular y estilizar la dirección de flujos (flow direction) a partir de capas ráster
- Calcular el orden Strahler
- Delinear corrientes (streams)
- Delinear cuencas (catchments)
- Estilizar capas para la visualización de las cuencas
Los pasos para delinear una Cuenca a partir de un DEM en el SIG tienen el siguiente orden (los detalles se tratan en esta lección):
1. Descarga los Tiles del DEM de tu área de estudio. Asegurate que los mosaicos (Tiles) cubren completamente el área de estudio y la cuenca que desea delinear. Se sugiere que sea un poco más grande que evite efectos del contorno.
2. Si el área de estudio se cubre por múltiples mosaicos (Tiles) DEM, requieres unir en una sola capa los mosaicos para crear una sola capa ráster DEM.
3. Los mosaicos DEM pueden estar en un Sistema de Coordenadas diferente al deseado, en este caso, puedes reproyectar la capa DEM a una proyección que se desees utilizar.
4. En el caso que los mosaicos fusionados sea mayor que la zona de estudio, puedes recortar con una operación Clip a un área más pequeña, con el fin de reducir el tiempo de cálculo.
5. Interpolar huecos en el DEM, en el caso de que existan pixeles con No Data.
6. Realizar un DEM hidrológicamente correcto con el llenado de huecos y eliminación de picos.
7. Calcular la dirección de flujo de cada celda.
8. Obtener la red de drenaje.
9. Delinear la Cuenca considerando el punto de salida de la Cuenca.
Cuando ya tienes un mapa con la red de corriente, el procedimiento puede mejorar con el quemado (burning) de la red fluvial al DEM. De esta forma, los valores de altura del DEM es menor a la altura de los ríos, la escorrentía seguirá el flujo de la actual. Este proceso se encuentra fuera del alcance de este capítulo.
Los materiales para esta lección pueden descargarse desde el menú principal. La siguiente sección muestra la forma de descargar los mosaicos (Tiles) SRTM desde USGS EarthExplorer, los cuales también se incluyen en los datos del curso.
2. Teoría
2.1. Los Modelos Digitales de Elevación en un SIG
Revisa el video acerca de los Modelos Digitales de Elevación en un SIG
2.2. Delineación de la Cuenca y corrientes
Estudia este video acerca de la teoría de la delineación de corrientes y cuencas
3. Descarga de mosaicos (Tiles) DEM
Para la zona de estudio Rur, descargaremos los mosaicos del Conjunto de Datos Globales SRTM 1 Arc-Second. Desde finales del año 2014, el Modelo de Elevación de 1-Arc second (aproximadamente 30 metros en el Ecuador) ha sido publicado como Datos Abiertos. Este Dataset cubre la mayor parte de la Tierra cubriendo desde 54 grados Latitud Sur y 60 grados latitud Norte. Una descripción de los productos SRTM se encuentran aquí.
Los siguientes pasos explican como descargar los mosaicos (Tiles) DEM SRTM para el área de estudio.
1. Abre un navegador y accede a la siguiente dirección http://earthexplorer.usgs.gov
La descarga de los mosaicos DEM STRM requiere una cuenta. En caso de fallar en la creación de una cuenta puedes usar los mosaicos SRTM proporcionados con los datos de la actual lección (4 GeoTiff files).
2. Haz Click Register y sigue el procedimiento. Después del registro, inicia sesión en el portal del USGS EarthExplorer.
Ahora, requieres indicar el parea de búsqueda. En EarthExplorer puede hacerse de formas distintas: (1) Búsqueda del nombre de lugares , (2) Acercarte en el mapa y usar las coordenas de la zona de acercamiento, (3) Cargar un archivo KML o Shapefile comprimido. En este caso, prepararemos un archivo Shapefile con un cuadro de delimitación (bounding box).
3. Escoge un archivo Shapefile y "sube" el archivo boundingbox.zip (requieres comprimir todos los archivos que pertenecen al Shapefile boundingbox).
4. Escoge Data Sets para ir a la siguiente pantalla donde seleccionaremos el Conjunto de datos (data set).
5. Da un Click sobre el signo de más (+) para elegir Digital Elevation | SRTM | SRTM 1-Arc Second Global. Asegurate de verificar la casilla.
6. Da Click en Results para navegar a la siguiente pantalla.
Ahora, observa los resultados de la búsqueda: cuatro Tiles. Antes de la descarga, puedes verificar los mosaicos (Tiles); con el icono de pie (foot), por ejemplo, puedes observar la miniatura del mosaico.
7. Descarga cada mosaico (Tile) con un click en el botón Download Options .
8. En la ventana de las opciones de descarga (Download Options) descarga el DEM en formato GeoTIFF en el folder del ejercicio en tu Disco duro.
Los mosaicos (tiles) pueden descargarse a tu proyecto QGIS.
9. Ejecuta QGIS Desktop.
Iniciaremos con un proyecto nuevo.
10. En el panel de navegador (Browser) agrega la carpeta con los datos para este capítulo a los Favoritos (Favorites).
11. Da un click derecho sobre uno de los mosaicos y selecciona las propiedades (Properties) para revisar los metadatos.
¿Cuál es la proyección?¿Cuales son las unidades del Mapa?
12. Selecciona y arrastra los cuatro mosaicos al Map Canvas.
4. Unir los mosaicos DEM
Antes de seguir, tenemos que fusionar los cuatro mosaicos DEM, el cual en la terminología SIG es llamado Mosaico. Hay dos formas en que se hace esto con los mosaicos (Tiles):
- Unir los mosaicos en un sólo archivo ráster (Por ejemplo GeoTIFF)
- unir los mosaicos en un archivo virtual (.vrt)
2. En la ventana de "Construye una Capa virtual" (Build Virtual Raster), selecciona cada archivo en forma individual o fusiona todos los archivos en un directorio. Es posible, la unión de los archivos visibles desde el Map Canvas. Utilizaremos la última opción:
- En las capas de entrada, haga click en .
- Usa el botón Select all que permite la selección de los cuatro mosaicos, después haga click en el botón OK (figura inferior).
- Navega a la ubicación donde se desees guardar los archivos de salida y nombra al archivo como
dem_mosaic.vrt
.
- Por default, la resolución (Resolution) se encuentra en el promedio (
average)
. En nuestro caso los archivos tienen la misma resolución(1 Arc Second).
- Deshabilita el checkbox en Place each input file into a separate band. Esta opción requiere verificarse sólo si requieres la creación de un mapstack, por ejemplo con bandas de las imágenes de Percepción Remota.
- Mantén el algoritmo Resampling algorithm como default a:
nearest
, el cual no tiene impacto porque el remuestreo no es requerido.
3. Haga Click Run que corre el algortimo. Después con un Click en Close regresa a la pantalla principal donde puedes observar el DEM unido.
Observa que en el Map Canvas, los bordes de los Tiles no son visibles más en el DEM unido porque QGIS "estira" la escala de grises utilizando el mínimo y máximo del DEM unido. Sólo es para la visualización, los valores en los Tiles son los mismos que en el Mosaico completo.
4. Ahora borra los Tiles individuales (not the dem_mosaic). Desde la lista de capas por medio de la selección mientras mantienes presionado el botón
Ctrl
. Así, con el Click-derecho sobre las capas seleccionadas, escoge la opción Remove Layer. Haz Click OK para confirmar. Esta acción borrará los Tiles en el monitor, pero permanecerá en el Disco Duro.5. Reproyecta el DEM
El DEM se encuentra en el sistema de coordenadas geográficas (GCS) Lat/Lon con datum WGS 84 (EPSG: 4326). No es recomendable el uso del GCS en un análisis DEM debido a que las unidades Z (por ejemplo metros) son diferentes que las unidades x y y (degrees). Se requiere elegir una proyección para nuestro proyecto, el cual, si cubre una proyección nacional. En nuestro caso, sin embargo, la proyección abarca múltiples países (Alemania, Holanda y Bélgica). De esta forma, utilizaremos una proyección global: UTM Zone 32 North, con WGS-84 como datum.
Buscaremos los códigos EPSG en http://www.spatialreference.org
1. Utilizaremos el Sitio web para buscar UTM 32N wgs 84
. Deja corriendo QGIS y abre un navegador.
La figura inferior muestra los resultados:
Revisemos los códigos EPSG. Utilizaremos el código EPSG: 32632 a lo largo de este proyecto; al hacer Click sobre él se proporcionarás más detalles.
Así, reproyectemos el DEM desde un Sistema Coordenado sin proyección GCS (Lat/Lon WGS 84 - EPSG: 4326) a UTM Zone 32 North / WGS 84 (EPSG: 32632).
2. From the main menu choose Raster | Projections | Warp (Reproject).
3. En la ventana de Reproyección Warp (reproject), Haz click para escoger Target CRS.
4. En la ventana de diálogo, teclea 32632
en el campo de Filter y selecciona WGS EPSG: 32632 en la sección de enmedio de la ventana, después da click en OK.
5. Ahora, completa en la ventana:
- En el método de muestreo (Resampling) selecciona
Nearest Neighbor
que preserva los valores de elevación de los archivos originales. - Coloca los valores de NoData (No data value for output bands) a
-9999
porque la capa ráster reproyectada tiene valores NoData en los bordes. De esta forma, definimos "no data" que tiene un valor -9999 y no se visualicen y procesen como "data". - Indica la resolución del archivo de salida (Output file resolution) a
30
m. - Navega a la carpeta de ejercicios y nombra el archivo de salida como
dem_reprojected.tif
La ventana de diálogo se muestra como la imagen inferior.
Observa que el comando gdalwarp
que se ejecutará en background.
6. Haz Click en Run que correo el algoritmo. Después de que termine el proceso da click en el botón Close que cierra la ventana.
El DEM reproyectado ahora aparece en la lista de capas. El DEM podría lucir distorcionado a pesar de que la Proyección en el Map Canvas aún se encuentra en Lat/Lon, como se indica en las coordenadas en la parte inferior del Map Canvas y el código EPSG del Map Canvas en la esquina inferior derecha. A pesar de la proyección al vuelo, lo cual, causa diferencias entre la proyección en el archivo y la capa visualizada.
7. Para completar esta operación y mostrar correctamente el nuevo conjunto de datos, es necesario, cambiar la proyección al vuelo. Haga click en EPSG:4326 en la parte inferior derecha de la pantalla.
8. En la ventada de diálogo selecciona la proyección EPSG:32632 que se encuentra en la lista de Sistemas de Coordenadas recientemente utilizados (Recently used coordinate reference systems), como se muestra en la figura inferior, y despúes haz un click en el botón OK.
En el Map Canvas, el DEM es ahora desplegada en la posición correcta( Con el Norte apuntando hacia arriba), como se muestra en la figura inferior. Las coordenadas del nuevo despliegue están en metros como se requiere.
9. Ahora, puedes borrar la capa dem_mosaic.
Es un buen momento, para salvar tu proyecto.
6. Creación de un recorte del DEM
Para reducir el cálculo de procesamiento de los algoritmos, realizaremos un subconjunto ( recorte) de la capa ráster, es importante realizar un límite un poco más allá de la zona de estudio porque esto nos ayuda en la prevención de efectos de borde.
Una forma fácil es la selección del límite de la zona de estudio es el uso de OpenStreetMap. OpenStreetMap contiene datos de fuentes múltiples. Si el complemento QuickMapServices se encuentra instalado, puedes agregar el mapa base OpenStreetMap a QGIS como sigue:
1. Desde el menú principal, navega Web | QuickMapServices | OSM | OSM Standard.
2. Así, la capa de OpenStreetMap se desplegará en el Map Canvas, temporalmente deshabilita el checkbox de la capa DEM o arrastra la capa OSM en la parte superior.
Para conocer la extensión aproximada de la Cuenca Rur, los datos del ejercicio contienen un archivo shapefile comprimido con el bounding box (el mismo que utilizamos para la descarga de los mosaicos DEM ). Considera que en la realidad no habrá un cuadro delimitador (bounding box), pero utiliza tu experiencia que ubiquen los límites aproximados de la zona de estudio.
3. Descomprime boundingbox.zip
en tu carpeta de ejercicios.
4. Arrastra la capa boundingbox.shp
en el Map Canvas.
5. Ahora, desde el menú principal selecciona Raster | Extraction | Clip Raster by Mask Layer.
6. En la ventana de Clip Raster by Mask Layer escoge dem_reprojected como la capa de entrada (Input layer). En el caso de la capa de mascara (Mask Layer), selecciona boundingbox. Mantén las demás opciones por default. Escoge la capa dem_subset
para el corte de salida Clipped (mask) y haz click en el botón Run. Despúes de hacer el proceso, cierra la ventana con un Click en Close.
Frecuentemente, no tendrás el shapefile del bounding. En este caso, puedes dirigirte a Raster | Extraction | Clip Raster by Extent. Entonces arrastra un cuadro en el Map Canvas y usalo para hacer un corte. Asegurate que la reproyección al vuelo es similar a la de la capa que quieres cortar porque las coordenadas del Map Canvas son utilizadas por el algoritmo.
7. Así puede borrar la capa dem_reprojected de la lista como hemos hecho anteriormentes desde otras capas que no se requieren más.
7. Rellenado de huecos (depresiones)
Los DEM crudos o sin procesamiento tienen artefactos como las depresiones. Los artefactos son resultado del proceso de adquisición y deben ser removidos antes de utilizar un DEM en el análisis hidrológico como la delineación de Cuencas y corrientes o el modelado hidrológico. Existen varios agoritmos para el llenado de huecos. Utilizaremos del algoritmo desarrollado por Wang and Liu (2006), el cual, es mucho más rápido que otros y funciona mejor en conjunto de datos de alta resolución. Este algortimo, implementado en las herramientas de SAGA, no sólo llena despresiones, también preserva la pendiente descendiente a lo largo de la trayectoria, logra esto preservando el gradiente de pendiente mínima entre celdas.
1. Usaremos la caja de herramientas (Processing Toolbox), en caso de que la sección no sea visible en la parte derecha de la ventana, habilitala desde el menú principal, seleccionando Processing | Toolbox.
2. En la caja de herramientas (Processing Toolbox) utiliza el campo de búsqueda Search y teclea fill sinks
.
Observarás que hay más algoritmos para el llenado de huecos, cada uno de ellos tienen sus ventajas y desventajas en términos de l tiempo de cálculo, uso de memoria y exactitud. En el caso de tu investigación requieres tomar decisiones basadas en la documentación o resultados de cada uno.
3. Sigue la siguiente ruta, ejecuta y haz click en SAGA | Terrain Analysis - Hydrology | Fill sinks (wang & liu).
4. En la ventana de diálogo, mantén los parámetros por default Minimum Slope [Degree] de 0.01. Asegúrate de seleccionar como entrada DEM dem_subset y como Filled DEM a dem_fill.sdat
. Sólo requerimos Filled DEM en este punto, así que puedes deshabilitar el checkbox para la apertura de la salida de la dirección de flujo (Flow Directions) y cuenca hidrográfica (Watershed Basins) (observa la figura inferior).
5. Después da Click en Run y cuando haya terminado el proceso, cierra la ventana con Close.
6. Ahora, borra la capa dem_subset de la lista de capas, porque ya no son necesarias.
7.1. Estilizado del DEM
Las leyendas por default hechas por el software GIS Desktop, por lo regular, son poco intuitivas porque (1) el softwtare desconoce que información es presentada (2) el tipo de datos ráster (ej. Boleano, discreto o datos continuos). Para apoyar la interpretación de los resultados es buena practica estilizar las capas intuitivamente.
La capa DEM es un ráster continuo que representa gradientes y por lo tanto puede contener números reales (punto flotante). Los ráster continuos son estilizados en QGIS utilizando rampas en la ventana de Singleband pseudocolor.
1. Selecciona la capa dem_filled y haz click en el icono para abrir el panel de estilizado de la capa (Layer Styling) o presiona F7
.
2. En el panel de estilizado (Layer Styling) escoge Singleband pseudocolor desde la lista desplegable.
La leyenda es automáticamente actualizada con la rampa default rojo a azul. Esto obviamente proporciona una impresión incorrecta de que existen lagos en el Sur, porque asociamos el color azul con las zonas bajas con agua. En el siguiente paso aplicaremos colores más intuitivos en el estilizado de los DEMs.
3. Con un click derecho sobre la rampa de color y selecciona Create New Color Ramp.
4. En la ventana emergente (popup) del tipo de rampa de color (Color ramp type) selecciona Catalog: cpt-city desde la lista desplegable.
El catálogo de cpt-city tiene muchas rampas de colores predefinidas que son muy útiles.
5. Selecciona Topography | Elevation. Observa que cd-a y sd-a son también opciones agradables.
Esto nos da colores más intuitivos en el DEM donde podemos distinguir claramente las áreas superiores e inferiores. Ahora mejoraremos aún más la visualización.
6. Da un click sobre la capa dem_fill en el panel Layers y selecciona Duplicate Layer.
Esto creará una copia de la capa dem_fill llamada dem_fill copy.
7. Desmarca la capa dem_fill y renombra la capa dem_fill copy a hillshade.
8. En el panel de estilizado (Layer Styling) panel, aún se encuentra abierto, asegurate de que la capa hillshade ahora este seleccionada. En la lista desplegable cambia Singleband pseudocolor a Hillshade.
Ahora la capa hillshade esta visualizada con un sombreado (shading).
¿De qué dirección proviene la iluminación? ¿Es posible en la realidad?
da los mejores resultados con una iluminación artificial en el noroeste, que en realidad no puede existir en el hemisferio norte.
Hillshade proporciona el mejor resultado con una iluminación artificial en el Noroeste, lo cual en realidad no existen en el Hemisferio Norte. Si mueves el marcador en el panel de estiizado (Layer Styling) al Sureste, observarás un relieve invertido. También observa la sección de muestreo (Resampling). El método de muestreo para el acercamiento (Zoomed in) y alejamiento (Zoomed out) es Vecinos cercanos (Nearest Neighbor). Este método es fino para la categorización de datos sin embargo, la elevación es considerada como datos continuos. Por lo tanto, elige en Zoomed in el méotodo de muestro Bilinear y en Zoomed out, el método de muestreo Average.
A continuación, mezclaremos el DEM con la capa hillshade.
9. Cambia a la capa dem_fill marcando el checkbox.
10. En el panel de estilizado (Layer Styling) asegurate que la capa dem_fill se encuentre seleccionada. En el bloque del panel Layer Rendering , cambia el modo de mezclado (Blending mode) a Multiply.
Como puedes observar, el mezclado produce un efecto mucho mas agradable que la transparencia, la cual, desvanecerá los colores. Ahora vemos claramente las diferencias de elevación: el gradiente de Sur a Norte y los valles donde se encuentran las corrientes.
Los modos de mezclado permiten una representación más elegante entre las capas SIG. Estos pueden ser más vistosos que simplemente ajustar la opacidad de la capa. Los modos de mezclado permiten efectos visuales donde la intensidad total de una capa subyacente es aún visible a través de la capa superior. Existen 13 modos de mezclado disponibles. Revisa aquí para mayor información.
8. Delineación de corrientes
Ahora que tenemos un DEM corregido hidrológicamente, gracias al algoritmo de llenado de huecos, procedemos con la delineación de las corrientes.
En primer lugar, calcularemos las ordenes de Strahler, después calibraremos el umbral de Strahler para determinar las corrientes. Finalmente, derivaremos la red de canales, la dirección de flujo y las cuencas.
8.1. Calibración del umbral de Strahler para determinar las corrientes
El siguiente paso es la aplicación del procedimiento de calibración para determinar el orden Strahler de las corrientes. Crearemos una capa boleana para el orden Strahler de las corrientes más largas, iguales que un valor umbral. Cada capa boleana se compara con una capa de referencia. En este punto, compararemos el orden Strahler de las corrientes con los ríos sobre OpenStreetMap. En áreas donde hay vacios de datos en OpenStreetMap, podrías utilizar Google Satellite; o ambas capas OpenStreetMap y Google Satellite que puedes encontrar desde el complemento QuickMapServices.2. Llena los parámetros de la ventana como se muestra en la figura inferior. Nombre al archivo de salida como
strahler5.tif
y haz un click en OK para ejecutar el cálculo.Es una buena practica estilizar la capa strahler5.
3. Dirígete al panel de estilizado (Layer Styling) y asegúrate que la capa strahler5 se encuentra seleccionada.
4. Selecciona Paletted/Unique values del menú desplegable.
Esta capa es boelana y por lo tanto contiene sólo pixeles con valores 0 y 1. Para la calibración es importante hacer los valores "cero" en transparentes y los valores "uno" en color azul así comparamos el ráster con las corrientes sobre la capa OpenStreetMap.
7. Hacer los pixeles con valores "cero" en transparentes, dirígete a la pestaña transparency y agrega el valor
0
como additional no data value.En el Map Canvas observa ahora todas las corrientes mayores o igual al orden Strahler 5 y podemos compararlas con los ríos en OpenStreetMap.
8. Repite los pasos en esta subsección para diferentes umbrales en los valores de la orden de las corrientes Strahler y determina aquella que corresponda mejor con los ríos sobre OpenStreetMap. Borra las otras capas boleanas.
Tip: Puedes copiar los estilos de las capas: click derecho sobre la capa, selecciona Styles | Copy Style entonces con un click derecho sonbre la capa destino y selecciona Styles | Paste Style. Aún requieres hacer los pixeles con valores zeros transparentes porque no son copiados.
8.2. Cálculo del orden Strahler de las corrientes
Antes de la obtención de las corrientes desde el DEM, requerimos determinar que consideramos como una corriente. Para este propósito utilizaremos el orden Strahler, el orden más alto tendrá la corriente más grande.
1. Teclea Strahler
en el cuadro de búsqueda de la caja de herramientas (Processing Toolbox) y selecciona SAGA | Terrain Analysis - Channels | Strahler order.
2. En la ventana de diálogo, selecciona la capa dem_fill para el campo Elevation, utiliza strahler.sdat
como el nombre del archivo de salida Strahler Order , y da un click en Run. Da un Click en Close cuando el algoritmo haya terminado.
Hacer que la capa strahler tenga más sentido, estilizaremos esta.
Es la capa de ordenes Strahler una capa boleana, discreta o un ráster continuo?
La capa de orden Strahler es una capa discreta, pero el orden de las clases es importante. Para ráster discretos en QGIS utilizaremos el estilizado Paletted/Unique values. El valor mayor del orden, la corriente más grande utilizará la rampa de color desde el color blanco al color azul.
3. Abre el panel de estilizado (Layer Styling) si es que cerraste este antes y asegurate que la capa strahler es seleccionada.
4. Escoge del menú desplegable Paletted/Unique values .
5. Da un Click en Classify.
Cada valor único en el ráster es un color aleatorio único que ha sido asignado.
6. Con un click derecho sobre Random colors y selecciona Blues.
La capa strahler ahora luce de una forma más entendible que las corrientes mayores serán mayores que las corrientes de menor orden.
8.3. Cálculo de la red de canales, dirección de flujo y cuencas.
El siguiente paso es calcular la red de canales basado en el umbral Strahler determinado en pasos anteriores. Con la misma herramienta podemos la capa de dirección de flujo y una capa de todas las cuencas.
En la caja de herramientas (Processing Toolbox) puedes encontrar herramientas relacionadas que utilizan diferentes algoritmos. Siempre, revisa la documentación del algoritmo que te permitirá conocer el uso adecuado de cada uno o hacer análisis comparativos.
1. Busca la palabra Channel
en la caja de herramientas y selecciona SAGA | Terrain Analysis - Channels | Channel network and drainage basin.
2. En la ventana de diálogo, asegurate que seleccionas la capa dem_fill para la elevación y coloca el valor del umbral que encontraste durante la calibración en la sección previa, ej. Significa que las corrientes con un orden Strahler igual o mayor a este valor se consideran como canales. El algoritmo calculará la dirección de flujo y el orden Strahler determinará los canales y cuencas de drenaje. Selecciona los checkboxes para las siguientes salidas:
- Flow direction, guarda como
flowdir.sdat
- Channels, guarda como
channels.shp
- Drainage basins, guarda como
basins.shp
(la otra es para el ráster)
3. Una vez que la ventana de diálogo se encuentre configurado en forma similar a la figura inferior, da click en el botón Run.
La capa de dirección de flujo (flow direction) muestra la dirección de cada pixel, que puede expresarse como la dirección de una brújula, sin embargo no es posible almacenar texto en una capa ráster. La dirección de la brújula puede expresarse en grados sobre un círculo donde el norte tiene 0 grados, el este tiene 90 grados, etc. Para almacenar los 360 grados requerimos más de 8 bits (0-255), que incrementará el tamaño del archivo. Además el método D8 utiliza direcciones discretas con los pixeles circundantes. Por lo tanto el software SIG decodifica las 8 direcciones. Cada software, sin embargo, lo hace de modo propio. El algoritmo de SAGA de la caja de herramientas en QGIS toma el valor 0 para el Norte, 1 para el Noreste, etc; lo cual significa que almacena valores del 0 al 7. El valor de 255 es utilizado en superficies planas.
8.4. Estilizar la capa de canales
A continuación, estilizarás los canales para obtener un mejor entendimiento de los resultados.
1. Abre el panel de estilizado (Layer Styling) y coloca la capa objetivo a channels.
2. Escoge la representación de Graduated con las siguientes configuraciones:
- Coloca Value a ORDER y Precision a 0 debido a que los valores del orden Strahler son enteros.
- Coloca Method a Size.
- Coloca valores Size from desde
0.3
a1.0
mm. - Coloca Mode a Equal Interval con el mismo número de clases de la cantidad de ordenes de Strahler que existen en tus datos. Revisa la tabla de atributos si se requiere. En este ejemplo se utilizan 3 clases. Da un Click en Classify.
- Coloca los valores de la leyenda a valores enteros si es necesario.
3. Haz un Click en el botón Change y coloca Color a un valor de RGB de 15 | 66 | 220.
4. Haz Click en el botón Go back para regresar al menú principal de la simbología.
Ahora, Tu mapa debería asemejarse a la figura inferior, que muestra las corrientes identificadas con ordenes de Strahler superiores son los canales principales y las corrientes con valores de orden Strahler menores son los tributarios.
8.5. Estilizar la capa de dirección de flujo
A continuación, trabajarás sobre el estilo de la capa de dirección de flujo, esto ayudará en la revisión de los resultados y a asegurarse que corresponden con lo esperado. Emplearás una rampa de color circular que mostrará de forma intuitiva los flujos en la dirección Sur con colores cálidos y los flujos orientados al Norte con colores fríos.
1. Utiliza el panel de estilizado (Layer Styling panel) e indica la capa objetivo a flowdir.
2. Coloca la representación a Paletted/Unique values y da un click en Classify.
3. Dado que el valor 255 representa superficies planas, selecciona el valor y haz un click en el botón Minus para borrar este valor de la clasificación.
4. Escoge la rampa de color Spectral como punto de partida. Con un click derecho sobre la rampa de color y escoge Edit Color Ramp.
5. La ventana Select color ramp se abrirá, en esta puedes editar los límites (stops) de color para la rampa. La idea es que el lado extremo derecho e izquierdo representen los flujos hacia el Norte y colocarás estos límites con el mismo color azul. Después coloca en el centro un límite en color amarillo brillante que represente los flujos del Sur. Otro límite en el medio hacia la derecha representará los flujos del Oeste y un límite en medio hacia la izquierda representará los flujos del Este.
6. Da un click en el menú desplegable para Color 1 y selecciona Pick color del menú contextual. Con el cursor de cuentagotas haz un click sobre Color 2.
7. Da un Click then el límite central que representará los flujos del Sur y asigna un amarillo brillante con los valores RGB : 255 | 255 | 0.
8. Para el límite Este (eastern) asigna el color verde, introduce los valores RGB: 0 | 255 | 0.
9. Para el límite Oeste (western) utiliza el color magenta con valores RGB: 214 | 60 | 170.
10. Da un Click OK para aceptar la nueva rampa de color.
11. Da un Click en el botón Plus que agrega un nuevo valor a la clasificación Value con 255. Este representa las zonas planas que fueron borradas antes de la creación de la rampa de color. Da un click en el círculo y selecciona white.
12. Para revisar los datos flowdir, ahora en el modo de mezclado (Blending mode) de Multiply y habilita la capa hillshade. Esta acción nos permite observar los detalles del hillshade a través de la capa flowdir . En este punto, ya contamos con las capas channels, flowdir y hillshade habilitadas. Acercate en algunas montañas al Sur de la zona de estudio donde hay canales. Inspecciona los datos.
Deberías observar las laderas orientadas hacia el Sur representadas por colores cálidos, los flujos del Este tienen color más verdes y los flujos del Oeste de un color púrpura.
9. Definir un punto de flujo de salida
Una cuenca es un área o extensión de superficie donde el agua superficial proveniente de la lluvia, derretimiento de la nieve o hielo, convergen en un punto único en una elevación baja, por lo regular la salida de la cuenca, donde las corrientes de agua se unen con otros cuerpos de agua como ríos, lagos, reservorios, estuarios, humedales, el mar u oceáno. Para delinear la cuenca requerimos:
- las coordenadas de nuestro punto de salida con el mismo sistema coordenado del mapa que estamos utilizando.
- la red de canales que coincide con las direcciones de flujo calculado a partir de un DEM hidrológicamente corregido.
El punto de salida en la Cuenca Rur esta en Roermond, donde el Río Rur encuentra al Río Meuse (Maas in Holandés). La red de canales que ha sido hecha en pasos previos se encuentra en la capa channels.
1. Asegurate de que la capa channels aparece en la parte superior, arriba de la capa OSM Standard agregada desde el complemento QuickMapServices.
2. Busca la ubicación donde el Río Rur desemboca en el Río Meuse. Considera que la capaNote de OpenStreetMap utliza los nombre en holandés porque se encuentran en los Países Bajos. Rur se escribe como Roer y Meuse se escribe comoMaas.
Considera que los canales delineados no corresponden totalmente con los canales sobre OpenStreetMap. Puede deberse a las siguientes razones: (1) Delineación automática incorrecta de las corrientes porque se hizo con un DEM con errores o áreas que son muy planas (2) Distorsión a pesar de la reproyección (al vuelo) y (3) la infuencia humana sobre el curso natural de los canales. La delineación de cuencas, sin embargo, sólo funciona cuando un punto de salida es definido sobre la red de canales porque corresponde con el DEM rellenado.
3. Selecciona un punto sobre el canal delineado que se encuentre cercano a un punto de salida real del Río Rur al Río Meuse (step 2).
4. Utiliza la herramienta Coordinate Capture para obtener las coordenadas del punto de salida sobre la capa channels . La herramienta Coordinate Capture la encuentras desde el menú principal Vector | Coordinate capture. Observa que el complemento requiere activarse (habilita el checkbox en el administrador de complementos - Plugins Manager). El panel de Coordinate Capture aparecerá.
5. Da un Click sobre Start capture y haz un click en un punto sobre el río delineado cercano a la desembocadura del Río Rur como se muestra en la figura inferior. Asegúrate de capturar el punto desde el Río en la capa channels, porque el mapa base (OpenStreetMap) no ha sido derivado desde el DEM y habrá errores en la delineación de la cuenca.
¿Qué proyección estamos utilizando?¿Cómo sabemos esto?
10. Delineación de la Cuenca
Ahora, utilizaremos estas coordenas para el cálculo de el área aguas arriba (ej. cuenca de captación) que produce la descarga en este punto.
1. En la caja de herramientas (Processing Toolbox) busca Upslope Area
y selecciona SAGA | Terrain Analysis - Hydrology | Upslope Area.
2. Llena la ventana de diálogo como se muestra en la figura inferior. Utiliza las coordenas que que capturaste, cercano al punto de salida. Para Elevation utiliz la capa dem_fill. Utiliza el método por dafault D8. Nombra al archivo de salida Upslope Area como catchment.sdat
. Da Click en Run. Después del procesamiento cierra la ventana con un Click en Close.
El resultado luce como la pantalla en la figura inferior, sí te acercas a la extensión completa de la capa.
En el caso dela superposición del límite de la cuenca con otros datos, es mejor convertirlo de ráster a vector (polígono).
3. Para convertir la capa ráster a una capa vectorial, dirígete desde el menú principal y selecciona Raster | Conversion | Polygonize (Raster to vector).
4. Asegurate de elegir la entrada adecuada y nombra la salida como Rur_catchment.shp
. Da Click en Run. Después del procesamiento, cierra la ventana principal con un Click en Close.
5. Observa el resultado. También verifica la tabla de atributos (con un click derecho sobre el nombre de la capa y selecciona Open attribute table).
¿Qué valor se asigna a la celdas dentro del polígono de la Cuenca de captación?¿Qué valores se encuentran fuera?
En el cálculo de la Upslope Area, las celdas dentro de la Cuenca con un valor de 100 mientras que las demás celdas contendrán un valor de 0. Durante la conversión de polígonos, puede haber errores de geometría (topología incorrecta) porque el límite del polígono hace un loop. Esto puede ocasionar errores cuando utilizamos el polígono para el procesamiento.Observa la figura inferior.
Por supuesto, estamos interesados sólo en el área de la cuenca, por lo que eliminaremos el polígono exterior.
6. En la tabla de atributos, In the attribute table, cambia al modo de edición utilizando el botón , después selecciona el registro que deseas borrar con un click sobre el renglón. La selección resaltará en amarillo sobre el mapa. Da un Click en el icono para borrar el objeto seleccionado, desactiva la edición con un click otra vez , y guarda los cambios.
7. Ahora borra todas las capas innecesarias de la lista de capas para que sólo tengamos las capas channels, Rur_catchment, dem_fill, y OSM Standard (en ese orden).
Para la visualización es más agradable recortar las capas al límite de la Cuenca.
8. Primero recortemos la capa vectorial channels para ver sólo las corrientes que se encuentran dentro de la Cuenca. Dirígete desde el menú principal y selecciona Vector | Geoprocessing Tools | Clip.
9. Llena la ventana de diálogo como en la figura inferior, utiliza la capa de la cuenca como "cookie cutter" que corta la capa de entrada Input layer channels al límite de la capa Overlay layer Rur_catchment. Nombra la capa Clipped como channels_clip.shp
. Da Click en Run para correr la herramienta. Da un Click en Close para regresar a la ventana principal.
Si obtienes un error relacionado con la geometría, puedes arreglarlo. Desde la caja de herramientas (Processing Toolbox) selecciona Vector geometry | Fix geometries. Otra solución es utilizar las herramientas desde la caja de herramientas que son menos sensibles a los errores geométricos.
10. Borra la capa channels desde la lista de capas y arrástra la capa channels_clip en la parte superior de la lista.
11. Corta el DEM en la misma forma como vimos en section 6 e a Subset of the DEM). La única diferencia es que indicaremos valores de NoData a aquellos valores fuera de rango de las elevaciones. Puedes utilizar el valor -9999
.
11. Almacenamiento de datos en un GeoPackage
Mantener los datos juntos y facilitar la distribución es recomendable guardar las capas como un GeoPackage.
1. En la caja de herramientas (Processing Toolbox) busca la herramienta Empaquetar capas (Package Layers).
2. En la ventana de diálogo de Package Layers haz un Click en el botón y selecciona todas las capas. Considera que estas únicamente son capas vectoriales.
3. Guarda con el nombre del archivos Rur_data.gpkg
y da click en Run y Close cuando haya terminado el proceso.
4. Podemos agregar las capas ráster desde el panel de navegación (Browser). Sólo arrastra las capas ráster (e.g. dem_catchment) a Rur_data.gpkg
. Requieres refrescar el panel de navegación para ver el archivo GeoPackage.
Finalmente, también podemos guardar el proyecto dentro el archivo GeoPackage. De esta forma los datos y configuraciones son almacenados en un archivo que puede compartirse con otros, en vez de archivos shapefiles separados y estilizado de archivos.
5. Desde el menú principal, elige Project | Save to | GeoPackage...
6. Escoge la conexión Rur_data.gpkg
y guarda el proyecto con un nombre ej. e.g. Rur
12. Estilizado de la Cuenca resultante
Mostrar los resultados de su análisis, puedes usar una técnica nombrada Inverted Polygon Shapeburst Fills que enfoca a la Cuenca Rur.
1. Abre el panel (Layer Styling panel) con un Click sobre el botón . Indica la capa objetivo (target) a Rur_Catchment.
2. Cambia la representación de Single symbol a Inverted polygons. Esta renderea los datos como una geometría inversa, ej. estiliza todos los objetos fuera del polígono. Se crea una mascara alrededor de el Valle Rur.
3. Después selecciona el componente Simple Fill . Así cambia el tipo de símbolo (Symbol layer type) a Shapeburst fill. En la sección de Gradient colors usa el método por default Two color . Cambia el primer color a un valaor de RGB de 135 | 135 | 135. El segundo color cambialo a blanco.
4. En la sección de Shading style, con un click en la opción Set distance y coloca la distancia a 4 mm.
5. Bajo la sección de Layer rendering indica la Opacity a 65%.
6. Finalmente en la parte superior del panel de estlizado Layer Styling, selecciona el componente Shapeburst fill.
7. Haz un Click en el botón Add symbol layer . Cambia la representación Simple fill a Symbol layer type of Outline simple line. Proporciona en color en negro y Stroke width de 0.46 mm.
Esto nos proporcionará un estilo agradable de la Cuenca como se muestra en la figura inferior.
13. Conclusiones
En esta lección haz aprendido a:
- encontrar y descargar SRTM DEM tiles desde el sitio web de EarthExplorer
- hacer un mosaico completo de capas ráster en una capa ráster virtual
- reproyectar ráster
- crear subconjuntos de capas ráster
- llenar huecos en un DEM
- calcular y estilizar capas ráster flow direction
- calcular orden Strahler de corrientes
- delinear corrientes
- delinear cuencas
- estilizar capas para la visualización de Cuencas