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

Last change on this file since 10386 was 8756, checked in by mervart, 5 years ago

Sort MP data (make the largest values always visible)

File size: 5.6 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 "ephemeris.h"
32#include "satObs.h"
33#include "polarplot.h"
34
35class t_plotData {
36 public:
37 struct t_hlpStatus {
38 QVector<double> _ok;
39 QVector<double> _slip;
40 QVector<double> _gap;
41 };
42 QVector<double> _mjdX24;
43 QVector<double> _numSat;
44 QVector<double> _PDOP;
45 QVector<double> _eleDeg;
46 QMap<char, t_hlpStatus> _status;
47};
48
49class t_reqcAnalyze : public QThread {
50Q_OBJECT
51
52 public:
53 t_reqcAnalyze(QObject* parent);
54 virtual void run();
55
56 protected:
57 ~t_reqcAnalyze();
58
59 class t_skyPlotData {
60 public:
61 t_skyPlotData() {
62 _sys = ' ';
63 _data = new QVector<t_polarPoint*>;
64 }
65 t_skyPlotData(char sys) {
66 _sys = sys;
67 _data = new QVector<t_polarPoint*>;
68 }
69 ~t_skyPlotData() {}
70 char _sys;
71 QString _title;
72 QVector<t_polarPoint*>* _data;
73 };
74
75 signals:
76 void finished();
77 void dspSkyPlot(const QString&, QVector<t_skyPlotData> skyPlotData,
78 const QByteArray&, double);
79 void dspAvailPlot(const QString&, const QByteArray&);
80
81 private:
82
83 class t_qcFrq {
84 public:
85 t_qcFrq() {
86 _phaseValid = false;
87 _codeValid = false;
88 _slip = false;
89 _gap = false;
90 _setMP = false;
91 _rawMP = 0.0;
92 _stdMP = 0.0;
93 _SNR = 0.0;
94 }
95 QString _rnxType2ch;
96 bool _phaseValid;
97 bool _codeValid;
98 bool _slip;
99 bool _gap;
100 bool _setMP;
101 double _rawMP;
102 double _stdMP;
103 double _SNR;
104 };
105
106 class t_qcSat {
107 public:
108 t_qcSat() {
109 _slotSet = false;
110 _eleSet = false;
111 _slotNum = 0;
112 _eleDeg = 0.0;
113 _azDeg = 0.0;
114 }
115 bool _slotSet;
116 bool _eleSet;
117 int _slotNum;
118 double _eleDeg;
119 double _azDeg;
120 QVector<t_qcFrq> _qcFrq;
121 };
122
123 class t_qcEpo {
124 public:
125 bncTime _epoTime;
126 double _PDOP;
127 QMap<t_prn, t_qcSat> _qcSat;
128 };
129
130 class t_qcFrqSum {
131 public:
132 t_qcFrqSum() {
133 _numObs = 0;
134 _numSlipsFlagged = 0;
135 _numSlipsFound = 0;
136 _numGaps = 0;
137 _numSNR = 0;
138 _sumSNR = 0.0;
139 _numMP = 0;
140 _sumMP = 0.0;
141 }
142 int _numObs;
143 int _numSlipsFlagged;
144 int _numSlipsFound;
145 int _numGaps;
146 int _numSNR;
147 double _sumSNR;
148 int _numMP;
149 double _sumMP;
150 };
151
152 class t_qcSatSum {
153 public:
154 QMap<QString, t_qcFrqSum> _qcFrqSum;
155 };
156
157 class t_qcFile {
158 public:
159 t_qcFile() {
160 clear();
161 _interval = 1.0;
162 }
163 void clear() {_qcSatSum.clear(); _qcEpo.clear();}
164 bncTime _startTime;
165 bncTime _endTime;
166 QString _antennaName;
167 QString _markerName;
168 QString _receiverType;
169 double _interval;
170 QMap<t_prn, t_qcSatSum> _qcSatSum;
171 QVector<t_qcEpo> _qcEpo;
172 };
173
174 private slots:
175 void slotDspSkyPlot(const QString& fileName, QVector<t_skyPlotData> skyPlotData,
176 const QByteArray& scaleTitle, double maxValue);
177
178 void slotDspAvailPlot(const QString& fileName, const QByteArray& title);
179
180 private:
181 void checkEphemerides();
182
183 void analyzePlotSignals();
184
185 void analyzeFile(t_rnxObsFile* obsFile);
186
187 void updateQcSat(const t_qcSat& qcSat, t_qcSatSum& qcSatSum);
188
189 void setQcObs(const bncTime& epoTime, const ColumnVector& xyzSta,
190 const t_satObs& satObs, QMap<QString, bncTime>& lastObsTime, t_qcSat& qcSat);
191
192 void setExpectedObs(const bncTime& startTime, const bncTime& endTime,
193 double interval, const ColumnVector& xyzSta);
194
195 void analyzeMultipath();
196
197 void preparePlotData(const t_rnxObsFile* obsFile);
198
199 double cmpDOP(const ColumnVector& xyzSta) const;
200
201 void printReport(const t_rnxObsFile* obsFile);
202
203 static bool mpLessThan(const t_polarPoint* p1, const t_polarPoint* p2);
204
205 QString _logFileName;
206 QFile* _logFile;
207 QTextStream* _log;
208 QStringList _obsFileNames;
209 QVector<t_rnxObsFile*> _rnxObsFiles;
210 QStringList _navFileNames;
211 QString _reqcPlotSignals;
212 QMap<char, QVector<char> > _signalTypes;
213 QMap<t_prn, int> _numExpObs;
214 QVector<char> _navFileIncomplete;
215 QStringList _defaultSignalTypes;
216 QVector<t_eph*> _ephs;
217 t_rnxObsFile::t_rnxEpo* _currEpo;
218 t_qcFile _qcFile;
219};
220
221#endif
Note: See TracBrowser for help on using the repository browser.