計算多邊形面積時出現“無效投影”

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()