Перетин двох полігонів
У цій вправі ми використаємо 2 окремі полігони: район Мюетт, отриманий з набору даних 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(____ / ____)