Een punt projecteren
In het vorige hoofdstuk werkten we met de locatie van de Eiffeltoren. We gaven je de coördinaten in een geprojecteerd coördinatensysteem, zodat je bijvoorbeeld afstanden kon berekenen. Laten we terugkeren naar dit iconische monument en de locatie uitdrukken in geografische coördinaten: 48°51′29.6″N, 2°17′40.2″E. Of, in decimalen: een breedtegraad van 48,8584 en een lengtegraad van 2,2945.
Shapely-geometrieobjecten hebben geen notie van een CRS en kunnen dus niet rechtstreeks naar een andere CRS worden geconverteerd. Daarom gaan we GeoPandas gebruiken om het punt van de Eiffeltoren naar een alternatieve CRS te transformeren. We zetten het enkele punt in een GeoSeries, gebruiken de methode to_crs() en halen het punt daarna weer op.
GeoPandas is al geïmporteerd.
Deze oefening maakt deel uit van de cursus
Werken met georuimtelijke data in Python
Oefeninstructies
- Maak een shapely Point-object met de coördinaten van de Eiffeltoren en sla dit op in een variabele
eiffel_tower. - Maak een GeoSeries (genaamd
s_eiffel_tower) met de Eiffeltoren als enig element en specificeer de CRS als EPSG:4326. - Converteer
s_eiffel_towernaar EPSG:2154 en noem het resultaats_eiffel_tower_projected
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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)