In [1]:
import pandas as pd
import geopandas as gpd

In [2]:
departamentos = {
    'Comuna 1':'02001',
    'Comuna 2':'02002',
    'Comuna 3':'02003',
    'Comuna 4':'02004',
    'Comuna 5':'02005',
    'Comuna 6':'02006',
    'Comuna 7':'02007',
    'Comuna 8':'02008',
    'Comuna 9':'02009',
    'Comuna 10':'02010',
    'Comuna 11':'02011',
    'Comuna 12':'02012',
    'Comuna 13':'02013',
    'Comuna 14':'02014',
    'Comuna 15':'02015',
    'ALMIRANTE BROWN':'06028',
    'AVELLANEDA':'06035',
    'BERAZATEGUI':'06091',
    'ESCOBAR':'06252',
    'ESTEBAN ECHEVERRIA':'06260',
    'EZEIZA':'06270',
    'FLORENCIO VARELA':'06274',
    'GENERAL RODRIGUEZ':'06364',
    'GENERAL SAN MARTIN':'06371',
    'HURLINGHAM':'06408',
    'ITUZAINGO':'06410',
    'JOSE C PAZ':'06412',
    'LA MATANZA':'06427',
    'LANUS':'06434',
    'LOMAS DE ZAMORA':'06490',
    'MALVINAS ARGENTINAS':'06515',
    'MARCOS PAZ':'06525',
    'MERLO':'06539',
    'MORENO':'06560',
    'MORON':'06568',
    'PILAR':'06638',
    'PRESIDENTE PERON':'06648',
    'QUILMES':'06658',
    'SAN FERNANDO':'06749',
    'SAN ISIDRO':'06756',
    'SAN MIGUEL':'06760',
    'SAN VICENTE':'06778',
    'TIGRE':'06805',
    'TRES DE FEBRERO':'06840',
    'VICENTE LOPEZ':'06861'
}

In [3]:
#CARGO CARTOGRAFIA NACIONAL NO OFICIAL
pais = gpd.GeoDataFrame.from_file('data/shapes/paisNoOficial/rad_paisp.shp')
print pais.shape
pais.head()


(52800, 11)
Out[3]:
AREA DEPTO FRAC LINK PERIMETER PROV RADIO RADPAIS_ RADPAIS_ID TIPO geometry
0 207407200.0 077 01 380770101 77083.800 38 01 2 2 M POLYGON ((3470906.5 7569104.5, 3470872.5 75691...
1 279095700.0 077 01 380770102 81703.910 38 02 3 3 R POLYGON ((3487339 7588166.5, 3492205 7586542, ...
2 110506900.0 077 02 380770202 48507.350 38 02 4 4 M POLYGON ((3495669.75 7577169.5, 3499717.5 7576...
3 151297.6 077 01 380770103 1736.481 38 03 5 5 U POLYGON ((3494430.25 7574442.5, 3494745 757444...
4 462011900.0 077 03 380770302 126718.300 38 02 6 6 M POLYGON ((3470906.5 7569104.5, 3471078.25 7569...

In [4]:
#CREO ID REDCODE PARA TODAS LAS PROVINCIAS PORQUE CABA NO TIENE LINK
pais['REDCODE'] = pais.PROV + pais.DEPTO + pais.FRAC + pais.RADIO
#CREO EL CODIGO DEPARTAMENTO PARA FILTRAR LOS PERTENECIENTES A PARTIDOS DEL AGBA
pais['DEPARTAMENTO'] = pais.REDCODE.map(lambda x: x[0:5])

In [5]:
#FILTRO LOS RADIOS QUE ESTAN EN DEPARTAMENTO
print pais.shape
pais = pais.loc[pais.DEPARTAMENTO.isin(departamentos.values()),:]
print pais.shape


(52800, 13)
(13550, 13)

In [6]:
#LOS RADIOS FALTANTES DE CABA SON URBANOS
pais.TIPO.fillna('U',inplace=True)

In [7]:
pais = pais.loc[pais.TIPO!='R',:]

In [8]:
pais.head()


Out[8]:
AREA DEPTO FRAC LINK PERIMETER PROV RADIO RADPAIS_ RADPAIS_ID TIPO geometry REDCODE DEPARTAMENTO
26688 16948780.0 252 01 062520101 26902.490 06 01 26690 26690 M POLYGON ((4158793.5 6185928, 4158824.25 618588... 062520101 06252
26704 12351870.0 252 01 062520102 30469.020 06 02 26706 26707 M POLYGON ((4159777 6185123.5, 4159847 6185160, ... 062520102 06252
26719 11063170.0 252 02 062520208 25350.040 06 08 26721 26723 U POLYGON ((4161392.5 6184536.5, 4161399.75 6184... 062520208 06252
26741 796720.3 252 02 062520201 6843.383 06 01 26743 26743 U POLYGON ((4164165.5 6183091, 4164223.5 6183021... 062520201 06252
26747 2996825.0 252 02 062520202 8393.760 06 02 26749 26749 U POLYGON ((4166774.25 6182343.5, 4166960.5 6182... 062520202 06252

In [11]:
departamentosName = {
    '02001':'Comuna 1',
    '02002':'Comuna 2',
    '02003':'Comuna 3',
    '02004':'Comuna 4',
    '02005':'Comuna 5',
    '02006':'Comuna 6',
    '02007':'Comuna 7',
    '02008':'Comuna 8',
    '02009':'Comuna 9',
    '02010':'Comuna 10',
    '02011':'Comuna 11',
    '02012':'Comuna 12',
    '02013':'Comuna 13',
    '02014':'Comuna 14',
    '02015':'Comuna 15',
    '06028':'ALMIRANTE BROWN',
    '06035':'AVELLANEDA',
    '06091':'BERAZATEGUI',
    '06252':'ESCOBAR',
    '06260':'ESTEBAN ECHEVERRIA',
    '06270':'EZEIZA',
    '06274':'FLORENCIO VARELA',
    '06364':'GENERAL RODRIGUEZ',
    '06371':'GENERAL SAN MARTIN',
    '06408':'HURLINGHAM',
    '06410':'ITUZAINGO',
    '06412':'JOSE C PAZ',
    '06427':'LA MATANZA',
    '06434':'LANUS',
    '06490':'LOMAS DE ZAMORA',
    '06515':'MALVINAS ARGENTINAS',
    '06525':'MARCOS PAZ',
    '06539':'MERLO',
    '06560':'MORENO',
    '06568':'MORON',
    '06638':'PILAR',
    '06648':'PRESIDENTE PERON',
    '06658':'QUILMES',
    '06749':'SAN FERNANDO',
    '06756':'SAN ISIDRO',
    '06760':'SAN MIGUEL',
    '06778':'SAN VICENTE',
    '06805':'TIGRE',
    '06840':'TRES DE FEBRERO',
    '06861':'VICENTE LOPEZ'
}

deptos = []
for i in range(pais.shape[0]):
    deptos.append(departamentosName[pais.DEPARTAMENTO.iloc[i]])
pais['DEPTO_NAME'] = deptos
pais.head()


Out[11]:
AREA DEPTO FRAC LINK PERIMETER PROV RADIO RADPAIS_ RADPAIS_ID TIPO geometry REDCODE DEPARTAMENTO DEPTO_NAME
26688 16948780.0 252 01 062520101 26902.490 06 01 26690 26690 M POLYGON ((4158793.5 6185928, 4158824.25 618588... 062520101 06252 ESCOBAR
26704 12351870.0 252 01 062520102 30469.020 06 02 26706 26707 M POLYGON ((4159777 6185123.5, 4159847 6185160, ... 062520102 06252 ESCOBAR
26719 11063170.0 252 02 062520208 25350.040 06 08 26721 26723 U POLYGON ((4161392.5 6184536.5, 4161399.75 6184... 062520208 06252 ESCOBAR
26741 796720.3 252 02 062520201 6843.383 06 01 26743 26743 U POLYGON ((4164165.5 6183091, 4164223.5 6183021... 062520201 06252 ESCOBAR
26747 2996825.0 252 02 062520202 8393.760 06 02 26749 26749 U POLYGON ((4166774.25 6182343.5, 4166960.5 6182... 062520202 06252 ESCOBAR

In [12]:
pais.to_file('data/shapes/agba')

In [ ]: