Extraer información de una capa ráster
Ahora vamos a extraer información de la capa ráster a partir de un archivo vectorial. Esta funcionalidad la ofrece el paquete rasterstats. En concreto, en este ejercicio determinaremos el tipo de vegetación en todos los sitios mineros obteniendo el valor del píxel ráster más cercano en cada punto del conjunto de datos de sitios mineros.
Ya tienes cargado un subconjunto del conjunto de datos de sitios mineros (mining_sites), y GeoPandas y matplotlib ya están importados.
Este ejercicio forma parte del curso
Trabajo con datos geoespaciales en Python
Instrucciones del ejercicio
- Importa el paquete
rasterstats. - Extrae el valor ráster más cercano para las ubicaciones de los sitios mineros (
mining_sites) usando la funciónrasterstats.point_query()y asígnalo a una nueva columna'vegetation'. - Crea una visualización de los datos de sitios mineros usando la columna
'vegetation'para colorear los puntos. Asegúrate de añadir una leyenda.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# 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()