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

Last change on this file since 4497 was 4403, checked in by mervart, 13 years ago
File size: 6.3 KB
RevLine 
[3615]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 *
[3626]29 * Class: t_postProcessing
[3615]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"
[3684]43#include "bncapp.h"
[3625]44#include "bncsettings.h"
[3642]45#include "pppopt.h"
[3643]46#include "bncpppclient.h"
[3717]47#include "rinex/rnxobsfile.h"
[3646]48#include "rnxnavfile.h"
[3688]49#include "corrfile.h"
[3695]50#include "bncsettings.h"
51#include "bncutils.h"
[3615]52
53using namespace std;
54
[3626]55// Constructor
[3615]56////////////////////////////////////////////////////////////////////////////
[3626]57t_postProcessing::t_postProcessing(QObject* parent) : QThread(parent) {
[3646]58 _opt = new t_pppOpt();
59 _rnxObsFile = 0;
60 _rnxNavFile = 0;
[3688]61 _corrFile = 0;
[3646]62 _pppClient = 0;
[3695]63
64 bncSettings settings;
65
66 QString outFileName = settings.value("postOutFile").toString();
67 if (outFileName.isEmpty()) {
68 _outFile = 0;
69 _outStream = 0;
70 }
71 else {
72 expandEnvVar(outFileName);
73 _outFile = new QFile(outFileName);
74 _outFile->open(QIODevice::WriteOnly | QIODevice::Text);
75 _outStream = new QTextStream();
76 _outStream->setDevice(_outFile);
77 }
[3625]78}
79
[3626]80// Destructor
81////////////////////////////////////////////////////////////////////////////
82t_postProcessing::~t_postProcessing() {
[3646]83 delete _pppClient;
84 delete _rnxNavFile;
85 delete _rnxObsFile;
[3688]86 delete _corrFile;
[3642]87 delete _opt;
[3695]88 delete _outStream;
89 delete _outFile;
[3626]90}
91
[3644]92// Write a Program Message
93////////////////////////////////////////////////////////////////////////////
94void t_postProcessing::slotMessage(QByteArray msg, bool /* showOnScreen */) {
[3695]95 if (_outStream) {
[3696]96 *_outStream << endl << msg;
[3695]97 _outStream->flush();
98 }
99 else {
100 ((bncApp*) qApp)->slotMessage(msg, false);
101 }
[3644]102}
103
[4261]104// Set Observations from RINEX File
105////////////////////////////////////////////////////////////////////////////
106void t_postProcessing::setObsFromRnx(const t_rnxObsFile* rnxObsFile,
107 const t_rnxObsFile::t_rnxEpo* epo,
108 const t_rnxObsFile::t_rnxSat& rnxSat,
109 t_obs& obs) {
110
111 strncpy(obs.StatID, rnxObsFile->markerName().toAscii().constData(),
112 sizeof(obs.StatID));
[4389]113
[4261]114 obs.satSys = rnxSat.satSys;
115 obs.satNum = rnxSat.satNum;
116 obs.GPSWeek = epo->tt.gpsw();
117 obs.GPSWeeks = epo->tt.gpssec();
[4389]118
[4261]119 for (int iType = 0; iType < rnxObsFile->nTypes(obs.satSys); iType++) {
[4391]120 QString type = rnxObsFile->obsType(obs.satSys,iType).toAscii();
[4403]121 obs.setMeasdata(type, rnxObsFile->version(), rnxSat.obs[iType]);
[4389]122 // TOOD: handle slip flags
[4261]123 }
124}
125
[3626]126//
127////////////////////////////////////////////////////////////////////////////
128void t_postProcessing::run() {
[3615]129
[3644]130 if (_pppClient) {
131 return;
132 }
133 else {
[4363]134 try {
135 _rnxObsFile = new t_rnxObsFile(_opt->obsFileName, t_rnxObsFile::input);
136 }
137 catch (...) {
[4364]138 delete _rnxObsFile; _rnxObsFile = 0;
[4363]139 emit finished();
140 deleteLater();
141 return;
142 }
[3999]143 _rnxNavFile = new t_rnxNavFile(_opt->navFileName, t_rnxNavFile::input);
[3646]144 _pppClient = new bncPPPclient("POST", _opt, false);
[3644]145
[3688]146 if (!_opt->corrFileName.isEmpty()) {
147 _corrFile = new t_corrFile(_opt->corrFileName);
148 connect(_corrFile, SIGNAL(newCorrections(QList<QString>)),
[3752]149 _pppClient, SLOT(slotNewCorrections(QList<QString>)),
150 Qt::DirectConnection);
[3688]151 }
[3644]152
153 connect(_pppClient, SIGNAL(newMessage(QByteArray,bool)),
154 this, SLOT(slotMessage(const QByteArray,bool)));
[3643]155 }
156
[3685]157 // Read/Process Observations
158 // -------------------------
159 int nEpo = 0;
[3717]160 const t_rnxObsFile::t_rnxEpo* epo = 0;
[3683]161 while ( (epo = _rnxObsFile->nextEpoch()) != 0 ) {
[3685]162 ++nEpo;
[3689]163
[3748]164 // Get Corrections
165 // ---------------
[3689]166 if (_corrFile) {
167 _corrFile->syncRead(epo->tt);
168 }
169
[3748]170 // Get Ephemerides
171 // ----------------
172 t_eph* eph = 0;
[3757]173 const QMap<QString, int>* corrIODs = _corrFile ? &_corrFile->corrIODs() : 0;
174 while ( (eph = _rnxNavFile->getNextEph(epo->tt, corrIODs)) != 0 ) {
[3748]175 if (_pppClient->putNewEph(eph) != success) {
176 delete eph; eph = 0;
177 }
178 }
179
[3717]180 for (unsigned iObs = 0; iObs < epo->rnxSat.size(); iObs++) {
[4261]181
[3717]182 const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iObs];
[4261]183
[3680]184 t_obs obs;
[4261]185 t_postProcessing::setObsFromRnx(_rnxObsFile, epo, rnxSat, obs);
[3762]186
187 // Get Glonass Channel Number
188 // -------------------------
[4261]189 bool obsOK = true;
[3762]190 if (obs.satSys == 'R') {
191 QString prn = QString("%1%2").arg(obs.satSys)
192 .arg(obs.satNum, 2, 10, QChar('0'));
193 const bncEphUser::t_ephPair* ephPair = _pppClient->ephPair(prn);
194 if (ephPair && ephPair->last) {
195 obs.slotNum = ((t_ephGlo*) ephPair->last)->slotNum();
196 }
[3763]197 else {
198 obsOK = false;
199 }
[3762]200 }
201
202 // Put the new Observation to the Client
203 // -------------------------------------
[3763]204 if (obsOK) {
205 _pppClient->putNewObs(obs);
206 }
[3680]207 }
[3685]208 if (nEpo % 10 == 0) {
209 emit progress(nEpo);
210 }
[3674]211 }
212
[3973]213 bncApp* app = (bncApp*) qApp;
214 if ( app->mode() != bncApp::interactive) {
215 app->exit(0);
216 }
217 else {
218 emit finished();
219 deleteLater();
220 }
[3615]221}
Note: See TracBrowser for help on using the repository browser.