1 |
|
---|
2 | #ifndef BNCCOMB_H
|
---|
3 | #define BNCCOMB_H
|
---|
4 |
|
---|
5 | #include <fstream>
|
---|
6 | #include <newmat.h>
|
---|
7 | #include "bncephuser.h"
|
---|
8 | #include "satObs.h"
|
---|
9 |
|
---|
10 | class bncRtnetDecoder;
|
---|
11 | class bncSP3;
|
---|
12 | class bncAntex;
|
---|
13 |
|
---|
14 | class bncComb : public bncEphUser {
|
---|
15 | Q_OBJECT
|
---|
16 | public:
|
---|
17 | bncComb();
|
---|
18 | virtual ~bncComb();
|
---|
19 | int nStreams() const {return _ACs.size();}
|
---|
20 |
|
---|
21 | public slots:
|
---|
22 | void slotProviderIDChanged(QString mountPoint);
|
---|
23 | void slotNewOrbCorrections(QList<t_orbCorr> orbCorr);
|
---|
24 | void slotNewClkCorrections(QList<t_clkCorr> clkCorr);
|
---|
25 |
|
---|
26 | signals:
|
---|
27 | void newMessage(QByteArray msg, bool showOnScreen);
|
---|
28 | void newClkCorrections(QList<t_clkCorr> clkCorr);
|
---|
29 |
|
---|
30 | private:
|
---|
31 | enum e_method{singleEpoch, filter};
|
---|
32 |
|
---|
33 | class cmbParam {
|
---|
34 | public:
|
---|
35 | enum parType {offACgps, offACglo, offACSat, clkSat};
|
---|
36 | cmbParam(parType type_, int index_, const QString& ac_, const QString& prn_);
|
---|
37 | ~cmbParam();
|
---|
38 | double partial(const QString& AC_, const QString& prn_);
|
---|
39 | QString toString() const;
|
---|
40 | parType type;
|
---|
41 | int index;
|
---|
42 | QString AC;
|
---|
43 | QString prn;
|
---|
44 | double xx;
|
---|
45 | double sig0;
|
---|
46 | double sigP;
|
---|
47 | bool epoSpec;
|
---|
48 | const t_eph* eph;
|
---|
49 | };
|
---|
50 |
|
---|
51 | class cmbAC {
|
---|
52 | public:
|
---|
53 | cmbAC() {
|
---|
54 | weight = 0.0;
|
---|
55 | numObs = 0;
|
---|
56 | }
|
---|
57 | ~cmbAC() {}
|
---|
58 | QString mountPoint;
|
---|
59 | QString name;
|
---|
60 | double weight;
|
---|
61 | unsigned numObs;
|
---|
62 | };
|
---|
63 |
|
---|
64 | class cmbCorr : public t_clkCorr {
|
---|
65 | public:
|
---|
66 | QString acName;
|
---|
67 | ColumnVector diffRao;
|
---|
68 | QString ID() {return acName + "_" + QString(_prn.toString().c_str());}
|
---|
69 | };
|
---|
70 |
|
---|
71 | class cmbEpoch {
|
---|
72 | public:
|
---|
73 | cmbEpoch() {}
|
---|
74 | ~cmbEpoch() {
|
---|
75 | QVectorIterator<cmbCorr*> it(corrs);
|
---|
76 | while (it.hasNext()) {
|
---|
77 | delete it.next();
|
---|
78 | }
|
---|
79 | }
|
---|
80 | QVector<cmbCorr*> corrs;
|
---|
81 | };
|
---|
82 |
|
---|
83 | void processEpoch();
|
---|
84 | t_irc processEpoch_filter(QTextStream& out, QMap<QString, t_clkCorr*>& resCorr,
|
---|
85 | ColumnVector& dx);
|
---|
86 | t_irc processEpoch_singleEpoch(QTextStream& out, QMap<QString, t_clkCorr*>& resCorr,
|
---|
87 | ColumnVector& dx);
|
---|
88 | t_irc createAmat(Matrix& AA, ColumnVector& ll, DiagonalMatrix& PP,
|
---|
89 | const ColumnVector& x0, QMap<QString, t_clkCorr*>& resCorr);
|
---|
90 | void dumpResults(const QMap<QString, t_clkCorr*>& resCorr);
|
---|
91 | void printResults(QTextStream& out, const QMap<QString, t_clkCorr*>& resCorr);
|
---|
92 | void switchToLastEph(const t_eph* lastEph, t_clkCorr* corr);
|
---|
93 | t_irc checkOrbits(QTextStream& out);
|
---|
94 | t_irc mergeOrbitCorr(const cmbCorr* orbitCorr, cmbCorr* clkCorr);
|
---|
95 |
|
---|
96 | QVector<cmbCorr*>& corrs() {return _buffer[_resTime].corrs;}
|
---|
97 |
|
---|
98 | QList<cmbAC*> _ACs;
|
---|
99 | bncTime _resTime;
|
---|
100 | QVector<cmbParam*> _params;
|
---|
101 | QMap<bncTime, cmbEpoch> _buffer;
|
---|
102 | bncRtnetDecoder* _rtnetDecoder;
|
---|
103 | SymmetricMatrix _QQ;
|
---|
104 | QByteArray _log;
|
---|
105 | bncAntex* _antex;
|
---|
106 | double _MAXRES;
|
---|
107 | QString _masterOrbitAC;
|
---|
108 | unsigned _masterMissingEpochs;
|
---|
109 | e_method _method;
|
---|
110 | bool _useGlonass;
|
---|
111 | int _cmbSampl;
|
---|
112 | QMap<QString, cmbCorr*> _orbitCorrs;
|
---|
113 | };
|
---|
114 |
|
---|
115 | #endif
|
---|