Layer: Percentage of households without sewage per AGEB "falta_dren_v1"

[marcoajh]

Currently available only for Mexico City. Version 1 Ubicación de la capa: https://www.dropbox.com/work/MEGADAPT_Integracion/Layers/Final

[paola-gp]

@marcoajh: Layer needs to cover the ZMCM, therefore the index needs to change.

[marcoajh]

Mapa de peticion usuarios (drenaje)

13 abril 2017

Datos generales del mapa.

Cita oficial:

Texto aquí

Resumen:

En este mapa se expresa la necesidad de infraestructura de drenaje en cada AGEB.

Nivel de avance:

En desarrollo

Ubicación geográfica.

Área geográfica:

Ciudad 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:

Falta de infraestructura - valor 1 en AGEB sin red de desalojo y % de viviendas sin drenaje en AGEB con red de desalojo

Proceso en GRASS

Desde Qgis, con la herramienta de consulta espacial, se seleccionan las AGEB que no tienen red primaria de desalojo y se importan a GRASS
v.import input=/home/marco/CVM_SIG/AGEB/AGEB_sin_desalojo.shp layer=AGEB_sin_desalojo output=AGEB_sin_desalojo
Se convierte a raster conservando el ID de AGEB
v.to.rast input=AGEB_sin_desalojo output=AGEB_sin_desalojo use=attr attribute_column=ageb_id label_column=cvgeo --o
Se rellenan los nulos con ceros
r.null map=AGEB_sin_desalojo null=0
r.mapcalc 'AGEB_sin_desalojo = float(AGEB_sin_desalojo)' --o

Se convierte a raster la capa AGEB_vivienda con el campo VIV23_R que es el porcentaje de las viviendas particulares que no disponern de drenaje
v.to.rast input=AGEB_vivienda output=viv_sdrenP use=attr attribute_column=VIV23_R_N_ label_column=ageb_id
Se normaliza

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

r.mapcalc 'Falta_dren = if(AGEB_sin_desalojo==0,viv_sdrenP_norm,1)' --o

Se normaliza y se crea la capa

MIN=`r.stats -1n Falta_dren | gawk 'NR==1 { MIN=$1; next }  $1 < MIN { MIN=$1 } END{ print MIN }'` 
MAX=`r.stats -1n Falta_dren | gawk 'NR==1 { MAX=$1; next } $1 > MAX { MAX=$1 } END{ print MAX }'`
r.mapcalc "falta_dren_v1 = float((Falta_dren - "$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	--

Nombre del campo	falta_dren
Tipo	Real
Unidades	--
Descripción	Viviendas sin drenaje / viviendas en el AGEB
Observaciones	--