大地坐标转换为空间直角坐标

import math
print（‘大地纬度按度分秒输入‘）
d1 = float（input（‘d1 = ‘））
f1 = float（input（‘f1 = ‘））
m1 = float（input（‘m1 = ‘））
print（‘大地经度按度分秒输入‘）
d2 = float（input（‘d2 = ‘））
f2 = float（input（‘f2 = ‘））
m2 = float（input（‘m2 = ‘））
H = float（input（‘H = ‘））
a = 6378137
b = 6356752.3142
e = pow（a*a - b*b1/2）/a
e1 = pow（a*a