Surface en km carrés, géométrie en degrés décimaux
Dans cet exercice, vous repartez du GeoDataFrame council_districts et du DataFrame permits. Vous allez passer council_districts au système de coordonnées EPSG 3857 avant de créer une colonne area. Une fois la colonne area créée, vous repasserez le CRS en EPSG 4326 afin que la géométrie soit en degrés décimaux.
Cet exercice fait partie du cours
Visualiser des données géospatiales en Python
Instructions
- Modifiez le système de référence de coordonnées de
council_districtsen EPSG 3857, puis affichez à nouveaucrset les deux premières lignes. - Créez une colonne nommée
area. Divisez l’aire de chaque polygone parsqm_to_sqkmpour obtenir la surface en kilomètres carrés. - Remettez le système de référence de coordonnées de
council_districtsà EPSG 4326. Affichezcrset les deux premières lignes.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de 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())