Buffer okolo bodu
Vezmi v úvahu město Goma, hlavní město provincie Severní Kivu v Kongu, které leží poblíž hranice s Rwandou. Jeho souřadnice jsou 1,66° j. š. a 29,22° v. d. (proměnná goma je již poskytnuta jako Point v souřadnicích UTM).
Kolik těžebních lokalit se nachází do 50 km od Gomy? A jak velká část národního parku? Zjistíme to pomocí operace buffer. Nezapomeň, že vzdálenosti musí být vyjádřeny v jednotce daného CRS (v tomto případě v metrech).
Poznámka: máš-li booleovskou sérii (například jako výsledek metody prostorového vztahu), spočítáš počet hodnot True (tedy počet geometrií, které podmínku splnily) tak, že hodnoty sečteš – True a False se totiž v takovém případě interpretují jako jedničky a nuly.
Toto cvičení je součástí kurzu
Práce s geoprostorovými daty v Pythonu
Pokyny k cvičení
- Vytvoř buffer o poloměru 50 km okolo bodu
gomaa ulož ho do proměnnégoma_buffer. - Zkontroluj, které body z
mining_sitesleží uvnitř tohoto bufferu (výsledek ulož jakomask). Spočítej tyto lokality součtem booleovské sériemaska výsledek vypiš. - Vypočítej průniky národních parků s bufferem. Sečti plochy těchto průniků a výsledek vypiš v km².
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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))