Taller No. 3 Amenaza ILWIS
Taller No. 3 Amenaza ILWIS
Taller No. 3 Amenaza ILWIS
Por:
C.J. van Westen
Department of Earth Resources Surveys,
International Institute for Aerospace Survey and Earth Sciences (ITC),
P.O. Box 6, 7500 AA Enschede, The Netherlands.
Tel: +31 53 4874263, Fax: +31 53 4874336, e-mail: WESTEN@ITC.NL
Resumen
Este caso de estudio pretende ilustrar el significado de amenaza, vulnerabilidad y
riesgo, para ello se utilizará un conjunto de datos referenciados a escala nacional de
Colombia.
La ocurrencia de un desastre depende de dos factores:
Amenaza: Es la probabilidad de ocurrencia de un fenómeno potencialmente dañino.
También se define como la probabilidad de ocurrencia de un evento potencialmente
desastroso durante cierto período de tiempo en un sitio dado.
Vulnerabilidad: Es el grado de pérdida de un elemento o grupo de elementos bajo
riesgo resultado de la probable ocurrencia de un evento desastroso, expresada en
una escala desde 0 o sin daño a 1 o pérdida total.
Primero se elaborará un mapa cualitativo de la amenaza, para ello cruzará varios mapas
cada uno con diferentes factores. Posteriormente hará un mapa de la vulnerabilidad, y
finalmente estos mapas se cruzarán y obtendrá un mapa de riesgo.
Nota: Es posible que ciertas instrucciones o comandos varíen dependiendo la versión de ILWIS que
utilice. Al final de la práctica debe entregar mapa de amenaza, vulnerabilidad y riesgo y la tabla
diligenciada
Iniciando
Ubíquese en el directorio de trabajo y/o cree uno, por ejemplo c:\AEA
Extraiga los datos del archivo: Datos amenaza vulnerabilidad y riesgo
(appDat01.RAR)
Haga doble clic en el icono del programa ILWIS y ubíquese en el
directorio de trabajo, creado c:\AEA
Doble clic en el mapa de polígonos Colombia. Se muestra el cuadro de
dialogo Display Options en el que puede modificar el mapa.
Seleccione Boundaries Only. Acepte todo por defecto, presione OK. El
mapa será desplegado en la ventana.
Lleve el mapa a un lado de la pantalla.
Arrastre y libere el mapa raster Colombia dentro de la ventana del
mapa. El cuadro de dialogo Display Options está abierto. Clic OK. El
mapa raster Colombia está ahora desplegado junto con los límites del
mapa de polígonos.
Seleccione con el cursor varias unidades del mapa para averiguar su
significado.
Ahora debe mirar los mapas que contienen información de los diferentes tipos de
amenaza en Colombia.
Arrastre y libere el mapa Seismic dentro de la ventana del mapa.
Acepte los defectos en el cuadro de dialogo Display Options y clic OK.
El mapa raster Colombia es ahora reemplazado por el mapa raster
Seismic. Verifique el significado de las unidades de mapeo.
Evalué los otros mapas (Landslid, Volcanic, Tsunami,
Inundat, Rivers, Beach y Topograp) con el mismo
procedimiento.
Abra Pixel information y arrastre y libere los mapas Colombia,
Seismic, Landslid, Volcanic, Tsunami, Inundat, Rivers,
Beach y Topograp. Verifique la información combinada de los mapas
mientras mueve el puntero del mouse en la ventana del mapa.
Arrastre y libere el mapa Seismic en la ventana del mapa. Acepte los
defectos en el cuadro de dialogo Display Options y haga clic OK. El
mapa raster Colombia ahora esta reemplazado por el mapa raster
Seismic.
Verifique el significado de las unidades de mapeo.
Evalúe los otros mapas (Industry, Infraseg y Concentr)
arrastrando y liberándolos en la ventana del mapa. Verifique su
contenido.
Agregue estos mapas también a la ventana de pixel information. Verifique
la información de los mapas combinados mientras mueve el puntero en la
ventana del mapa.
Cierre las ventanas del mapa y de pixel information.
Cuando tenga una idea clara de los datos de entrada continué con el próximo tema.
En este ejercicio estudiará el primer aspecto: amenaza. En los ejercicios posteriores evaluará la
vulnerabilidad y finalmente el riesgo. El grado de la amenaza de una cierta área está determinado
por una combinación de factores. Los diferentes factores, que influyen en el grado de amenaza, se
pueden observar separadamente, aunque ellos se influyan uno al otro (p.e. en un área con una
amenaza alta de sismicidad habrá más deslizamientos de tierras). Los factores se proporcionan en
forma de mapas de parámetro, cada uno de ellos describen un fenómeno potencialmente dañino.
Los fenómenos siguientes son tomados en cuenta:
Deslizamientos
Amenaza de sismicidad
Amenaza tsumani: terremoto induciendo olas de inundación
Amenaza volcánica
Amenaza de inundación
Amenaza de inundación por ríos torrenciales
Erosión de playa y aluvión
Regiones topográficas
Se seguirá el segundo enfoque. El impacto de las diferentes amenazas a la actividad humana no es igual
para todas las clases de amenazas. Los terremotos fuertes tendrán un mayor efecto devastador en un área
que deslizamientos de tierras. Por lo tanto, se tienen que asignar valores de ponderación a cada una de las
clases dentro de los mapas individuales de la amenaza, teniendo en cuenta su importancia a producir
daños.
El valor de ponderación que será dado a un cierto factor y la manera en que este factor se clasificará es
altamente subjetivo. Este método se refiere también en la literatura como ponderación ciega, cuando no
se dispone de datos cuantitativos para decidir si la relación entre la amenaza alto de sismicidad y amenaza
alto de deslizamientos de tierras debe ser o tomar valores entre 1 o 10 o 100. La opinión de expertos es
utilizada para definir estos valores de ponderación. Esto tiene como consecuencia que cada experto puede
asignar valores diferentes.
Posteriormente los valores de ponderación se usarán en vez de los valores del mapa. De esta manera, los
mapas de parámetro separados de los factores implicados en el análisis de la amenaza llegan a ser mapas
de ponderación.
Paso 1: Asignar los valores de ponderación a las clases de los mapas de parámetro.
Los valores de ponderación serán asignados en las tablas conectadas a los mapas raster. Debe hacer una
tabla para cada mapa y luego adicionará una columna weight en cada tabla, en la cual editará los
valores de ponderación para las diferentes clases.
Sismicidad
Leyenda Ponderación
Alto riesgo de sismicidad 10
Moderado riesgo de sismicidad 5
Bajo riesgo de sismicidad 0
Volcánico
Leyenda Ponderación
Sin amenaza volcánico 0
Amenaza volcánica 10
Deslizamiento
Leyenda Ponderación
Nariño 2
Región Huila 2
Valle del Cauca 3
Zona Cafetera 4
Manizales y alrededores 4
Valle de Aburrá 4
Cundinamarca 3
Boyacá y Santander 2
Bucaramanga 2
Sin amenaza de deslizamiento 0
Tsunami
Leyenda Ponderación
Amenaza tsunami 10
Sin amenaza de tsunami 0
Playa
Leyenda Ponderación
Acumulación 0
Acumulación y erosión 1
Erosión 2
Sin acumulación o erosión 0
Todas las clases en el mapa Inundat recibirán un valor de ponderación de 5. La clase “Sin amenaza
de inundación " (no inundation hazard) recibirá un valor de ponderación de 0.
Todas las clases en el mapa River recibirán un valor de ponderación de clase 4. "No río" (no river)
recibirá un valor de ponderación de 0.
Todas las clases en el mapa Topograp recibirán un valor de ponderación de clase 2, y solo la clase
"Altitude menos que 1000m" (altitude less than 1000m) un valor de ponderación de 0.
Aquí solo el primer procedimiento del mapa (Seismic) es explicado. Los otros mapas deben ser
elaborados con el mismo procedimiento.
Seleccione del menú options: File, Create, Create Table.
Escriba el siguiente nombre de archivo: seismic. Seleccione el
dominio Seismic y clic en OK. Un dominio en ILWIS determina el
contenido de un mapa, o una tabla. Si tenemos clases, el dominio es una
lista de posibles nombres de las clases. Una tabla que será unida a un
mapa deberá por eso tener el mismo dominio. La tabla se muestra ahora.
Puede ver los ítems del dominio en la columna gris al lado izquierdo.
Seleccione: Columns, Add Column. Escriba el nombre de la columna:
Weight, y seleccione el dominio: Value. Escriba los valores de 0 y
10 para los value range, y 1.0 para la precisión. Clic OK.
Llene los valores de ponderación de las tres clases, como se indica en la
tabla de la página anterior.
Cierre la tabla.
Ahora ha hecho una tabla con los valores de ponderación para el mapa Seismic.
Repita el mismo procedimiento para los mapas Landslid,
Volcanic, Tsunami y Beach. Creando una tabla para cada mapa
usando su propio dominio.
Antes de pasar al próximo paso tiene que tener en cuenta que las tablas deben ser unidas a los mapas. Esto
se hace cambiando las propiedades de un mapa.
Haga clic con el botón del mouse derecho en el mapa seismic y
escoja desde el menú la opción Properties. Ahora en el cuadro de
dialogo las propiedades del mapa se muestran, tales como el dominio, el
número de pixeles, etc.
Haga un clic sobre la opción Attribute Table y seleccione la tabla
Seismic.
Presione OK. Ahora la tabla seismic está unida al mapa seismic.
Siga el mismo procedimiento para los otros mapas (Landslid,
Volcanic, Tsunami y Beach).
La creación de valores de ponderación para los otros tres mapas (Inundat, River, y Topograp) se
harán de forma diferente. Cuando prácticamente todas las unidades reciben el mismo valor de
ponderación es mejor usar otro método de reclasificación: a partir de la opción Map Calculation o sea
utilizando una fórmula de cálculo de mapa.
Ahora que tiene una tabla con valores de ponderación para los mapas de parámetros, puede renombrar los
mapas a mapas de ponderación. Para esta operación use la operación Attribute map (Attribras).
Desde el menú de la ventana principal, seleccione: Operations, Raster
operations, Attribute Map.
Seleccione Seismic como un mapa raster; seleccione Seismic como la
tabla y Weight como el atributo. Escriba Wseismic como el nombre
de salida del mapa.
Escriba: Mapa Seismic reclasificado con valores de
ponderación como Descripción.
Clic Show y OK.
En el cuadro de dialogo Display Options, seleccione la representación
Pseudo. Clic OK. El mapa de ponderación se muestra.
Si hace un clic sobre una unidad de mapa, no verá la descripción, pero sí
un número, el cual es el valor de ponderación asignado anteriormente.
Para verificar con el mapa original Seismic use el PixelInfo. Puede
hacer esto seleccionando File, Open Pixel Information. En la ventana
del pixel information, seleccione File, Add Map y seleccione el mapa
Seismic. Cuando mueve el puntero sobre el mapa, verá el nombre de
la unidad original en el mapa Seismic y los valores de ponderación.
Cierre la ventana del mapa y la de pixel information.
Repita el mismo procedimiento para los mapas Landslid,
Volcanic, Tsunami y Beach. Creando mapas de ponderación
llamados: Wlandsli, Wvolcani, Wtsunami y Wbeach.
Ahora otro método es usado para reclasificar los mapas Inundat, River y Topograp. Para ello
creará directamente los mapas de ponderación, usando una fórmula, sin almacenar primero los valores de
ponderación en una tabla.
Escriba la siguiente fórmula sobre la línea de comando de la ventana del
ILWIS: Winundat=iff(Inundat=”No inundation hazard”,0,5)
En el cuadro de dialogo Raster Map Definition, seleccione el dominio
Value con un rango de valores entre 0 y 5 y una precisión de 1.0.
Clic OK. Ningún cálculo sucede todavía. Solo la definición se almacena
y el icono del mapa raster para el mapa Winundat se hace. El cálculo
realmente sucede solo después de abrir el mapa.
Haga doble clic en el mapa raster Winundat. El mapa es calculado,
antes el cuadro de dialogo Display Options estaba abierta. Seleccionar la
representación Pseudo, y haga clic en OK.
Diseñe la fórmula para crear el mapa Wriver. Todas las clases del
mapa de entrada River obtendrán el valor de 4, con excepción de la
clase “No river” que tiene 0.
Haga también un mapa Wtopogra, en que todas las clases del mapa
El próximo paso en la creación del mapa de amenaza es la combinación de los mapas de ponderación
individuales.
Escriba la siguiente fórmula en la línea de comando:
Hazard=Wseismic+Wlandsli+Wvolcani+Wbeach+
Wtsunami+Wriver+Winundat+Wtopogra
Clic OK en el cuadro de dialogo Raster Map Definition.
Despliegue el mapa Hazard. Use la representación Pseudo.
Use la ventana de pixel information para leer los valores de los mapas
junto con el mapa de la amenaza. Verifique si las sumas de los valores se
han realizado correctamente.
Cierre las ventanas del mapa y también la de pixel information.
El mapa de amenaza elaborado en la sección anterior tiene muchos valores diferentes. El rango de estos
valores puede ser evaluado calculando un histograma.
En la ventana principal; seleccione del menú: Operations, Statistics,
Histogram.
Seleccione el mapa Hazard y clic OK. El histograma es calculado y los
resultados se muestran en la tabla.
Evalúe los valores en el histograma. Con base a estos valores puede
hacer una subdivisión (marcas de clase) en clases y/o categorías de
amenaza.
Cierre la gráfica y la tabla del histograma.
Desde el rango de los valores puede ver que la mejor manera de clasificar el mapa de amenaza es
dividirlo en cuatro (4) clases.
En la ventana principal, seleccione File, Create, Create Domain. Crear
un nuevo dominio Hazclas. Seleccione Class y Group. Clic OK.
El Domain Editor está abierto.
Seleccione Edit, Add Item y llene las siguientes clases: clases de límites:
0 sin Amenaza
5 Amenaza baja
12 Amenaza media
50 Amenaza alta
- Primero calcular el histograma del mapa Colombia, y averiguar el área (en metros cuadrados)
para cada departamento.
Segundo hacer una tabla Colombia y unir con la tabla del histograma para obtener el área de cada
departamento.
Tercero crear una nueva columna Population en la tabla Colombia, la cual debe llenar con
la población de cada departamento
Luego dividir la columna Population con la columna Area para averiguar la densidad de población.
Finalmente reclasificar el mapa Colombia con la columna Popdens y obtener el mapa de
densidad poblacional.
Calcule el histograma para el mapa raster Colombia.
Cree una nueva tabla, Colombia, con el dominio Colombia.
En la ventana de la tabla, lea las áreas de los departamentos ejecutando
la fórmula siguiente:
Area=Colombia.his.Area
Acepte los valores por defecto.
Cree una columna Population (dominio= Value, con un rango de 0
a 10000000, y una precisión de 1).
Llene la tabla con los siguientes valores.
Nota: Solo debe llenar los datos de población de los departamentos relacionados en el dominio Colombia
Después de ingresar los valores, calcular una columna Popdens
(Densidad de población se define como número de personas por km 2).
Note que los valores en la columna Área están en m2.
Cierre la ventana de tabla.
De la misma forma como preparó los mapas de ponderación, usados para el mapa de la amenaza, varios
mapas de ponderación se deben elaborar para producir un mapa de vulnerabilidad. Los valores de
ponderación, que se emplearán en este ejercicio, están basados sobre la importancia relativa de cada uno
de los elementos de riesgo respecto al daño causado por un desastre. En la práctica estos valores de
vulnerabilidad tendrán diferente ponderación para disímiles tipos de amenaza.
Por ejemplo: la vulnerabilidad en los caminos es menor para inundaciones que para deslizamiento de
tierras o terremotos.
Por ejemplo, si ocurre una inundación, los caminos no se pueden usar temporeramente, mientras que, en
un deslizamiento de tierras o terremoto, los caminos pueden ser (en parte) destruidos y no podrán ser
utilizados. Este criterio de asignación de pesos no será tenido en cuenta.
Los siguientes valores de ponderación son usados:
Industria
Leyenda Ponderación
Oleoductos principales (pipeline) 8
Sin vulnerabilidad (No vulnerability) 0
Centro industrial de 1 orden (1 st order industrial center) 10
Centro industrial de 2 orden (2 nd order industrial center) 7
Concentración
Leyenda Ponderación
Alta (high) 5
Baja (low) 0
Infraestructura
Leyenda Ponderación
Ciudades principales (city) 10
Caminos principales (main road) 8
Ferrovías (railway) 8
Caminos (road) 4
Sin vulnerabilidad (No 0
vulnerability)
popdens=Iff(densidad=0,?,iff(densidad<1,0,Iff(densidad<20,3,Iff(densidad<50,7,10
))))
Densidad poblacional
Leyenda Ponderación
2 2
< 1 persona/km (< 1 person/km ) 0
1-20 personas/km2 (1-20 persons/km2) 3
20-50 personas/km2 7
>50 personas/km2 10
Puesto que el cálculo del mapa de vulnerabilidad es muy semejante al del mapa de amenaza, no todos los
pasos se explican con detalle, tenga en cuenta este aspecto.
Cree las tablas para los mapas Industry, Concentr, e Infraseg,
y llene los valores de ponderación.
Cree mapas de ponderación utilizando los mapas Industry,
Concentr, e Infraseg, y las tablas que ya ha creado recientemente.
Los mapas resultantes deben llamarse Windustr, Wconcent, y
Winfras.
Cree los valores de ponderación para densidad de población en la tabla
Colombia usando los valores de la anterior página.
Cree el mapa de ponderación Wpopdens para reclasificar el mapa
Colombia con la densidad poblacional reclasificada.
Combine los cuatro mapas de ponderación Windustr, Wconcent,
Winfrase, y Wpopdens mediante suma. El mapa resultante es llamado
Vulnerab.
Cree un dominio Vulclas tipo class/group con 4 clases (Muy baja
vulnerabilidad, Baja vulnerabilidad, Moderada
vulnerabilidad, Alta vulnerabilidad).
Clasificar el mapa Vulnerab en 4 clases usando la operación slicing.
El resultado se llamará Vulclas.
Hasta aquí tiene un mapa de la amenaza y otro de vulnerabilidad, ahora combínelos para obtener el mapa
final de riesgo.
La etapa final en un análisis de la amenaza es la creación de un mapa de amenaza. Según Varnes (1984),
la amenaza se puede definir como "el grado esperado de pérdida debido a fenómenos naturales
particulares". En un proyecto real esto sería la multiplicación de costos* vulnerabilidad* período de
retorno del fenómeno natural que daña. En esta aplicación no se disponen de datos de costo ni del
período de retorno.
En la ventana principal, seleccione File, Create, Create 2
Dimensional Table.
Ingrese Table Name Risk. Seleccione Hazclas como el Primary
domain, y Vulclas como el Secondary domain.
Presione el botón Create del cuadro de la lista del dominio (Botón
amarillo). El Create domain se abre. Cree un dominio tipo clase: Risk,
con tres clases (Riesgo Bajo, Moderado, y alto riesgo).
Cierre el dominio, y haga clic OK en el cuadre de dialogo Create de la
tabla de 2 Dimensiones.
La tabla de 2 dimensiones está desplegada. Contiene los valores indefinidos para todas las
combinaciones. Abajo se han llenado ya algunos de ellos, por ejemplo:
Cuando la amenaza es muy bajo, no importa si la vulnerabilidad es baja o alta; el riesgo será bajo
en todos los casos
Cuando la vulnerabilidad es muy baja (significa que el área no contiene ningún elemento
importante en el riesgo), el riesgo es siempre bajo.
Llene las clases con que tiene el signo? en la tabla de dos dimensiones
arriba. Justifique la asignación de las clases del riesgo alto, medio, bajo o
sin riesgo.
Haga clic en el campo izquierdo superior de la tabla de 2 dimensiones.
Edite los campos usando la flecha hacia la izquierda (). Oprima la
flecha hacia abajo () para ir al próximo campo. Complete toda la tabla.
Es más rápido trabajar con la tecla de la flecha, que hacer clic en cada
campo, seleccione los nombres del lado derecho de la tabla, y luego haga
clic en otro campo. Sólo tiene que hacer clic en el campo superior de la
próxima columna cuando ha terminado con el último campo de la
columna anterior.
Cuando haya terminado cierre la ventana.
La tabla de dos dimensiones debe estar completa. Si prefiere puede hacer los cálculos en la línea de
comando. La fórmula del cálculo debe tener la sintaxis siguiente:
Output map = Two-dim table [Map1, Map2]
Escriba la siguiente fórmula en la línea de comando de la ventana
principal:
Risk=Risk[Hazclas,Vulclas]
Presione Enter. El cuadro de dialogo del raster Map Definition está
abierto. Clic OK.
Abra el mapa Risk y verifique su contenido.
Cierre la ventana del mapa.
Elabore la leyenda para el mapa Risk e imprímalo si tiene una
impresora en línea.
La creación de un mapa de riesgo no es el paso final en un manejo del desastre. Es simplemente una
herramienta que debería ser usada para tomar decisiones en el proceso de planificación.
Uno de los aspectos importantes es evaluar y determinar el porcentaje de riesgo ya sea alto, medio o bajo
en cada departamento. Estos porcentajes se pueden calcular, así:
Cruce los mapas Colombia y Risk y cree una tabla de cruce.
Abra la tabla de cruce y escriba la siguiente fórmula en la línea de
comando:
High:=iff(Risk=“High”,Npix,0)
También calcule la columna Medium y Low de la misma manera.
Cierre la tabla de cruce y abra la tabla Colombia. Use la operación join
para leer el valor máximo de las columnas High, Medium, y Low de la
tabla de cruce. Use la tabla join otra vez para obtener la suma de la
columna Npix de la tabla de cruce. Calcule los porcentajes de riesgo alto,
medio y bajo para cada departamento.
Encuentre los cinco departamentos de Colombia con mayor amenaza (donde el
mayor porcentaje de su área es clasificada como alto riesgo).
Referencias
Romero, J. A., Florez, A. and Sanchez, H.A. (1989). Inventario inicial de riesgos
Naturales en Colombia. Análisis Geográficos No. 16. Instituto Geográfico Augustín
Codazzi, Bogotá, Colombia. 56 pp.
Varnes, D.J. (1984). Landslide Hazard Zonation: a review of principles and practice.
Commission on landslides of the IAEG, UNESCO. Natural Hazards No. 3, 61 pp.
Van Westen, C.J., Van Duren, I, Kruse, H.M.G. and Terlien, M.T.J. (1993). GISSIZ:
training package for Geographic Information Systems in Slope Instability Zonation.
ITC-Publication Number 15, ITC, Enschede, The Netherlands.Volume 1: Theory,
245 pp. Volume 2: Exercises, 359 pp. 10 diskettes.
Velasquez, A. and Meyer, H.J. (1990). Un ensayo de evaluación de las amenazas, de los
riesgos y de los desastres en Colombia. Proceedings 1er Seminario Andino de
Geología Ambiental, 1er Conferencia Colombiana de Geología Ambiental, 3er
Conferencia de Riesgos Geológicos del Valle de Aburra. AGID Report No.13, paper
42, pp 547-580.