python编写一个程序，使得根据输入的起点和终点坐标值计算出坐标方位角

在Python中，我们可以使用Haversine公式来计算两个经纬度之间的方位角。首先，你需要安装`geopy`库，它提供了一个方便的方式来处理地理位置相关的操作。如果你还没有安装，可以使用`pip install geopy`命令。

下面是一个简单的示例程序，用于计算两个点之间的方位角：

from math import radians, cos, sin, asin, sqrt
from geopy.distance import distance

def calculate_bearing(start_point, end_point):
    # 将坐标转换为弧度
    lat1, lon1 = map(radians, start_point.split(','))
    lat2, lon2 = map(radians, end_point.split(','))

    # Haversine formula to get the difference in longitude and latitude
    dlon = lon2 - lon1
    dlat = lat2 - lat1

    # Calculate bearing using arctangent
    bearing = atan2(dsin(dlat), dcot(lat1) * cos(lat2) - tan(lon2) * sin(lat2)) * 180 / pi

    # Adjust for compass direction (0-360 degrees)
    if bearing < 0:
        bearing += 360

    return bearing

# 测试函数
start_point = "40.7128, -74.0060"  # 纽约坐标
end_point = "51.5074, -0.1278"   # 伦敦坐标

bearing = calculate_bearing(start_point, end_point)
print(f"从{start_point}到{end_point}的方位角大约是 {bearing} 度.")