Interseccionando um GeoDataFrame com um polígono
Combinando o conjunto de dados de uso do solo e o de distritos, agora podemos investigar qual é o uso do solo em um determinado distrito.
Para isso, primeiro precisamos determinar a interseção do conjunto de dados de uso do solo com um distrito específico. Vamos usar novamente o distrito Muette como exemplo.
Os conjuntos de dados de uso do solo e de distritos já foram importados como land_use e districts, e o distrito Muette foi extraído para o polígono shapely muette. Além disso, GeoPandas e matplotlib estão importados.
Este exercício faz parte do curso
Trabalhando com Dados Geoespaciais em Python
Instruções do exercício
- Calcule a interseção dos polígonos de
land_usecom o polígono únicomuette. Chame o resultado deland_use_muette. - Faça um gráfico rápido dessa interseção e passe
edgecolor='black'para ver mais claramente os limites dos diferentes polígonos. - Imprima as cinco primeiras linhas de
land_use_muette.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Print the land use datset and Notre-Dame district polygon
print(land_use.head())
print(type(muette))
# Calculate the intersection of the land use polygons with Notre Dame
land_use_muette = ____
# Plot the intersection
land_use_muette.____
plt.show()
# Print the first five rows of the intersection
print(land_use_muette.____)