Изучение результата наложения
Теперь, когда мы создали наложение наборов данных о землепользовании и административных районах, анализировать землепользование по отдельным районам стало значительно удобнее. Вернёмся к примеру с районом Мюэт и изучим структуру землепользования в нём.
Библиотеки GeoPandas и Matplotlib уже импортированы. Результат функции overlay() из предыдущих упражнений доступен как combined.
Это упражнение является частью курса
Работа с геопространственными данными в Python
Инструкции к упражнению
- Добавьте в GeoDataFrame
combinedновый столбец'area'с площадью каждого полигона. - Создайте подмножество
land_use_muette, отфильтровав строки, в которых'district_name'равно "Muette". - Постройте график
land_use_muette, используя столбец'class'для раскраски полигонов. - Вычислите суммарную площадь для каждого значения
'class'вland_use_muetteс помощью методаgroupby()и выведите результат.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Print the first rows of the overlay result
print(combined.head())
# Add the area as a column
____ = ____
# Take a subset for the Muette district
land_use_muette = combined[____]
# Visualize the land use of the Muette district
land_use_muette.____(____)
plt.show()
# Calculate the total area for each land use class
print(land_use_muette.____(____)['area'].____() / 1000**2)