Buffer alrededor de un punto
Considera la ciudad de Goma, capital de la provincia de Kivu del Norte en el Congo, cerca de la frontera con Ruanda. Sus coordenadas son 1,66°S 29,22°E (el Point ya está proporcionado en coordenadas UTM como la variable goma).
¿Cuántos yacimientos mineros están situados a menos de 50 km de Goma? ¿Y cuánta superficie de parque nacional? Vamos a determinarlo usando la operación de buffer. Recuerda que las distancias deben expresarse en la unidad del CRS (es decir, en metros en este caso).
Nota: si tienes una Serie booleana (por ejemplo, como resultado de un método de relación espacial), puedes calcular cuántos valores True hay (es decir, cuántas geometrías superaron la comprobación) sumando esos booleanos, porque en ese caso los valores True y False se consideran unos y ceros.
Este ejercicio forma parte del curso
Trabajo con datos geoespaciales en Python
Instrucciones del ejercicio
- Crea un buffer de 50 km alrededor de
gomay asígnalo a una variablegoma_buffer. - Comprueba cuáles de los puntos de
mining_sitesestán dentro de este buffer (llama a estomask). Cuenta esos sitios sumando la máscara booleana y muestra este resultado por pantalla. - Calcula las intersecciones de los parques nacionales con el buffer. Suma el área de esas intersecciones y muestra el resultado en km².
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# goma is a Point
print(type(goma))
# Create a buffer of 50km around Goma
goma_buffer = ____
# The buffer is a polygon
print(type(goma_buffer))
# Check how many sites are located within the buffer
mask = ____
print(____)
# Calculate the area of national park within the buffer
print(national_parks.____ / (1000**2))