Travailler avec les systèmes de coordonnées dans GeoPandas

Travailler avec des données géospatiales en Python

Joris Van den Bossche

Open source software developer and teacher, GeoPandas maintainer

Renseignements CRS dans GeoPandas

L'attribut .crs d'un GeoDataFrame/GeoSeries :

import geopandas
gdf = geopandas.read_file("countries.shp")
print(gdf.crs)

{'init': 'epsg:4326'}

gdf_noCRS = geopandas.read_file("countries_noCRS.shp")
print(gdf_noCRS.crs)

{}

Ajouter les renseignements CRS à crs :

# Option 1
gdf.crs = {'init': 'epsg:4326'}

# Option 2
gdf.crs = {'proj': 'longlat', 'datum': 'WGS84', 'no_defs': True}

import geopandas
gdf = geopandas.read_file("countries_web_mercator.shp")
print(gdf.crs)

{'init': 'epsg:3857', 'no_defs': True}

La méthode to_crs() :

# Option 1
gdf2 = gdf.to_crs({'proj': 'longlat', 'datum': 'WGS84', 'no_defs': True})

# Option 2
gdf2 = gdf.to_crs(epsg=4326)

1) Sources avec un CRS différent

df1 = geopandas.read_file(...)
df2 = geopandas.read_file(...)

df2 = df2.to_crs(df1.crs)

1) Sources avec un CRS différent

2) Cartographie (distorsion des formes et des distances)

1) Sources avec un CRS différent

2) Cartographie (distorsion des formes et des distances)

3) Calculs de distance/surface

Conseils :

Sites utiles :

Sites utiles :

Travailler avec des données géospatiales en Python