python - 如何计算经纬度多边形的面积？

Question

我有以下方法检查点是否是多边形

from shapely import wkt

def validate_polygon(points):
    try:
        wkt.loads("POLYGON((%s))" % points)
    except Exception as ex:
        raise WrongRequestDataError("Incorrect points format. " + str(ex))

如何计算多边形的面积并在大于 400 平方英里时抛出异常？

我试过检查这个:

polygon = wkt.loads("POLYGON((%s))" % "34.093523 -118.274893,34.091414 -118.275887,34.092082 -118.278062,34.093867 -118.276609,34.093523 -118.274893")
print(polygon.area)

4.406979500001112e-06

但似乎答案不正确或者它使用什么类型的单位来表示这个值以及如何将其转换为平方英里或 km^2？

Answer 1

通常你会像这样计算面积:

from shapely.geometry import Polygon
points = [(34.093523, -118.274893), (34.091414, -118.275887), (34.092082, -118.278062), (34.093867, -118.276609), (34.093523, -118.274893)]
polygon = Polygon(points)
# the area in square degrees
area_sdeg = polygon.area

注意:对于以平方米为单位的面积计算，您必须使用 https://gist.github.com/robinkraft/c6de2f988c9d3f01af3c 中描述的投影