给定两点在特定时间找到gps位置？

Question

如果我有两个已知位置和一个已知速度，我如何计算距离 d（以公里为单位）的当前位置？

例如，给定：

ES4236 中的两个 gps 位置：

37.783333, -122.416667   # San Francisco
32.715, -117.1625        # San Diego

以1km/min直线行驶（忽略高度）

如何在一定距离找到gps坐标？类似的SO question使用geopy中的VincentyDistance根据方位和距离计算下一个点。

我猜，更具体地说：

• 如何使用geopy 计算两个gps 点之间的方位角？

• 使用 VincentyDistance 通过方位和距离获取下一个 gps 点，我如何知道我是否已到达目的地，或者是否应该继续前进？它不需要完全在目的地上才能被视为到达。也许任何距离目的地 0.5 公里半径的点都被认为是“到达”。

即，

import geopy

POS1 = (37.783333, -122.416667)     # origin
POS2 = (32.715, -117.1625)          # dest

def get_current_position(d):
    # use geopy to calculate bearing between POS1 and POS2
    # then use VincentyDistance to get next coord

    return gps_coord_at_distance_d

# If current position is within .5 km of destination, consider it 'arrived'
def has_arrived(curr_pos):
    return True/False

d = 50     # 50 km
print get_current_position(d)
print has_arrived(get_current_position(d))

Answer 1

好的，我想我会回到这个问题，并尽我所能，因为它还没有看到任何其他解决方案。不幸的是，我现在无法测试代码，但我相信使用 geopy 和geographiclib 可以解决您的问题。来了。

从终端（可能使用 sudo）

pip install geographiclib
pip install geopy

现在使用 Python

获取当前位置

import geographiclib
from geopy import geopy.distance

# Get the first azimuth, which should be the bearing
bearing = geographiclib.WGS84.Inverse(37.783333, -122.416667, 32.715, -117.1625)[2]


# Now we use geopy to calculate the distance over time
dist = geopy.distance.VincentyDistance(kilometers = 1)
san_fran = geopy.Point(37.783333, -122.416667)

print dist.destination(point=san_fran, bearing=bearing)

到货了

def has_arrived(d):
    return geopy.distance.vincenty(curr_pos, (32.715, -117.1625)).kilometers < .5

就像我说的，很遗憾我无法对此进行测试，但我相信这是正确的。轴承计算可能会有一些单位差异：it calculates bearing off of North as seen here。抱歉，如果这不完全正确，但就像我说的那样，因为我想我不妨把我知道的东西扔掉，所以还没有收到回复。