sostenibilidad-unam / layers

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

Layer: Age of the construction of the AGEB. "antiguedad_v1" #2

Open marcoajh opened 7 years ago

marcoajh commented 7 years ago

Based on urban data coverage from the 1950´s to 2000´s. Seven categories. Version 1 https://www.dropbox.com/work/MEGADAPT_Integracion/Layers/Final

marcoajh commented 7 years ago

Mapa de Antigüedad

13 abril 2017

Datos generales del mapa.

Cita oficial:
Resumen:

Antigüedad de cada AGEB en años

Nivel de avance:

Terminada

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 antigüedad del crecimiento de la macha urbana fueron proporcionado por el equipo de urbasnismo a cargo del Mtro. Sergio Flores.

Proceso en GRASS
Se importa la capa de antigüedad de la mancha urbana
v.import input=/home/marco/CVM_SIG/Urbano_Crecim_Hist/UTM/Crec_Urb.shp layer=Crec_Urb output=Crec_Urb
Se convierte a raster
v.to.rast input=Crec_Urb@Marco output=Crec_Urb_hist use=attr attribute_column=CLAVE_EDAD label_column=EDAD

Los datos de antigüedad se presentan en siete intervalos
1 - 4 a 6 años
2 - 6 a 9 años
3 - 9 a 14 años
4 - 14 a 18 años
5 - 18 a 33 años
6 - 33 a 64 años
7 - Mas de 64 años

Se cruzan los datos para saber a qué categoría de antigüedad y en cuánta superficie corresponde cada AGEB.
r.stats -lan AGEB_edos,Crec_Urb_hist separator=, > AGEB_ant.csv

Con esos datos y usando una hoja de cálculo, se asigna el limite inferior del intervalo de antigüedad correspondiente a la mayor superficie en el AGEB, se crea un archivo de reclasificación antiguedad.rec

Se reclasifica la capa de AGEB
r.reclass i=AGEB_edos o=antiguedad rules=antiguedad.rec --o

Se normaliza y se crea la capa

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

Con el archivo antiguedad.rec se crea un archivo vectorial desde QGis

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