source: ntrip/trunk/BNC/src/rinex/reqcanalyze.h@ 4704

Last change on this file since 4704 was 4704, checked in by mervart, 10 years ago
File size: 4.9 KB
Line 
1// Part of BNC, a utility for retrieving decoding and
2// converting GNSS data streams from NTRIP broadcasters.
3//
4// Copyright (C) 2007
5// German Federal Agency for Cartography and Geodesy (BKG)
6// http://www.bkg.bund.de
7// Czech Technical University Prague, Department of Geodesy
8// http://www.fsv.cvut.cz
9//
10// Email: euref-ip@bkg.bund.de
11//
12// This program is free software; you can redistribute it and/or
13// modify it under the terms of the GNU General Public License
14// as published by the Free Software Foundation, version 2.
15//
16// This program is distributed in the hope that it will be useful,
17// but WITHOUT ANY WARRANTY; without even the implied warranty of
18// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19// GNU General Public License for more details.
20//
21// You should have received a copy of the GNU General Public License
22// along with this program; if not, write to the Free Software
23// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
24
25#ifndef REQCANALYZE_H
26#define REQCANALYZE_H
27
28#include <QtCore>
29#include "rnxobsfile.h"
30#include "rnxnavfile.h"
31#include "RTCM/GPSDecoder.h"
32#include "RTCM3/ephemeris.h"
33
34class t_polarPoint;
35
36class t_availData {
37 public:
38 QVector<double> _L1ok;
39 QVector<double> _L2ok;
40 QVector<double> _L1slip;
41 QVector<double> _L2slip;
42 QVector<double> _L1gap;
43 QVector<double> _L2gap;
44 QVector<double> _eleDeg;
45 QVector<double> _eleTim;
46};
47
48class t_prnStat {
49 public:
50 t_prnStat() {
51 _numObs = 0;
52 _numSlipsFlagged = 0;
53 _numSlipsFound = 0;
54 _numGaps = 0;
55 }
56 int _numObs;
57 int _numSlipsFlagged;
58 int _numSlipsFound;
59 int _numGaps;
60};
61
62class t_obsStat {
63 public:
64 void reset() {
65 _mjdX24.clear();
66 _numSat.clear();
67 _PDOP.clear();
68 _prnStat.clear();
69 }
70 QVector<double> _mjdX24;
71 QVector<double> _numSat;
72 QVector<double> _PDOP;
73 bncTime _startTime;
74 bncTime _endTime;
75 QString _antennaName;
76 QString _markerName;
77 QString _receiverType;
78 double _interval;
79 QMap<QString, t_prnStat> _prnStat;
80};
81
82class t_reqcAnalyze : public QThread {
83Q_OBJECT
84
85 public:
86 t_reqcAnalyze(QObject* parent);
87
88 protected:
89 ~t_reqcAnalyze();
90
91 signals:
92 void finished();
93 void dspSkyPlot(const QString&,
94 const QByteArray&,
95 QVector<t_polarPoint*>*,
96 const QByteArray&,
97 QVector<t_polarPoint*>*,
98 const QByteArray&, double);
99
100 void dspAvailPlot(const QString&, const QByteArray&);
101
102 private slots:
103 void slotDspSkyPlot(const QString& fileName,
104 const QByteArray& title1,
105 QVector<t_polarPoint*>* data1,
106 const QByteArray& title2,
107 QVector<t_polarPoint*>* data2,
108 const QByteArray& scaleTitle, double maxValue);
109
110 void slotDspAvailPlot(const QString& fileName, const QByteArray& title);
111
112 public:
113 virtual void run();
114
115 private:
116 class t_oneObs {
117 public:
118 t_oneObs(int GPSWeek, double GPSWeeks) {
119 _GPSWeek = GPSWeek;
120 _GPSWeeks = GPSWeeks;
121 _hasL1 = false;
122 _hasL2 = false;
123 _slipL1 = false;
124 _slipL2 = false;
125 _MP1 = 0.0;
126 _MP2 = 0.0;
127 _SNR1 = 0.0;
128 _SNR2 = 0.0;
129 }
130 int _GPSWeek;
131 double _GPSWeeks;
132 bool _hasL1;
133 bool _hasL2;
134 bool _slipL1;
135 bool _slipL2;
136 double _MP1;
137 double _MP2;
138 double _SNR1;
139 double _SNR2;
140 };
141
142 class t_allObs {
143 public:
144 t_allObs() {}
145 ~t_allObs() {
146 for (int ii = 0; ii < _oneObsVec.size(); ii++) {
147 delete _oneObsVec[ii];
148 }
149 }
150 t_irc addObs(const t_obs& obs);
151 QVector<t_oneObs*> _oneObsVec;
152 };
153
154 void analyzeFile(t_rnxObsFile* obsFile);
155 void preparePlotData(const QString& prn, const ColumnVector& xyzSta,
156 double obsInterval,
157 QVector<t_polarPoint*>* dataMP1,
158 QVector<t_polarPoint*>* dataMP2,
159 QVector<t_polarPoint*>* dataSNR1,
160 QVector<t_polarPoint*>* dataSNR2);
161 void prepareObsStat(unsigned iEpo, double obsInterval,
162 const ColumnVector& xyzSta);
163 double cmpDOP(const ColumnVector& xyzSta) const;
164 void printReport(QVector<t_polarPoint*>* dataMP1,
165 QVector<t_polarPoint*>* dataMP2,
166 QVector<t_polarPoint*>* dataSNR1,
167 QVector<t_polarPoint*>* dataSNR2);
168
169 QString _logFileName;
170 QFile* _logFile;
171 QTextStream* _log;
172 QStringList _obsFileNames;
173 QVector<t_rnxObsFile*> _rnxObsFiles;
174 QStringList _navFileNames;
175 QVector<t_eph*> _ephs;
176 t_rnxObsFile::t_rnxEpo* _currEpo;
177 QMap<QString, t_allObs> _allObsMap;
178 QMap<QString, t_availData> _availDataMap;
179 t_obsStat _obsStat;
180 QMutex _mutex;
181};
182
183#endif
Note: See TracBrowser for help on using the repository browser.