search

lancis-apc / espejos-lancis

GNU General Public License v3.0
0 stars 1 forks source link

Generar el modelo de vulnerabilidad costera para BCS #papiit #75

Open edithvilla opened 5 years ago

edithvilla commented 5 years ago

Modelo de Vulnerabilidad Costera para el Pacífico Norte

ruta: \PAPIIT\pap_insumos\integracion\PVuln\Zofemat\Modelos - Susceptibilidad, resiliencia, exposición

Exposición

Resiliencia

Susceptibilidad

edithvilla commented 5 years ago

Información geográfica del modelo de Vulnerabilidad costera para el Pacífico Norte

ruta: \PAPIIT\pap_insumos\integracion\zofemat\II-UNAM\Datos

Criterio de Susceptibilidad

\PAPIIT\pap_insumos\integracion\zofemat\II-UNAM\Datos\Hidrodinámica

image

image

 `Geomorfología`

\PAPIIT\pap_insumos\integracion\zofemat\II-UNAM\Datos\Geomorfología\Clasificación\1_TIPO_DE_COSTA

image

image

image

edithvilla commented 5 years ago

Criterio de Susceptibilidad

\PAPIIT\pap_insumos\integracion\zofemat\II-UNAM\Datos\Ecología

image

image

image

edithvilla commented 5 years ago

Criterio de Resiliencia

ruta: \PAPIIT\pap_insumos\integracion\zofemat\II-UNAM\Datos\Ecología

image

edithvilla commented 5 years ago

Criterio de Exposición

image

image

image

image

image

image

image

image

image

edithvilla commented 5 years ago

Funciones de valor:

Ruta: \PAPIIT\pap_insumos\integracion\PVuln\Zofemat\FV Zofemat

FuncionesValor_ZOFEMAT-PN.xls

en las siguientes carpetas se encuentran los archivos de illustrator (en la misma ruta)

edithvilla commented 5 years ago

Modelo multicriterio

ruta: \PAPIIT\pap_insumos\integracion\PVuln\Zofemat\Modelos - Susceptibilidad, resiliencia, exposición

Modelos&Pesos_Enero2014.xlsx

en este archivo viene la imagen del superdecisions y los pesos locales de los criterios

edithvilla commented 5 years ago

Shells

Ruta: \PAPIIT\pap_insumos\PVuln\SIG_PacificoNorte_Grass

en esta carpeta se encuentran varios archivos de tipo txt, así también se encuentra la carpeta "shells"

Dentro de la carpeta "shells" se encuentran las carpetas:

Exposicion, Resiliencia, Suceptibilidad, en cada carpeta se encuentran archivos de tipo .sh y también una bitácora

edithvilla commented 5 years ago

@yosunemiquela @paola-gp @vichdzgeo

en este issue encontrarán la documentación del modelo de vulnerabilidad costera de pacífico norte

edithvilla commented 5 years ago

@yosunemiquela @vichdzgeo

Informe del modelo de vulnerabilidad costera en La Paz

ruta: \PAPIIT\pap_informes\informe_año2\vulnerabilidad_lapaz.docx

edithvilla commented 5 years ago

@CorinaCantero

Información para generar los metadatos de las capas geográficas para el modelo de vulnerabilidad de La Paz

Para el criterio Exposición:

Marino: 

Terrestre: 

Para el criterio susceptibilidad

Físico:   

CEPAL. (2011). Efectos del cambio climático en la costa de América Latina y el Caribe. Dinámicas, tendencias y variabilidad climática. ONU - Universidad de Cantabria (Vol. I). Santiago de Chile. Retrieved from http://www.cepal.org/id.asp?id=48025

Silva-Casarín, R., Lacouture, M., Durón, F., Pedroza-Paez, D., Pérez, M. A., Baldwin, E., … Cid, A. (2011). Caracterización de la zona costera y planteamiento de elementos técnicos para la elaboración de criterios de regulación y manejo sustentable.

ruta: \PAPIIT\pap_insumos\integracion\zofemat\II-UNAM\II-UNAM_ManualCFE.pdf

edithvilla commented 5 years ago

@CorinaCantero

Información para generar los metadatos de las capas geográficas para el modelo de vulnerabilidad de La Paz

Archivo: RelaciónCapasSIG_EcorregiónMarina.xlsx

ruta: /sig/desarrollo/sig_papiit/procesamiento/lapaz

En esta misma ruta voy a colocar las capas que se utilizan en cada criterio:

edithvilla commented 5 years ago

@CorinaCantero

Exposición

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\terrestre liga de descarga: https://www.inegi.org.mx/app/geo2/elevacionesmex/

Susceptibilidad

Resiliencia

ruta de las capas: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\resiliencia

edithvilla commented 5 years ago

@CorinaCantero ya está el proyecto en QGIS de las capas que son insumos para el modelo de vulnerabilidad costera de La Paz, BCS

RUTA: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\lapaz.qgs

Dentro del proyecto en qgis, las capas están organizadas por criterio: Resiliencia, Susceptibilidad y me hace falta exposición, la cual mañana la agrego

edithvilla commented 5 years ago

Documentación del mapa de oleaje

esta documentación la encontré en el archivo Suceptibilidad.sh ruta: \PAPIIT\pap_insumos\PVuln\SIG_PacificoNorte_Grass\shells\Suceptibilidad y por lo que entiendo ocupa el shape: H_media.shp ruta: \PAPIIT\pap_insumos\zofemat\II-UNAM\Datos\Hidrodinámica\Oleaje

~~~~~~~~~

Hidrodinámica

~~~~~~~~~

REGIMEN MAREAL

La region del Pacifico Norte presenta un régimen mesomareal y es constante a lo largo de la costa

En este caso en particular, el valor de "Regimen mareal" está dado por la función de valor de "Oleaje"

# ##############################

OLEAJE

Función de valor: Coseno

min=0.48

max=1.18

Parámetro a.... 0.0

Parámetro b.... 0.4571428571

Parámetro c.... 1.0

Parámetro d.... 1.0

#

Se normaliza el mapa al que se le aplicará la función de valor coseno

r.mapcalc 'oleaje_norm=((H_media-0.48)/(1.18-0.48))'

#

Se crea el mapa de la sección creciente de la curva (en radianes)

r.mapcalc 'alpha1=(1-((oleaje_norm-0.0)/(0.4571428571-0.0)))*(3.14159265358979/2)'

Se convierte a grados (la función coseno de GRASS necesita que x esté dado en grados)

r.mapcalc 'alpha1gra=(alpha1*180)/3.14159265358979'

Se crea el mapa de la sección decreciente de la curva (en radianes)

r.mapcalc 'alpha2=((oleaje_norm-1.0)/(1.0-1.0))*(3.14159265358979/2)'

Se convierte a grados (la función coseno de GRASS necesita que x esté dado en grados)

r.mapcalc 'alpha2gra=(alpha2*180)/3.14159265358979'

Se crea el mapa normalizado de la función de valor coseno

r.mapcalc 'fv_oleaje_2=if(oleaje_norm_2 < Parametro_a,0,if(oleaje_norm_2 > Parametro_d,0,if(oleaje_norm_2 >= Parametro_a && oleaje_norm_2 <= Parametro_b,(exp(cos(alpha1gra_2),2)),if(oleaje_norm_2 > Parametro_b && oleaje_norm_2 < Parametro_c,1,if(oleaje_norm_2 >= Parametro_c && oleaje_norm_2 <= Parametro_d,(exp(cos(alpha2gra_2),2)))))))'

r.mapcalc 'fv_oleaje=if(oleaje_norm < 0.0,0,if(oleaje_norm > 1.0,0,if(oleaje_norm >= 0.0 && oleaje_norm <= 0.4571428571,(exp(cos(alpha1gra),2)),if(oleaje_norm > 0.4571428571 && oleaje_norm < 1.0,1,if(oleaje_norm >= 1.0 && oleaje_norm <= 1.0,(exp(cos(alpha2gra),2)))))))'

Multiplicación por restricciones (cuenca baja)

r.mapcalc 'fv_oleaje=fv_oleaje*cuenca_baja'

#

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_oleaje_5c=if(fv_oleaje<0.2,1,if(fv_oleaje<0.4,2,if(fv_oleaje<0.6,3,if(fv_oleaje<0.8,4,5))))'

Se establece tabla de colores

r.colors m=fv_oleaje_5c rules=color_fv.txt

# ############################################################################################################### ##############################

OLEAJE

Función de valor: Coseno

min=0.705

max=0.92

Parámetro a.... -0.023256

Parámetro b.... 1.0

Parámetro c.... 1.0

Parámetro d.... 1.0

#

Se normaliza el mapa al que se le aplicará la función de valor coseno

r.mapcalc 'oleaje_norm=((H_media-0.705)/(0.92-0.705))' #

Se crea el mapa de la sección creciente de la curva (en radianes)

r.mapcalc 'alpha1=(1-((oleaje_norm-(-0.023256))/(1.0-(-0.023256))))*(3.14159265358979/2)'

Se convierte a grados (la función coseno de GRASS necesita que x esté dado en grados)

r.mapcalc 'alpha1gra=(alpha1*180)/3.14159265358979'

Se crea el mapa de la sección decreciente de la curva (en radianes)

r.mapcalc 'alpha2=((oleaje_norm-1.0)/(1.0-1.0))*(3.14159265358979/2)'

Se convierte a grados (la función coseno de GRASS necesita que x esté dado en grados)

r.mapcalc 'alpha2gra=(alpha2*180)/3.14159265358979'

Se crea el mapa normalizado de la función de valor coseno

r.mapcalc 'fv_oleaje=if(oleaje_norm < Parametro_a,0,if(oleaje_norm > Parametro_d,0,if(oleaje_norm >= Parametro_a && oleaje_norm <= Parametro_b,(exp(cos(alpha1gra),2)),if(oleaje_norm > Parametro_b && oleaje_norm < Parametro_c,1,if(oleaje_norm >= Parametro_c && oleaje_norm <= Parametro_d,(exp(cos(alpha2gra),2)))))))'

r.mapcalc 'fv_oleaje=if(oleaje_norm < -0.023256,0,if(oleaje_norm > 1.0,0,if(oleaje_norm >= -0.023256 && oleaje_norm <= 1.0,(exp(cos(alpha1gra),2)),if(oleaje_norm > 1.0 && oleaje_norm < 1.0,1,if(oleaje_norm >= 1.0 && oleaje_norm <= 1.0,(exp(cos(alpha2gra),2)))))))'

Multiplicación por restricciones (cuenca baja)

r.mapcalc 'fv_oleaje=if(cuenca_baja>=0,fv_oleaje,null())'

Se aplica restricción de elevación a 20 metros sobre el nivel del mar

r.mapcalc 'elev_20=if(Elevacion<=20,1,null())' r.mapcalc 'fv_oleaje=elev_20*fv_oleaje' #

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_oleaje_5c=if(fv_oleaje<0.2,1,if(fv_oleaje<0.4,2,if(fv_oleaje<0.6,3,if(fv_oleaje<0.8,4,5))))'

Se establece tabla de colores

r.colors m=fv_oleaje_5c rules=color_fv.txt #

edithvilla commented 5 years ago

Informe de vulnerabilidad costera para BCS

Para el criterio de Susceptibilidad con respecto a Oleaje y Régimen mareal dice lo siguiente:

Los procesos de Hidrodinámica se componen de Oleaje y Régimen mareal. La región de Pacífico Norte presenta un Régimen mareal de tipo mesomareal (mayor a 1 metro). Por consiguiente, el cúmulo Hidrodinámica se representó sólo con la función de valor de Oleaje.

edithvilla commented 5 years ago

Documentación del Criterio exposición / Marino / Aumento del nivel del mar

/bin.sh

Shell para evaluar la Exposición

Se establece la máscara general

r.mask -r r.mask i=mascara5@marco maskcats=1 # ###################################

MEDIO MARINO

###################################

Aumento del nivel del mar

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

Elevación SLR

~~~~~~~~~

Elevación del nivel del mar

Se usarán las elevaciones de 0 a 5 msnm como restricción para este atributo

Se crea el mapa de restriccion por elevación

r.mapcalc 'restr_elev=if(Elevacion < 5,1.0,0)'

Función de valor para elevación del nivel del mar

r.mapcalc 'rslr=float(RSLR_2040_20k)' r.mapcalc 'fv_rslr=if(rslr<63.88,0.06,if(rslr<64.75,0.125,if(rslr<66.5,0.25,if(rslr<70,0.5,1.0))))'

Multiplicación por restricciones

r.mapcalc 'fv_rslr=fv_rslr*restr_elev'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_rslr_5c=if(fv_rslr==0,0,if(fv_rslr<0.2,1,if(fv_rslr<0.4,2,if(fv_rslr<0.6,3,if(fv_rslr<0.8,4,5)))))'

Se establece tabla de colores

r.colors m=fv_rslr_5c rules=color_fv.txt

Pendiente SLR

~~~~~~~~~

Función de valor de Pendente

# r.mapcalc 'fv_pend_rslr=if(pend_gra<=2,1.00,if(pend_gra<=6,0.33,0))'

Multiplicación por restricciones

r.mapcalc 'fv_pend_rslr=fv_pend_rslr*restr_elev'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_pend_rslr_5c=if(fv_pend_rslr==0,0,if(fv_pend_rslr<0.2,1,if(fv_pend_rslr<0.4,2,if(fv_pend_rslr<0.6,3,if(fv_pend_rslr<0.8,4,5)))))'

Se establece tabla de colores

r.colors m=fv_pend_rslr_5c rules=color_fv.txt #

Cálculo del atributo Elevación del nivel del mar

r.mapcalc 'fv_aum_nm=(fv_rslr0.5)+(fv_pend_slr0.5)'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_aum_nm_5c=if(fv_aum_nm==0,0,if(fv_aum_nm<0.2,1,if(fv_aum_nm<0.4,2,if(fv_aum_nm<0.6,3,if(fv_aum_nm<0.8,4,5)))))' r.null map=fv_aum_nm_5c setnull=0

Se establece tabla de colores

r.colors m=fv_aum_nm_5c rules=color_fv.txt

edithvilla commented 5 years ago

Documentación del Criterio exposición / Marino / Ciclones tropicales

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

Ciclones tropicales

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

#

El atributo se califica por la altura de la ola asiciada a tormentas y la geomofologia de la costa (sistemas abiertos o sistemas semicerrados)

Se aplica como restricción la cuenca baja (zonas funcionales)

Función de valor de la Altura significativa de la ola

# r.mapcalc 'fv_Hs=if(Ciclones_Hs<4,0.06,if(Ciclones_Hs<5,0.125,if(Ciclones_Hs<6,0.25,if(Ciclones_Hs<8,0.5,1.00))))' # #

Función de valor de Ciclones tropicales

La capa Geo_Costa corresponde a la geomorfología de la costa

#

cat Etiqueta valor

1.... sistemas abiertos........... 1.00

2.... sistemas cerrados........... 0.33

# r.mapcalc 'fv_cic_trop=if(Geo_Costa==1,(fv_Hs1.00),if(Geo_Costa==2,(fv_Hs0.33)))' #

Se aplica restricción de elevación a 20 metros sobre el nivel del mar

r.mapcalc 'elev_20=if(Elevacion<=20,1,null())' r.mapcalc 'fv_cic_trop=fv_cic_trop*elev_20'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_cic_trop_5c=if(fv_cic_trop<0.2,1,if(fv_cic_trop<0.4,2,if(fv_cic_trop<0.6,3,if(fv_cic_trop<0.8,4,5))))'

Se establece tabla de colores

r.colors m=fv_cic_trop_5c rules=color_fv.txt # #

Cálculo de Medio Marino

# r.mapcalc 'Medio_marino=(fv_aum_nm0.5)+(fv_cic_trop0.5)'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'Medio_marino_5c=if(Medio_marino<0.2,1,if(Medio_marino<0.4,2,if(Medio_marino<0.6,3,if(Medio_marino<0.8,4,5))))'

Se establece tabla de colores

r.colors m=Medio_marino_5c rules=color_fv.txt

edithvilla commented 5 years ago

Documentación del Criterio exposición / Terrestre

###################################

MEDIO TERRESTRE

################################### #

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

Elevación

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

#

Se restringe la elevacion a sólo 20 msnm

r.mapcalc 'elev20=if(Elevacion<=20,Elevacion,-9999)' r.null m=elev20 setnull=-9999

Se aplica la función de valor a la elevacion:

funcion de valor:Campana

Escala mínima................... 0

Escala máxima................... 20

Valor máximo de Y en X.......... 0

Amplitud........................ 10

Unidades........................ metros

r.mapcalc 'fv_elev20=exp ((-exp((elev20-0)/10, 2))-(0.0183156388887342))/((1)-(0.0183156388887342))' #

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_elev20_5c=if(fv_elev20<0.2,1,if(fv_elev20<0.4,2,if(fv_elev20<0.6,3,if(fv_elev20<0.8,4,5))))'

Se establece tabla de colores

r.colors m=fv_elev20_5c rules=color_fv.txt # #

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

Pendiente

»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»»

#

Se reclasifica el mapa de pendiente en grados para asignar la función de valor

r.mapcalc 'Pendiente_costa=if(pend_gra<=2,1,if(pend_gra>2 && pend_gra<=6,2,0))' r.mapcalc 'fv_pend_costa=if(Pendiente_costa==1,1.00,if(Pendiente_costa==2,0.33,0))'

Se aplica restricción de elevación a 20 metros sobre el nivel del mar

r.mapcalc 'elev_20=if(Elevacion<=20,1,null())'

r.mapcalc 'fv_pend_costa=fv_pend_costa*elev_20'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'fv_pend_costa_5c=if(fv_pend_costa<0.2,1,if(fv_pend_costa<0.4,2,if(fv_pend_costa<0.6,3,if(fv_pend_costa<0.8,4,5))))'

Se establece tabla de colores

r.colors m=fv_pend_costa_5c rules=color_fv.txt # #

Funcion de valor de Medio Terrestre

# r.mapcalc 'Medio_terrestre=(fv_pend_costa0.50)+(fv_elev200.5)'

Multiplicación por restricciones (continente)

r.mapcalc 'Medio_terrestre=if(uga_ter>0,Medio_terrestre,null())'

Multiplicación por restricciones (cuenca baja)

r.mapcalc 'Medio_terrestre=Medio_terrestre*cuenca_baja'

Reclasificación de la función de valor en cinco categorías

r.mapcalc 'Medio_terrestre_5c=if(Medio_terrestre<0.2,1,if(Medio_terrestre<0.4,2,if(Medio_terrestre<0.6,3,if(Medio_terrestre<0.8,4,5))))'

Se establece tabla de colores

r.colors m=Medio_terrestre_5c rules=color_fv.txt

edithvilla commented 5 years ago

@CorinaCantero estoy colocando las capas geográficas en el proyecto de qgis del criterio Exposicón

Aumento del nivel del mar / Elevación y Pendiente

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\marino

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\terrestre

edithvilla commented 5 years ago

Para el criterio Terrestre se ocupan las siguiente información geográfica:

Elevación de la costas y Pendiente de la costa

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\marino

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\terrestre

edithvilla commented 5 years ago

@CorinaCantero

Exposición / Marino / Aumento del nivel del mar

El insumo que se requiere para el criterio de Pendiente y Elevación es:

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\marino

también están las siguientes capas que se indican en el shell:

ruta: \SIG\desarrollo\sig_papiit\procesamiento\lapaz\exposicion\marino

CorinaCantero commented 5 years ago

Susceptibilidad/ Biótico / Acuáticos

El insumo que se ocupa para Marismas, es el correspondiente a la vegetación halófila en la Serie IV Uso de Suelo y Vegetación de INEGI. Esto está indicado en el archivo "RelaciónCapasSIG_EcorregiónMarina", ubicado en:

SIG\desarrollo\sig_papiit\procesamiento\lapaz

edithvilla commented 5 years ago

Tipo de proyección para Pacífico Norte

image