source: ntrip/trunk/BNC/src/combination/bnccomb.h@ 6159

Last change on this file since 6159 was 6159, checked in by mervart, 8 years ago
File size: 3.5 KB
Line 
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
10class bncRtnetDecoder;
11class bncSP3;
12class bncAntex;
13
14class 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> orbCorrections);
24 void slotNewClkCorrections(QList<t_clkCorr> clkCorrections);
25
26 signals:
27 void newMessage(QByteArray msg, bool showOnScreen);
28 void newClkCorrections(QList<t_clkCorr> clkCorrections);
29 void newOrbCorrections(QList<t_orbCorr> orbCorrections);
30
31 private:
32 enum e_method{singleEpoch, filter};
33
34 class cmbParam {
35 public:
36 enum parType {offACgps, offACglo, offACSat, clkSat};
37 cmbParam(parType type_, int index_, const QString& ac_, const QString& prn_);
38 ~cmbParam();
39 double partial(const QString& AC_, const QString& prn_);
40 QString toString() const;
41 parType type;
42 int index;
43 QString AC;
44 QString prn;
45 double xx;
46 double sig0;
47 double sigP;
48 bool epoSpec;
49 const t_eph* eph;
50 };
51
52 class cmbAC {
53 public:
54 cmbAC() {
55 weight = 0.0;
56 numObs = 0;
57 }
58 ~cmbAC() {}
59 QString mountPoint;
60 QString name;
61 double weight;
62 unsigned numObs;
63 };
64
65 class cmbCorr {
66 public:
67 cmbCorr() {
68 _eph = 0;
69 _iod = 0;
70 }
71 ~cmbCorr() {}
72 QString _prn;
73 bncTime _time;
74 int _iod;
75 const t_eph* _eph;
76 t_orbCorr _orbCorr; // used for input
77 t_clkCorr _clkCorr; // used for input
78 QString _acName;
79 double _dClk; // used for output
80 ColumnVector _diffRao;
81 QString ID() {return _acName + "_" + _prn;}
82 };
83
84 class cmbEpoch {
85 public:
86 cmbEpoch() {}
87 ~cmbEpoch() {
88 QVectorIterator<cmbCorr*> it(corrs);
89 while (it.hasNext()) {
90 delete it.next();
91 }
92 }
93 QVector<cmbCorr*> corrs;
94 };
95
96 void processEpoch();
97 t_irc processEpoch_filter(QTextStream& out, QMap<QString, cmbCorr*>& resCorr,
98 ColumnVector& dx);
99 t_irc processEpoch_singleEpoch(QTextStream& out, QMap<QString, cmbCorr*>& resCorr,
100 ColumnVector& dx);
101 t_irc createAmat(Matrix& AA, ColumnVector& ll, DiagonalMatrix& PP,
102 const ColumnVector& x0, QMap<QString, cmbCorr*>& resCorr);
103 void dumpResults(const QMap<QString, cmbCorr*>& resCorr);
104 void printResults(QTextStream& out, const QMap<QString, cmbCorr*>& resCorr);
105 void switchToLastEph(const t_eph* lastEph, cmbCorr* corr);
106 t_irc checkOrbits(QTextStream& out);
107 QVector<cmbCorr*>& corrs() {return _buffer[_resTime].corrs;}
108
109 QList<cmbAC*> _ACs;
110 bncTime _resTime;
111 QVector<cmbParam*> _params;
112 QMap<bncTime, cmbEpoch> _buffer;
113 bncRtnetDecoder* _rtnetDecoder;
114 SymmetricMatrix _QQ;
115 QByteArray _log;
116 bncAntex* _antex;
117 double _MAXRES;
118 QString _masterOrbitAC;
119 unsigned _masterMissingEpochs;
120 e_method _method;
121 bool _useGlonass;
122 int _cmbSampl;
123 QMap<QString, QMap<t_prn, t_orbCorr> > _orbCorrections;
124};
125
126#endif
Note: See TracBrowser for help on using the repository browser.