Projeter un point
Au chapitre précédent, nous avons travaillé avec l’emplacement de la tour Eiffel. Une fois encore, nous vous avons fourni les coordonnées dans un système de coordonnées projetées, afin que vous puissiez, par exemple, calculer des distances. Revenons à ce monument emblématique et exprimons sa position en coordonnées géographiques : 48°51′29,6″N, 2°17′40,2″E. Ou, en décimal : latitude 48,8584 et longitude 2,2945.
Les objets géométriques Shapely n’intègrent pas la notion de SCR et ne peuvent donc pas être convertis directement vers un autre SCR. Nous allons donc utiliser GeoPandas pour transformer le point de la tour Eiffel vers un autre SCR. Nous placerons ce point unique dans une GeoSeries, utiliserons la méthode to_crs(), puis extrairons de nouveau le point.
GeoPandas est déjà importé.
Cet exercice fait partie du cours
Travailler avec des données géospatiales en Python
Instructions
- Créez un objet point Shapely avec les coordonnées de la tour Eiffel et assignez-le à une variable appelée
eiffel_tower. - Créez une GeoSeries (appelée
s_eiffel_tower) avec la tour Eiffel comme unique élément et indiquez le SCR EPSG:4326. - Convertissez
s_eiffel_toweren EPSG:2154 et nommez le résultats_eiffel_tower_projected
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Construct a Point object for the Eiffel Tower
from shapely.geometry import Point
eiffel_tower = ____
# Put the point in a GeoSeries with the correct CRS
s_eiffel_tower = geopandas.GeoSeries([____], crs={'init': '____'})
# Convert to other CRS
s_eiffel_tower_projected = s_eiffel_tower.____
# Print the projected point
print(s_eiffel_tower_projected)