Applicare operazioni spaziali personalizzate

Lavorare con i dati geospaziali in Python

Joris Van den Bossche

Open source software developer and teacher, GeoPandas maintainer

Lunghezza totale dei fiumi entro 50 km da ogni città?

Per un singolo punto (cairo):

area = cairo.buffer(50000)

rivers_within_area = rivers.intersection(area)

print(rivers_within_area.length.sum() / 1000)

186.397219642

Il metodo apply()

Series.apply(): esegue una funzione su ciascun valore della Series

Series.apply(function, **kwargs)

function: la funzione chiamata su ogni valore; il valore è passato come primo argomento
**kwargs: argomenti aggiuntivi passati alla funzione

Per una GeoSeries, la funzione è chiamata come function(geom, **kwargs) per ogni geom nella GeoSeries

Applicare un'operazione spaziale personalizzata

La funzione da applicare:

def river_length(geom, rivers):
    area = geom.buffer(50000)
    rivers_within_area = rivers.intersection(area)
    return rivers_within_area.length.sum() / 1000

Chiamare la funzione sulla singola geometria:

river_length(cairo, rivers=rivers)

186.3972196423455

Applicarla a tutte le città:

cities.geometry.apply(river_length, rivers=rivers)

Applicare un'operazione spaziale personalizzata

Applicarla a tutte le città:

cities.geometry.apply(river_length, rivers=rivers)

0        0.000000
1        0.000000
2      106.072198
          ...

Applicare un'operazione spaziale personalizzata

Applicarla a tutte le città e assegnare il risultato a una nuova colonna:

cities['river_length'] = cities.geometry.apply(river_length, rivers=rivers)
cities.head()

             name                       geometry  river_length
0    Vatican City    POINT (1386304.6 5146502.5)      0.000000
1      San Marino    POINT (1385011.5 5455558.1)      0.000000
2           Vaduz    POINT (1059390.7 5963928.5)    106.072198
..            ...                            ...           ...

Esercitiamoci!

Lavorare con i dati geospaziali in Python