Superficie 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 d'abord passer council_districts au système de référence de coordonnées EPSG 3857 avant de créer une colonne area. Une fois la colonne area créée, vous remettrez le SRC à EPSG 4326 afin que la géométrie soit en degrés décimaux.
Cette activité fait partie du cours
Visualisation de données géospatiales avec Python
Instructions de l’exercice
- Changez le système de référence de coordonnées de
council_districtsà EPSG 3857 et affichez de nouveaucrsainsi que les deux premières lignes. - Créez une colonne appelée
area. Divisez l'aire de chaque polygone parsqm_to_sqkmpour obtenir la superficie 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 ce code d’exemple.
# 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())