在Python中按纬度计算地球半径-复制公式

import math

def calculate_earth_radius(latitude):
    equatorial_radius = 6378.137  # 赤道半径，单位：千米
    polar_radius = 6356.7523  # 极半径，单位：千米

    latitude_rad = math.radians(latitude)  # 将纬度转换为弧度

    numerator = math.pow(math.pow(equatorial_radius, 2) * math.cos(latitude_rad), 2) + math.pow(math.pow(polar_radius, 2) * math.sin(latitude_rad), 2)
    denominator = math.pow(equatorial_radius * math.cos(latitude_rad), 2) + math.pow(polar_radius * math.sin(latitude_rad), 2)

    radius = math.sqrt(numerator / denominator)  # 地球半径，单位：千米

    return radius