"Invalid projection" when calculating area of a polygon

Question

I try to calculate the area of a polygon in km2 and projected in EPSG: 3857, but it doesn't recognize my crs

raise GeodError("Invalid geometry provided.")
pyproj.exceptions.GeodError: Invalid geometry provided.

should I use another one? is the function wrong?

def area(polygon):
        geod = Geod('EPSG: 3857')
        x, y = polygon.exterior.coords.xy
        area, perimeter = geod.geometry_area_perimeter(x,y)
        return area

What is the CRS of the points in the polygon? – snowman2 Oct 23 '21 at 18:47 — snowman2, Oct 23 '21 at 18:47

score 0 · Accepted Answer · answered Oct 25 '21 at 15:08

Geodesic calculations expect geographic data. Here is an example for your data in WGS84/EPSG:4326:

https://pyproj4.github.io/pyproj/stable/examples.html#creating-geod-class

from pyproj import Geod, CRS
geod = Geod(ellps='WGS84')
# or
geod = CRS("EPSG:4326").get_geod()

"Invalid projection" when calculating area of a polygon

1 Answers1