Python：计算两个纬度/经度之间的方位角-解网

问：

我正在尝试计算两个纬度/经度之间的方位角。

我对函数/公式本身没有疑问，

提供：

def get_bearing(lat1, long1, lat2, long2):
    dLon = (long2 - long1)

    y = math.sin(dLon) * math.cos(lat2)
    x = math.cos(lat1) * math.sin(lat2) - math.sin(lat1) * math.cos(lat2) * math.cos(dLon)

    brng = math.atan2(y, x)

    brng = np.rad2deg(brng)

    return brng

问题是结果不是预期的。

该函数的预期用法返回（非常长）列表中两个纬度/经度对之间的方位角，即

    lat1 = path[int(len(path) * location / 1000)][0]
    lat2 = path[int(len(path) * location / 1000) + 1][0]
    lng1 = path[int(len(path) * location / 1000)][1]
    lng2 = path[int(len(path) * location / 1000) + 1][1]

然后，方位角结果会改变图的视图方向，其中方位角可以假定 [-180， 180] 范围内的值。理想情况下，结果应使得 lat1、lng1 和 lat2、lng2 之间形成的线在图中完全“垂直”（lat/lon 注释在图中切换），见下文

我希望有人能够从函数返回的轴承以及预期的轴承应该是什么来推断出问题。以下是几个实例：

Current Location: 30.07134 -97.23076
Next in path: 30.0709 -97.22907
Calculated Bearing: 88.39967863143139
Expected Bearing: ~-70.67

Current Location: 29.91581 -96.85068
Next in path: 29.91556 -96.85021
Calculated Bearing: 118.9170342272798
Expected Bearing: ~122.67

Current Location: 29.69419 -96.53487
Next in path: 29.69432 -96.53466
Calculated Bearing 141.0271357781952
Expected Bearing: ~56

Current Location: 29.77357 -96.07924
Next in path: 29.77349 -96.07876
Calculated Bearing 165.24612555483893
Expected Bearing: ~104

很乐意提供更多信息，提前感谢您的任何/所有帮助。

python 函数 latitude-longitude plotly-dash

你有错误的纬度和长度，请参阅其中的 pyproj4.github.io/pyproj/stable/api/geod.html，例如： >>> from pyproj import Geod >>> g = Geod（ellps='clrk66'） # 使用 Clarke 1866 椭圆体。>>> boston_lat = 42.+（15./60.）;boston_lon = -71.-（7./60.） >>> portland_lat = 45.+（31./60.）;portland_lon = -123.-（41./60.） >>> # 计算前后方位角，加上波士顿和波特兰之间的距离 >>> #。>>> az12，az21，dist = g.inv（boston_lon，boston_lat，portland_lon，portland_lat） >>> “%7.3f %6.3f %12.3f” % （az12，az21，dist） '-66.531 75.654 4164192.708'

10赞 Sterling Butters 2/26/2019 #2

最终更改了功能：

from geographiclib.geodesic import Geodesic
...
def get_bearing(lat1, lat2, long1, long2):
    brng = Geodesic.WGS84.Inverse(lat1, long1, lat2, long2)['azi1']
    return brng

import numpy
import math

def get_bearing(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))
    brng = numpy.arctan2(x,y)
    brng = numpy.degrees(brng)

    return brng

Python：计算两个纬度/经度之间的方位角

Python: Calculate bearing between two lat/long

评论

评论

评论

评论