Beräkna skärningspunkten mellan en GeoDataFrame och en polygon
Genom att kombinera markanvändningsdatasetet och distriksdatasetet kan vi undersöka markanvändningen i ett visst distrikt.
För det behöver vi först beräkna skärningspunkten mellan markanvändningsdatasetet och ett givet distrikt. Vi använder distriktet Muette som exempel.
Markanvändnings- och distriktsdataseten har redan importerats som land_use respektive districts, och Muette-distriktet har extraherats till Shapely-polygonen muette. Dessutom är GeoPandas och matplotlib importerade.
Den här övningen är en del av kursen
Arbeta med geospatial data i Python
Övningsinstruktioner
- Beräkna skärningspunkten mellan polygonerna i
land_useoch den enskilda polygonenmuette. Lagra resultatet i variabelnland_use_muette. - Skapa en enkel visualisering av skärningspunkten och ange
edgecolor='black'för att tydligare se gränserna mellan de olika polygonerna. - Skriv ut de fem första raderna i
land_use_muette.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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.____)