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