Interseção de dois polígonos
Neste exercício, vamos usar 2 polígonos individuais: o distrito de Muette, extraído do conjunto de dados districts, e a área urbana verde de Boulogne, um grande parque público na zona oeste de Paris, extraída do conjunto de dados land_use. Os dois polígonos já foram atribuídos às variáveis muette e park_boulogne.
Vamos primeiro visualizar os dois polígonos. Você vai ver que eles se sobrepõem, mas o parque não está totalmente dentro do distrito de Muette. Vamos determinar a parte sobreposta.
GeoPandas e matplotlib já estão importados.
Este exercício faz parte do curso
Trabalhando com Dados Geoespaciais em Python
Instruções do exercício
- Plote os dois polígonos em um único mapa para examinar visualmente o grau de sobreposição
- Calcule a interseção dos polígonos
park_boulogneemuette. - Imprima a proporção da área do distrito que é ocupada pelo parque.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Plot the two polygons
geopandas.GeoSeries([park_boulogne, muette]).plot(alpha=0.5, color=['green', 'blue'])
plt.show()
# Calculate the intersection of both polygons
intersection = ____
# Plot the intersection
geopandas.GeoSeries([intersection]).plot()
plt.show()
# Print proportion of district area that occupied park
print(____ / ____)