Area i km², geometri i decimalgrader
I den här övningen börjar du igen med GeoDataFrame:en council_districts och DataFrame:en permits. Du ändrar council_districts till koordinatreferenssystemet EPSG 3857 innan du skapar en kolumn för area. När area-kolumnen är skapad ändrar du CRS tillbaka till EPSG 4326 så att geometrin uttrycks i decimalgrader.
Den här övningen är en del av kursen
Visualisering av geospatial data i Python
Övningsinstruktioner
- Ändra koordinatreferenssystemet för
council_districtstill EPSG 3857 och skriv utcrssamt de två första raderna igen. - Skapa en kolumn som heter
area. Dela arean för varje polygon medsqm_to_sqkmför att få arean i kvadratkilometer. - Ändra koordinatreferenssystemet för
council_districtstillbaka till EPSG 4326. Skriv utcrsoch de två första raderna.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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())