Projekce bodu
V předchozí kapitole jsme pracovali s polohou Eiffelovy věže. Tehdy jsme ti souřadnice poskytli v projected souřadnicovém systému, aby sis mohl/a například počítat vzdálenosti. Vraťme se teď k této ikonické stavbě a vyjádřeme její polohu v geografických souřadnicích: 48°51′29,6″ s. š., 2°17′40,2″ v. d. Nebo v desetinném zápisu: zeměpisná šířka 48,8584 a zeměpisná délka 2,2945.
Geometrické objekty Shapely nemají žádnou informaci o CRS, a proto je nelze přímo převádět do jiného souřadnicového systému. Použijeme tedy GeoPandas, abychom transformovali polohu Eiffelovy věže do alternativního CRS. Vložíme bod do GeoSeries, použijeme metodu to_crs() a bod pak z výsledku opět získáme.
GeoPandas je již importovaný.
Toto cvičení je součástí kurzu
Práce s geoprostorovými daty v Pythonu
Pokyny k cvičení
- Vytvoř shapely objekt typu bod se souřadnicemi Eiffelovy věže a přiřaď ho do proměnné
eiffel_tower. - Vytvoř GeoSeries (s názvem
s_eiffel_tower), která bude obsahovat Eiffelovu věž jako jediný prvek, a nastav CRS na EPSG:4326. - Převeď
s_eiffel_towerdo EPSG:2154 a výsledek ulož do proměnnés_eiffel_tower_projected.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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)