5

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

谁能指出我做错了什么？

4

2 回答 2

9

谁能指出我做错了什么？

  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
于 2011-02-20T21:20:12.003 回答
2

我不太确定你想在你的代码中做什么，但我确实看到了一些可能需要清理的奇怪之处。

  1. 您在条件之前测试dy<=0之后进行测试。dy>=0如果dy==0和dx==0. 你的代码应该怎么做？
  2. 您的测试((dy>=0)and((dx>0)or(dx<0)))等效于 (dy>=0 and dx!=0)，这是您的意图吗？
  3. 您在所有条件句中基本上都在做同样的事情。不能return math.degrees(math.atan2(dy,dx))+360)%360在所有情况下工作？在这种情况下，无论如何您都不需要使用 if 语句。
于 2011-02-20T17:33:43.520 回答