Agregación de vialidades

En este tutorial vamos a usar el módulo coberturas para agregar la red vial de la CDMX en una retícula.

Importamos los que vamos a necesitar

import matplotlib.pyplot as plt
import geopandas as gpd
from pyLandUseMX.descargas import *
from pyLandUseMX.usodesuelo import *

red = gpd.read_file("datos/descargas/red_zmvm.gpkg")
# zmvm = gpd.read_file("datos/descargas/silhuetas.shp.zip")
# zmvm = zmvm.loc[zmvm.CVE_SUN == '13']
# red = red.sjoin(zmvm.to_crs(red.crs))

id_acceso_c = [100, 101, 102, 103, 104, 105]
id_primarias = [106, 107]
id_secundarias = 108
id_terciarias = 109
id_residencial = 110

acceso_controlado = red.loc[red.tag_id.isin(id_acceso_c)]
primarias = red.loc[red.tag_id.isin(id_primarias)]
secundarias = red.loc[red.tag_id == id_secundarias]
terciarias = red.loc[red.tag_id == id_terciarias]
residenciales = red.loc[red.tag_id == id_residencial]

mapas = {
    "Vialidades de acceso controlado": {
        "datos":acceso_controlado,
        "color": "darkred",
        "peso": 1},
    "Vialidades primarias": {
        "datos":primarias,
        "color": "tomato",
        "peso": 0.8},
    "Vialidades secundaris": {
        "datos":secundarias,
        "color": "chocolate",
        "peso": 0.6},
    "Vialidades terciarias": {
        "datos":terciarias,
        "color": "tan",
        "peso": 0.5},
    "Calles residenciales": {
        "datos":residenciales,
        "color": "black",
        "peso": 0.02},
}

fig, axs = plt.subplots(3,2,figsize=(15, 15))
axs = axs.ravel()
for i, (k,v) in enumerate(mapas.items()):
    ax = v['datos'].plot(ax=axs[i], linewidth=v['peso'], color=v['color'])
    ax.set_title(k)
    ax.axis('off')

fig.delaxes(axs[-1])
plt.subplots_adjust(wspace=0, hspace=0)
# plt.tight_layout(pad=0.05)

malla = (Malla
         .desde_capa(red, 1000)
         .agrega_lineas(acceso_controlado, campo='Acceso controlado')
         .agrega_lineas(primarias, campo='Vialidades primarias')
         .agrega_lineas(secundarias, campo='Vialidades secundarias')
         .agrega_lineas(terciarias, campo='Vialidades terciarias')
         .agrega_lineas(residenciales, campo='Calles residenciales'))

malla.datos

	grid_id	geometry	Acceso controlado	Vialidades primarias	Vialidades secundarias	Vialidades terciarias	Calles residenciales
0	0	POLYGON ((455987.604 2101522.619, 456987.604 2...	0.0	0.0	0.0	0.0	0.0
1	1	POLYGON ((455987.604 2102522.619, 456987.604 2...	0.0	0.0	0.0	0.0	0.0
2	2	POLYGON ((455987.604 2103522.619, 456987.604 2...	0.0	0.0	0.0	0.0	0.0
3	3	POLYGON ((455987.604 2104522.619, 456987.604 2...	0.0	0.0	0.0	0.0	0.0
4	4	POLYGON ((455987.604 2105522.619, 456987.604 2...	0.0	0.0	0.0	0.0	0.0
...	...	...	...	...	...	...	...
7609	7609	POLYGON ((535987.604 2190522.619, 536987.604 2...	0.0	0.0	0.0	0.0	0.0
7610	7610	POLYGON ((535987.604 2191522.619, 536987.604 2...	0.0	0.0	0.0	0.0	0.0
7611	7611	POLYGON ((535987.604 2192522.619, 536987.604 2...	0.0	0.0	0.0	0.0	0.0
7612	7612	POLYGON ((535987.604 2193522.619, 536987.604 2...	0.0	0.0	0.0	0.0	0.0
7613	7613	POLYGON ((535987.604 2194522.619, 536987.604 2...	0.0	0.0	0.0	0.0	0.0

7614 rows × 7 columns

rasters = malla.to_xarray()

rasters

<xarray.Dataset>
Dimensions:                 (y: 95, x: 82)
Coordinates:
  * y                       (y) float64 2.102e+06 2.102e+06 ... 2.196e+06
  * x                       (x) float64 5.365e+05 5.355e+05 ... 4.555e+05
    spatial_ref             int64 0
Data variables:
    Acceso controlado       (y, x) float64 nan nan nan nan ... 0.0 0.0 0.0 nan
    Vialidades primarias    (y, x) float64 nan nan nan nan ... 0.0 0.0 0.0 nan
    Vialidades secundarias  (y, x) float64 nan nan nan nan ... 0.0 0.0 0.0 nan
    Vialidades terciarias   (y, x) float64 nan nan nan nan ... 0.0 0.0 0.0 nan
    Calles residenciales    (y, x) float64 nan nan nan nan ... 0.0 0.0 0.0 nan

xarray.Dataset

• Dimensions:
  • y: 95
  • x: 82
• Coordinates: (3)
  • y
    (y)
    float64
    2.102e+06 2.102e+06 ... 2.196e+06

    axis :

    Y

    long_name :

    y coordinate of projection

    standard_name :

    projection_y_coordinate

    units :

    metre

    array([2101500., 2102500., 2103500., 2104500., 2105500., 2106500., 2107500.,
           2108500., 2109500., 2110500., 2111500., 2112500., 2113500., 2114500.,
           2115500., 2116500., 2117500., 2118500., 2119500., 2120500., 2121500.,
           2122500., 2123500., 2124500., 2125500., 2126500., 2127500., 2128500.,
           2129500., 2130500., 2131500., 2132500., 2133500., 2134500., 2135500.,
           2136500., 2137500., 2138500., 2139500., 2140500., 2141500., 2142500.,
           2143500., 2144500., 2145500., 2146500., 2147500., 2148500., 2149500.,
           2150500., 2151500., 2152500., 2153500., 2154500., 2155500., 2156500.,
           2157500., 2158500., 2159500., 2160500., 2161500., 2162500., 2163500.,
           2164500., 2165500., 2166500., 2167500., 2168500., 2169500., 2170500.,
           2171500., 2172500., 2173500., 2174500., 2175500., 2176500., 2177500.,
           2178500., 2179500., 2180500., 2181500., 2182500., 2183500., 2184500.,
           2185500., 2186500., 2187500., 2188500., 2189500., 2190500., 2191500.,
           2192500., 2193500., 2194500., 2195500.])

  • x
    (x)
    float64
    5.365e+05 5.355e+05 ... 4.555e+05

    axis :

    X

    long_name :

    x coordinate of projection

    standard_name :

    projection_x_coordinate

    units :

    metre

    array([536500., 535500., 534500., 533500., 532500., 531500., 530500., 529500.,
           528500., 527500., 526500., 525500., 524500., 523500., 522500., 521500.,
           520500., 519500., 518500., 517500., 516500., 515500., 514500., 513500.,
           512500., 511500., 510500., 509500., 508500., 507500., 506500., 505500.,
           504500., 503500., 502500., 501500., 500500., 499500., 498500., 497500.,
           496500., 495500., 494500., 493500., 492500., 491500., 490500., 489500.,
           488500., 487500., 486500., 485500., 484500., 483500., 482500., 481500.,
           480500., 479500., 478500., 477500., 476500., 475500., 474500., 473500.,
           472500., 471500., 470500., 469500., 468500., 467500., 466500., 465500.,
           464500., 463500., 462500., 461500., 460500., 459500., 458500., 457500.,
           456500., 455500.])

  • spatial_ref
    ()
    int64
    0

    crs_wkt :

    PROJCS["WGS 84 / UTM zone 14N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-99],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32614"]]

    semi_major_axis :

    6378137.0

    semi_minor_axis :

    6356752.314245179

    inverse_flattening :

    298.257223563

    reference_ellipsoid_name :

    WGS 84

    longitude_of_prime_meridian :

    0.0

    prime_meridian_name :

    Greenwich

    geographic_crs_name :

    WGS 84

    horizontal_datum_name :

    World Geodetic System 1984

    projected_crs_name :

    WGS 84 / UTM zone 14N

    grid_mapping_name :

    transverse_mercator

    latitude_of_projection_origin :

    0.0

    longitude_of_central_meridian :

    -99.0

    false_easting :

    500000.0

    false_northing :

    0.0

    scale_factor_at_central_meridian :

    0.9996

    spatial_ref :

    PROJCS["WGS 84 / UTM zone 14N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-99],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32614"]]

    GeoTransform :

    537000.0 -1000.0 0.0 2101000.0 0.0 1000.0

    array(0)

• Data variables: (5)
  • Acceso controlado
    (y, x)
    float64
    nan nan nan nan ... 0.0 0.0 0.0 nan

    name :

    Acceso controlado

    long_name :

    Acceso controlado

    _FillValue :

    nan

    array([[nan, nan, nan, ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan]])

  • Vialidades primarias
    (y, x)
    float64
    nan nan nan nan ... 0.0 0.0 0.0 nan

    name :

    Vialidades primarias

    long_name :

    Vialidades primarias

    _FillValue :

    nan

    array([[nan, nan, nan, ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan]])

  • Vialidades secundarias
    (y, x)
    float64
    nan nan nan nan ... 0.0 0.0 0.0 nan

    name :

    Vialidades secundarias

    long_name :

    Vialidades secundarias

    _FillValue :

    nan

    array([[nan, nan, nan, ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan]])

  • Vialidades terciarias
    (y, x)
    float64
    nan nan nan nan ... 0.0 0.0 0.0 nan

    name :

    Vialidades terciarias

    long_name :

    Vialidades terciarias

    _FillValue :

    nan

    array([[nan, nan, nan, ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan]])

  • Calles residenciales
    (y, x)
    float64
    nan nan nan nan ... 0.0 0.0 0.0 nan

    name :

    Calles residenciales

    long_name :

    Calles residenciales

    _FillValue :

    nan

    array([[nan, nan, nan, ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0.,  0., nan]])

• Attributes: (0)

rasters['Acceso controlado'].plot()

<matplotlib.collections.QuadMesh>

(rasters
 .to_array()
 .plot(col = 'variable', col_wrap=3, figsize=(15, 10)))

<xarray.plot.facetgrid.FacetGrid>