Delimitar una cuenca hidrográfica en ArcGIS con Hydrology de Spatial Analyst
En este tutorial de ArcGis aprenderás a delimitar una cuenca hidrográfica utilizando la herramienta Hydrology de Spatial Analyst. En primer lugar, debemos contar un modelo digital de elevación del área de estudio, el cual se puede obtener de algún servidor gratuito o interpolar a partir coordenadas XYZ o curvas de nivel con la ayuda de la herramienta 3D Analyst.
Para ello utilizaremos la herramienta Hydrology contenida en Arctoolbox-Spatial Analyst Tools (ver artículo relacionado).
![Arctoolbox de ArcGis Arctoolbox de ArcGis](https://images-blogger-opensocial.googleusercontent.com/gadgets/proxy?url=http%3A%2F%2F1.bp.blogspot.com%2F-66obqgjPZvQ%2FUl9wNQwwulI%2FAAAAAAAADfA%2FxrGDV7JOOUU%2FHydrology-de-arcgis_thumb%2525255B1%2525255D.png%3Fimgmax%3D800&container=blogger&gadget=a&rewriteMime=image%2F*)
Paso 1. Fill Sinks. Con esta herramienta se rellenan las imperfecciones existentes en la superficie del modelo digital de elevaciones, de tal forma que las celdas en depresión alcancen el nivel del terreno de alrededor, con el objetivo de poder determinar de forma adecuada la dirección del flujo. Para ello a partir de Hydrology se da clic en Fill, se abre una ventana donde se debe rellenar la siguiente información.
Input surface raster: se selecciona el dtm que vamos a utilizar para el procesamiento, en este caso es dtm_cuencas.
Output surface raster: aquí seleccionamos la ruta y el nombre del archivo de salida, por defecto le colocará el nombre “Fill_dtm_cue1”.
Z limit (optional): indica la máxima profundidad de los sumideros que queremos rellenar. Las profundidades de sumideros o imperfecciones mayores al valor colocado en este campo no se rellenaran. En caso de dejar el campo en blanco, el programa tomará por defecto rellenar todos los sumideros, independientemente de la profundidad. Para el caso del ejemplo la dejaremos en blanco.
![Rellenando sumideros con arcgis Rellenando sumideros con arcgis](https://4.bp.blogspot.com/-rXdb6y6qAXA/Um1KQGV011I/AAAAAAAADmg/JgsUprEwyQs/s1600/rellenando-sumideros-con-arcgis.png)
Como resultado se obtiene el raster denominado fill_dtm_cue1
Paso 2. Flow direction. Se define aquí la dirección del flujo buscando el camino descendente de una celda a otra. A partir de Hydrology se da clic en Flow direction, se abre una ventana donde se debe rellenar la siguiente información.
Input surface raster: se selecciona el raster creado en el paso anterior que se denomina Fill_dtm_cue1
Output surface raster: aquí seleccionamos la ruta y el nombre del archivo de salida, por defecto le colocará el nombre “FlowDir_fill1”.
Output drop raster (optional): Es una salida opcional. El drop raster muestra la relación entre el cambio máximo en la elevación de cada celda a lo largo de la dirección del flujo, expresada en porcentajes. No escribimos nada en este campo.
![Ventana direccion de flujo con arcgis Ventana direccion de flujo con arcgis](https://3.bp.blogspot.com/-EbJTyzpTehQ/Um1KcgW5BsI/AAAAAAAADm8/F2tPKwheBBA/s1600/flow-direction-con-arcgis.png)
Como resultado se obtiene el raster denominado FlowDir_fill1 (dirección de flujo)
![Mapa de direccion flujo con arcgis Mapa de direccion flujo](https://4.bp.blogspot.com/-uL4W3LG2S3A/Um1Kdfk29tI/AAAAAAAADns/GNyHCPkGmjU/s320/mapa-flow-direction-con-arcgis.png)
Paso 3. Flow accumulation. Crea el raster de acumulación de flujo en cada celda. Se determina el número de celdas de aguas arriba que vierten sobre cada una de las celdas inmediatamente aguas abajo de ella. A partir de Hydrology se da clic en Flow accumulation, se abre una ventana donde se debe rellenar la siguiente información.
Input direction raster: se selecciona el raster creado en el paso anterior que se denomina FlowDir_fill1
Output accumulation raster: aquí seleccionamos la ruta y el nombre del archivo de salida, por defecto le colocará el nombre “FlowAcc_flow1”.
Input weight raster (optional): Es una salida opcional.
Output data type raster (optional): Es una salida opcional. Por defecto dejamos FLOAT
![Ventana acumulacion de flujo Ventana acumulacion de flujo](https://4.bp.blogspot.com/-DBZb5ZR5Bic/Um1KfnotkSI/AAAAAAAADoU/pffL0Q5cEjM/s1600/ventana-flow-accumulation.png)
Como resultado se obtiene el raster denominado FlowAcc_flow1
![Mapa de acumulacion de flujo con arcgis Mapa de acumulacion de flujo con arcgis](https://4.bp.blogspot.com/-1fc3MaB7ITs/Um1KcXJhdGI/AAAAAAAADnQ/e892ceM4A2k/s1600/flow-accumulation-con-arcgis.png)
Paso 4. Stream definition. En esta fase se clasifican las celdas con acumulación de flujo superior a un umbral especificado por el usuario como celdas pertenecientes a la red de flujo. El umbral debe ser especificado como el número de celdas vertientes a la que se está clasificando en cada momento. Aquí se debe entrar a sopesar que valor sería el más indicado, ya que si el valor de acumulación es muy bajo muchos pixeles serán seleccionados como pertenecientes a la red hídrica, si por lo contrario, el valor del pixel es muy alto solo aquellos drenajes de orden alto serían definidos como red hídrica. En otras palabras, seleccionar un valor bajo del umbral significa que obtendremos afluentes pequeños en nuestra red de drenajes, en cambio un valor alto, modela los drenajes de mayor tamaño.
Para crear una red de corriente a partir del raster de acumulación de flujo, se procede de la siguiente forma.
a. En el menú de la barra de herramienta desplegable de Spatial Analyst selecciona la Calculadora Raster.
b. Con el fin de crear una red de drenajes, tendrá que especificar un umbral para la cantidad de píxeles adyacentes que constituyen una corriente. Aquí vamos a especificar un valor límite de 20000 píxeles de acumulación (si hay más de 20000 píxeles que desembocan en él se parte de la red de corriente).
c. En la calculadora raster escribimos la expresión siguiente:
stream_Red = con([FlowAcc_flow1] > 20000, 1)
![Calculadora raster Calculadora raster](https://3.bp.blogspot.com/-1837HtMeJ5k/Um1KbyU6kRI/AAAAAAAADm0/0SvqParzcjc/s1600/calculadora-raster-de-arcgis.png)
Clic en Evaluate y como resultado se obtiene el raster stream_Red(red de flujo) que se muestra a continuación.
![Red de rios con arcgis Red de rios con arcgis](https://4.bp.blogspot.com/-KAK2lhWLR_A/Um1KdBEseBI/AAAAAAAADnc/KXmHKpWEnxQ/s1600/flujo-raster-con-arcgis.png)
Paso 5. Stream Link. Divide el cauce en segmentos no interrumpidos. Es decir, que dichas secciones en las que se divide el recorrido del flujo serán segmentos que conectan dos uniones sucesivas, una unión y un punto de desague o una unión y una división del área de drenaje. Para ello se procede de la siguiente forma.
Clic en Hydrology, luego en Stream Link, en la ventana que aparece se debe rellenar los siguientes campos.
Input stream raster: se selecciona el raster streamRed creado en el paso anterior.
Input flow direction raster: Se selecciona el raster de dirección de flujo, denominado flowdir_fill1 creado en el paso 2.
Output raster: aquí seleccionamos la ruta y el nombre del archivo de salida, le colocaremos el nombre “StreamLink1”.
![Segmentación de ríos con Arcgis Segmentación de ríos con Arcgis](https://2.bp.blogspot.com/-Ia0rvQZlFMo/Um1KfDvUhwI/AAAAAAAADoE/LDsfXcuKxms/s1600/ventan-stream-link.png)
Después de dar clic en OK, se obtiene el siguiente mapa.
![Ríos segmentados Ríos segmentados](https://1.bp.blogspot.com/-lIbNmyRaJyc/Um1KecEO3VI/AAAAAAAADog/6jI2x4F5gdQ/s1600/stream-link-con-arcgis.png)
Paso 6. Stream Order. Crea un raster del orden de las corrientes. Para ello se usan dos métodos:
En el método Strahler, el orden de la corriente se incrementa cuando se cruzando dos drenajes del mismo orden. Dos drenajes de diferentes ordenes no se traducirá en un aumento del orden de la siguiente corriente.
En el método Shreve los órdenes de corrientes son aditivos.
El procedimiento es el siguiente: Clic en Hydrology, luego en Stream Order, en la ventana que aparece se debe rellenar los siguientes campos.
Input stream raster: se selecciona el raster stream_red creado el paso 4.
Input flow direction raster: Se selecciona el raster de dirección de flujo, denominado flowdir_fill1 creado en el paso 2.
Output raster: aquí seleccionamos la ruta y el nombre del archivo de salida, le colocaremos el nombre “Stream_order”.
Method of stream ordering (optional). Se selecciona el método STRAHLER.
![Orden de las corrientes con arcgis Orden de las corrientes con arcgis](https://4.bp.blogspot.com/-gpgfBJDNg5Y/Um1KgLSqUGI/AAAAAAAADos/fwrjc7Un3a0/s1600/ventana-stream-order.png)
El resultado es el siguiente.
![corrientes con orden corrientes con orden](https://2.bp.blogspot.com/-MZc3BhWjjY8/Um1KedN_nBI/AAAAAAAADn8/_9Kvsnw3Lv8/s1600/stream-order-con-arcgis.png)
Paso 7. Stream Feature. Crea un shape de drenajes. El procedimiento es el siguiente: Clic en Hydrology, luego en Stream Feature, en la ventana que aparece se debe rellenar los siguientes campos.
Input stream raster: se selecciona el raster stream_red creado el paso 4.
Input flow direction raster: Se selecciona el raster de dirección de flujo, denominado flowdir_fill1 creado en el paso 2.
Output polyline feature: aquí seleccionamos la ruta y el nombre del archivo de salida, le colocaremos el nombre “Stream20000”.
![Ventana-para-convertir-stream Ventana-para-convertir-stream](https://3.bp.blogspot.com/-EJEyhZSbJAM/Um1KgbNZaXI/AAAAAAAADoo/8VKK-Y4kOK8/s1600/ventana-stream-to-feature.png)
Clic en OK, Obtenemos o siguiente.
![](https://2.bp.blogspot.com/-7XQXFw72N1A/Um1Kd5BZARI/AAAAAAAADno/fgjv030ZpJQ/s1600/shape-de-rio-con-arcgis.png)
Paso 8. Feature Vertice To Point. Esta herramienta permite determinar los puntos donde se cortan cada uno de los drenajes, es decir convierte los vértices a punto. Podemos determinar un punto al inicio, la mitad o al final de cada tramo de corriente, para este caso nos interesan los puntos finales que es donde hay acumulación de flujo y es el punto importante para determinación de las cuencas. Para ello seguimos los siguientes pasos.
Vamos a ArcToolsbox, luego Data Management tools, seguido de Feature y finalmente Feature Vertice to Point.
![Feature vertices to point Feature vertices to point](https://1.bp.blogspot.com/-KBxL_gYlgvc/Um1KcMCE1aI/AAAAAAAADnA/wkBwz5fM5Kc/s1600/convertir-vertices-a-puntos-con-arcgis.png)
Aparece una ventana donde se debe rellenar la siguiente información.
Input Feature: introducimos el shape de la red de drenaje creada en el paso 7 (Stream20000).
Output Feature Class: aquí seleccionamos la ruta y el nombre del archivo de salida, le colocaremos el nombre “vertToPoint.shape”.
Point Type (optional): Seleccionamos la opción END que agregará los puntos al final de cada tramo de corriente.
![convertir vértices a punto con arcgis convertir vértices a punto con arcgis](https://2.bp.blogspot.com/--RGcGqtK78k/Um1KfaKU0YI/AAAAAAAADoc/ZxtGjWqOP-A/s1600/ventana-convertir-vertices-a-puntos.png)
Como resultado obtenemos lo siguiente.
![Vértices convertidos en puntos Vértices convertidos en puntos](https://2.bp.blogspot.com/-SLKPDWG4NU4/Um1KdzWkTlI/AAAAAAAADnk/9Wc-_9vroqc/s1600/puntos-de-cierre-microcuencas.png)
Paso 9. Watershed Delineation. Delinea una subcuenca por cada uno de los segmentos de cauce definidos en el paso anterior. Se procede de la siguiente forma: Clic en Hydrology, luego en Watershed Delineation, en la ventana que aparece se debe rellenar los siguientes campos.
Input flow direction raster: se intruduce el raster de dirección de flujo creado en el paso 2, denominado FlowDir_fill1.
Input raster or feature pour point data: se introduce el shape de puntos creado en el paso anterior “vertToPoint”.
Pour point fiel (optional): Dejamos la opción por defecto (ARCID).
Output raster: La ruta y el nombre del archivo salida, le damos el nombre Watershed…
![Delimitación de cuencas con arcgis Delimitación de cuencas con arcgis](https://3.bp.blogspot.com/-BaQIQEy0ks8/Um1Ke1kJvOI/AAAAAAAADoA/99VJKTS5ns4/s1600/delimitacion-de-cuencas-con-arcgis.png)
Finalmente obtenemos la delimitación de las cuencas. Le cambiamos la simbología y se obtiene lo siguiente. En total salieron 341 cuencas, si no queremos tantas o queremos modelar cuencas principales o más grandes, en el paso 5 colocamos un número mayor de celdas.
![Cuencas delimitadas con arcgis Cuencas delimitadas con arcgis](https://4.bp.blogspot.com/-2-7g30J1cYw/Um1KcDTwL6I/AAAAAAAADm4/bvLPOFYN4jE/s1600/cuencas-delimitadas-con-arcgis.png)
La creación de la red de drenaje y delimitación de cuencas, también se puede implementar a través del Modelbuider que ha explicado con anterioridad, mas adelante lo veremos.