Proyectar un punto
En el capítulo anterior trabajamos con la ubicación de la Torre Eiffel. De nuevo, te damos las coordenadas en un sistema de coordenadas proyectadas para que, por ejemplo, puedas calcular distancias. Volvamos a este monumento emblemático y expresemos su ubicación en coordenadas geográficas: 48°51′29.6″N, 2°17′40.2″E. O, en decimales: latitud 48.8584 y longitud 2.2945.
Los objetos de geometría de Shapely no tienen noción de CRS y, por tanto, no se pueden convertir directamente a otro CRS. Por eso vamos a usar GeoPandas para transformar el punto de la Torre Eiffel a un CRS alternativo. Pondremos el punto en una GeoSeries, usaremos el método to_crs() y volveremos a extraer el punto.
GeoPandas ya está importado.
Este ejercicio forma parte del curso
Trabajo con datos geoespaciales en Python
Instrucciones del ejercicio
- Crea un objeto punto de Shapely con las coordenadas de la Torre Eiffel y asígnalo a una variable llamada
eiffel_tower. - Crea una GeoSeries (llamada
s_eiffel_tower) con la Torre Eiffel como único elemento y especifica el CRS EPSG:4326. - Convierte
s_eiffel_towera EPSG:2154 y llama al resultados_eiffel_tower_projected.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# 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)