ПочатиПочніть безкоштовно

Добування інформації з растрового шару

Тепер добудьмо інформацію з растрового шару на основі векторного файла. Цю можливість надає пакет 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()
Редагувати та запускати код