Generar capa de elevación para el criterio de exposición #PAPIIT

[vichdzgeo]

Objetivo

Ruta de trabajo

sig_papiit\procesamiento\exposicion\fisica\elev_yuc\

Insumos

• cem_yuc_30m.tif

[vichdzgeo]

Procesamiento

1. descargar el insumo cem_yuc_30m.tif INEGI
2. Alinear insumo con el script alinear_cem_extent_1.py -> cem_yuc_30m_v2.tif
3. Cortar el insumo con la capa buffer_yuca_terrestre_3k.shp ubicada en : sig_papiit/procesamiento/procesamiento_inicial/buffer_yucatan/buffer_yuca_terrestre_3k/

Elevación mínima del raster: -16 Elevación máxima del raster: 31

Aplicar función de valor

La función de valor acordada en términos de exposición se comporta como una función logística y la expresión es la siguiente:

Para la construcción de la función de valor los parámetros son Saturación = 3 Centro = 7

[CarlosCruzVz]

Insumo

cem_yuc_30m_clip.tif

Procesamiento

1. Se añadió el insumo a QGIS
2. Del menu Raster, se utilizó la calcuadora raster
3. En la zona de Expresión de la calculadora de campos se escribio el siguiente comando

   1 /  ( 1+ ( 2.71828 ^  (  ( -0.083499* ( -  (  ( 100*7 ) - ( 100 *  ( -16 )  )  / ( 31 -  (  - 16)  )  ) +  ( (  ( 100* ( "cem_yuc_30m_clip@1" )  ) - ( 100 *  ( -16 )  )  / ( 31 -  (  - 16)  )  )   )  )   )  )  )  )  

4. En la zona de Capa de resultado, se coloco el nombre de la capa fv_elev_yuc.tif
5. Se guardo en la siguiente ruta Dropbox\sig_papiit\procesamiento\exposicion\fisica\elev_yuc

Resultado

Como resultado se genero la siguiente imagen

[vichdzgeo]

Corrección de capa

En la ubicacion del puerto de abrigo de san crisanto se localiza que la capa de elevación no contiene datos, sin embargo esta dentro de la zona de estudio, causando conflicto en la integración de las capas.

Se procede a corregir la capa de exposición física de la siguiente manera

( "exp_fisica@1" > 2.9 ) ( 0.13 "fv_distancia_playa@1")+ ( "exp_fisica@1"<2.9) * "exp_fisica@1"