Projeter un point
Dans le chapitre précédent, nous avons travaillé avec l'emplacement de la tour Eiffel. Nous vous avons de nouveau 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 son emplacement en coordonnées géographiques : 48°51′29,6″N, 2°17′40,2″E. Ou, en décimales : latitude de 48,8584 et longitude de 2,2945.
Les objets géométriques de Shapely n'ont pas de notion de SCR et ne peuvent donc pas être convertis directement vers un autre SCR. Nous allons donc utiliser GeoPandas pour transformer le point représentant 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é.
Cette activité fait partie du cours
Travailler avec des données géospatiales en Python
Instructions de l’exercice
- 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 seul élément et précisez le SCR EPSG:4326. - Convertissez
s_eiffel_towervers EPSG:2154 et nommez le résultats_eiffel_tower_projected
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# 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)