Shapely / Pyproj查找由纬度和经度创建的多边形的面积(以m^2为单位)

pip install shapely pyproj

from shapely.geometry import Polygon
from pyproj import Proj, transform

def calculate_polygon_area(latitudes, longitudes):
    # 创建多边形对象
    polygon = Polygon(zip(longitudes, latitudes))

    # 定义投影坐标系
    in_proj = Proj(init='epsg:4326')  # WGS84坐标系
    out_proj = Proj(init='epsg:3857')  # Web Mercator坐标系

    # 将经纬度坐标转换为投影坐标
    projected_polygon = transform(in_proj, out_proj, polygon)

    # 计算多边形的面积
    area = projected_polygon.area

    return area