НачатьНачать бесплатно

Пересечение GeoDataFrame с полигоном

Объединив набор данных о землепользовании и набор данных о районах, можно исследовать, как используются земли в конкретном районе.

Для этого сначала нужно найти пересечение набора данных о землепользовании с заданным районом. Возьмём в качестве примера район Muette.

Наборы данных о землепользовании и районах уже импортированы как land_use и districts, а район Muette извлечён в полигон muette библиотеки Shapely. GeoPandas и matplotlib также импортированы.

Это упражнение является частью курса

Работа с геопространственными данными в Python

Посмотреть курс

Инструкции к упражнению

  • Вычислите пересечение полигонов land_use с отдельным полигоном muette. Сохраните результат в переменную land_use_muette.
  • Постройте быстрый график этого пересечения, передав параметр edgecolor='black', чтобы границы отдельных полигонов были чётко видны.
  • Выведите первые пять строк land_use_muette.

Интерактивное практическое упражнение

Попробуйте выполнить это упражнение, дополнив этот пример кода.

# 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.____)
Редактировать и запускать код