Delimitar una cuenca hidrológica con GRASS en QGIS

En el artículo anterior, escribimos sobre cómo crear una localización y conjunto de mapas para GRASS en QGIS. El ejemplo que se tomó fue con el objetivo de delimitar una cuenca hidrólógica en esa área. Allí se consideró el uso de QGIS en entorno Linux pero, para variar, se va a usar Windows XP (se usó inicialmente Seven pero, a pesar de que las herramientas de GRASS funcionan adecuadamente, los DEM no pueden visualizarse por ensanchamiento de contraste y “desaparecen” de la vista). En el caso de XP, la localización se creará en el directorio C:\Document and Settings\nombre_usuario\grassdata y llevará el mismo nombre: tiznados. El DEM importado fue el mismo del artículo anterior y se renombró para el complemento de GRASS como N09W068_dem_canoaUTM19. Tiene una superficie de algo más de un millón de hectáreas (1.225.469,61586 Ha); obtenida de la tabla atributiva del polígono de extensión de la capa (Vectorial -> Herramientas de investigación -> Polígono de extensión de la capa).

La delimitación se hará con base en el DEM y usando la herramienta de GRASS r.watershed (Complementos -> GRASS -> Abrir herramientas de GRASS). Cuando se abre la ventana de diálogo de la herramienta aparece en uno de los campos el “tamaño mínimo para cada cuenca (número de celdas)”. El ráster es de 30×30 por lo que una celda ocupa una superficie de 900 m2. Para considerar tamaños mínimos de 10.000, 100.000 y 1.000.000 Has para cada cuenca entonces los números de celdas serán aproximadamente de 111111, 1111111 y 11111111; respectivamente.

Cuando se ejecuta la herramienta (cuidando, en cada caso, de clickar también en el botón que circunscribe el primer cuadrado rojo de la imagen anterior), el resultado es “vacío” para una cuenca de un millón de hectáreas (lo cual era previsible porque el DEM refleja una superficie similar) y produce lo siguiente:

para las superfices de 10.000 y 100.000 Has de tamaño mínimo para cada cuenca; respectivamente. En cada situación se han sobrepuesto los mismos shapefiles que representan los ríos y cuerpos de agua que se encuentran en la región. En el segundo caso se observa que la zona de color azúl oscuro central es la que contiene el río Tiznados; aunque este es un afluente del Río Portuguesa y se une con él aproximadamente a unos 87 Km en línea recta, hacia el SW, desde el borde inferior de la cuenca que lo incluye en mayor extensión.

Para vectorizar el ráster como un shapefile tipo polígono se utiliza la herramienta r.to.vect.area. El resultado es el siguiente:

después de colorizar las áreas por un categorizado del campo value (generado por GRASS). Se puede observar que es posible pasar de partes sencillas a multiparte para eliminar los registros repetidos y luego agregar columnas de geometría para estimar el área y el perímetro. Esto lo hice en el ámbito de QGIS; no de GRASS.

En el próximo artículo se ejecutarán y analizarán las otras opciones de r.watershed.

Esta entrada fue publicada en GRASS, QGIS, SIG, Software Libre. Guarda el enlace permanente.

18 respuestas a Delimitar una cuenca hidrológica con GRASS en QGIS

  1. Pingback: Delimitar una cuenca hidrológica con GRASS en QGIS | Geoprocessing - Geoprocessamento | Scoop.it

  2. Pingback: Delimitar una cuenca hidrológica con GRASS en QGIS | geoprocessamento.org

  3. Pingback: Generar escurrimiento superficial (“ríos”) con GRASS en QGIS |

  4. Matías dijo:

    Hola José,has usado el r.water.outlet? porque estoy tratando de generar una cuenca para un punto de cierre determinado a partir del modelo digital de terreno, pero no me funciona.Le ingreso el archivo de drenaje generado con el r.watershed, pero no me genera nada coherente, solamente unos puntos negros cerca del pto de cierre.
    Desde ya muchas gracias
    Matías desde Montevideo,Uruguay

  5. yam dijo:

    que tal josé, de que manera puedo delimitarla a partir de un punto (boquilla de la cuenca) como lo realizaría con ArcGIS, por favor si pudieras ilustrarme…

  6. Estimado, eso que mencionas sobre el “ensanchamiento” del DEM y que “desaparecen” te refieres a la visualización? porque yo estoy usando 7 con home basic y 32 bits, logro ver el dem adecuadamente si utilizo opciones/estilo/mejora de contraste, usando “estirar a min/max”. No se si es válido mi comentario pero ahi está a ver a quien le sirve.

  7. Aqui el mensaje de error es:
    r.watershed elevation=ASTGM2_UTM_CP3@CRIVAS threshold=111111 basin=Cuenca10k
    SECCI�N 1a (de 5): iniciando memoria.
    El mapa raster se encuentra en una zona diferente (16) a la de la regi�n actual (0).
    No se puede abrir el mapa raster .
    no se puede abrir capa de mapa de elevaci�n
    Ha fallado el subproceso con c�digo de salida 1
    La informaci�n de categor�a para [Cuenca10k] en [CRIVAS] no se encuentra o no es v�lida.
    Finalizado con error

    • Este:

      “El mapa raster se encuentra en una zona diferente (16) a la de la regi�n actual (0).
      No se puede abrir el mapa raster .”

      es el mismo error que en tu primer post (donde aseguraste que ya habías resuelto el problema). Asegúrate ahora de tener en los settings de QGIS la misma proyección que el ráster.

  8. Estimado, es al aplicar el r.watershed que me da el error:

    r.watershed elevation=ASTGM2_UTM_CP3@CRIVAS threshold=111111 basin=Cuenca1
    SECCI�N 1a (de 5): iniciando memoria.
    El mapa raster se encuentra en una zona diferente (16) a la de la regi�n actual (0).
    No se puede abrir el mapa raster .
    no se puede abrir capa de mapa de elevaci�n
    Ha fallado el subproceso con c�digo de salida 1
    La informaci�n de categor�a para [Cuenca1] en [CRIVAS] no se encuentra o no es v�lida.
    Finalizado con error

    Aqui puse la captura de pantalla del error, como ves todo el DEM esta correcto, el area roja indica el área de trabajo de Grass (habia puesto Nic en el paso que indicaste), pero aun asi como ves no hay manera de que me funcione. Que podrá pasar?

  9. Lo logré!, aqui esta el resultado:
    r.watershed elevation=ASTGM2_UTM_CP3@CRIVAS threshold=111111 basin=W1
    SECCI�N 1a (de 5): iniciando memoria.
    SECCI�N 1b (de 5): Determinando flujo Offmap.
    SECCI�N 2: A * Buscar.
    GRASS_INFO_PERCENT: -61
    SECCI�N 3: Acumulando Flujo de Superficie con SFD.
    GRASS_INFO_PERCENT: -61
    SECCI�N 4: determinaci�n de cuenca hidrogr�fica.
    SECCI�N 5: cerrando mapas.
    Finalizado correctamente

    Está interesante el resultado, debo empezar ahora a analizar según los datos que tengo de hojas a escala 50,000

    Saludos!
    y para los que quieran información relacionada los amigos de USGSS y WWF han hecho esto:
    http://worldwildlife.org/pages/hydrosheds

  10. Pingback: Delimitar una cuenca hidrológica con GRA...

  11. Gracias José por tu tremendo aporte en mi trabajo, estoy comparando el uso de software libres con otro que no lo es. Aprendiendo sobre como exportar e importar. gracias nuevamente

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s