Attributs et méthodes de GeoSeries I

Visualisation de données géospatiales avec Python

Mary van Valkenburg

Data Science Program Manager, Nashville Software School

Attributs et méthodes de Shapely

# la colonne geometry est une GeoSeries
type(school_districts.geometry)

geopandas.geoseries.GeoSeries
  • GeoSeries.area - renvoie l'aire de chaque géométrie d'une GeoSeries
  • GeoSeries.centroid - renvoie le centroïde de chaque géométrie d'une GeoSeries
  • GeoSeries.distance(other) - renvoie la distance minimale à other
Visualisation de données géospatiales avec Python

GeoSeries.area

  • renvoie l'aire de chaque géométrie d'une GeoSeries
# aire du premier polygone dans districts 
print(districts.geometry[0].area)

325.78

polygone

Visualisation de données géospatiales avec Python

Aires des commissions scolaires

GeoSeries.area
# afficher les 4 premières lignes des commissions scolaires et le nombre total de lignes
print(school_districts.head(4))
print('There are ', school_districts.shape[0], ' school districts.' )

first_name  last_name    position  district  geometry
Sharon      Gentry       Member       1      (POLYGON ((-86.771 36.383...)))
Jill        Speering     Vice-Chair   3      (POLYGON ((-86.753 36.404...)))
Jo Ann      Brannon      Member       2      (POLYGON ((-86.766 36.083...)))
Anna        Shepherd     Chair        4      (POLYGON ((-86.580 36.209...)))
There are  9  school districts.
Visualisation de données géospatiales avec Python
# calculer l'aire de chaque commission scolaire
district_area = school_districts.area
# afficher les aires et le SCR utilisé
print(district_area.sort_values(ascending = False))
print(school_districts.crs)

0    0.036641
4    0.023030
8    0.015004
1    0.014205
3    0.014123
5    0.010704
2    0.008328
7    0.007813
6    0.006415
dtype: float64
epsg:4326
Visualisation de données géospatiales avec Python
# créer une copie de school_districts qui utilise EPSG:3857
school_districts_3857 = school_districts.to_crs(epsg = 3857)

# définir une variable pour m^2 en km^2 et obtenir l'aire en kilomètres carrés
sqm_to_sqkm = 10**6
district_area_km = school_districts_3857.area / sqkm_to_sqm
print(district_area_km.sort_values(ascending = False))
print(school_districts_3857.crs)

0    563.134380
4    353.232132
8    230.135653
1    218.369949
3    216.871511
5    164.137548
2    127.615396
7    119.742279
6     98.469632
dtype: float64
epsg:3857
Visualisation de données géospatiales avec Python

Passons à la pratique !

Visualisation de données géospatiales avec Python

Preparing Video For Download...