Area in km squared, geometry in decimal degrees

In this exercise, you'll start again with the council_districts GeoDataFrame and the permits DataFrame. You will change the council_districts to use the EPSG 3857 coordinate reference system before creating a column for area. Once the area column has been created, you will change the CRS back to EPSG 4326 so that the geometry is in decimal degrees.

This exercise is part of the course

Visualizing Geospatial Data in Python

View Course

Exercise instructions

Change the coordinate reference system for the council_districts to EPSG 3857, and print the crs and first two rows again.
Create a column called area. Divide the area of each polygon by sqm_to_sqkm to get the area in kilometers squared.
Change the coordinate reference system for the council_districts back to EPSG 4326. Print the crs and first two rows.

Hands-on interactive exercise

Have a go at this exercise by completing this sample code.

# 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())

Edit and Run Code