source: ntrip/trunk/BNC/src/rinex/rnxobsfile.h@ 5908

Last change on this file since 5908 was 5883, checked in by mervart, 10 years ago
File size: 6.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 RNXOBSFILE_H
26#define RNXOBSFILE_H
27
28#include <QtCore>
29
30#include <fstream>
31#include <vector>
32#include <map>
33
34#include "newmat.h"
35#include "bncconst.h"
36#include "bnctime.h"
37#include "t_prn.h"
38#include "GPSDecoder.h"
39
40class t_rnxObsHeader {
41 public:
42 static const double defaultRnxObsVersion2 = 2.11;
43 static const double defaultRnxObsVersion3 = 3.02;
44
45 t_rnxObsHeader();
46 ~t_rnxObsHeader();
47
48 t_irc read(QTextStream* stream, int maxLines = 0);
49 int nTypes(char sys) const;
50 const QString& obsType(char sys, int index) const;
51 QStringList obsTypesStrings() const;
52 void write(QTextStream* stream,
53 const QMap<QString, QString>* txtMap = 0) const;
54
55 static const QString _emptyStr;
56 float _version;
57 double _interval;
58 QString _antennaNumber;
59 QString _antennaName;
60 QString _markerName;
61 QString _markerNumber;
62 QString _observer;
63 QString _agency;
64 QString _receiverNumber;
65 QString _receiverType;
66 QString _receiverVersion;
67 QStringList _comments;
68 ColumnVector _antNEU;
69 ColumnVector _antXYZ;
70 ColumnVector _antBSG;
71 ColumnVector _xyz;
72 QVector<QString> _obsTypesV2;
73 QMap<char, QVector<QString> > _obsTypesV3;
74 int _wlFactorsL1[t_prn::MAXPRN_GPS+1];
75 int _wlFactorsL2[t_prn::MAXPRN_GPS+1];
76 bncTime _startTime;
77};
78
79class t_rnxObsFile {
80 public:
81
82 static bool earlierStartTime(const t_rnxObsFile* file1, const t_rnxObsFile* file2) {
83 return file1->startTime() < file2->startTime();
84 }
85
86 class t_rnxSat {
87 public:
88 char satSys;
89 int satNum;
90 std::vector<double> obs;
91 std::vector<int> lli;
92 std::vector<int> snr;
93 };
94
95 class t_rnxEpo {
96 public:
97 void clear() {
98 rnxSat.clear();
99 }
100 bncTime tt;
101 std::vector<t_rnxSat> rnxSat;
102 };
103
104 enum e_inpOut {input, output};
105
106 t_rnxObsFile(const QString& fileName, e_inpOut inpOut);
107 ~t_rnxObsFile();
108
109 float version() const {return _header._version;}
110 double interval() const {return _header._interval;}
111 int nTypes(char sys) const {return _header.nTypes(sys);}
112 const QString& fileName() const {return _fileName;}
113 const QString& obsType(char sys, int index) const {return _header.obsType(sys, index);}
114
115 const QString& antennaName() const {return _header._antennaName;}
116 const QString& markerName() const {return _header._markerName;}
117 const QString& receiverType() const {return _header._receiverType;}
118
119 void setInterval(double interval) {_header._interval = interval;}
120 void setAntennaName(const QString& antennaName) {_header._antennaName = antennaName;}
121 void setMarkerName(const QString& markerName) {_header._markerName = markerName;}
122 void setReceiverType(const QString& receiverType) {_header._receiverType = receiverType;}
123
124 const ColumnVector& xyz() const {return _header._xyz;}
125 const ColumnVector& antNEU() const {return _header._antNEU;}
126 const ColumnVector& antXYZ() const {return _header._antXYZ;}
127 const ColumnVector& antBSG() const {return _header._antBSG;}
128
129 const bncTime& startTime() const {return _header._startTime;}
130 void setStartTime(const bncTime& startTime) {_header._startTime = startTime;}
131
132 t_rnxEpo* nextEpoch();
133 int wlFactorL1(unsigned iPrn) {
134 return iPrn <= t_prn::MAXPRN_GPS ? _header._wlFactorsL1[iPrn] : 1;
135 }
136 int wlFactorL2(unsigned iPrn) {
137 return iPrn <= t_prn::MAXPRN_GPS ? _header._wlFactorsL2[iPrn] : 1;
138 }
139
140 const t_rnxObsHeader& header() const {return _header;}
141 void setHeader(const t_rnxObsHeader& header, double version);
142 void checkNewHeader(const t_rnxObsHeader& header);
143 void writeEpoch(const t_rnxEpo* epo);
144
145 QTextStream* stream() {return _stream;}
146
147 static void setObsFromRnx(const t_rnxObsFile* rnxObsFile,
148 const t_rnxObsFile::t_rnxEpo* epo,
149 const t_rnxObsFile::t_rnxSat& rnxSat,
150 t_obs& obs);
151
152 private:
153 enum e_trafo {trafoNone, trafo2to3, trafo3to2};
154
155 t_rnxObsFile() {};
156 void openRead(const QString& fileName);
157 void openWrite(const QString& fileName);
158 void close();
159 void writeEpochV2(const t_rnxEpo* epo);
160 void writeEpochV3(const t_rnxEpo* epo);
161 t_rnxEpo* nextEpochV2();
162 t_rnxEpo* nextEpochV3();
163 void handleEpochFlag(int flag, const QString& line, bool& headerReRead);
164
165 QString type2to3(char sys, const QString& typeV2);
166 QString type3to2(const QString& typeV3);
167
168 QMap<char, QMap<int, int> > _indexMap2to3;
169 QMap<char, QMap<int, int> > _indexMap3to2;
170
171 e_inpOut _inpOut;
172 QFile* _file;
173 QString _fileName;
174 QTextStream* _stream;
175 t_rnxObsHeader _header;
176 t_rnxEpo _currEpo;
177 bool _flgPowerFail;
178 e_trafo _trafo;
179};
180
181#endif
Note: See TracBrowser for help on using the repository browser.