Utiliser les systèmes de coordonnées avec GeoPandas

Travailler avec des données géospatiales en Python

Joris Van den Bossche

Open source software developer and teacher, GeoPandas maintainer

Informations 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 informations 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 (déformation des formes et des distances)

1) Sources avec un CRS différent

2) Cartographie (déformation des formes et des distances)

3) Calculs de distance/surface

Conseils :

Sites utiles :

Sites utiles :

Travailler avec des données géospatiales en Python