Tampon autour d'un point
Considérez la ville de Goma, la capitale de la province du Nord-Kivu en République démocratique du Congo, près de la frontière avec le Rwanda. Ses coordonnées sont 1,66°S 29,22°E (le Point est déjà fourni en coordonnées UTM dans la variable goma).
Combien de sites miniers se trouvent dans un rayon de 50 km autour de Goma? Et quelle superficie de parc national? Déterminons cela à l'aide de l'opération de tampon. N'oubliez pas que les distances doivent être exprimées dans l'unité du SCR (c.-à-d. en mètre dans ce cas-ci).
Remarque : si vous avez une série booléenne (par exemple le résultat d'une méthode de relation spatiale), vous pouvez calculer combien de valeurs True (c.-à-d. combien de géométries ont satisfait la condition) en faisant la somme de ces booléens, puisque dans ce cas les valeurs True et False sont respectivement interprétées comme des uns et des zéros.
Cette activité fait partie du cours
Travailler avec des données géospatiales en Python
Instructions de l’exercice
- Créez un tampon de 50 km autour de
gomaet assignez-le à une variablegoma_buffer. - Vérifiez lesquels des points de
mining_sitesse trouvent à l'intérieur de ce tampon (appelez cecimask). Comptez ces sites en faisant la somme du booléenmask, puis affichez ce résultat. - Calculez les intersections des parcs nationaux avec le tampon. Faites la somme des superficies de ces intersections et affichez le résultat en km².
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# 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))