source: ntrip/branches/BNC_2.12/src/ephemeris.h@ 8167

Last change on this file since 8167 was 8167, checked in by stuerze, 6 years ago

IRNSS support is added in RINEX QC
File size: 12.5 KB
Line 
1#ifndef EPHEMERIS_H
2#define EPHEMERIS_H
3
4#include <newmat.h>
5#include <QtCore>
6#include <stdio.h>
7#include <string>
8#include "bnctime.h"
9#include "bncconst.h"
10#include "t_prn.h"
11#include "gnss.h"
12
13
14class t_orbCorr;
15class t_clkCorr;
16
17class t_eph {
18 public:
19 enum e_type {unknown, GPS, QZSS, GLONASS, Galileo, SBAS, BDS, IRNSS};
20 enum e_checkState {unchecked, ok, bad, outdated};
21
22 t_eph();
23 virtual ~t_eph();
24
25 virtual e_type type() const = 0;
26 virtual QString toString(double version) const = 0;
27 virtual unsigned int IOD() const = 0;
28 virtual int slotNum() const {return 0;}
29 bncTime TOC() const {return _TOC;}
30 bool isNewerThan(const t_eph* eph) const {return earlierTime(eph, this);}
31 e_checkState checkState() const {return _checkState;}
32 void setCheckState(e_checkState checkState) {_checkState = checkState;}
33 t_prn prn() const {return _prn;}
34 t_irc getCrd(const bncTime& tt, ColumnVector& xc, ColumnVector& vv, bool useCorr) const;
35 void setOrbCorr(const t_orbCorr* orbCorr);
36 void setClkCorr(const t_clkCorr* clkCorr);
37 const QDateTime& receptDateTime() const {return _receptDateTime;}
38 static QString rinexDateStr(const bncTime& tt, const t_prn& prn, double version);
39 static QString rinexDateStr(const bncTime& tt, const QString& prnStr, double version);
40 static bool earlierTime(const t_eph* eph1, const t_eph* eph2) {return eph1->_TOC < eph2->_TOC;}
41
42 protected:
43 virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const = 0;
44 t_prn _prn;
45 bncTime _TOC;
46 QDateTime _receptDateTime;
47 e_checkState _checkState;
48 t_orbCorr* _orbCorr;
49 t_clkCorr* _clkCorr;
50};
51
52
53class t_ephGPS : public t_eph {
54 friend class t_ephEncoder;
55 friend class RTCM3Decoder;
56 public:
57 t_ephGPS() {
58 _clock_bias = 0.0;
59 _clock_drift = 0.0;
60 _clock_driftrate = 0.0;
61 _IODE = 0.0;
62 _Crs = 0.0;
63 _Delta_n = 0.0;
64 _M0 = 0.0;
65 _Cuc = 0.0;
66 _e = 0.0;
67 _Cus = 0.0;
68 _sqrt_A = 0.0;
69 _TOEsec = 0.0;
70 _Cic = 0.0;
71 _OMEGA0 = 0.0;
72 _Cis = 0.0;
73 _i0 = 0.0;
74 _Crc = 0.0;
75 _omega = 0.0;
76 _OMEGADOT = 0.0;
77 _IDOT = 0.0;
78 _L2Codes = 0.0;
79 _TOEweek = 0.0;
80 _L2PFlag = 0.0;
81 _ura = 0.0;
82 _health = 0.0;
83 _TGD = 0.0;
84 _IODC = 0.0;
85 _TOT = 0.0;
86 _fitInterval = 0.0;
87 }
88 t_ephGPS(float rnxVersion, const QStringList& lines);
89 virtual ~t_ephGPS() {}
90
91 virtual e_type type() const {
92 switch (_prn.system()) {
93 case 'J':
94 return t_eph::QZSS;
95 case 'I':
96 return t_eph::IRNSS;
97 };
98 return t_eph::GPS;
99 }
100 virtual QString toString(double version) const;
101 virtual unsigned int IOD() const { return static_cast<unsigned int>(_IODE); }
102 double TGD() const {return _TGD;} // Timing Group Delay (P1-P2 DCB)
103
104 private:
105 virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;
106
107 double _clock_bias; // [s]
108 double _clock_drift; // [s/s]
109 double _clock_driftrate; // [s/s^2]
110
111 double _IODE; // IODEC in case of IRNSS
112 double _Crs; // [m]
113 double _Delta_n; // [rad/s]
114 double _M0; // [rad]
115
116 double _Cuc; // [rad]
117 double _e; //
118 double _Cus; // [rad]
119 double _sqrt_A; // [m^0.5]
120
121 double _TOEsec; // [s]
122 double _Cic; // [rad]
123 double _OMEGA0; // [rad]
124 double _Cis; // [rad]
125
126 double _i0; // [rad]
127 double _Crc; // [m]
128 double _omega; // [rad]
129 double _OMEGADOT; // [rad/s]
130
131 double _IDOT; // [rad/s]
132 double _L2Codes; // Codes on L2 channel (not valid for IRNSS)
133 double _TOEweek;
134 double _L2PFlag; // L2 P data flag (not valid for IRNSS)
135
136 mutable double _ura; // SV accuracy
137 double _health; // SV health
138 double _TGD; // [s]
139 double _IODC; // (not valid for IRNSS)
140
141 double _TOT; // Transmisstion time
142 double _fitInterval; // Fit interval (not valid for IRNSS)
143
144};
145
146class t_ephGlo : public t_eph {
147 friend class t_ephEncoder;
148 friend class RTCM3Decoder;
149 public:
150 t_ephGlo() {
151 _xv.ReSize(6);
152 _gps_utc = 0.0;
153 _tau = 0.0;
154 _gamma = 0.0;
155 _tki = 0.0;
156 _x_pos = 0.0;
157 _x_velocity = 0.0;
158 _x_acceleration = 0.0;
159 _health = 0.0;
160 _y_pos = 0.0;
161 _y_velocity = 0.0;
162 _y_acceleration = 0.0;
163 _frequency_number = 0.0;
164 _z_pos = 0.0;
165 _z_velocity = 0.0;
166 _z_acceleration = 0.0;
167 _E = 0.0;
168 }
169 t_ephGlo(float rnxVersion, const QStringList& lines);
170 virtual ~t_ephGlo() {}
171
172 virtual e_type type() const {return t_eph::GLONASS;}
173 virtual QString toString(double version) const;
174 virtual unsigned int IOD() const;
175 virtual int slotNum() const {return int(_frequency_number);}
176
177 private:
178 virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;
179 static ColumnVector glo_deriv(double /* tt */, const ColumnVector& xv, double* acc);
180
181 mutable bncTime _tt; // time
182 mutable ColumnVector _xv; // status vector (position, velocity) at time _tt
183
184 double _gps_utc;
185 double _tau; // [s]
186 double _gamma; //
187 mutable double _tki; // message frame time
188
189 double _x_pos; // [km]
190 double _x_velocity; // [km/s]
191 double _x_acceleration; // [km/s^2]
192 double _health; // 0 = O.K.
193
194 double _y_pos; // [km]
195 double _y_velocity; // [km/s]
196 double _y_acceleration; // [km/s^2]
197 double _frequency_number; // ICD-GLONASS data position
198
199 double _z_pos; // [km]
200 double _z_velocity; // [km/s]
201 double _z_acceleration; // [km/s^2]
202 double _E; // Age of Information [days]
203};
204
205class t_ephGal : public t_eph {
206 friend class t_ephEncoder;
207 friend class RTCM3Decoder;
208 public:
209 t_ephGal() {
210 _clock_bias = 0.0;
211 _clock_drift = 0.0;
212 _clock_driftrate = 0.0;
213 _IODnav = 0.0;
214 _Crs = 0.0;
215 _Delta_n = 0.0;
216 _M0 = 0.0;
217 _Cuc = 0.0;
218 _e = 0.0;
219 _Cus = 0.0;
220 _sqrt_A = 0.0;
221 _TOEsec = 0.0;
222 _Cic = 0.0;
223 _OMEGA0 = 0.0;
224 _Cis = 0.0;
225 _i0 = 0.0;
226 _Crc = 0.0;
227 _omega = 0.0;
228 _OMEGADOT = 0.0;
229 _IDOT = 0.0;
230 _TOEweek = 0.0;
231 _SISA = 0.0;
232 _E5aHS = 0.0;
233 _E5bHS = 0.0;
234 _E1_bHS = 0.0;
235 _BGD_1_5A = 0.0;
236 _BGD_1_5B = 0.0;
237 _TOT = 0.0;
238 };
239 t_ephGal(float rnxVersion, const QStringList& lines);
240 virtual ~t_ephGal() {}
241
242 virtual QString toString(double version) const;
243 virtual e_type type() const {return t_eph::Galileo;}
244 virtual unsigned int IOD() const { return static_cast<unsigned long>(_IODnav); }
245
246 private:
247 virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;
248
249 double _clock_bias; // [s]
250 double _clock_drift; // [s/s]
251 double _clock_driftrate; // [s/s^2]
252
253 double _IODnav;
254 double _Crs; // [m]
255 double _Delta_n; // [rad/s]
256 double _M0; // [rad]
257
258 double _Cuc; // [rad]
259 double _e; //
260 double _Cus; // [rad]
261 double _sqrt_A; // [m^0.5]
262
263 double _TOEsec; // [s]
264 double _Cic; // [rad]
265 double _OMEGA0; // [rad]
266 double _Cis; // [rad]
267
268 double _i0; // [rad]
269 double _Crc; // [m]
270 double _omega; // [rad]
271 double _OMEGADOT; // [rad/s]
272
273 double _IDOT; // [rad/s]
274 double _TOEweek;
275 // spare
276
277 mutable double _SISA; // Signal In Space Accuracy
278 double _E5aHS; // [0..3] E5a Health Status
279 double _E5bHS; // [0..3] E5b Health Status
280 double _E1_bHS; // [0..3] E1-b Health Status
281 double _BGD_1_5A; // group delay [s]
282 double _BGD_1_5B; // group delay [s]
283
284 double _TOT; // [s]
285 /** Data comes from I/NAV when <code>true</code> */
286 bool _inav;
287 /** Data comes from F/NAV when <code>true</code> */
288 bool _fnav;
289 /** EE Data is not valid */
290 bool _e1DataInValid;
291 /** E5A Data is not valid */
292 bool _e5aDataInValid;
293 /** E5B Data is not valid */
294 bool _e5bDataInValid;
295};
296
297class t_ephSBAS : public t_eph {
298 friend class t_ephEncoder;
299 friend class RTCM3Decoder;
300 public:
301 t_ephSBAS() {
302 _IODN = 0;
303 _TOW = 0.0;
304 _agf0 = 0.0;
305 _agf1 = 0.0;
306 _x_pos = 0.0;
307 _x_velocity = 0.0;
308 _x_acceleration = 0.0;
309 _y_pos = 0.0;
310 _y_velocity = 0.0;
311 _y_acceleration = 0.0;
312 _z_pos = 0.0;
313 _z_velocity = 0.0;
314 _z_acceleration = 0.0;
315 _ura = 0.0;
316 _health = 0.0;
317 }
318 t_ephSBAS(float rnxVersion, const QStringList& lines);
319 virtual ~t_ephSBAS() {}
320
321 virtual e_type type() const {return t_eph::SBAS;}
322 virtual unsigned int IOD() const;
323 virtual QString toString(double version) const;
324
325 private:
326 virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;
327
328 int _IODN;
329 double _TOW; // not used (set to 0.9999e9)
330 double _agf0; // [s] clock correction
331 double _agf1; // [s/s] clock correction drift
332
333 double _x_pos; // [m]
334 double _x_velocity; // [m/s]
335 double _x_acceleration; // [m/s^2]
336
337 double _y_pos; // [m]
338 double _y_velocity; // [m/s]
339 double _y_acceleration; // [m/s^2]
340
341 double _z_pos; // [m]
342 double _z_velocity; // [m/s]
343 double _z_acceleration; // [m/s^2]
344
345 mutable double _ura;
346 double _health;
347};
348
349class t_ephBDS : public t_eph {
350 friend class t_ephEncoder;
351 friend class RTCM3Decoder;
352 public:
353 t_ephBDS() : _TOEweek(-1.0) {
354 _TOT = 0.0;
355 _AODE = 0;
356 _AODC = 0;
357 _URAI = 0;
358 _URA = 0.0;
359 _clock_bias = 0.0;
360 _clock_drift = 0.0;
361 _clock_driftrate = 0.0;
362 _Crs = 0.0;
363 _Delta_n = 0.0;
364 _M0 = 0.0;
365 _Cuc = 0.0;
366 _e = 0.0;
367 _Cus = 0.0;
368 _sqrt_A = 0.0;
369 _Cic = 0.0;
370 _OMEGA0 = 0.0;
371 _Cis = 0.0;
372 _i0 = 0.0;
373 _Crc = 0.0;
374 _omega = 0.0;
375 _OMEGADOT = 0.0;
376 _IDOT = 0.0;
377 _TGD1 = 0.0;
378 _TGD2 = 0.0;
379 _SatH1 = 0.0;
380 _TOW = 0.0;
381 _TOEsec = 0.0;
382 _TOEweek = 0.0;
383 }
384 t_ephBDS(float rnxVersion, const QStringList& lines);
385 virtual ~t_ephBDS() {}
386
387 virtual e_type type() const {return t_eph::BDS;}
388 virtual unsigned int IOD() const;
389 virtual QString toString(double version) const;
390
391 private:
392 virtual t_irc position(int GPSweek, double GPSweeks, double* xc, double* vv) const;
393
394 double _TOT;
395 bncTime _TOE;
396 int _AODE;
397 int _AODC;
398 int _URAI; // [0..15] index from RTCM stream
399 mutable double _URA; // user range accuracy
400 double _clock_bias; // [s]
401 double _clock_drift; // [s/s]
402 double _clock_driftrate; // [s/s^2]
403 double _Crs; // [m]
404 double _Delta_n; // [rad/s]
405 double _M0; // [rad]
406 double _Cuc; // [rad]
407 double _e; //
408 double _Cus; // [rad]
409 double _sqrt_A; // [m^0.5]
410 double _Cic; // [rad]
411 double _OMEGA0; // [rad]
412 double _Cis; // [rad]
413 double _i0; // [rad]
414 double _Crc; // [m]
415 double _omega; // [rad]
416 double _OMEGADOT; // [rad/s]
417 double _IDOT; // [rad/s]
418 double _TGD1; // [s]
419 double _TGD2; // [s]
420 int _SatH1; //
421 double _TOW; // [s] of BDT week
422 double _TOEsec; // [s] of BDT week
423 double _TOEweek; // BDT week will be set only in case of RINEX file input
424};
425
426#endif
Note: See TracBrowser for help on using the repository browser.