Given a GeoPandas's GeoDataFrame
containing a series of polygons, I would like to get the area in km sq of each feature in my list.
This is a pretty common problem, and the usual suggested solution in the past has been to use shapely
and pyproj
directly (e.g. here and here).
Is there a way to do this in pure GeoPandas?
Best Answer
If the crs of the GeoDataFrame is known (EPSG:4326 unit=degree, here), you don't need Shapely, nor pyproj in your script because GeoPandas uses them).
Now copy your GeoDataFrame and change the projection to a Cartesian system (EPSG:3857, unit= m as in the answer of ResMar)
Now the area in square kilometers
But the surfaces in the Mercator projection are not correct, so with other projection in meters.