Пересечение двух полигонов
В этом упражнении мы будем работать с двумя отдельными полигонами: районом Мюэт, извлечённым из набора данных districts, и зелёной зоной Булонского леса — большого общественного парка на западе Парижа, извлечённого из набора данных land_use. Оба полигона уже присвоены переменным muette и park_boulogne.
Сначала мы визуализируем эти два полигона. Вы увидите, что они перекрываются, однако парк расположен не целиком в границах района Мюэт. Давайте определим область их пересечения.
GeoPandas и matplotlib уже импортированы.
Это упражнение является частью курса
Работа с геопространственными данными в Python
Инструкции к упражнению
- Отобразите оба полигона на одной карте, чтобы визуально оценить степень их перекрытия.
- Вычислите пересечение полигонов
park_boulogneиmuette. - Выведите долю площади района, которую занимает парк.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# 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(____ / ____)