source: ntrip/trunk/BNC/rnxobsfile.cpp@ 3682

Last change on this file since 3682 was 3682, checked in by mervart, 13 years ago
File size: 6.0 KB
RevLine 
[3645]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/* -------------------------------------------------------------------------
26 * BKG NTRIP Client
27 * -------------------------------------------------------------------------
28 *
29 * Class: t_rnxObsFile
30 *
31 * Purpose: Reads RINEX Observation File
32 *
33 * Author: L. Mervart
34 *
35 * Created: 24-Jan-2012
36 *
37 * Changes:
38 *
39 * -----------------------------------------------------------------------*/
40
41#include <iostream>
[3677]42#include <iomanip>
[3645]43#include "rnxobsfile.h"
[3646]44#include "bncutils.h"
[3645]45
46using namespace std;
47
48// Constructor
49////////////////////////////////////////////////////////////////////////////
[3646]50t_rnxObsFile::t_rnxObsHeader::t_rnxObsHeader() {
[3672]51 _antNEU.ReSize(3);
52 _xyz.ReSize(3);
53 _antNEU = 0.0;
54 _xyz = 0.0;
55 _version = 0.0;
[3645]56}
57
58// Destructor
59////////////////////////////////////////////////////////////////////////////
[3646]60t_rnxObsFile::t_rnxObsHeader::~t_rnxObsHeader() {
61}
62
63// Read Header
64////////////////////////////////////////////////////////////////////////////
65t_irc t_rnxObsFile::t_rnxObsHeader::read(QTextStream* stream) {
[3649]66 while (stream->status() == QTextStream::Ok && !stream->atEnd()) {
[3646]67 QString line = stream->readLine();
[3656]68 if (line.isEmpty()) {
69 continue;
70 }
[3651]71 QString value = line.left(60).trimmed();
72 QString key = line.mid(60).trimmed();
73 if (key == "END OF HEADER") {
[3650]74 break;
75 }
[3652]76 else if (key == "RINEX VERSION / TYPE") {
77 QTextStream in(value.toAscii(), QIODevice::ReadOnly);
[3651]78 in >> _version;
79 }
[3652]80 else if (key == "MARKER NAME") {
81 _markerName = value;
82 }
83 else if (key == "ANT # / TYPE") {
84 _antennaName = value.mid(20);
85 }
[3656]86 else if (key == "APPROX POSITION XYZ") {
87 QTextStream in(value.toAscii(), QIODevice::ReadOnly);
88 in >> _xyz[0] >> _xyz[1] >> _xyz[2];
89 }
90 else if (key == "ANTENNA: DELTA H/E/N") {
91 QTextStream in(value.toAscii(), QIODevice::ReadOnly);
92 in >> _antNEU[2] >> _antNEU[1] >> _antNEU[0];
93 }
[3672]94 else if (key == "# / TYPES OF OBSERV") {
95 QTextStream in(value.toAscii(), QIODevice::ReadOnly);
96 int nTypes;
97 in >> nTypes;
98 for (int ii = 0; ii < nTypes; ii++) {
99 QString hlp;
100 in >> hlp;
101 _obsTypes << hlp;
102 }
103 }
[3646]104 }
105
106 return success;
107}
108
109// Constructor
110////////////////////////////////////////////////////////////////////////////
111t_rnxObsFile::t_rnxObsFile(QString fileName) {
112 expandEnvVar(fileName);
113 _file = new QFile(fileName);
114 _file->open(QIODevice::ReadOnly | QIODevice::Text);
115 _stream = new QTextStream();
116 _stream->setDevice(_file);
117 _header.read(_stream);
118}
119
120// Destructor
121////////////////////////////////////////////////////////////////////////////
[3645]122t_rnxObsFile::~t_rnxObsFile() {
[3646]123 delete _stream;
124 delete _file;
[3645]125}
126
127// Retrieve single Epoch
128////////////////////////////////////////////////////////////////////////////
[3675]129const t_rnxObsFile::t_epo* t_rnxObsFile::nextEpoch() {
[3676]130
131 _currEpo.clear();
132
[3673]133 if (version() < 3.0) {
[3675]134 return nextEpochV2();
[3673]135 }
136 else {
[3675]137 return nextEpochV3();
[3673]138 }
139}
140
141// Retrieve single Epoch (RINEX Version 3)
142////////////////////////////////////////////////////////////////////////////
[3675]143const t_rnxObsFile::t_epo* t_rnxObsFile::nextEpochV3() {
144 return 0; // TODO
[3673]145}
146
147// Retrieve single Epoch (RINEX Version 2)
148////////////////////////////////////////////////////////////////////////////
[3675]149const t_rnxObsFile::t_epo* t_rnxObsFile::nextEpochV2() {
[3656]150 while (_stream->status() == QTextStream::Ok && !_stream->atEnd()) {
[3680]151
[3656]152 QString line = _stream->readLine();
[3680]153
[3656]154 if (line.isEmpty()) {
155 continue;
156 }
[3680]157
[3673]158 QTextStream in(line.toAscii());
[3680]159
160 // Epoch Time
161 // ----------
[3673]162 int year, month, day, hour, min, flag;
163 double sec;
164 in >> year >> month >> day >> hour >> min >> sec >> flag;
165 if (year < 80) {
166 year += 2000;
167 }
168 else if (year < 100) {
169 year += 1900;
170 }
[3680]171 _currEpo.tt.set(year, month, day, hour, min, sec);
[3673]172
[3680]173 // Number of Satellites
174 // --------------------
[3673]175 int numSat;
176 readInt(line, 29, 3, numSat);
[3674]177
[3677]178 _currEpo.satObs.resize(numSat);
179
[3676]180 // Read Satellite Numbers
181 // ----------------------
[3674]182 int pos = 32;
183 for (int iSat = 0; iSat < numSat; iSat++) {
184 if (iSat > 0 && iSat % 12 == 0) {
185 line = _stream->readLine();
186 pos = 32;
187 }
188 QString prn = line.mid(pos, 3);
189 pos += 3;
[3676]190
[3677]191 _currEpo.satObs[iSat].prn = prn;
192 _currEpo.satObs[iSat].ReSize(_header.nTypes());
[3674]193 }
194
[3676]195 // Read Observation Records
196 // ------------------------
197 for (int iSat = 0; iSat < numSat; iSat++) {
198 line = _stream->readLine();
199 pos = 0;
[3677]200 QString prn = _currEpo.satObs[iSat].prn;
[3676]201 for (int iType = 0; iType < _header.nTypes(); iType++) {
202 if (iType > 0 && iType % 5 == 0) {
203 line = _stream->readLine();
204 pos = 0;
205 }
[3677]206 readDbl(line, pos, 14, _currEpo.satObs[iSat][iType]);
207 readInt(line, pos + 14, 1, _currEpo.satObs[iSat].lli);
208 readInt(line, pos + 15, 1, _currEpo.satObs[iSat].snr);
[3676]209 pos += 16;
210 }
211 }
212
[3678]213 return &_currEpo;
[3656]214 }
[3674]215
[3678]216 return 0;
[3645]217}
Note: See TracBrowser for help on using the repository browser.