-
Notifications
You must be signed in to change notification settings - Fork 9
/
Ephemeris.h
103 lines (97 loc) · 3.42 KB
/
Ephemeris.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
#ifndef EPHEMERIS_H_
#define EPHEMERIS_H_
//#include "Time_utility.h"
#include "Rinex.h"
#include <vector>
//#include <map>
#include <set>
#include <string>
const double ura_value[] = { /* ura max values */
2.4,3.4,4.85,6.85,9.65,13.65,24.0,48.0,96.0,192.0,384.0,768.0,1536.0,
3072.0,6144.0
};
int uraindex(double value);
class Sat_EPH {
public:
Sat_EPH(int SatID, vector<double> ::iterator Satnavmessage);
//virtual void Sat_EPH_GBparaset(const double * Satnavmessage) = 0;
void calculate_Satpos_inorbit(gtime_t time);//传入此处的time是补偿了伪距时间后的发射时间
int getCRS(double & CRS);
int getSVH() { return svh; }
double cal_POSVAR();
gtime_t gettoe() { return toe; }
~Sat_EPH() {};
protected:
/*计算卫星位置有关参数*/
double mu; //为了计算卫星位置的参数
double omge;// RTKILIB中用宏表示，GPSBDS 不一样
double sinE;//中间参数
double tk;//发射时间与toc之差
double x_inorbit;//轨道平面位置，可通过它们进一步计算卫星位置
double y_inorbit;
double i_inorbit;
/*星历参数*/
int sat; /* satellite number */
int iode, iodc; /* IODE,IODC */
int sva; /* SV accuracy (URA index) */
int svh; /* SV health (0:ok) */
int week; /* GPS/QZS: gps week, GAL: galileo week */
//int code; /* GPS/QZS: code on L2, GAL/CMP: data sources */
gtime_t toe, toc, ttr; /* Toe,Toc,T_trans */
/* SV orbit parameters */
double A, e, i0, OMG0, omg, M0, deln, OMGd, idot;
double crc, crs, cuc, cus, cic, cis;
double toes; /* Toe (s) in week */
double f0, f1, f2; /* SV clock parameters (af0,af1,af2) */
//double tgd[2]; /* group delay parameters */
/* GPS/QZS:tgd[0]=TGD */
/* GAL :tgd[0]=BGD E5a/E1,tgd[1]=BGD E5b/E1 */
/* CMP :tgd[0]=BGD1,tgd[1]=BGD2 */
//double Adot, ndot; /* Adot,ndot for CNAV */
};
class GPS_Sat_EPH :public Sat_EPH {
public:
GPS_Sat_EPH(int SatID, vector<double>::iterator Satnavmessage, int SatToc);
int get_iode() { return iode; }
~GPS_Sat_EPH() {}
void Cali_GPS_Satinfo(gtime_t obstime, vector<double > &SAT_Pos, double &CLK_BIAS);
double getTGD() { return tgd; }
private:
int flag;
double fit;
int code;
double tgd;
};
class BDS_Sat_EPH :public Sat_EPH {
public:
BDS_Sat_EPH(int SatID, vector<double>::iterator Satnavmessage, int SatToc);
double getdeln() { return deln; }
~BDS_Sat_EPH() {}
void Cali_BDS_Satinfo(int ID,gtime_t obstime, vector<double > &SAT_Pos, double &CLK_BIAS);
double getTGD() { return tgd[0]; }
//void Cali_BDS_Satinfo(double obstime, vector<double > &SAT_Pos_Vel, vector<double> &CLK_BIAS_DRIFT);
private:
double tgd[2];
};
typedef std::pair<int, GPS_Sat_EPH> gepair;
typedef std::pair<int, BDS_Sat_EPH> bepair;
typedef std::multimap <int, GPS_Sat_EPH> Gsatephmmp;
typedef std::multimap <int, BDS_Sat_EPH> Bsatephmmp;
class All_Sats_EPH {
public:
All_Sats_EPH() {};
All_Sats_EPH(const NAV_RINEX &nav);
void satephsolve(NAV_RINEX & nav);
vector<double > getGPSIONCORR() { return GPS_ION_CORR; }
vector<double > getBDSIONCORR() { return BDS_ION_CORR; }
int Cal_OneSatinfo(int SatID,char mode, gtime_t obs_time,vector<double> & pos_vel, vector<double> &bias_drift, double &var,int &svh,double &TGD);
private:
int LEAP_SECONDS;
vector<double > GPS_ION_CORR;
vector<double > BDS_ION_CORR;
vector<double > GPS_TIME_CORR;
vector<double > BDS_TIME_CORR;
Gsatephmmp Gmmp;
Bsatephmmp Bmmp;
};
#endif