地球表面位置解算python代码实现，要求正反解算两地的球面距离和方位角

在Python中，我们可以使用地理信息系统的库，如`geopy`或`pyproj`，来进行地球表面两点间的球面距离和方位角计算。这里以`geopy.distance`为例，展示如何实现正反解算：

首先，确保已经安装了`geopy`库，如果没有，可以使用`pip install geopy`命令进行安装。

**正解算（计算两点之间的球面距离）**:

from geopy.distance import great_circle

def calculate_distance(lat1, lon1, lat2, lon2):
    # 使用经纬度坐标
    point1 = (lat1, lon1)
    point2 = (lat2, lon2)

    # 计算大圆距离（弧度）
    distance_in_radians = great_circle(point1, point2).kilometers
    distance_km = round(distance_in_radians * 6371)  # 地球平均半径约为6371公里

    return distance_km

# 示例
distance = calculate_distance(51.5074, -0.1278, 40.7128, -74.0060)  # 北京到纽约的距离

**反解算（计算方位角，通常指两点间的大圆航向角度）**:

由于`geopy`库本身并不直接提供方位角的功能，我们需要借助其他数学工具。假设我们已经有了一个辅助函数`calculate_bearing()`来计算两点间的方位角：

import math

def calculate_bearing(lat1, lon1, lat2, lon2):
    """
    计算两个点之间的大圆航向角
    """
    dLon = math.radians(lon2 - lon1)
    y = math.sin(dLon) * math.cos(math.radians(lat2))
    x = math.cos(math.radians(lat1)) * math.sin(math.radians(lat2)) - math.sin(math.radians(lat1)) * math.cos(math.radians(lat2)) * math.cos(dLon)
    bearing = math.degrees(math.atan2(y, x))

    # 如果结果在南、北半球相反，加上180°或减去180°
    if bearing < 0:
        bearing += 360

    return bearing

# 示例
bearing = calculate_bearing(51.5074, -0.1278, 40.7128, -74.0060)  # 北京到纽约的方位角