source: ntrip/trunk/BNC/src/PPP/pppSatObs.h@ 5830

Last change on this file since 5830 was 5830, checked in by mervart, 10 years ago
File size: 3.7 KB
Line 
1#ifndef SATOBS_H
2#define SATOBS_H
3
4#include <string>
5#include <map>
6#include <newmat.h>
7#include "pppInclude.h"
8#include "bnctime.h"
9
10namespace BNC_PPP {
11
12class t_pppStation;
13
14class t_pppSatObs {
15 public:
16 t_pppSatObs(const t_satObs& satObs);
17 ~t_pppSatObs();
18 bool isValid() const {return _valid;}
19 void prepareObs();
20 const t_prn& prn() const {return _prn;}
21 const ColumnVector& xc() const {return _xcSat;}
22 const bncTime& time() const {return _time;}
23 t_irc cmpModel(const t_pppStation* station);
24 double obsValue(t_lc::type tLC) const;
25 double cmpValue(t_lc::type tLC) const;
26 double cmpValueForBanc(t_lc::type tLC) const;
27 double rho() const {return _model._rho;}
28 double sagnac() const {return _model._sagnac;}
29 double eleSat() const {return _model._eleSat;}
30 bool modelSet() const {return _model._set;}
31 void printModel() const;
32 t_irc createDifference(const t_pppSatObs* obsBase);
33 double lc(t_lc::type tLC, double L1, double L2,
34 double C1, double C2, ColumnVector* coeff = 0) const;
35 double lambda(t_lc::type tLC) const;
36 double sigma(t_lc::type tLC) const;
37 bool outlier() const {return _outlier;}
38 void setOutlier() {_outlier = true;}
39 int channel() const {return _channel;}
40 void setRes(t_lc::type tLC, double res);
41 double getRes(t_lc::type tLC) const;
42
43 bool slip() const {
44 for (unsigned ii = 0; ii < _allObs.size(); ii++) {
45 if (_allObs[ii]->_slip) {
46 return true;
47 }
48 }
49 return false;
50 }
51
52 int slipCounter() const {
53 int cnt = -1;
54 for (unsigned ii = 0; ii < _allObs.size(); ii++) {
55 if (_allObs[ii]->_slipCounter > cnt) {
56 cnt = _allObs[ii]->_slipCounter;
57 }
58 }
59 return cnt;
60 }
61
62 int biasJumpCounter() const {
63 int jmp = -1;
64 for (unsigned ii = 0; ii < _allObs.size(); ii++) {
65 if (_allObs[ii]->_biasJumpCounter > jmp) {
66 jmp = _allObs[ii]->_biasJumpCounter;
67 }
68 }
69 return jmp;
70 }
71
72 private:
73 class t_model {
74 public:
75 t_model() {reset();}
76 ~t_model() {}
77 void reset() {
78 _set = false;
79 _rho = 0.0;
80 _eleSat = 0.0;
81 _azSat = 0.0;
82 _recClkM = 0.0;
83 _satClkM = 0.0;
84 _sagnac = 0.0;
85 _antEcc = 0.0;
86 _tropo = 0.0;
87 _tide = 0.0;
88 _windUp = 0.0;
89 _antPco1 = 0.0;
90 _antPco2 = 0.0;
91 _biasC1 = 0.0;
92 _biasC2 = 0.0;
93 _biasL1 = 0.0;
94 _biasL2 = 0.0;
95 }
96 bool _set;
97 double _rho;
98 double _eleSat;
99 double _azSat;
100 double _recClkM;
101 double _satClkM;
102 double _sagnac;
103 double _antEcc;
104 double _tropo;
105 double _tide;
106 double _windUp;
107 double _antPco1;
108 double _antPco2;
109 double _biasC1;
110 double _biasC2;
111 double _biasL1;
112 double _biasL2;
113 };
114
115 t_prn _prn;
116 bncTime _time;
117 int _channel;
118 std::vector<t_frqObs*> _allObs;
119 bool _valid;
120 t_frqObs* _validObs1;
121 t_frqObs* _validObs2;
122 double _f1;
123 double _f2;
124 double _rawC1;
125 double _rawC2;
126 double _rawL1;
127 double _rawL2;
128 ColumnVector _xcSat;
129 ColumnVector _vvSat;
130 t_model _model;
131 bool _outlier;
132 std::map<t_lc::type, double> _res;
133};
134
135}
136
137#endif
Note: See TracBrowser for help on using the repository browser.