给定两个经纬度，求方位算法

给定两个经纬度计算方位通常涉及到地理坐标系统中的方向判断，比如两点间的大致朝向。这可以使用Haversine公式或者更简单的卡西尼法则来进行计算。以下是基本步骤：

1. **角度转换**：将经度从弧度转换成角度。地球半径（例如6371公里）乘以弧度差得到平面上的距离。

2. **确定大圆距离**：如果两点位于北纬45度以上的范围，需要考虑地球是个椭球体，使用反余弦函数（如acos）找出两点之间的大圆航向角。

3. **确定偏航角度**：如果第一个点在东边，则航向角加180度；如果在南边，则航向角取负值。这是因为地球自转方向导致的结果。

4. **结果处理**：对于最终的方向，通常会转换为常见的四个象限（如东北、东南等），或者直接给出角度数值。

这里是一个简化版的Python示例，假设经纬度是以浮点数形式表示的：

from math import radians, cos, sin, asin, sqrt

def calculate_direction(lat1, lon1, lat2, lon2):
    R = 6371  # 地球平均半径，单位为公里
    
    dlat = radians(lat2 - lat1)
    dlon = radians(lon2 - lon1)
    
    a = sin(dlat / 2)**2 + cos(radians(lat1)) * cos(radians(lat2)) * sin(dlon / 2)**2
    c = 2 * asin(sqrt(a))

    direction_radians = atan2(sin(dlon) * cos(radians(lat2)), cos(radians(lat1)) * sin(radians(lat2)) - sin(radians(lat1)) * cos(radians(lat2)) * cos(dlon))
    
    if lon1 > lon2:  # 如果第一个点在东边
        direction_degrees = (direction_radians + pi) % (2*pi) - pi + 180
    else:
        direction_degrees = direction_radians
    
    return round(direction_degrees, 2)

# 示例
direction = calculate_direction(51.5074, -0.1278, 40.7128, -74.0060)  # 假设伦敦到纽约
print(f"两地之间的大致方向是: {direction}°")