我一直在使用geopy包,它做得很好,但是我得到的一些结果不一致或者有一个相对大的位移,我怀疑问题在于我的轴承计算:
def gb(x,y,center_x,center_y):
dx=x-center_x
dy=y-center_y
if ((dy>=0)and((dx>0)or(dx<0))):
return math.degrees(math.atan2(dy,dx))
elif (dy<=0)and((dx>0)or (dx<0)):
return (math.degrees(math.atan2(dy,dx))+360)
else:
return (math.degrees(math.atan2(dy,dx))+360)%360
Run Code Online (Sandbox Code Playgroud)
我需要计算轴承,st center_x和center_y是枢轴.之后我使用geopy对gps坐标进行逆向工程:
latlon = VincentyDistance(miles=dist).destination(Point(lat1, lon1), bearing)
Run Code Online (Sandbox Code Playgroud)
谁能指出我可能做错了什么?