sostenibilidad-unam / layers

data layers
GNU General Public License v3.0
0 stars 0 forks source link

Layer: Media pressure "pres_med_v1" #29

Open marcoajh opened 7 years ago

marcoajh commented 7 years ago

Media pressure, based on the number of Economic Units in each AGEB (census block)

Version 1

https://www.dropbox.com/work/MEGADAPT_Integracion/Layers/Final

marcoajh commented 7 years ago

Mapa de Presion de medios

10 febrero 2017

Datos generales del mapa.

Cita oficial:

-

Resumen:

En este mapa se expresa la presión en medios de comunicación basada en la cantidad de Unidades Económicas en cada AGEB

Nivel de avance:

Aún en procesamiento

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:

Con datos del sistema SINCE de INEGI, se creó una capa de puntos combinando las unidades económicas de DF, Edo Mex, Hidalgo y Tlaxcala. Con ésta capa y la capa de AGEB, se usó el geoalgoritmo de qgis "Contar puntos en polígono" para obtener un campo con el número total de unidades económicas en cada AGEB. Posteriormente se importó la capa vectorial a GRASS

Proceso en GRASS
Se importa la capa
v.import input=/home/marco/CVM_SIG/AGEB_UE/AGEB_UE.shp layer=AGEB_UE output=AGEB_UE
Se convirete a raster

v.to.rast --overwrite input=AGEB_UE@Marco output=AGEB_UE use=attr attribute_column=NUMPOINTS

r.mapcalc 'AGEB_UE = float(AGEB_UE)' --o

Se normaliza y se crea la capa

MIN=`r.stats -1n AGEB_UE | gawk 'NR==1 { MIN=$1; next } $1 < MIN { MIN=$1 } END{ print MIN }'` 
MAX=`r.stats -1n AGEB_UE | gawk 'NR==1 { MAX=$1; next } $1 > MAX { MAX=$1 } END{ print MAX }'` 
r.mapcalc "pres_med_v1 = ((AGEB_UE - "$MIN")/("$MAX" - "$MIN"))" --o

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 Texto
Nombre del campo pres_med
Tipo Real
Unidades --
Descripción Valor normalizado
Observaciones --