X = (N+H)cosBcosL
Y = (N+H)cosBsinL
Z = (N(1-e
2
)+H)sinB
获取大地坐标
B, L, H
计算第一辅助函数
W = sqrt(1-e
2
sin
2
B)
N = a/W
开始
已知
a, b, f, e
输出 X, Y, Z
开始
已知
a, b, f, e
计算大地经度
L = arctan(Y/X)
计算大地纬度
tanB
0
= Z/(sqrt(X
2
+Y
2
))
tanB
1
= 1/(sqrt(X
2
+Y
2
))*
(Z+ae
2
tanB
0
/sqrt(1+tan
2
B-e
2
tan
2
B
0
))
abs(tanB
1
-tanB
0
) <= 5e-10
获取空间直角坐标
X, Y, Z
否
tanB
0
= tanB
1
是
B = arctan(tanB
1
)
求出大地高
H = sqrt(X
2
+Y
2
)/cosB-N
输出 B, L, H