Extrair informações da camada raster
Agora vamos extrair informações da camada raster com base em um arquivo vetorial. Essa funcionalidade é fornecida pelo pacote rasterstats. Especificamente neste exercício, vamos determinar o tipo de vegetação em todos os locais de mineração, obtendo o valor do pixel raster mais próximo em cada ponto do conjunto de dados de locais de mineração.
Um subconjunto do conjunto de dados de locais de mineração (mining_sites) já está carregado, e o GeoPandas e o matplotlib já foram importados.
Este exercício faz parte do curso
Trabalhando com Dados Geoespaciais em Python
Instruções do exercício
- Importe o pacote
rasterstats. - Extraia o valor raster mais próximo para as localizações dos locais de mineração (
mining_sites) usando a funçãorasterstats.point_query()e atribua o resultado a uma nova coluna'vegetation'. - Faça um gráfico dos dados de locais de mineração usando a coluna
'vegetation'para colorir os pontos. Não esqueça de adicionar uma legenda.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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()