"Invalid projection" 计算多边形面积时

Question

我尝试计算以 km2 为单位并以 EPSG 投影的多边形面积：3857，但它无法识别我的 crs

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

我应该使用另一个吗？函数有错吗？

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

Answer 1

测地线计算需要地理数据。这是 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" 计算多边形面积时

"Invalid projection" when calculating area of a polygon

python

area

geode

pyproj