Proiectarea unui punct
În capitolul anterior, am lucrat cu locația Turnului Eiffel. Ți-am furnizat coordonatele într-un sistem de coordonate proiectat, pentru a putea calcula, de exemplu, distanțe. Să revenim la acest reper iconic și să exprimăm locația sa în coordonate geografice: 48°51′29,6″N, 2°17′40,2″E. Sau, în format zecimal: latitudine 48,8584 și longitudine 2,2945.
Obiectele de tip geometrie Shapely nu au noțiunea unui CRS și, prin urmare, nu pot fi convertite direct la un alt CRS. De aceea, vom folosi GeoPandas pentru a transforma locația punctului Turnului Eiffel într-un CRS alternativ. Vom plasa punctul într-un GeoSeries, vom aplica metoda to_crs() și vom extrage punctul din nou.
GeoPandas este deja importat.
Acest exercițiu face parte din cursul
Lucrul cu date geospațiale în Python
Instrucțiuni pentru exercițiu
- Creează un obiect de tip punct Shapely cu coordonatele Turnului Eiffel și atribuie-l unei variabile numite
eiffel_tower. - Creează un GeoSeries (numit
s_eiffel_tower) cu Turnul Eiffel ca unic element și specifică CRS-ul ca EPSG:4326. - Convertește
s_eiffel_towerla EPSG:2154 și numește rezultatuls_eiffel_tower_projected.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# 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)