python作业怎么做求多边形面积_如何使用python计算地球表面多边形的面积？

也许有点晚了，但这里有一个不同的方法，使用吉拉德定理。它指出一个大圆圈的多边形的面积是多边形之间的角度之和减去(N-2)*pi的R**2倍，其中N是角的数目。

我认为这是值得发表的，因为它不依赖于任何其他的库，除了numpy，这是一个非常不同的方法。当然，这只适用于一个球体，所以当把它应用到地球上时会有一些不准确的地方。

首先，我定义了一个函数来计算从点1沿大圆到点2的方位角：import numpy as np

from numpy import cos, sin, arctan2

d2r = np.pi/180

def greatCircleBearing(lon1, lat1, lon2, lat2):

dLong = lon1 - lon2

s = cos(d2r*lat2)*sin(d2r*dLong)

c = cos(d2r*lat1)*sin(d2r*lat2) - sin(lat1*d2r)*cos(d2r*lat2)*cos(d2r*dLong)

return np.arctan2(s, c)

现在我可以用这个来找到角度，然后是面积(在下面，lons和lats当然应该被指定，并且它们应该是正确的顺序)。此外，还应指定球体的半径。)N = len(lons)

angles = np.empty(N)

for i in range(N):

phiB1, phiA, phiB2 = np.roll(lats, i)[:3]

LB1, LA, LB2 = np.roll(lons, i)[:3]

# calculate angle with north (eastward)

beta1 = greatCircleBearing(LA, phiA, LB1, phiB1)

beta2 = greatCircleBearing(LA, phiA, LB2, phiB2)

# calculate angle between the polygons and add to angle array

angles[i] = np.arccos(cos(-beta1)*cos(-beta2) + sin(-beta1)*sin(-beta2))

area = (sum(angles) - (N-2)*np.pi)*R**2

根据另一份回复中给出的科罗拉多坐标，地球半径6371公里，我得到的面积是268930758560.74808