UTILIZACION SAGA GIS

1. LOS ARCHIVOS DEBEN ESTAR EN UTM

2. IMPORTAMOS DESDE GRID EL MDT DESCARGADO – Y EL SHAPE DEL PUNTO DE CONTROL

3. CORREGIMOS LAS DEPRESIONES

4. OBTENEMOS FLUJO ACUMULADO

TERRAIN ANALISIS  HYDROLOGY  FLOW ACUMULATION

5. DETERMINAMOS LA CUENCA A UN PUNTO SELECCIONADO

TERRAIN ANALISIS  HYDROLOGY  FLOW UPSLOPE AREA INTERACTIVE

6. CONVERTIR EN POLIGONO

SHAPE GRID TOOLS  VECTORISING GRID CLASSES  BORRAMOS EL AREA QUE NO NOS
INTERESA (FUERA DE CUENCA)

7. OBTENER ARCHIVO RASTER SOLO DE LA CUENCA

SHAPE GRID TOOLS  CLIP GRID WITH POLYGON  BORRAMOS EL AREA QUE NO NOS INTERESA
(FUERA DE CUENCA).
8. GENERAR LA RED HIDRICA DESDE TERRAIN ANALYSIS

9. DETERMINAR LAS SUBCUENCAS DESDE

1. Iniciamos Google Earth y nos posicionamos en la zona de trabajo. Seleccionamos el punto
de control que hace referencia a la ubicación de nuestro proyecto. Luego seleccionamos el
área aproximada donde podría estar la cuenca de aporte. Guardamos cada uno en una
carpeta especial.
2. Iniciamos Global Mapper y abrimos los archivos guardados anteriormente.
Vamos a File-Download online imagenery- Worldwide data- SRTM v3.
Luego vamos a Tool-Projection- Selecciono coordenadas UTM (Esto nos permite pasar las
coordenadas geográficas de longitud y latitud obtenidas en el google earth, a coordenadas
métricas).
Para curvas de nivel: Análisis- Generar contorno- Selecciono el intervalo de las curvas de
nivel de s.
Para generar cursos de ríos: Análisis- Generar Watershed.
Guardados:
Guardamos el área de trabajo (File- Export- Export elevation grid format- GEOTIFF).
Guardamos el punto de Control ((File- Export- Export vector/ lisar data- SHAPEFILE).
3. Iniciamos Saga Gis.
 Abrimos los archivos que guardamos anteriormente. Para el área de trabajo
(GEOTIFF), vamos a: File- Grid- All files- Abrimos el área de la cuenca guardado en
formato geotiff- Data: para visualizar el área. Para el punto de control (SHAPEFILE)
vamos a: File- Shape- Abro punto de control. Después abriremos el punto de
control del el archivo (no sinks) para visualizarse.
 Corrección de Depresiones (ver Nota 1): Tool- Terrain Analysis- Preprocessing- Fill
Sink (Rellenar sumideros o depresiones) y elegimos el algoritmo de Wang and Liu,
que es el más recomendado para utilizar.
Se ingresa:
Grid system- Selecciono el archivo con las coordenadas
>>DEM (Modelo de elevación digital)- Selecciono el área que debo corregir.
Se crea un archivo con (No sinks) y también Flow Directions (direcciones de flujo) y
Watershed Busins (Cuenca hidrográfica). Pero estos 2 últimos no son de nuestro
interés, por lo tanto lo eliminamos (no se borra el archivo, solo lo estamos
sacando del panel de visualización).
 Obtención de flujo acumulado: Vamos a Tools- Terrain analysis- Hydrology- Flow
Acumulation Flujo acumulado (Flow Tracing que significa trazado de flujo).
Se ingresa:
Grid system- Selecciono el archivo con las coordenadas
>>DEM (Modelo de elevación digital)- Selecciono el área ya corregida (No sinks).
En data se crea un archivo que dice Flow Acumulation.
Ver Nota 2
También: El método de trazado de escurrimiento es el algoritmo de ruteo de onda
cinemática, es decir la ecuación de Manning que utiliza únicamente la variable
dependiente, pendiente y rugosidad.
 Abrimos el punto de control.
 Delimitación de la cuenca con respecto al punto de control: Vamos a Tools-
Terrain Analysis- Hidrology- Flow uslope área [interactive] que sería área de
deslizamiento de flujo interactivo.
Se ingresa:
Grid system (sistema de red)- Selecciono el archivo con las coordenadas
>>DEM (Modelo de elevación digital)- Selecciono el área ya corregida (No sinks).
Options- Convergence- Determinisic 8- Ok (Esto es un algoritmo que el programa
aplica internamente para hacer millones de cálculos para poder definir por donde
va el escurrimiento, y éste es el más recomendado para este tipo de hidrología).
 Luego, vamos a Data, selecciono la flecha de la parte superior y hago click sobre la
línea azul cercana al punto de control de tal forma de la cuenca contenga dicho
punto. Esto se sigue ejecutando, por lo tanto para terminar con la ejecución voy a
Tool- Terrain Analysis- Hidrology- Flow uslope área [interactive] y hago click de
vuelta.
Esto permite definir que el 100% de la cuenca delimitada está aportando al punto
seleccionado que representa el punto de control en estudio o el emplazamiento
de la obra propiamente dicha.
Ya tenemos la cuenca total determinada, la de estudio que nos importa.
 Convertir en polígono la cuenca: Vamos a Tools- Shapes- Grid Tools- Vectoring grid
classes (verctorizar clases de grilla), lo que estamos haciendo es convertir las
grillas en polígonos
Se ingresa:
Grid system- Selecciono el archivo con las coordenadas
>>Grid- Usplope área
Polygons<<créate>>- Ok
Data- polygon- 01. Upslope área- New- Selecciono la flecha de arriba y hago click
en el polígono. Esto permite mostrarme que hay 2 polígonos y se representan con
diferentes colores. El que importa es el interior (cuenca), por lo tanto procedo a
eliminar el otro con click derecho y pongo Delete Selection.
Con esto eliminamos todos aquellos puntos que no aportan al punto de control de
nuestro interés, es decir, se eliminan aquellos que aportan 0%.
 Obtener el archivo raster de la cuenca: Con esto lo que hacemos es obtener un
MDT (modelo digital del terreno) con la forma del polígono que creamos
anteriormente, es decir, el polígono permite cortar toda aquella zona por fuera del
mismo que no es de nuestro interés, ya que no pertenece a la cuenca. Vamos a
Tools- Shapes- Grid Tools- Clip grid with polygon (Copiar la grilla con un polígono).
Grid system- Selecciono el archivo con las coordenadas
>>Imput- Selecciono el archivo de la cuenca corregido (No sinks)
>>Polygons- 01. Uslope área- OK
Si voy a data puedo visualizar que se creó otro archivo (No sinks) pero
definitivamente de la cuenca de interés.
Aquí borro todos los demás archivos, menos el nuevo (No sinks), que sería el
polígono de la cuenca ya delimitado.
 Generar la red hídrica: Vamos a Tools- Terrain analysis- Chanels- Chanel network
and drainage basins (canales y áreas de drenaje)
>>Grid sistem- Selecciono las coordenadas
>>DEM- MDT (no sinks) es el corregido
>>Elevation- Cuenca (no sinks)
Options (el número de orden de los canales, por lo general pueden ir de 3 a 5,
pero seleccionamos el orden 3, ya que es más conveniente).
 Guardamos:
Punto de control en formato shape
El archivo MDT (No sinks) en formato grid
La red de escurrimiento en formato shape
El área de la cuenca (Upslope área).
 Creación de Subcuencas (Opcional): Vamos a Tools- Geoprocessing- Terrain
analysis- Chanels- Watershape Basins (Extended)
Se ingresa:
>>Grid sistem- Selecciono las coordenadas
>>DEM- MDT (no sinks) es el corregido
>>Drainage Network- MDT (no sinks) es el corregido
Estas Subcuencas las podemos trabajar por separado dentro de la cuenca mayor,
para conocer los aportes de cada una. Haríamos como una red de aporte que
luego se cargaran en QGIS para ser estudiados. Debo definir un punto de control
para cada subcuenca.
Las Subcuencas se dividen de acuerdo al número de orden que dimos
anteriormente.
 Iniciamos QGIS. Abrimos los siguientes archivos extraídos del Saga Gis, que son los
siguientes:
Punto de control
MDT (No sinks)
Red de escurrimiento
Área de la cuenca
 Procedimientos para obtener los parámetros geomorfológicos de la cuenca:

 Área: En la solapa de capas, seleccionamos el área de la cuenca y luego en

la barra de tareas que se encuentra en la parte superior, abrimos tabla de
atributos (F6). Editar (Ctrl E)- Calculadora de campos (Ctrl M)- Creamos un
campo nuevo y le asignamos un nombre al campo de Salida “Área [km2]”.-
Asignamos un tipo de campo de salida- Nro. decimal (real) con una
precisión de 2.- Busco Geometría y selecciono $área- Luego tengo que
modificar esta expresión por conversión de unidades, ya que este
programa trabaja en metros y la información que necesito es en km2. Por
lo tanto coloco”$área/1000000”. Aceptar!
Una vez realizado esto, debo cerrar el modo de edición para que se
guarden los cambios realizado.
  Perímetro: Se procede análogamente al inciso anterior, pero con la
diferencia que en este caso, la conversión de unidades es diferente, de
modo tal que resulta, de dividir el perímetro en 1000!
 Ancho medio: Se define como la relación entre el área y la longitud de la
cuenca. Esta longitud es la que se toma a partir, del aforo o punto de
control aguas abajo y el punto más alejado aguas arriba. La longitud se
debe medir manualmente con un comando desde la barra de tareas
(medir línea)
 Pendiente y relieve: Para obtener estos datos, previamente se deben crear
los mapas de relieve y pendiente para que estos sean analizados. Para ello
vamos a barra de tareas- Raster- Análisis de terreno- Pendiente- Capa de
entrada: MDT (No sinks) y en capa de salida: Se asigna nombre y luego
guardamos. Análogamente se realiza los mismos pasos para el mapa de
relieve.
En caso de no encontrar este comando en Raster, abro la caja de
herramientas (Ctrl+Alt+T) y procedo como se mencionó anteriormente.
Una vez realizado este paso previo, lo siguiente es obtener los datos de
pendiente y relieve. Para esto vamos a Caja de Herramientas- Análisis
raster- Estadística de la zona. En capa raster cargamos los mapas que
fueron creados. En el caso de pendiente se coloca en “estadística a
calcular” se colocan la información que requerimos que por lo general
resulta ser Pendiente Media, Mínima y Máxima. Por otro lado, para el
relieve solo es de nuestro interés el Relieve Medio!.

PARAMETROS DE LA RED DE DRENAJE:

 PENDIENTE MEDIA DEL CURSO PRINCIPAL SIGUIENDO LOS CRITERIOS DE TAYLOR

 NUMERO DE ORDEN
 DENSIDAD DE DRENAJE
 PATRON DE ALINEACION

PARAMETRO GEOMORFOLOGICOS DE LA CUENCA:

 AREA
 PERIMETRO
 ANCHO MEDIO
 RELIEVE
 FORMA: INDICE DE COMPACIDAD O DE GLAVELIUS – FACTOR DE FORMA DE HORTON
 PENDIENTE MEDIA: CALCULAR SIGUIENDO LOS CRITERIOS ALBOR Y DE HORTON.
 CURVA HIPSOMETRICA: COMPARAR CON LA CURVA TIPO Y DEFINIR SI SE TRATA DE UN
RIO MEDIANO O VIEJO
 ALTURA MEDIA DE LA CUENCA
  ALTURA MEDIANA DE LA CUENCA

ANALIZAR LOS ALGORITMOS!!

SABER COMO SE DELIMITAN SUBCUENCAS DENTRO DE NUESTRA CUENCA DE IMPORTANCIA

GLOSARIO:

Grid: Un Esri grid es un formato GIS ráster de almacenamiento de datos raster nativo de ESRI. ...
Utilizamos grids de tipo entero para representar datos discretos y grids de punto flotante para
representar datos continuos. Los datos de elevación son un ejemplo de un grid de punto flotante.

Shape: El formato ESRI Shapefile (SHP) es un formato de archivo informático propietario de datos

espaciales desarrollado por la compañía ESRI, quien crea y comercializa software para Sistemas de
Información Geográfica como Arc/Info o ArcGIS. Originalmente se creó para la utilización con su
producto ArcView GIS, pero actualmente se ha convertido en formato estándar de facto para el
intercambio de información geográfica entre Sistemas de Información Geográfica por la
importancia que los productos ESRI tienen en el mercado SIG y por estar muy bien documentado.
Un shapefile es un formato vectorial de almacenamiento digital donde se guarda la localización de
los elementos geográficos y los atributos asociados a ellos. No obstante carece de capacidad para
almacenar información topológica. Es un formato multiarchivo, es decir está generado por varios
ficheros informáticos. El número mínimo requerido es de tres y tienen las extensiones siguientes:

 .shp - es el archivo que almacena las entidades geométricas de los objetos.

 .shx - es el archivo que almacena el índice de las entidades geométricas.
 .dbf - es la base de datos, en formato dBASE, donde se almacena la información de los
atributos de los objetos.
Además de estos tres archivos requeridos, opcionalmente se pueden utilizar otros para mejorar el
funcionamiento en las operaciones de consulta a la base de datos, información sobre
la proyección cartográfica, o almacenamiento de metadatos. Estos archivos son:

 .prj - Es el archivo que guarda la información referida al sistema de coordenadas en

formato WKT.
 .sbn y .sbx - Almacena el índice espacial de las entidades.
 .fbn y .fbx - Almacena el índice espacial de las entidades para los shapefiles que son
inalterables (solo lectura).
 .ain y .aih - Almacena el índice de atributo de los campos activos en una tabla o el tema de
la tabla de atributos.
 .xml - Almacena los metadatos del shapefile.

Raster: En su forma más simple, un ráster consta de una matriz de celdas (o píxeles) organizadas
en filas y columnas (o una cuadrícula) en la que cada celda contiene un valor que representa
información, como la temperatura.

Raster de imagen digital: Una imagen en mapa de bits, imagen ráster (calcos del inglés) o imagen
de píxeles es una estructura o fichero de datos que representa una rejilla rectangular de píxeles o
puntos de color, denominada matriz, que se puede visualizar en un monitor, papel u otro
dispositivo de representación.

Lidar: El LIDAR (de light detection and ranging) es una técnica de teledetección óptica que utiliza la
luz de láser para obtener una muestra densa de la superficie de la tierra produciendo mediciones
exactas de x, y y z. LIDAR, que se utiliza principalmente en aplicaciones de representación
cartográfica láser aéreas, está surgiendo como una alternativa rentable para las técnicas de
topografía tradicionales como una fotogrametría. LIDAR produce datasets de nube de puntos
masivos que se pueden administrar, visualizar, analizar y compartir usando ArcGIS.

Los componentes de hardware principales de un sistema lidar incluyen un vehículo de recolección

(avión, helicóptero, vehículo y trípode), sistema de escáner láser, GPS (Sistema de
posicionamiento global) e INS (sistema de navegación por inercia). Un sistema INS mide la
rotación, inclinación y encabezamiento del sistema lidar.

LIDAR es un sensor óptico activo que transmite rayos láser hacia un objetivo mientras se mueve a
través de rutas de topografía específicas. El reflejo del láser del objetivo lo detectan y analizan los
receptores en el sensor lidar. Estos receptores registran el tiempo preciso desde que el pulso láser
dejó el sistema hasta cuando regresó para calcular la distancia límite entre el sensor y el objetivo.
Combinado con la información posicional (GPS e INS), estas medidas de distancia se transforman
en medidas de puntos tridimensionales reales del objetivo reflector en el espacio del objeto.

Los datos de punto se procesan posteriormente después de que la recopilación de datos lidar se
reconocen dentro de las coordenadas x,y,z georreferenciadas con alta precisión al analizar el
rango de tiempo láser, ángulo de escaneo láser, posición del GPS e información del INS.
Determinisic 8: Módulo de acumulación de flujo (recursivo). Procesamiento ascendente recursivo
de celdas para el cálculo de la acumulación de flujo y parámetros relacionados. Este conjunto de
algoritmos procesa recursivamente todas las celdas conectadas hacia arriba hasta que cada celda
del DEM haya sido procesada.

MDT: El tipo de Modelo Digital del Terreno (MDT) más conocido es el Modelo Digital de
Elevaciones (MDE), un caso particular de aquel, en el que la variable representada es la cota del
terreno en relación a un sistema de referencia concreto. 1 No obstante no hay un uso normalizado
en la literatura científica de los términos Modelo Digital de Elevaciones (MDE), Modelo Digital del
Terreno (MDT) y Modelo Digital de Superficie (MDS). En la mayoría de los casos, el término
Modelo Digital de Superficie se refiere a la superficie de la tierra e incluye todos los objetos que
esta contiene. En cambio un MDT representa la superficie de suelo desnudo y sin ningún objeto,
como la vegetación o los edificios

Nota 1:

Rellenar los vacíos y depresiones

Para que SAGAGIS pueda interpretar y delimitar la cuenca y la red hídrica, el raster debe ser
"hidrológicamente correcto", es decir que todas elevaciones sean continuamente decrecientes y
que no hayan sumideros. Nuestro raster se tiene que corregir para poder delimitar las cuentas,
para lo cual usamos el comando:

Geoproccessing / Terrain Analisys / Preprocessing / Fill Sinks (Wang Liu)

Se selecciona el Grid System que es único disponible y el DEM que es el nombre del archivo como
se ve en la figura.

Nota 2:

La herramienta Acumulación de flujo calcula el flujo acumulado como el peso acumulado de todas

las celdas que fluyen en cada celda de pendiente descendente en el ráster de salida. Si no se
proporciona un ráster de peso, se aplica un peso de 1 a cada celda, y el valor de celdas en el ráster
de salida es el número de celdas que fluye en cada celda.

En el siguiente gráfico siguiente, la imagen superior izquierda muestra la dirección de viaje desde
cada celda y la superior derecha el número de celdas que fluyen hacia cada celda.

Determinar la acumulación de
flujo

Las celdas con una acumulación de flujo alta son áreas de flujo concentrado y pueden ser útiles
para identificar canales de arroyos. Esto se analiza en Identificación de redes de arroyos. Las celdas
con una acumulación de flujo de 0 son alturas topográficas locales y se pueden utilizar para
identificar crestas.

Ejemplo

Uno de los usos de muestra de la herramienta Acumulación de flujo con un ráster de peso de

entrada puede ser determinar cuánta lluvia ha caído en una cuenca hidrográfica determinada. En
ese caso, el ráster de peso puede ser un ráster continuo que representa las precipitaciones
promedio durante una tormenta determinada. La salida desde la herramienta puede representar
la cantidad de lluvia que puede fluir por cada celda, asumiendo que toda la lluvia se convirtió en
escorrentía y que no hubo intercepción, evapotranspiración o pérdida al agua subterránea.
También se puede visualizar como la cantidad de lluvia que cayó en la superficie, a una pendiente
ascendente desde cada celda.

Los resultados de Acumulación de flujo se pueden utilizar para crear una red de arroyos al aplicar
un valor de umbral para seleccionar celdas con un flujo acumulado alto.

Por ejemplo, el procedimiento para crear un ráster en el que el valor 1 represente una red de
arroyos en un fondo de NoData, puede utilizar la siguiente manera:

 Realizar una operación condicional con la herramienta Evaluación condicional con la

siguiente configuración:

o Ráster condicional de entrada: Flowacc

Expresión: Valor > 100

Ráster verdadero o constante de entrada: 1

 Alternativamente, ejecute la herramienta Establecer nulos con las siguientes

configuraciones:

o Ráster condicional de entrada: : Flowacc

Expresión: Valor <= 100

Ráster falso o constante de entrada: : 1

En los dos ejemplos, a todas las celdas que tienen más de 100 celdas que fluyen hacia ellas se les
asignan el valor 1 y a todas las otras celdas se les asigna NoData. Para procesos futuros, es
importante que la red de arroyos, un conjunto de entidades lineales de ráster, se represente como
valores en un fondo de NoData.

Este método de derivar un flujo acumulado desde un DEM se presenta en Jenson and Domingue
(1988). En Tarboton et al, se presenta un método analítico para determinar un valor de umbral
apropiado para definir una red de arroyos. (1991).