pythongps距离_关于python：在给定点，方位角和距离的情况下计算gps坐标

我有一个问题，现在我在某个项目上退缩了一段时间。

我基本上是想用一些脚本编写的x，y点捕获一个多边形。

lat1，lon1是多边形的中心gps线，我正在寻找其周围的多边形。

这是我在python中的代码的一部分：

def getcords(lat1,lon1,dr,bearing):

lat2=asin(sin(lat1)*cos(dr)+cos(lat1)*sin(dr)*cos(bearing))

lon2=lon1+atan2(sin(bearing)*sin(dr)*cos(lat1),cos(dr)-sin(lat1)*sin(lat2))

return [lat2,lon2]

我的输入是这样的：

lat1，lon1-以十进制度给出。

-dr是通过以英里为单位的距离除以地球的-raiuds(= 3958.82)而计算出的角度

-介于0-360度之间。

但是对于输入getcorsds1(42.189275，-76.85823,0.5 / 3958.82,30)

我得到[-1.3485899508698462，-76.8576637627568]，但是[42.2516666666667，-76.8097222222222]是正确的答案。

至于角距离，我只需将距离(以英里为单位)除以地球的突袭来计算即可(= 3958.82)。

有人吗

您为什么不使用漂亮的库？

from geopy import Point

from geopy.distance import distance, VincentyDistance

# given: