@@ -684,6 +684,9 @@ l'oeil.
684
684
avec ` cartiflette ` . Pour cela, utiliser ` download_vectorfile_url_all `
685
685
depuis ` cartiflette.s3 ` en fixant l'option ` level ` à ` COMMUNE_ARRONDISSEMENT ` .
686
686
Nommer cet objet ` df ` .
687
+ 2 . Afin que les calculs ultérieurs de surface ne soient pas faussés,
688
+ assurez-vous que les données sont en Lambert 93 en reprojetant
689
+ nos contours (code EPSG: 2154).
687
690
2 . Créer un objet ` departements ` avec ` dissolve ` pour également disposer
688
691
d'un fond de carte des départements
689
692
3 . Créer une variable ` surface ` et utilisant la méthode ` area ` . L'unité
@@ -712,16 +715,18 @@ df = download_vectorfile_url_all(
712
715
values = ["75", "92", "93", "94"],
713
716
level="COMMUNE_ARRONDISSEMENT",
714
717
decoupage="departement")
715
- # 2/ Départements
718
+ # 2/ Lambert 93
719
+ df = df.to_crs(2154)
720
+ # 3/ Départements
716
721
departements = df.dissolve("INSEE_DEP")
717
722
```
718
723
719
724
``` {python}
720
725
#| echo: false
721
- # 3 / et 4 /
726
+ # 4 / et 5 /
722
727
df['surface'] = df.area.div(10**6)
723
728
df['densite'] = (df['POPULATION']/df['surface'])
724
- #5 /
729
+ #6 /
725
730
df['markersize'] = 8*df['densite'].div(df.area.div(10**6).sum())
726
731
df['markercolor'] = pd.cut(
727
732
df['densite'],
@@ -730,7 +735,7 @@ df['markercolor'] = pd.cut(
730
735
731
736
``` {python}
732
737
#| echo: false
733
- # 6 /
738
+ # 7 /
734
739
df_points = df.copy()
735
740
df_points["geometry"] = df_points["geometry"].centroid
736
741
```
@@ -739,7 +744,7 @@ La carte obtenue devrait ressembler à celle-ci:
739
744
740
745
``` {python}
741
746
#| echo: false
742
- # 7 /
747
+ # 8 /
743
748
ax = df.plot(
744
749
color="lightgray", edgecolor="grey",
745
750
figsize=(7,7),
0 commit comments