Análisis de sensibilidad

La prueba de sensibilidad por remoción de capas mide la importancia de cada mapa que se utiliza en un índice cartográfico como el que resulta de la aplicación de la combinación lineal ponderada.

Descargar el código de ejemplo

sensibilidad.py.

Requerimientos generales

Para asegurar la ejecución correcta del código es importante verificar la instalación y funcionamiento de los siguientes elementos:

• Qgis 3.4 o superior y librerías de Osgeo4W
• Librerías python:

• copy
• pprint
• string
• osgeo/gdal
• gdal_calc
• os

Requerimientos generales de los insumos

Es importante que todas las capas raster cumplan con las siguientes condiciones:

• Misma proyección cartográfica
• Mismo tamaño de pixel
• Misma extensión de capa
• Mismo valor de NoData

Ejemplo

Insumos

Crear una (1)carpeta en el directorio raíz o en la unidad C que se llame analisis_sensibilidad, para descargar los insumos hacer clic aqui (2)guarde el archivo insumos.zip en la carpeta analisis_sensibilidad, posteriormente hacer clic derecho sobre el archivo y elegir la opción (3)Extract to insumos

Una vez terminado el proceso, crear en la carpeta analisis_sensibilidad una (4) carpeta con el nombre salida

Procedimiento

1. Abrir el código

Abrir el código sensibilidad.py en Qgis 3.4 o superior, Para resolver cualquier duda al respecto, consultar la guia

2. Actualizar el diccionario

Ingresar la (1) ponderación del compontente según corresponda (Exposición, Susceptibilidad, Resiliencia), posteriormente ingresar la (2)ponderación del subcomponente (biológico,físico), Ingresar el (3)nombre de la capa raster de entrada con su respectiva (4)ponderación y su (5)ruta repita los pasos siguiendo la estructura y hasta ingresar cada una de las capas.

3. Indicar el direcctorio de salida

Indicar el directorio donde guardarán los archivos necesarios para realizar el análisis de sensibilidad y el archivo analisis_sensibilidad.csv que contendrá los resultados.

p_procesamiento = 'C:/analisis_sensibilidad/salida/'

4. Ejecutar el código

hacer clic en el (1) botón de ejecutar código, puede demorar 10 minutos o más dependiendo el procesador y memoria RAM que tenga el equipo en donde se ejecute, al concluir aparecerá en la (2) consola una lista que indica que ha procesado cada una de las capas.

Bibliografía

Documentación dentro del código

#Autores Fidel Serrano, Víctor Hernández

Afiliación

Área de Planeación Colaborativa, Laboratorio Nacional de Ciencias de la Sostenibilidad, Instituto de Ecología, Universidad Nacional Autónoma de México

Reconocimiento

Este trabajo fue realizado con financiamiento del Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (PAPIIT IV100118) del proyecto Análisis integrado de sistemas socio-ambientales acoplados: desarrollo de capacidades para la evaluación de la vulnerabilidad costera.

sensibilidad_por_remocion_capas.crea_capa(ecuacion, rasters_input, salida)

Esta función crea una capa mediante la calculadora raster de GDAL, esta función esta limitada hasta 14 variables en la ecuación.

Parámetros:

• ecuacion (String) – ecuación expresada en formato gdal, es este caso es la salida de la funcion ecuacion_clp
• rasters_input (lista) – lista de los paths de los archivos rasters, salida de la función separa_ruta_pesos
• salida (String) – ruta con extensión tiff de la salida

sensibilidad_por_remocion_capas.ecuacion_clp(pesos)

Esta función recibe una lista de pesos para regresar la ecuación en la estructura requerida por gdal para la combinación lineal ponderada.

Parámetros:pesos (lista) – lista de los pesos de las capas, salida de la función separa_ruta_pesos

sensibilidad_por_remocion_capas.ecuacion_vulnerabilidad(n)

Esta función expresa la ecuación para el cálculo de la vulnerabilidad

\[ \begin{align}\begin{aligned}vulnerabilidad = \exp^{( 1 - sus)^{(1 + ca)}}\\| exp = Exposición | sus = Susceptibilidad | ca = Capacidad adaptativa\end{aligned}\end{align} \]

Devuelve:str ecuacion

sensibilidad_por_remocion_capas.get_region(path_layer)

Esta función regresa en forma de cadena de texto las coordenadas de la extensión de una capa raster

param path_layer: ruta de la capa raster type path_layer: str

sensibilidad_por_remocion_capas.lista_criterios(dicc)

Esta función regresa una lista de los criterios de un diccionario

Parámetros:dicc – Diccionario que contiene nombres, rutas y pesos para el

análisis de vulnerabilidad / sensibilidad :type dicc: diccionario python

sensibilidad_por_remocion_capas.lista_pesos_ruta(dicc)

Funcion para sacar listas por subcriterio

sensibilidad_por_remocion_capas.media_raster(path_raster)

Esta función regresa el promedio de todos los pixeles válidos de un archivo raster

Parámetros:path_raster (String) – Ruta del archivo raster

sensibilidad_por_remocion_capas.nombre_capa(path_capa)

Esta función regresa el nombre de una capa sin extensión

Parámetros:path_capa (str) – ruta de la capa

sensibilidad_por_remocion_capas.quita(dicc, key)

Esta función retira un elemento del diccionario y regresa un nuevo diccionario sin dicho elemento <<dicc_q>>.

Parámetros:

• dicc (diccionario) – Diccionario con la estructura requerida
• key (String) – nombre de la variable a quitar

sensibilidad_por_remocion_capas.quita_reescala(dicc, key)

Función que integra las funciones quita y reescala y regresa un diccionario sin la variable y con los pesos reescalados.

Parámetros:

• dicc (diccionario) – Diccionario con la estructura requerida
• key (String) – nombre de la variable a quitar

sensibilidad_por_remocion_capas.raster_min_max(path_raster)

Esta funcion regresa los valores maximos y minimos de una capa raster

Parámetros:path_raster (str) – ruta de la capa raster

sensibilidad_por_remocion_capas.reescala(dicc_q)

Esta función rescala un diccionario que se le a quitado un criterio y regresa el diccionario con los pesos rescalados

Parámetros:dicc_q – salida de la función quita