Function calc_central_angle

backend/orbit.py:43–51

计算球面大圆地心角 (Degrees)

(lat1, lon1, lat2, lon2)

def calc_central_angle(lat1, lon1, lat2, lon2):
    """计算球面大圆地心角 (Degrees)"""
    rad = math.pi / 180.0
    phi1, phi2 = lat1 * rad, lat2 * rad
    dphi = (lat2 - lat1) * rad
    dlam = (lon2 - lon1) * rad
    a = math.sin(dphi / 2) ** 2 + math.cos(phi1) * math.cos(phi2) * math.sin(dlam / 2) ** 2
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))
    return c * 180.0 / math.pi


class OrbitalElements: