marcoajh
commented
7 years ago Mapa de Subsidencia
13 abril 2017
Datos generales del mapa.
Cita oficial:
Texto aquí
Resumen:
En este mapa se expresa la velocidad de subsidencia en centímetros por año en cada AGEB para el periodo 2002 - 2010
Nivel de avance:
Terminado
Ubicación geográfica.
Área geográfica:
Cuenca del valle de México
Coordenadas extremas:
xMín 348621.68
yMín 2054434.61
xMáx 644398.35
yMáx 2344228.00
Calidad de los datos.
Metodología:
Los datos de subsidencia porvienen de un estudio contratado por la Subdirección General Técnica, Gerencia de Aguas Subterráneas de la CONAGUA (Contrato: CNA-SGT-GAS-19-2011 "Estudio de impacto por la subsidencia ocasionada por la extracción, uso y manejo de los recursos hídricos en los acuíferos Texcoco y Chalco-Amecameca, Distrito Federal y Estado de México”.) Consultor: Centro GEO (CONACYT) Año: 2011
CONAGUA, Impacto por la subsidencia ocasionada por la extracción, uso y manejo de los recursos hídricos en los acuíferos Texcoco y Chalco-Amecameca, DF y Estado de México, 2011, 14 pp.
Proceso en GRASS
El raster proporcionado por la CNA, tenía huecos de informacion, para rellenarlos se usó la herramienta Rellenar sin datos de Qgis
Se importa a GRASS
r.import input=/home/marco/CVM_SIG/Subsidencia/Subsidencia_utm.tif output=Subsidencia
Se convierte a raster (era sólo un reclasificado)
r.mapcalc 'Subsidencia = Subsidencia*1' --o
r.mapcalc 'Subsidencia = Subsidencia*100000'
r.mapcalc 'Subsidencia = int(Subsidencia)' --o
r.null map=Subsidencia null=0
Se obtiene la subsidencia promedio en cada AGEB
r.statistics base=AGEB_edos cover=Subsidencia method=average output=AGEB_subs --o
Se obtiene la lista de reclasificación
r.stats -ln AGEB_subs | gawk '{printf("%i %s %i\n", $1, "=", $2)}' > AGEB_subs_rec
Se crea la capa
r.reclass i=AGEB_edos o=AGEB_subsidencia rules=AGEB_subs_rec
Se restablecen los valores decimales
r.mapcalc 'subsidencia = AGEB_subsidencia*0.00001' --o
Se normaliza y se crea la capa
MIN=`r.stats -1n subsidencia | gawk 'NR==1 { MIN=$1; next } $1 < MIN { MIN=$1 } END{ print MIN }'`
MAX=`r.stats -1n subsidencia | gawk 'NR==1 { MAX=$1; next } $1 > MAX { MAX=$1 } END{ print MAX }'`
r.mapcalc "subsidencia_v1 = ((subsidencia - "$MIN")/("$MAX" - "$MIN"))" --o
r.null map=subsidencia_v1 null=0
r.mapcalc 'subsidencia_v1 = if(Pob_saent>=0,subsidencia_v1,null())' --oDatos espaciales.
Estructura del dato:
Vector
Tipo de datos:
Polígonos
Proyección geográfica.
Sistema de coordenadas:
Universal Transversa de Mercator
Proyección:
WGS 84 / UTM zone 14N
Paralelos estándar:
Greenwich
Esferoide:
--
Datum:
World Geodetic System 1984
Elipsoide:
WGS84
Atributos del mapa.
| Nombre del campo | AGEB_ID |
|---|---|
| Tipo | Entero |
| Unidades | -- |
| Descripción | Identificador numérico para cada AGEB |
| Observaciones | -- |
| Nombre del campo | subsidencia |
|---|---|
| Tipo | Real |
| Unidades | -- |
| Descripción | Valor normalizado |
| Observaciones | -- |
Subsidence rate in centimeters per year in each AGEB for the period 2002-2010 Version 1 Ubicación de la capa: https://www.dropbox.com/work/MEGADAPT_Integracion/Layers/Final