清华大佬耗费三个月吐血整理的几百G的资源，免费分享！....>>>

import math
wgs84={"a":6378137,"b":6356752.3142,}
def num2blh(num):
    '''
    23.434=>23:
    '''
    dd=int(num)
    mm=int((num-dd)*60)
    ss=(num*60-int(num*60))*60
    return ":".join([str(x) for x in (dd,mm,ss)])
def XYZ2BLH(x,y,z,epi):
    '''
    x,y,z:分别对应空间直角坐标系坐标
    epi:椭球长短轴
    '''
    x,y,z=float(x),float(y),float(z)
    a,b=epi["a"],epi["b"]
    e=math.sqrt(a**2-b**2)/a
    ep=math.sqrt(a**2-b**2)/b
    L=math.atan(y/x)*180/math.pi+180
    r=math.sqrt(x**2+y**2)
    R=math.sqrt(r**2+z**2)
    u=math.atan(b*z*(1+ep**2*b/R)/(a*r))
    B=math.atan((z+ep**2*b*math.sin(u)**3)/(r-e**2*a*math.cos(u)**3))
    BB=B*180/math.pi
    H=r*math.cos(B)+z*math.sin(B)-a*math.sqrt(1-e**2*math.sin(B)**2)
    return (BB,L,H)