Unidad5 Pract TelAmb
Unidad5 Pract TelAmb
Unidad5 Pract TelAmb
Índice
1. Introducción 1
1.3. Actividades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.4. Requerimientos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
2. Anexos 3
1. Introducción
El uso de un software que, mediante el análisis de diferencias de altura permite la extracción de los límites
de una cuenca y el diagrama de la redes de drenaje, es un salto hacia hacia la automatización para un
fotointerprete.
¿Por qué es posible diagramar con un software una cuenca hidrográfica, tomando cómo datos de entrada
diferencias de altura, propias de un Modelo de Elevación Digital (MDE)?
3. Implementar complementos GRASS gis (GRASS) para delimitar una cuenca hidrográfica a partir de
un MDE.
1.3. Actividades
1.4. Requerimientos
QGIS
Datos oficiales de: Instituto Geográfico Nacional (IGN), Infraestructura de Datos Espaciales de la
Provincia de Córdoba (IDECOR), Shuttle Radar Topography Mission (srtm) y adquisiciones LANDSAT
8 OLI.
Para ser evaluado, envíanos el mapa con todas las actividades plasmadas, de la cuenca del dique La
Quebrada (Córdoba, Argentina), el escenario es a criterio del autor respetando las consignas, con la
memoria respaldo (la misma no debe superar una carilla).
Cabe destacar que la memoria corresponde al escrito que justifica la elaboración del mapa, las técnicas
aplicadas y el problema ambiental que estaríamos resolviendo. La memoria, es igual de importante que el
mapa en sí, ya que ambos productos se complementan.
Los anexos siguientes constituyen la guía para generar paso a paso todas las actividades propuestas.
¡¡Que se disfrute el proceso!! estamos a disposición.
2. Anexos
Pasos y resultados en el programa QGIS (instalar la versión más actualizada a la fecha), para la delimitación
de la cuenca, utilizando el material preparado para la unidad.
Lo vimos en la primer práctica, va el repaso. QGIS por defecto cuando se abre trae una tipología de
coordenadas, ajustar las adecuadas para efectuar el práctico (LAT/LONG y POSGAR son las más comunes).
Es necesario que las coordenadas de todas las capas se presenten en la misma georreferenciación, para
controlar y cambiar se recurre, en el programa QGIS, en un proyecto EN BLANCO las coordenadas de
trabajo de interés y ajustar cada capa que se agregue.
Primero SIEMPRE, cargar datos raster, para fijar las coordenadas del proyecto, revisar el formato y la
compatibilidad, ésta tarea permite que los datos que estamos manejando se ubiquen específicamente en
la región del globo terráqueo donde estamos trabajando. Un tip clave es revisar el header de los datos
raster y saber en que tipología de coordenadas vienen.
Importante, cargar primero los datos tipo raster y posteriormente los datos vectoriales (ver Figura 1).
Revisar que todos los datos necesarios se encuentren en el panel de tareas del programa QGIS. Detalle no
menor, es que posteriormente a la carga de datos se ubiquen por encima los datos vectoriales en el panel
de capas, caso contrario, no se visualizan.
¿Que es GRASS?
GRASS (acrónimo inglés de Geographic Resources Analysis Support System) es un software para tratar
datos provenientes de: Sistema de Información Geográfica (SIG) bajo licencia GPL (software libre);
básicamente soporta y procesa información tanto raster como vectorial.
Para crear un directorio de mapas en el complemento GRASS de QGIS, elegir una denominación adecuada
para la “localización de mapas” y finalmente otra para el directorio de trabajo. La secuencia sería:
La Figura 2, complementa los pasos para crear un directorio en GRASS, dicha función es necesaria, ya
que la propuesta es generar a partir de un MDE (el cual posee básicamente posee variaciones de altura)
una cuenca. Eso se logra con una sola función, que es el procesamiento.
(a) (b)
(c) (d)
(e)
Figura 2: Crear directorio en GRASS a) Acceso y b) Nombre del directorio c) Selección de coordenadas y d) Crear nueva
directorio de mapas e) Directorio creado
Propiamente en el MDE, previo a procesar, es necesario analizar el histograma de los datos. Ésta operación
básicamente es para certificar la coherencia de la altitud media sobre el nivel del mar de la zona de estudio.
De acuerdo a la finalidad del ejercicio, cargamos el MDE en el directorio de GRASS que hemos generado,
para ello la secuencia sería (ver Figura 4 ):
Importante en éste apartado, saber cual es la nomenclatura de análisis para la cuenca seleccionada. De
aquí se desprende la necesidad de recopilar datos referentes al área de interés y construir la memoria.
Concepto que hemos presentado en capítulos anteriores.
La Figura 5 indica uno de los procesos más importantes para poder procesar el MDE con GRASS,
recordemos que, hay versiones de QGIS que tienen resuelta la interrelación con GRASS; y la secuencia
que se detalla está directamente implementada.
Las coordenadas de la sección de interés, constituyen la limitante a considerar más importante, porque
básicamente corresponde al escenario de trabajo. Para ello, en el apartado región, en el directorio abierto,
se delimita la zona de interés en el MDE.
ATENCIÓN: Hemos notado que dependiendo de la versión de QGIS, hay formatos del complemento
GRASS diferentes. Importante saber que, para poder hacer procesamientos en GRASS, debo tener
cargado el MDE propiamente en GRASS (recordemos que hay versiones automáticas) y darle coordenadas
a la capa raster de interés, cómo en la Figura 6:
Nuevamente consideramos el directorio abierto (ver Figura 7). A partir del MDE, que constituye un mapa
de elevaciones, se genera el primer mapa de flujo (ver Figura 8).
El procedimiento consiste en rellenar todas las depresiones tratando de encontrar una dirección única para
cada celda (unidad más pequeña: píxel). Si el programa detecta áreas con depresiones delimita esta zona,
del resto del área y rellena nuevamente las depresiones usando la técnica de los vecinos más próximos1 .
La salida final será una capa sin depresiones y con direcciones de flujo única.
1
Algoritmo computacional de seleccionar la menor distancia entre un punto y otro.
La Figura 9 resume la tendencia principal de escorrentía dependiendo del relieve y la división de áreas
en la sección por altura. Constituyen los datos de entrada para empezar a pensar en cómo se forma una
cuenca a partir del relieve (representado por el MDE).
Resumiendo para los resultados una división del área de interés por topografía (zonas altas y bajas)
y detalle de escurrimiento de flujo de acuerdo a la distribución de elevaciones ya enunciada y litología
presente.
(a) (b)
Figura 9: Productos de salida a) Mapa relieve de la dirección de flujo sin depresiones y b) Topografía
Aplicación de GRASS que delimita las subcuencas de cada curso de agua, a lo largo de la red de drenaje.
Se extrae a partir del mapa de flujos generado en el paso anterior, en la Figura 10 detalles del procedimiento.
Cabe destacar que, dentro de las opciones de la solapa, el número de celdas alude a la cantidad de píxeles
que se agrupan según las características del relieve.
Figura 10: Apartado para elaborar características de la cuenca a partir del mapa de flujos
El mapa de la cuenca generado a partir del mapa de flujos suavizado, depende directamente del MDE que
contiene las características de altura y orientación de laderas correspondiente a la topografía.
La Figura 11 resume los límites de cuenca posibles, y se distingue en todos los productos, la necesidad de
acotar las posibilidades de distribución de la red de drenaje para captar un área más exacta. Eso lo vemos
en el apartado siguiente; antes te proponemos probar opciones de número de celda (50, 100, 200, 350,
500, etc.) fijándote cómo varía la delimitación de la cuenca.
(a) (b)
(c) (d)
Figura 11: Mapa de cuencas generado a partir de a) 50 unidades, b) 250 unidades, c) 350 unidades y d) 500 unidades
El MDE que tenemos a disposición tiene cómo unidad mínima un píxel de 30×30m por lo que una celda
ocupa una superficie de 900 m2 equivalente a 0,09 Ha. Recordemos que en QGIS, podemos medirlo todo
con la función:
El área de la cuenca es inferior a 600 Has. para cada subcuenca entonces, utilizando la famosa regla de 32
los números de celdas adecuados, corresponderían aproximadamente a 6666 unidades.
Para certificar lo adecuado del cálculo, probamos el análisis y efectivamente la Figura 24, muestra una
distribución de subcuencas conforme al relieve; recordemos que trabajamos con datos de altura y el relieve,
responde a las características de la roca, las condiciones del medio y por ende define la disposición general
de la/las cuenca/s.
c) Subcuencas.
2
Regla de tres simple: Operación matemática que nos ayuda a resolver rápidamente problemas de proporcionalidad.
(a) (b)
(c) (d)
Figura 12: Detalle de los productos resultantes a) Red de drenaje principal por segmentos, b) Relieve por el cual discurre la red
de drenaje conformando la cuenca, c) Subcuencas y d) Áreas de captación y acumulación de la red de drenaje
Ésta función, permite supervisar la generación del área de la cuenca, según convergencia de redes de
drenaje, detalle importante para modelar éste tipo de geomorfología, demostrado en el apartado anterior.
Para ello cargamos en el programa QGIS la red de drenaje del punto anterior, definiendo el punto según
coordenadas geográficas que se ajuste al cierre de la cuenca de interés, y copiamos la ubicación del punto
en LATITUD y LONGITUD:
Se despliega una ventana cómo se muestra en la Figura 13; donde se carga el relieve por el cual discurre
la red de drenaje y las coordenadas de cierre de la cuenca.
Figura 13: Relieve por el cual discurre la red de drenaje conformando la cuenca
En el teórico detallamos el porque de la distribución de una red de drenaje, en éste apartado hacemos
foco en el punto geométrico de convergencia de todos los causes (o la gran mayoría). Cobra especial
importancia para el análisis, el mapa de segmentos de la red de drenaje.
La figura 14, resume cómo hacer la extracción de las coordenadas de la geolocalización del punto de unión
de todos los cauces que pertenecen a la cuenca que estamos tratando de delimitar.
Una vez que obtengamos el píxel de cierre de la cuenca, ubicamos el cursor en el centro de la celda
y copiamos las coordenadas que especifica el directorio de QGIS en la parte inferior, esas son las
coordenadas que se cargan en la solapa presentada en el principio de ésta sección.
Figura 14: Pasos para extraer las coordenadas del punto de cierre de una cuenca
¿Y sí no aparece?
Puede que el raster resultante no distinga la información, ese arreglo corresponde sólo al formato visual,
para corregirlo, recurrimos a los paneles detallados más abajo y completamos las características requeridas
en estilo (ver Figura 16).
La cuenca resultante para ser utilizada en correlación con otros productos, requiere ser pasada a formato
vectorial, ya que cómo raster no suma ni resta información por píxel, sólo delimita la presencia/ausencia de
la geoforma.
En formato vectorial básicamente funciona cómo herramienta de corte ó máscara, fundamentalmente para
datos raster, que sí contienen información variable por píxel (ver Figura 18).
ATENCIÓN: Los productos obtenidos del complemento GRASS se tienen que guardar, eso se efectúa en
las propiedades de la capa. Acordate de tener una carpeta específica para archivos raster (geoTIFF) y otra
para vectoriales (.shp o .kml). Una vez guardados, el directorio de GRASS puede ser cerrado, ya que para
ésta práctica, no lo seguiremos indagando.
A partir del vector cuenca obtenido, recortar el MDE cargado en QGIS (no en el directorio GRASS) y las
bandas 5, 4 y 3 de las imágenes LANDSAT 8 OLI.
No es el objetivo de ésta práctica, sin embargo a disposición para profundizar, ¡¡CONSULTANOS!! Hay
muchas imágenes satelitales y MDE disponibles para descargar en forma gratuita.
En el programa QGIS, cargar: Las bandas de LANDSAT 8 (543), el MDE, ubicar el vector cuenca generado
en el apartado anterior, cómo molde. La secuencia y la figura 18 detallan los pasos a seguir para encontrar
la función y hacer el recorte.
Una vez seleccionados los datos, dejar correr la función con la rueda verde (extremo inferior derecho).
¡¡Esperá!! puede demorar un poco en recortar todas las capas.
Una vez obtenidos los recortes, eliminá del panel de capas las bandas y MDE sin recortar (sí está en la
carpeta de tu directorio adecuadamente guardado, podés volver a cargarlo las veces que lo necesites.).
El consejo es para que el programa no funcione lento en tú computadora, al tener info de más cargada
(primer motivación para recortar datos raster).
Sí demora demasiado en recortar, podés seleccionar las capas raster de a una y repetir el procedimiento.
Para hacer recortes, combinaciones RGB, operaciones de banda, recortar, entre muchas otras opciones,
QGIS ofrece un serie de complementos útiles de instalar en el software previo a trabajar, se sugiere instalar
SCP (Semi Automatic Classification), OpenLayers (para visualizar Google Earth, Bing Maps, entre otros, en
la pantalla), se supone lo usamos en otros cursos, por las dudas, refrescamos:
Para hacer un análisis mineralógico y civil de la zona de estudio, combinar éstas bandas resulta muy útil
cómo primer aproximación. Tiene relación directa con la reflectancia de los objetos.
Las bandas propuestas corresponden al infrarrojo y al rojo del ya presentado en módulos anteriores
espectro electromagnético, por lo que las firmas de suelo alterado, agua estancada, población y
vegetación vigorosa, resaltarán automáticamente con ésta combinación.
La secuencia en INFRARROJO, se resumen en la Figura 19, cuidar el orden en la cual se disponen las
bandas y la configuración de la longitud de onda.
Figura 20: Combinación INFRARROJO (543). Corresponde a bandas del infrarrojo y rojo de LANDSAT 8. En celeste canteras de
áridos perfectamente delimitadas, en rojo vegetación vigorosa correspondiente a una adquisición en temporada estival, en gris
ciudad y en azul oscuro dique La Quebrada
El apartado anterior, necesita ser validado con datos oficiales, es por ello que dada la instancia, recurrimos
a nuestra carpeta de datos para realizar la práctica y cargamos los vectores disponibles, detallados en el
enunciado.
La figura 21 resume toda la información, recordar que dentro de las propiedades de cada capa, en el
programa QGIS, podemos modificar el estilo.
Figura 21: Datos vectoriales disponibles: En celeste red de drenaje, en blanco con contorno rojo canteras, los puntos
corresponden a localidades y en celeste con contorno blanco al dique "La Quebrada"
Delimitar zonas de influencia y la correlación de las mismas, paso importante en el marco de un análisis
ambiental. A modo de ejemplo en éste apartado aprendemos a preparar un buffer e intersección para la
cuenca.
Recurrimos al marco teórico de la unidad para entender el ¿Por qué?; y en la secuencia detallada a
continuación (ver Figura 22) cómo prepararlo en el programa QGIS.
ATENCIÓN: Nuestros datos vectoriales están expuestos en grados, sí cambiamos esa disposición de
coordenadas, debemos hacerlo para todo el proyecto, o convertir el valor en grados a unidades en metros.
Figura 22: Ventana para generar buffer del vector que se desee
Para generar la intersección es importante seleccionar el buffer o zona de influencia de 2.5Km de diámetro
correspondiente a las canteras y a la red de drenaje, y correlacionar para obtener las secciones de cruce
de las áreas.
Figura 23: Ventana para generar la intersección del buffer que se desee
ATENCIÓN: Repetimos, las unidades para éste proyecto están en grados, la distancia 0.01 corresponde a
2.5Km de diámetro de la zona de influencia a partir de un vector. A veces cuando nos perdemos haciendo
cálculos y no sabemos cuanto exactamente mide el resultado, lo podemos respaldar con la escala o
comparar el equivalente en distancia (sí conocemos la región), de ese modo certificamos que el resultado
sea coherente.
La Figura 24 resume las zonas de influencia elegidas para las canteras y la red de drenaje y su
correspondiente correlación a través de la función intersección.
Los detalles de estilo se pueden cambiar en el apartado propiedades, para presentar resultados en un
formato más prolijo, la escala de colores (va de rojo a amarillo) para las canteras, tiene que ver con la
cercanía a la zona de explotación.
(a) (b)
Figura 24: Detalle de los productos resultantes a) Buffer en diámetro: Red de drenaje de 2.5Km (0.01°) en celeste, canteras a
2.5Km en naranja y 5Km en amarillo respectivamente, b) Intersección: En rojo correlación entre el buffer de 2.5Km de la red de
drenaje y el de 2.5Km de las canteras
En un estudio de impacto ambiental, analizar la cantera, el tipo de productos que requiere para la extracción
de áridos y, en caso de que sea oportuno, una posible contaminación de la red de drenaje, parte de un
análisis de zonas de influencia.
Recordemos que la extracción de áridos no requiere de químicos fuertes, sí hay una modificación del
paisaje y sí hay mucho polvo en suspensión en caso de movimientos grandes, voladuras, etc.
Recordando el concepto de punto de cierre para una red de drenaje, reforcemos la dirección de
escurrimiento del flujo, a través de las curvas de nivel, las mismas permiten al fotointerprete ubicarse en las
zonas de mayor y menor altura con facilidad, y conocer específicamente la equidistancia entre una sección
y otra, obteniendo un mayor detalle de lo abrupto o plano que se presenta el relieve.
Las curvas de nivel, se extraen de un MDE, porque básicamente representan la variación de alturas. La
secuencia siguiente resume cómo extraerlas del programa QGIS y cerrar el análisis de la distribución de la
red de drenaje, según alturas y tipo de litología (ver Figura 27).
Hemos visto ese concepto en la introducción a la práctica, una forma de fotointerpretar isopiezas depende
de la distribución de las curvas de nivel, dejamos ese apartado de tarea para curiosos.
Volviendo al eje, la Figura 26 muestra cómo de la selección de diferentes equidistancias, de las que
obtenemos mayor o menor densidad de curvas de nivel, entrenando el ojo, podemos enunciar que la Figura
a, expresa que la zona escarpada y alta se encuentra en la mitad izquierda de la misma (geográficamente
zona Oeste).
(a) (b)
Figura 26: Detalle de los productos resultantes a) Equidistancia 10-50-100m, b) Equidistancia 100m
A criterio del autor. Acordate que en la unidad I especificamos todo lo que debe poseer un mapa. Dada la
instancia, y los múltiples productos a mostrar, es momento de ponerse creativo y pensar cómo un monitor
ambiental que necesita mostrar resultados, de la forma más clara y concisa que considere.
¡MUCHA SUERTE!
Un modelo . . .
Lista de Acrónimos