Извлечение информации из растрового слоя
Теперь извлечём информацию из растрового слоя на основе векторного файла. Эту возможность предоставляет пакет rasterstats. В данном упражнении мы определим тип растительности во всех местах добычи, получив ближайшее значение растрового пикселя для каждой точки набора данных о горнодобывающих объектах.
Подмножество набора данных о горнодобывающих объектах (mining_sites) уже загружено, GeoPandas и matplotlib импортированы.
Это упражнение является частью курса
Работа с геопространственными данными в Python
Инструкции к упражнению
- Импортируйте пакет
rasterstats. - Извлеките ближайшее растровое значение для местоположений горнодобывающих объектов (
mining_sites) с помощью функцииrasterstats.point_query()и сохраните результат в новый столбец'vegetation'. - Постройте график данных о горнодобывающих объектах, используя столбец
'vegetation'для окраски точек. Не забудьте добавить легенду.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Import the rasterstats package
____
# Extract the nearest value in the raster for all mining sites
vegetation_raster = "central_africa_vegetation_map_foraf.tif"
mining_sites['vegetation'] = ____(____, ____, interpolate='nearest')
print(mining_sites.head())
# Replace numeric vegation types codes with description
mining_sites['vegetation'] = mining_sites['vegetation'].replace(vegetation_types)
# Make a plot indicating the vegetation type
____
plt.show()