Tegucigalpa
Tegucigalpa
Tegucigalpa
Presentado por:
SILVANA PAOLA MORANTES OCHO
CODIGO: 2080857
Presentado a:
JUAN DIEGO COLEGIAL GUTIERREZ
Introducción…………………………………………………………………………Pág. 3
Objetivo ……………………………………………………………………….........Pág. 4
Marco Teórico……………………………………………………………...............Pág. 4
Metodología y resultados………………………………………………………….Pág. 5
Conclusiones………………………………………………………………………..Pág. 22
INTRODUCCION
Tegucigalpa está asentada en un valle profundo localizado al interior de la zona central montañosa de
Honduras, aproximadamente 15°06’N/87°11’W, lo que corresponde con la latitud de Bangkok o de
África Sub-Sahariana. El fondo del valle se localiza sobre los 930 m. s.n.m (sobre el nivel del mar)
con picos (llamados cerros) alrededor de la cuenca que alcanzan hasta 1300m s.n.m, tales como
Cerro Grande y Cerro El Berrinche en el Nor-occidente o El Picacho al norte de la ciudad.
Los ríos que cruzan la cuenca juegan un importante papel en la caracterización morfológica de la
zona. El Río Grande Choluteca separa Tegucigalpa de la ciudad hermana de Comayagüela. El Río
Chiquito entra en la cuenca por el costado oriental, dividiendo la capital entre la ciudad vieja, en la
ribera norte, y el moderno distrito comercial en la ribera sur (figura 1).
Generar un mapa de zonificación cualitativa del riesgo por inundación y deslizamientos para un
sector de la ciudad de Tegucigalpa, Honduras utilizando las herramientas proporcionadas por los
sistemas de información geográfica, que para para lo cual se implementó el software ILWIS.
MARCO TEORICO
Preparación de datos
En esta sección la preparación de datos es general, es decir no solo aplicable para la evaluación de
riego o amenaza específica. Se realizó con diferentes herramientas que permite el software Ilwis las
cuales son útiles en la ejecución los estudios relacionados con el manejo de desastres. Se
necesitaron los siguientes pasos:
Se generó un estéreo par el cual requiere de dos imágenes raster, las cuales se traslapan parcial o
totalmente, por ejemplo dos fotografías aéreas en formato digital que presentan algún porcentaje de
traslape. En el estéreo par de salida se podrá ver el área de traslape en 3 dimensiones o un Modelo
Digital del Terreno (DTM) genera un par estereoscópico “artificial”. En el estéreo par de salida, usted
podrá ver el área del mapa de entrada que coincide con el DTM utilizado, como una imagen en 3-D.
METODOLOGIA Y RESULTADOS
Lidar_ tegu→ pseudo → zoom in → display options → 1maplidar_tegu→Strecth 920-1000 (figura 2).
Operation tree → imagen processing→ filter→ filtro lineal shadow→ mapa shadow (figura 3)
Figura 3: mapa shadow
Se hace clic del 1° punto en la imagen shadow y luego el mismo punto en la imagen ortho_tegu. Se
repitió este procedimiento 5 veces más para otros puntos. Los puntos se seleccionaron en diferentes
partes del mapa Shadow (distribuidos en toda la imagen); en otras palabras, puntos que no se
encuentran agrupados.
Antes del paso siguiente se asignó primero la georeferencia New al mapaLidar_tegu. Luego se realizó
un re-muestreo (resampling) de la imagen Lidar_tegu. Donde “resampling” (re-muestreo) significa
describir el proceso por medio del cual se ajusta el tamaño de píxel de una imagen con respecto a la
otra, de tal manera que el tamaño de píxel de ambos mapas sea igual. Para esto el “re-muestreo”
(ajustar el tamaño de píxel) se hizo utilizando la georeferencia del mapa Ortho_tegu.
Generación del estéreo par: se creó un estéreo par artificial a partir de la imagen LIDAR y la imagen
ortofoto detallada.
Operation tree → imaging processing→ stero par from DTM → raster: Ortho_tegu → DTM:
Lidar_dem → Look modus: left →out put map: Ortho_stereo (figura 11).
Figura 11. Imagen del estero par
Operation tree → imagen processing → color separation → raster: Ortho_tegu → out put map:
Ortho_tegu_BW (figura 12).
Para el estereo par: Operation tree → imaging processing→ stero par from DTM → raster:
Ortho_tegu_BW → DTM: Lidar_dem → Look modus: left →out put map: Stereo_BW (figura 12).
Figura 12: mapa Ortho_tegu_BW
En la primera parte del ejercicio se utilizó una imagen estereoscópica digital para la interpretación de
deslizamientos generada anteriormente y su digitalización en pantalla. Para el proceso de
interpretación usted puede utilizar una imagen anáglifo o utilizar el estereoscopio de pantalla (screen-
stereoscope). La interpretación de deslizamientos se concentrara en la delimitación de los grandes
deslizamientos inactivos que se observan en el área. El anáglifo se generó de la siguiente forma:
Para adicionar la capa de slope es necesario calsificarla, por lo cual se crea el dominio
Slopec mediante la herramienta Slicing (figura 18).
Figura 18: mapa Slopec con intervalos de 10.
Para el análisis de los datos es necesario volver raster los mapas slides y lithology (figura 19)
con la siguiente operación:
Operations →I mage Processing → Slicing.
Para finalizar se genero el mapa ACTIVE, el cual representa a el mapa slides donde las
áreas de no deslizamientos tiene un valor de 0 y las áreas de código Br y Sr tiene un valor de
1. Esto se obtiene mediante la siguiente formula:
Active:=iff(isundef(slides),0,iff((slides="br")or(Slides="sr"),1,0))
El resultado obtenido es el siguiente:
Para este paso se seleccionaron los mapas Slopec y lithology. El mapa Slopec como el primer mapa,
el mapa active como el segundo mapa, y se denominó la tabla de salida Actslope Figura 21). Seguido
se repite el procedimiento pero el primer mapa será el lithology y la tabla de salida se denominó
Actlithology (figura 22).
Una vez cruzados los anteriores mapas, lo siguiente es calcular los valores de densidad. Para lo
siguiente se realizan los siguientes sub pasos:
Crear una columna donde únicamente los deslizamientos activos son señalados. Se utiliza
la siguiente formula: AreaAct=iff(Active=1,area,0).
Calculo del área total en cada clase de pendiente.
Cálculo del área con deslizamientos activos para cada clase de pendiente.
Cálculo del área total del mapa.
Cálculo del área total ocupada por los deslizamientos en la zona de estudio.
Cálculo de la densidad de deslizamientos para cada clase de pendiente.
Cálculo de la densidad de deslizamientos para la totalidad del mapa. Se utiliza la siguiente
formula: Densmap=Areamapact/Areamaptot.
Figura 23: Tabla Actslope corregida.
Dado que el rango de valores obtenidos en WEIGHT es muy amplio y no puede ser presentando
como un mapa de susceptibilidad a deslizamientos es necesario clasificar dichos valores.
Para esto creamos un dominio llamado HAZARD donde seleccionamos límites para 3 clases
diferentes:
LOW HAZARD : 0 – 1.5 .
MODERATE HAZARD: 1.5 – 3.
HIGH HAZARD: 3 – 4.5
Se usó la herramienta “Slicing” para crear un mapa HAZARD donde se muestran las diferentes zonas
de amenaza de deslizamiento en intervalos según los criterios que cada estudiante determino. El
resultado se muestra en la siguiente imagen:
El resultado final obtenido y que es analizado en las conclusiones como mapa final de
amenaza de deslizamiento es el siguiente.
Para los resultados obtenidos, la zona de mayor impacto crítico por deslizamiento se
encuentra en la parte Suroeste del sector central del municipio según los datos de litología,
debido a acumulaciones de limolita y capas rojas de arenas y shale las cuales tienen una
pendiente que oscila entre los 30 – 41 grados representa un alto riesgo de deslizamiento.
Respecto a litologías presentes, las zonas con mayor riesgo por deslizamientos corresponden
a las formaciones del grupo Padre Miguel del terciario y el grupo cretácico valle de ángeles.
Las zonas con mayor probabilidad de presentar movimientos en masa de gran impacto se
encuentran en un intervalo de 38 – 51 grados, según la ponderación en densidades de
deslizamientos, respecto a las pendientes.