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