2

如何在 Python 中计算第二个地理坐标点的基本方向？我需要知道两个不同位置的距离和方向。

示例：工作：坐标 41.4107628,2.1745004 家庭：坐标 41.4126728,2.1704725

from geopy.distance import vincenty
work = geopy.Point(41.4107628,2.1745004)
home = geopy.Point(41.4126728,2.1704725)
print(vincenty(home, work))
0.398011015257 km

我想知道第二点在哪个方向上，（例如：北方，西北等）相对于第一个点......

提前谢谢了

4

3 回答 3

4

使用我的 python pacakge geolib。

后

pip install geographiclib

做

$ python
Python 2.7.14 (default, Nov  2 2017, 18:42:05) 
>>> from geographiclib.geodesic import Geodesic
>>> geod = Geodesic.WGS84
>>> g = geod.Inverse(41.4107628,2.1745004, 41.4126728,2.1704725)
>>> print "The initial direction is {:.3f} degrees.".format(g['azi1'])
The initial direction is -57.792 degrees.

方向是从北顺时针测量的。所以 -57.8 度 = 302.2 度 = 西北偏西；请参阅指南针的要点。

于 2017-12-12T19:36:10.157 回答
1

为清楚起见进行编辑：这是一个近似值，可以说明大致方向。使用下面的链接获取简单的在线计算器。否则请遵循提供的其他答案。

坐标系中的第一个值代表北/南方向，第二个值代表东/西。简单的减法将提供一个大致的方向。例如，从 A 中减去 B，我们得到：

41.4126728 - 41.4107628 = 0.00191

2.1704725 - 2.1745004 = - 0.0040279

这表明要从 A 点到达 B 点，您需要沿北（正值）西（负值）方向行驶。使用三角函数可以找到精确的角度（将每个值视为直角三角形的边）。正如评论中提到的，三角函数不足以进行准确的计算。

您可能会觉得这个网站很有趣：https ://www.movable-type.co.uk/scripts/latlong.html

于 2017-12-05T17:45:57.320 回答
0

FWIW，我需要北、南、东、西等作为文本。这是我的（旧式程序员）代码：

import math

def calcBearing (lat1, long1, lat2, long2):
    dLon = (long2 - long1)
    x = math.cos(math.radians(lat2)) * math.sin(math.radians(dLon))
    y = math.cos(math.radians(lat1)) * math.sin(math.radians(lat2)) - math.sin(math.radians(lat1)) * math.cos(math.radians(lat2)) * math.cos(math.radians(dLon))
    bearing = math.atan2(x,y)   # use atan2 to determine the quadrant
    bearing = math.degrees(bearing)

    return bearing

def calcNSEW(lat1, long1, lat2, long2):
    points = ["north", "north east", "east", "south east", "south", "south west", "west", "north west"]
    bearing = calcBearing(lat1, long1, lat2, long2)
    bearing += 22.5
    bearing = bearing % 360
    bearing = int(bearing / 45) # values 0 to 7
    NSEW = points [bearing]

    return NSEW

# White house 38.8977° N, 77.0365° W
lat1 = 38.8976763
long1 = -77.0365298
# Lincoln memorial 38.8893° N, 77.0506° W
lat2 = 38.8893
long2 = -77.0506

points = calcNSEW(lat1, long1, lat2, long2)
print ("The Lincoln memorial is " + points + " of the White House")
print ("Actually bearing of 231.88 degrees")

print ("Output says: The Lincoln memorial is south west of the White House ")
于 2021-03-19T20:43:50.660 回答