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 m². 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.