#?-*-?coding:?utf-8?-*-
“““
Created?on?Thu?May?23?15:22:07?2019

@author:?72987
“““
import?math

dFlattenting=0.00673949674227
dSemiMajorAxis=6378137.0

def?CARTtoGEO（stationpos）:
????
????x=stationpos[0]
????y=stationpos[1]
????z=stationpos[2]
????BLH=[]
????
????dFlattenting=0.00673949674227
????dSemiMajorAxis=6378137.0
????R=math.sqrt（x**2+y**2+z**2）
????e2=dFlattenting*（2-dFlattenting）
????L=math.atan2（yx）
????B0=B=math.atan2（zmath.sqrt（x**2+y**2））
????
????for?i?in?range（int（1e+6））:
????????B=B0
????????N=dSemiMajorAxis/math.sqrt（1-e2*math.sin（B0）**2）
????????sinB=math.sin（B0）
????????B=math.atan2（（z+N*e2*sinB）math.sqrt（x**2+y**2））
????????deltaB=B-B0
????????
????????if?deltaB>1e-6:
????????????break
????
????BLH.append（B）
????BLH.append（L）
???

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2351??2019-06-04?20:32??coorsystem.py
?????文件???????15456??2019-06-09?15:05??GPSmain.py
?????文件???????11481??2019-06-04?20:34??GPSreadfile.py
?????文件????????1100??2019-06-06?21:30??timetranslate.py