Open marcoajh opened 6 years ago
Precipitación anual total por AGEB para el año 2010
Terminado
Cuenca del valle de México
xMín 348621.68
yMín 2054434.61
xMáx 644398.35
yMáx 2344228.00
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
Vector
Polígonos
Universal Transversa de Mercator
WGS 84 / UTM zone 14N
Greenwich
--
World Geodetic System 1984
WGS84
Nombre del campo | AGEB_ID |
---|---|
Tipo | Entero |
Unidades | -- |
Descripción | Identificador numérico para cada AGEB |
Observaciones | -- |
Nombre del campo | pr_2010 |
---|---|
Tipo | Real |
Unidades | -- |
Descripción | Valor normalizado |
Observaciones | -- |
Total annual rainfall for the year 2010 by AGEB Version 1 Ubicación de la capa: https://www.dropbox.com/work/MEGADAPT_Integracion/Layers/Final