Bufor wokół punktu

Przyjrzyjmy się miastu Goma, stolicy prowincji Kiwu Północne w Kongu, położonej blisko granicy z Rwandą. Jej współrzędne to 1,66°S 29,22°E (zmienna goma zawiera już punkt Point w układzie współrzędnych UTM).

Ile miejsc wydobycia znajduje się w odległości do 50 km od Gomy? Jak dużą powierzchnię parków narodowych obejmuje ten obszar? Odpowiemy na te pytania za pomocą operacji bufora. Pamiętaj, że odległości należy podawać w jednostkach układu CRS (w tym przypadku – w metrach).

Uwaga: jeśli masz serię wartości logicznych (np. jako wynik metody sprawdzającej relację przestrzenną), możesz obliczyć liczbę wartości True (czyli geometrii, które spełniły warunek), sumując te wartości – True i False są traktowane odpowiednio jako jedynki i zera.

To ćwiczenie jest częścią kursu

Praca z danymi geoprzestrzennymi w Pythonie

Zobacz kurs

Instrukcje do ćwiczenia

Utwórz bufor o promieniu 50 km wokół punktu goma i przypisz go do zmiennej goma_buffer.
Sprawdź, które punkty ze zbioru mining_sites znajdują się wewnątrz tego bufora (wynik zapisz jako mask). Policz te miejsca, sumując wartości logiczne ze zmiennej mask, i wyświetl wynik.
Oblicz części wspólne parków narodowych z buforem. Zsumuj powierzchnie tych części wspólnych i wyświetl wynik w km².

Interaktywne ćwiczenie praktyczne

Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.

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

Edytuj i uruchom kod