Mapa de Precipitación anual total por AGEB 2013

25 de septiembre de 2017

Datos generales del mapa.

Cita oficial:

Resumen:

Precipitación anual total por AGEB para el año 2013

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:

Se interpolaron datos de precipitación de 49 estaciones hidrometricas de SACMEX correspondientes al periodo 2000-2014. Posteriormente se calculó la precipitación promedio para cada AGEB

Proceso en GRASS

Interpolación de datos de precipitación para Ciudad de México Datos proporcionados por SACMEX

Hofierka J., Parajka J., Mitasova H., Mitas L., 2002, Multivariate Interpolation of Precipitation Using Regularized Spline with Tension. Transactions in GIS 6, pp. 135-150. http://onlinelibrary.wiley.com/doi/10.1111/1467-9671.00101/epdf

Con los datos de la tabla de precipitacion anual total se creó un vectorial de puntos con las estaciones, se agregó el campo elev

Se importa a GRASS
v.import input=/home/marco/CVM_SIG/Precipitacion/Interp_SACMEX/Total_Annual_Rainfall.shp layer=Total_Annual_Rainfall output=Tot_An_Rain --o

Se cruza con el mapa de elevacion para obtener la elevación de cada estación.
v.what.rast map=Tot_An_Rain raster=CVM_DEM_15m column=ELEV

Se convierte a vector 3D
v.to.3d input=Tot_An_Rain out=Tot_An_Rain_3d column=ELEV

Se crea una máscara con un buffer a 8.5 km de Ciudad de México (la estación más lejana está a poco más de 8 km)
r.buffer input=Delegaciones output=mask_interp distances=8500 units=meters --o
r.mapcalc 'mask_interp = if(mask_interp>=1,1,0)' --o
Se establece como máscara
r.mask raster=mask_interp

Se establece la resolución Cima-Fondo de malla 3d
g.region t=5500 b=0 tbres=5500

Se interpolan los datos para los años 2000-2014

for ((year=2000; year<=2014; year++))
do
echo $year
v.vol.rst input=Tot_An_Rain_3d wcolumn=RAIN_$year tension=5 smooth=0.1 cross_input=CVM_DEM_15m elevation=RAIN_$year cross_output=rain_out_$year maskmap=mask_interp --o
r.colors map=rain_out_$year color=precipitation
done

Los mapas interpolados se convierten a enteros

for ((year=2000; year<=2014; year++))
do
echo $year
r.mapcalc 'rain_out_'$year' = int(rain_out_'$year')' --o
done

Se calculan cruzas de precipitación por AGEB para 2000-2014

for ((year=2000; year<=2014; year++))
do
echo $year
r.statistics base=AGEB_edos cover=rain_out_$year method=average output=precip_AGEB_$year --o
r.stats -ln precip_AGEB_$year | gawk '{printf "%i %s %i\n", $1,"=",$2}' > /home/marco/CVM_SIG/Precipitacion/Interp_SACMEX/precip_ageb_$year
done

Con los archivos obtenidos, se crearon las tablas de los archivos vectoriales usando QGis

Se crean las capas para cada año

for ((year=2000; year<=2014; year++))   
do   
echo $year   
r.reclass i=AGEB_edos o=Pr_"$year"_v1 rules=/home/marco/CVM_SIG/Precipitacion/Interp_SACMEX/precip_ageb_$year --o   
done

Se normaliza usando el mínimo y máximo del total de los datos (2000-2014)

for ((year=2000; year<=2014; year++))   
do   
echo $year
r.mapcalc 'Pr_'$year'_v1 = float(Pr_'$year'_v1)' --o
done

for ((year=2000; year<=2014; year++))   
do   
echo $year      
r.mapcalc 'Pr_'$year'_v1 = ((Pr_'$year'_v1 - 241)/(1387 - 241))' --o   
done

Datos 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	pr_2013
Tipo	Real
Unidades	--
Descripción	Valor normalizado
Observaciones	--

sostenibilidad-unam / layers

Layer: Pr_2013_v1 #53