Buffer rond een punt

Kijk naar de stad Goma, de hoofdstad van de provincie Noord-Kivu in Congo, dicht bij de grens met Rwanda. De coördinaten zijn 1,66°Z 29,22°O (het Point is al in UTM-coördinaten voorzien als de variabele goma).

Hoeveel mijnlocaties liggen binnen 50 km van Goma? En hoeveel oppervlakte van nationaal park? Dat bepalen we met de bufferoperatie. Onthoud dat afstanden moeten worden uitgedrukt in de eenheid van het CRS (dus in dit geval in meter).

Opmerking: als je een booleaanse Series hebt (bijvoorbeeld als resultaat van een methode voor ruimtelijke relaties), kun je tellen hoeveel True-waarden er zijn (dus hoeveel geometrieën aan de voorwaarde voldoen) door de som te nemen van die booleans, omdat True en False dan als enen en nullen worden gezien.

Deze oefening maakt deel uit van de cursus

Werken met georuimtelijke data in Python

Cursus bekijken

Oefeninstructies

Maak een buffer van 50 km rond goma en wijs die toe aan een variabele goma_buffer.
Controleer welke punten van mining_sites binnen deze buffer liggen (noem dit mask). Tel die locaties door de som te nemen van de booleaanse mask, en print dit resultaat.
Bereken de intersecties van de nationale parken met de buffer. Neem de som van de oppervlakte van die intersecties en print het resultaat in km².

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

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

Code bewerken en uitvoeren