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

Last change on this file since 6630 was 6563, checked in by stuerze, 10 years ago

try to catch incorrect user input

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