Sobrepondo conjuntos de dados espaciais

Trabalhando com Dados Geoespaciais em Python

Joris Van den Bossche

Open source software developer and teacher, GeoPandas maintainer

Interseção com um polígono

countries.intersection(circle)

Interseção com um polígono

Limitações de countries.intersection(circle):

Só cruza uma GeoSeries com um único polígono
Não preserva os atributos

Sobrepondo dois conjuntos de dados

countries.plot()

geologic_regions.plot()

Sobrepondo dois conjuntos de dados

geopandas.overlay(countries, geologic_regions, 
                  how='intersection')

Sobreposição vs. interseção

Método de interseção (com um único polígono)

countries.intersection(geologic_region_A)

0                            ()
1    POLYGON ((-1.661 48.803...
2    POLYGON ((1.201 51.145,...
dtype: object

Método de sobreposição

geopandas.overlay(countries, geologic_regions,
                  how='intersection')

     name geologic_region              geometry
1  France               C  POLYGON ((2.5 51....
2      UK               C  POLYGON ((0.7 52 ...
3  France               B  POLYGON ((-1.7 46...
..    ...             ...                   ...

Vamos praticar!

Trabalhando com Dados Geoespaciais em Python