source: ntrip/trunk/BNC/src/rinex/bncpostprocess.cpp@ 5262

Last change on this file since 5262 was 5262, checked in by mervart, 9 years ago
File size: 7.3 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/* -------------------------------------------------------------------------
26 * BKG NTRIP Client
27 * -------------------------------------------------------------------------
28 *
29 * Class: t_postProcessing
30 *
31 * Purpose: Precise Point Positioning in Post-Processing Mode
32 *
33 * Author: L. Mervart
34 *
35 * Created: 22-Jan-2012
36 *
37 * Changes:
38 *
39 * -----------------------------------------------------------------------*/
40
41#include <iostream>
42#include "bncpostprocess.h"
43#include "bnccore.h"
44#include "bncsettings.h"
45#include "pppopt.h"
46#include "bncpppclient.h"
47#include "rinex/rnxobsfile.h"
48#include "rnxnavfile.h"
49#include "corrfile.h"
50#include "bncsettings.h"
51#include "bncutils.h"
52
53using namespace std;
54
55// Constructor
56////////////////////////////////////////////////////////////////////////////
57t_postProcessing::t_postProcessing(QObject* parent, int maxSpeed) : QThread(parent) {
58
59 _maxSpeed = maxSpeed;
60 _opt = new t_pppOpt();
61 _rnxObsFile = 0;
62 _rnxNavFile = 0;
63 _corrFile = 0;
64 _pppClient = 0;
65 _isToBeDeleted = false;
66
67 bncSettings settings;
68
69 QString outFileName = settings.value("postOutFile").toString();
70 if (outFileName.isEmpty()) {
71 _outFile = 0;
72 _outStream = 0;
73 }
74 else {
75 expandEnvVar(outFileName);
76 _outFile = new QFile(outFileName);
77 _outFile->open(QIODevice::WriteOnly | QIODevice::Text);
78 _outStream = new QTextStream();
79 _outStream->setDevice(_outFile);
80 }
81}
82
83// Destructor
84////////////////////////////////////////////////////////////////////////////
85t_postProcessing::~t_postProcessing() {
86 delete _pppClient;
87 delete _rnxNavFile;
88 delete _rnxObsFile;
89 delete _corrFile;
90 delete _opt;
91 delete _outStream;
92 delete _outFile;
93}
94
95// Write a Program Message
96////////////////////////////////////////////////////////////////////////////
97void t_postProcessing::slotMessage(QByteArray msg, bool /* showOnScreen */) {
98 if (_outStream) {
99 *_outStream << endl << msg;
100 _outStream->flush();
101 }
102 else {
103 BNC_CORE->slotMessage(msg, false);
104 }
105}
106
107// Set Observations from RINEX File
108////////////////////////////////////////////////////////////////////////////
109void t_postProcessing::setObsFromRnx(const t_rnxObsFile* rnxObsFile,
110 const t_rnxObsFile::t_rnxEpo* epo,
111 const t_rnxObsFile::t_rnxSat& rnxSat,
112 t_obs& obs) {
113
114 strncpy(obs.StatID, rnxObsFile->markerName().toAscii().constData(),
115 sizeof(obs.StatID));
116
117 obs.satSys = rnxSat.satSys;
118 obs.satNum = rnxSat.satNum;
119 obs.GPSWeek = epo->tt.gpsw();
120 obs.GPSWeeks = epo->tt.gpssec();
121
122 for (int iType = 0; iType < rnxObsFile->nTypes(obs.satSys); iType++) {
123 QString type = rnxObsFile->obsType(obs.satSys,iType).toAscii();
124 obs.setMeasdata(type, rnxObsFile->version(), rnxSat.obs[iType]);
125 if (type.indexOf("L1") == 0) {
126 obs.snrL1 = rnxSat.snr[iType];
127 obs.slipL1 = (rnxSat.lli[iType] & 1);
128 }
129 else if (type.indexOf("L2") == 0) {
130 obs.snrL2 = rnxSat.snr[iType];
131 obs.slipL2 = (rnxSat.lli[iType] & 1);
132 }
133 else if (type.indexOf("L5") == 0) {
134 obs.snrL5 = rnxSat.snr[iType];
135 obs.slipL5 = (rnxSat.lli[iType] & 1);
136 }
137 }
138}
139
140//
141////////////////////////////////////////////////////////////////////////////
142void t_postProcessing::terminate() {
143 _isToBeDeleted = true;
144 if (!isRunning()) {
145 delete this;
146 }
147}
148
149//
150////////////////////////////////////////////////////////////////////////////
151void t_postProcessing::run() {
152
153 if (_pppClient) {
154 return;
155 }
156 else {
157 try {
158 _rnxObsFile = new t_rnxObsFile(_opt->obsFileName, t_rnxObsFile::input);
159 }
160 catch (...) {
161 delete _rnxObsFile; _rnxObsFile = 0;
162 emit finished();
163 return;
164 }
165 _rnxNavFile = new t_rnxNavFile(_opt->navFileName, t_rnxNavFile::input);
166 _pppClient = new bncPPPclient("POST", _opt, false);
167
168 if (!_opt->corrFileName.isEmpty()) {
169 _corrFile = new t_corrFile(_opt->corrFileName);
170 connect(_corrFile, SIGNAL(newCorrections(QList<QString>)),
171 _pppClient, SLOT(slotNewCorrections(QList<QString>)),
172 Qt::DirectConnection);
173 }
174
175 connect(_pppClient, SIGNAL(newMessage(QByteArray,bool)),
176 this, SLOT(slotMessage(const QByteArray,bool)));
177 qRegisterMetaType<bncTime>("bncTime");
178 connect(_pppClient, SIGNAL(newPosition(bncTime, double, double, double)),
179 this, SIGNAL(newPosition(bncTime, double, double, double)));
180 }
181
182 // Read/Process Observations
183 // -------------------------
184 int nEpo = 0;
185 const t_rnxObsFile::t_rnxEpo* epo = 0;
186 while ( (epo = _rnxObsFile->nextEpoch()) != 0 ) {
187 ++nEpo;
188
189 // Get Corrections
190 // ---------------
191 if (_corrFile) {
192 _corrFile->syncRead(epo->tt);
193 }
194
195 // Get Ephemerides
196 // ----------------
197 t_eph* eph = 0;
198 const QMap<QString, int>* corrIODs = _corrFile ? &_corrFile->corrIODs() : 0;
199 while ( (eph = _rnxNavFile->getNextEph(epo->tt, corrIODs)) != 0 ) {
200 if (_pppClient->putNewEph(eph) != success) {
201 delete eph; eph = 0;
202 }
203 }
204
205 for (unsigned iObs = 0; iObs < epo->rnxSat.size(); iObs++) {
206 if (_isToBeDeleted) {
207 QThread::exit(0);
208 return;
209 }
210
211 const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iObs];
212
213 t_obs obs;
214 t_postProcessing::setObsFromRnx(_rnxObsFile, epo, rnxSat, obs);
215
216 // Get Glonass Channel Number
217 // -------------------------
218 bool obsOK = true;
219 if (obs.satSys == 'R') {
220 QString prn = QString("%1%2").arg(obs.satSys)
221 .arg(obs.satNum, 2, 10, QChar('0'));
222 const bncEphUser::t_ephPair* ephPair = _pppClient->ephPair(prn);
223 if (ephPair && ephPair->last) {
224 obs.slotNum = ((t_ephGlo*) ephPair->last)->slotNum();
225 }
226 else {
227 obsOK = false;
228 }
229 }
230
231 // Put the new Observation to the Client
232 // -------------------------------------
233 if (obsOK) {
234 _pppClient->putNewObs(obs);
235 }
236 }
237 if (nEpo % 10 == 0) {
238 emit progress(nEpo);
239 }
240 }
241
242 if (BNC_CORE->mode() != t_bncCore::interactive) {
243 qApp->exit(0);
244 }
245 else {
246 emit finished();
247 deleteLater();
248 }
249}
250
251//
252////////////////////////////////////////////////////////////////////////////
253void t_postProcessing::slotSetSpeed(int speed) {
254 qDebug() << "speed" << speed;
255}
Note: See TracBrowser for help on using the repository browser.