Oppervlakte in km², geometrie in decimale graden

In deze oefening begin je opnieuw met de GeoDataFrame council_districts en de DataFrame permits. Je zet council_districts eerst om naar het coördinatenreferentiesysteem EPSG 3857 voordat je een kolom area aanmaakt. Zodra de kolom area is aangemaakt, zet je het CRS terug naar EPSG 4326 zodat de geometrie in decimale graden staat.

Deze oefening maakt deel uit van de cursus

Geospatiale data visualiseren in Python

Bekijk cursus

Oefeninstructies

Wijzig het coördinatenreferentiesysteem van council_districts naar EPSG 3857 en print opnieuw de crs en de eerste twee rijen.
Maak een kolom area. Deel de oppervlakte van elke polygoon door sqm_to_sqkm om de oppervlakte in vierkante kilometers te krijgen.
Zet het coördinatenreferentiesysteem van council_districts terug naar EPSG 4326. Print de crs en de eerste twee rijen.

Interactieve oefening met praktijkervaring

Probeer deze oefening door deze voorbeeldcode aan te vullen.

# Change council_districts crs to epsg 3857
council_districts = council_districts.____(____ = ____)
print(council_districts.crs)
print(council_districts.head())

# Create area in square km
sqm_to_sqkm = 10**6
council_districts['area'] = council_districts.____ / sqm_to_sqkm

# Change council_districts crs back to epsg 4326
council_districts = council_districts.____(____ = ____)
print(council_districts.crs)
print(council_districts.head())

Code bewerken en uitvoeren