两点之间的轴承

Question

我一直在使用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

我需要计算轴承,st center_x和center_y是枢轴.之后我使用geopy对gps坐标进行逆向工程:

latlon = VincentyDistance(miles=dist).destination(Point(lat1, lon1), bearing)

谁能指出我可能做错了什么？

Answer 1

谁能指出我可能做错了什么？

1. 没有显示您的"不一致或来自相对大的位移"结果的示例,也没有显示您的预期结果; 因此,回答者必须依靠猜测.

2. 没有说明您的输入(x,y等)的单位是什么,以及如何dist在destination计算中使用.我假设(在bearing2下面的计算中)正x以英里为东,正y以英里为北.如果您要修改问题以修复(1)和(2),这将有很大帮助.

3. 一种不太有利于让民众想要阅读它的编码风格......看看这个.

4. 在学校三角学中,从X轴(东)逆时针测量角度.在导航中,轴承从Y轴(北)顺时针测量.见下面的代码.有关轴承使用的示例,请点击此链接,向下滚动到"目的地点给定距离并从起点承载"部分,注意该示例正在谈论约96或97度的轴承,然后单击"查看地图" "你会注意到航向略微偏向东方(东方是90度).

码:

from math import degrees, atan2
    def gb(x, y, center_x, center_y):
        angle = degrees(atan2(y - center_y, x - center_x))
        bearing1 = (angle + 360) % 360
        bearing2 = (90 - angle) % 360
        print "gb: x=%2d y=%2d angle=%6.1f bearing1=%5.1f bearing2=%5.1f" % (x, y, angle, bearing1, bearing2)

    for pt in ((0, 1),(1,1),(1,0),(1,-1),(0,-1),(-1,-1),(-1, 0),(-1,1)):
        gb(pt[0], pt[1], 0, 0)

输出:

gb: x= 0 y= 1 angle=  90.0 bearing1= 90.0 bearing2=  0.0
gb: x= 1 y= 1 angle=  45.0 bearing1= 45.0 bearing2= 45.0
gb: x= 1 y= 0 angle=   0.0 bearing1=  0.0 bearing2= 90.0
gb: x= 1 y=-1 angle= -45.0 bearing1=315.0 bearing2=135.0
gb: x= 0 y=-1 angle= -90.0 bearing1=270.0 bearing2=180.0
gb: x=-1 y=-1 angle=-135.0 bearing1=225.0 bearing2=225.0
gb: x=-1 y= 0 angle= 180.0 bearing1=180.0 bearing2=270.0
gb: x=-1 y= 1 angle= 135.0 bearing1=135.0 bearing2=315.0