[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"
|
---|
[5070] | 43 | #include "bnccore.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 |
|
---|
| 53 | using namespace std;
|
---|
| 54 |
|
---|
[3626] | 55 | // Constructor
|
---|
[3615] | 56 | ////////////////////////////////////////////////////////////////////////////
|
---|
[5262] | 57 | t_postProcessing::t_postProcessing(QObject* parent, int maxSpeed) : QThread(parent) {
|
---|
| 58 |
|
---|
| 59 | _maxSpeed = maxSpeed;
|
---|
[3646] | 60 | _opt = new t_pppOpt();
|
---|
| 61 | _rnxObsFile = 0;
|
---|
| 62 | _rnxNavFile = 0;
|
---|
[3688] | 63 | _corrFile = 0;
|
---|
[3646] | 64 | _pppClient = 0;
|
---|
[5238] | 65 | _isToBeDeleted = false;
|
---|
[3695] | 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 | }
|
---|
[3625] | 81 | }
|
---|
| 82 |
|
---|
[3626] | 83 | // Destructor
|
---|
| 84 | ////////////////////////////////////////////////////////////////////////////
|
---|
| 85 | t_postProcessing::~t_postProcessing() {
|
---|
[3646] | 86 | delete _pppClient;
|
---|
| 87 | delete _rnxNavFile;
|
---|
| 88 | delete _rnxObsFile;
|
---|
[3688] | 89 | delete _corrFile;
|
---|
[3642] | 90 | delete _opt;
|
---|
[3695] | 91 | delete _outStream;
|
---|
| 92 | delete _outFile;
|
---|
[3626] | 93 | }
|
---|
| 94 |
|
---|
[3644] | 95 | // Write a Program Message
|
---|
| 96 | ////////////////////////////////////////////////////////////////////////////
|
---|
| 97 | void t_postProcessing::slotMessage(QByteArray msg, bool /* showOnScreen */) {
|
---|
[3695] | 98 | if (_outStream) {
|
---|
[3696] | 99 | *_outStream << endl << msg;
|
---|
[3695] | 100 | _outStream->flush();
|
---|
| 101 | }
|
---|
| 102 | else {
|
---|
[5068] | 103 | BNC_CORE->slotMessage(msg, false);
|
---|
[3695] | 104 | }
|
---|
[3644] | 105 | }
|
---|
| 106 |
|
---|
[4261] | 107 | // Set Observations from RINEX File
|
---|
| 108 | ////////////////////////////////////////////////////////////////////////////
|
---|
| 109 | void 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));
|
---|
[4389] | 116 |
|
---|
[4261] | 117 | obs.satSys = rnxSat.satSys;
|
---|
| 118 | obs.satNum = rnxSat.satNum;
|
---|
| 119 | obs.GPSWeek = epo->tt.gpsw();
|
---|
| 120 | obs.GPSWeeks = epo->tt.gpssec();
|
---|
[4389] | 121 |
|
---|
[4261] | 122 | for (int iType = 0; iType < rnxObsFile->nTypes(obs.satSys); iType++) {
|
---|
[4391] | 123 | QString type = rnxObsFile->obsType(obs.satSys,iType).toAscii();
|
---|
[4403] | 124 | obs.setMeasdata(type, rnxObsFile->version(), rnxSat.obs[iType]);
|
---|
[4561] | 125 | if (type.indexOf("L1") == 0) {
|
---|
[4609] | 126 | obs.snrL1 = rnxSat.snr[iType];
|
---|
[4619] | 127 | obs.slipL1 = (rnxSat.lli[iType] & 1);
|
---|
[4561] | 128 | }
|
---|
| 129 | else if (type.indexOf("L2") == 0) {
|
---|
[4609] | 130 | obs.snrL2 = rnxSat.snr[iType];
|
---|
[4619] | 131 | obs.slipL2 = (rnxSat.lli[iType] & 1);
|
---|
[4561] | 132 | }
|
---|
| 133 | else if (type.indexOf("L5") == 0) {
|
---|
[4609] | 134 | obs.snrL5 = rnxSat.snr[iType];
|
---|
[4619] | 135 | obs.slipL5 = (rnxSat.lli[iType] & 1);
|
---|
[4561] | 136 | }
|
---|
[4261] | 137 | }
|
---|
| 138 | }
|
---|
| 139 |
|
---|
[5238] | 140 | //
|
---|
| 141 | ////////////////////////////////////////////////////////////////////////////
|
---|
| 142 | void t_postProcessing::terminate() {
|
---|
| 143 | _isToBeDeleted = true;
|
---|
| 144 | if (!isRunning()) {
|
---|
| 145 | delete this;
|
---|
| 146 | }
|
---|
| 147 | }
|
---|
| 148 |
|
---|
[3626] | 149 | //
|
---|
| 150 | ////////////////////////////////////////////////////////////////////////////
|
---|
| 151 | void t_postProcessing::run() {
|
---|
[3615] | 152 |
|
---|
[3644] | 153 | if (_pppClient) {
|
---|
| 154 | return;
|
---|
| 155 | }
|
---|
| 156 | else {
|
---|
[4363] | 157 | try {
|
---|
| 158 | _rnxObsFile = new t_rnxObsFile(_opt->obsFileName, t_rnxObsFile::input);
|
---|
| 159 | }
|
---|
| 160 | catch (...) {
|
---|
[4364] | 161 | delete _rnxObsFile; _rnxObsFile = 0;
|
---|
[4363] | 162 | emit finished();
|
---|
| 163 | return;
|
---|
| 164 | }
|
---|
[3999] | 165 | _rnxNavFile = new t_rnxNavFile(_opt->navFileName, t_rnxNavFile::input);
|
---|
[3646] | 166 | _pppClient = new bncPPPclient("POST", _opt, false);
|
---|
[3644] | 167 |
|
---|
[3688] | 168 | if (!_opt->corrFileName.isEmpty()) {
|
---|
| 169 | _corrFile = new t_corrFile(_opt->corrFileName);
|
---|
| 170 | connect(_corrFile, SIGNAL(newCorrections(QList<QString>)),
|
---|
[3752] | 171 | _pppClient, SLOT(slotNewCorrections(QList<QString>)),
|
---|
| 172 | Qt::DirectConnection);
|
---|
[3688] | 173 | }
|
---|
[3644] | 174 |
|
---|
| 175 | connect(_pppClient, SIGNAL(newMessage(QByteArray,bool)),
|
---|
| 176 | this, SLOT(slotMessage(const QByteArray,bool)));
|
---|
[5215] | 177 | qRegisterMetaType<bncTime>("bncTime");
|
---|
[5214] | 178 | connect(_pppClient, SIGNAL(newPosition(bncTime, double, double, double)),
|
---|
| 179 | this, SIGNAL(newPosition(bncTime, double, double, double)));
|
---|
[3643] | 180 | }
|
---|
| 181 |
|
---|
[3685] | 182 | // Read/Process Observations
|
---|
| 183 | // -------------------------
|
---|
| 184 | int nEpo = 0;
|
---|
[3717] | 185 | const t_rnxObsFile::t_rnxEpo* epo = 0;
|
---|
[3683] | 186 | while ( (epo = _rnxObsFile->nextEpoch()) != 0 ) {
|
---|
[3685] | 187 | ++nEpo;
|
---|
[3689] | 188 |
|
---|
[3748] | 189 | // Get Corrections
|
---|
| 190 | // ---------------
|
---|
[3689] | 191 | if (_corrFile) {
|
---|
| 192 | _corrFile->syncRead(epo->tt);
|
---|
| 193 | }
|
---|
| 194 |
|
---|
[3748] | 195 | // Get Ephemerides
|
---|
| 196 | // ----------------
|
---|
| 197 | t_eph* eph = 0;
|
---|
[3757] | 198 | const QMap<QString, int>* corrIODs = _corrFile ? &_corrFile->corrIODs() : 0;
|
---|
| 199 | while ( (eph = _rnxNavFile->getNextEph(epo->tt, corrIODs)) != 0 ) {
|
---|
[3748] | 200 | if (_pppClient->putNewEph(eph) != success) {
|
---|
| 201 | delete eph; eph = 0;
|
---|
| 202 | }
|
---|
| 203 | }
|
---|
| 204 |
|
---|
[3717] | 205 | for (unsigned iObs = 0; iObs < epo->rnxSat.size(); iObs++) {
|
---|
[5238] | 206 | if (_isToBeDeleted) {
|
---|
| 207 | QThread::exit(0);
|
---|
| 208 | return;
|
---|
| 209 | }
|
---|
[4261] | 210 |
|
---|
[3717] | 211 | const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iObs];
|
---|
[4261] | 212 |
|
---|
[3680] | 213 | t_obs obs;
|
---|
[4261] | 214 | t_postProcessing::setObsFromRnx(_rnxObsFile, epo, rnxSat, obs);
|
---|
[3762] | 215 |
|
---|
| 216 | // Get Glonass Channel Number
|
---|
| 217 | // -------------------------
|
---|
[4261] | 218 | bool obsOK = true;
|
---|
[3762] | 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 | }
|
---|
[3763] | 226 | else {
|
---|
| 227 | obsOK = false;
|
---|
| 228 | }
|
---|
[3762] | 229 | }
|
---|
| 230 |
|
---|
| 231 | // Put the new Observation to the Client
|
---|
| 232 | // -------------------------------------
|
---|
[3763] | 233 | if (obsOK) {
|
---|
| 234 | _pppClient->putNewObs(obs);
|
---|
| 235 | }
|
---|
[3680] | 236 | }
|
---|
[3685] | 237 | if (nEpo % 10 == 0) {
|
---|
| 238 | emit progress(nEpo);
|
---|
| 239 | }
|
---|
[3674] | 240 | }
|
---|
| 241 |
|
---|
[5072] | 242 | if (BNC_CORE->mode() != t_bncCore::interactive) {
|
---|
[5066] | 243 | qApp->exit(0);
|
---|
[3973] | 244 | }
|
---|
| 245 | else {
|
---|
| 246 | emit finished();
|
---|
| 247 | deleteLater();
|
---|
| 248 | }
|
---|
[3615] | 249 | }
|
---|
[5262] | 250 |
|
---|
| 251 | //
|
---|
| 252 | ////////////////////////////////////////////////////////////////////////////
|
---|
| 253 | void t_postProcessing::slotSetSpeed(int speed) {
|
---|
| 254 | qDebug() << "speed" << speed;
|
---|
| 255 | }
|
---|