Przecinanie GeoDataFrame z wielokątem
Łącząc zbiór danych o użytkowaniu gruntów ze zbiorem danych o dzielnicach, możemy teraz zbadać, jakie jest użytkowanie gruntów w danej dzielnicy.
Aby to zrobić, musimy najpierw wyznaczyć przecięcie zbioru danych o użytkowaniu gruntów z wybraną dzielnicą. Weźmy ponownie dzielnicę Muette jako przykład.
Zbiory danych o użytkowaniu gruntów i dzielnicach zostały już zaimportowane jako land_use i districts, a dzielnica Muette została wyodrębniona do wielokąta shapely o nazwie muette. Biblioteki GeoPandas i matplotlib są również zaimportowane.
To ćwiczenie jest częścią kursu
Praca z danymi geoprzestrzennymi w Pythonie
Instrukcje do ćwiczenia
- Oblicz przecięcie wielokątów z
land_usez pojedynczym wielokątemmuette. Wynik nazwijland_use_muette. - Wykonaj szybki wykres tego przecięcia i przekaż argument
edgecolor='black', aby wyraźniej zobaczyć granice poszczególnych wielokątów. - Wyświetl pierwsze pięć wierszy
land_use_muette.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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.____)