Buffer în jurul unui punct

Să luăm în considerare orașul Goma, capitala provinciei Kivu de Nord din Congo, aflat aproape de granița cu Rwanda. Coordonatele sale sunt 1,66°S 29,22°E (variabila Point este deja furnizată în coordonate UTM sub numele goma).

Câte situri miniere se află în raza de 50 km față de Goma? Și ce suprafață din parcul național se încadrează în această zonă? Hai să determinăm asta folosind operațiunea de buffer. Reține că distanțele trebuie exprimate în unitatea CRS-ului (adică în metri, în acest caz).

Notă: dacă ai o Serie de valori booleene (de exemplu, ca rezultat al unei metode de relație spațială), poți calcula câte valori True există (adică câte geometrii au trecut verificarea) luând suma acelor valori booleene, deoarece valorile True și False vor fi tratate ca unu și zero.

Acest exercițiu face parte din cursul

Lucrul cu date geospațiale în Python

Vezi cursul

Instrucțiuni pentru exercițiu

Creează un buffer de 50 km în jurul variabilei goma și atribuie-l variabilei goma_buffer.
Verifică care dintre punctele din mining_sites se află în interiorul acestui buffer (și numește această variabilă mask). Numără siturile respective luând suma valorilor booleene din mask și afișează rezultatul.
Calculează intersecțiile parcurilor naționale cu bufferul. Calculează suma ariilor acestor intersecții și afișează rezultatul în km².

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# 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))

Editează și rulează codul