用编程求平行圈曲率半径

import math as m

# 建立类
class Radius():
    # 创建析构函数
    def __init__(self):
        self.a = 6378140.0   # 某椭球体的长轴，根据自己的需要选择设置长轴
        self.alpha = 1/298.257  # 扁率

    def Curvature_Radius(self, B):

        ee = 2 * self.alpha - m.pow(self.alpha, 2)
        R = self.a * m.cos(m.radians(B)) / m.sqrt(1-ee*m.pow(m.sin(m.radians(B)), 2))
        # m.radians(B)将纬度B的角度转换为弧度

        return R

B = eval(input("请输入纬度B(例如：30.7)："))
# 创建对象
radius = Radius()
T = radius.Curvature_Radius(B)
print(T)

运行结果示例：

请输入纬度B(例如：30.7)：30
5528259.243013336

Process finished with exit code 0