Buffer attorno a un punto
Considera la città di Goma, capoluogo della provincia del Nord Kivu, in Congo, vicino al confine con il Ruanda. Le sue coordinate sono 1,66°S 29,22°E (il Point è già fornito in coordinate UTM come variabile goma).
Quanti siti minerari si trovano entro 50 km da Goma? E quanta area di parco nazionale cade in questo raggio? Determiniamolo usando l’operazione di buffer. Ricorda che le distanze devono essere espresse nell’unità del CRS (cioè in metri in questo caso).
Nota: se hai una Series booleana (per esempio come risultato di un metodo di relazione spaziale), puoi calcolare quanti valori True ci sono (cioè quante geometrie hanno superato il controllo) facendo la somma di quei booleani, perché in questo caso i valori True e False sono interpretati come uno e zero.
Questo esercizio fa parte del corso
Lavorare con i dati geospaziali in Python
Istruzioni dell'esercizio
- Crea un buffer di 50 km attorno a
gomae assegnalo alla variabilegoma_buffer. - Verifica quali dei punti in
mining_sitessi trovano all’interno di questo buffer (chiama questo risultatomask). Conta questi siti facendo la somma dellamaskbooleana e stampa il risultato. - Calcola le intersezioni dei parchi nazionali con il buffer. Somma l’area di queste intersezioni e stampa il risultato in km².
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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))