同椭球下空间直角坐标转换为大地坐标，椭球大地,import mathw

同椭球下空间直角坐标转换为大地坐标，椭球大地,import mathw

import mathwgs84={"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)