Tampon autour d’un point
Considérez la ville de Goma, 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 à moins de 50 km de Goma ? Et quelle surface de parc national cela représente-t-il ? Déterminons-le à l’aide de l’opération de tampon. N’oubliez pas que les distances doivent être exprimées dans l’unité du SCR (c’est-à-dire en mètres dans ce cas).
Remarque : si vous avez une Series booléenne (par exemple résultat d’une méthode de relation spatiale), vous pouvez calculer combien de valeurs True vous avez (c.-à-d. combien de géométries satisfont le test) en faisant la somme de ces booléens, car dans ce cas les valeurs True et False sont respectivement interprétées comme 1 et 0.
Cet exercice fait partie du cours
Travailler avec des données géospatiales en Python
Instructions
- Créez un tampon de 50 km autour de
goma, et affectez-le à une variablegoma_buffer. - Vérifiez quels points de
mining_sitesse situent à l’intérieur de ce tampon (appelez celamask). Comptez ces sites en faisant la somme du masque booléen, puis affichez ce résultat. - Calculez les intersections des parcs nationaux avec le tampon. Faites la somme des surfaces de ces intersections et affichez le résultat en km².
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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))