GPS标准单点定位代码.zip

基于python编写的GPS单点定位,包含时间转换坐标系统转换，RINEX3.0 文件读取

# -*- 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）
   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        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