source: ntrip/trunk/BNC/src/pppRun.cpp@ 11179

// Part of BNC, a utility for retrieving decoding and
// converting GNSS data streams from NTRIP broadcasters.
//
// Copyright (C) 2007
// German Federal Agency for Cartography and Geodesy (BKG)
// http://bkg.bund.de
// Czech Technical University Prague, Department of Geodesy
// http://www.fsv.cvut.cz
//
// Email: euref-ip@bkg.bund.de
//
// This program is free software; you can redistribute it and/or
// modify it under the terms of the GNU General Public License
// as published by the Free Software Foundation, version 2.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.

/* -------------------------------------------------------------------------
 * BKG NTRIP Client
 * -------------------------------------------------------------------------
 *
 * Class:      t_pppRun
 *
 * Purpose:    Single Real-Time PPP Client
 *
 * Author:     L. Mervart
 *
 * Created:    29-Jul-2014
 *
 * Changes:
 *
 * -----------------------------------------------------------------------*/


#include <iostream>
#include <iomanip>
#include <sstream>
#include <string.h>
#include <map>

#include "pppRun.h"
#include "pppThread.h"
#include "bnccore.h"
#include "bncephuser.h"
#include "bncsettings.h"
#include "bncoutf.h"
#include "bncsinextro.h"
#include "rinex/rnxobsfile.h"
#include "rinex/rnxnavfile.h"
#include "rinex/corrfile.h"
#include "combination/bnccomb.h"

using namespace BNC_PPP;
using namespace std;

// Constructor
////////////////////////////////////////////////////////////////////////////
t_pppRun::t_pppRun(const t_pppOptions* opt) {

  _opt = opt;

  connect(this, SIGNAL(newMessage(QByteArray,bool)),
          BNC_CORE, SLOT(slotMessage(const QByteArray,bool)));
  connect(this, SIGNAL(newLogRecord(t_logRecord,bool)),
          BNC_CORE, SLOT(slotLogRecord(t_logRecord,bool)));
  t_logRecord::registerStation(QByteArray(_opt->_roverName.c_str()));

  connect(this,     SIGNAL(newPosition(QByteArray, bncTime, QVector<double>)),
          BNC_CORE, SIGNAL(newPosition(QByteArray, bncTime, QVector<double>)));

  connect(this,     SIGNAL(newNMEAstr(QByteArray, QByteArray)),
          BNC_CORE, SIGNAL(newNMEAstr(QByteArray, QByteArray)));

  _pppClient = new t_pppClient(_opt);

  bncSettings settings;

  if (_opt->_realTime) {
    _ephMounts = QString::fromStdString(_opt->_ephMount).split(QRegularExpression("[,\\s]+"),
                                                                Qt::SkipEmptyParts);
    if (!_ephMounts.isEmpty()) {
      emit newMessage(QByteArray(_opt->_roverName.c_str()) + ": PPP uses the ephemerides of "
                      + _ephMounts.join(", ").toLatin1() + " only", true);
    }

    Qt::ConnectionType conType = Qt::AutoConnection;
    if (BNC_CORE->mode() == t_bncCore::batchPostProcessing) {
      conType = Qt::BlockingQueuedConnection;
    }

    connect(BNC_CORE->caster(), SIGNAL(newObs(QByteArray, QList<t_satObs>)),
            this, SLOT(slotNewObs(QByteArray, QList<t_satObs>)),conType);

    connect(BNC_CORE, SIGNAL(newGPSEph(t_ephGPS)),
            this, SLOT(slotNewGPSEph(t_ephGPS)),conType);

    connect(BNC_CORE, SIGNAL(newGlonassEph(t_ephGlo)),
            this, SLOT(slotNewGlonassEph(t_ephGlo)),conType);

    connect(BNC_CORE, SIGNAL(newGalileoEph(t_ephGal)),
            this, SLOT(slotNewGalileoEph(t_ephGal)),conType);

    connect(BNC_CORE, SIGNAL(newBDSEph(t_ephBDS)),
            this, SLOT(slotNewBDSEph(t_ephBDS)),conType);

    connect(BNC_CORE, SIGNAL(newTec(t_vTec)),
            this, SLOT(slotNewTec(t_vTec)),conType);

    connect(BNC_CORE, SIGNAL(newMetaData(t_metaData)),
            this, SLOT(slotNewMetaData(t_metaData)),conType);

    connect(BNC_CORE, SIGNAL(newOrbCorrections(QList<t_orbCorr>)),
            this, SLOT(slotNewOrbCorrections(QList<t_orbCorr>)),conType);

    connect(BNC_CORE, SIGNAL(newClkCorrections(QList<t_clkCorr>)),
            this, SLOT(slotNewClkCorrections(QList<t_clkCorr>)),conType);

    connect(BNC_CORE, SIGNAL(newCodeBiases(QList<t_satCodeBias>)),
            this, SLOT(slotNewCodeBiases(QList<t_satCodeBias>)),conType);

    connect(BNC_CORE, SIGNAL(newPhaseBiases(QList<t_satPhaseBias>)),
            this, SLOT(slotNewPhaseBiases(QList<t_satPhaseBias>)),conType);

    connect(BNC_CORE, SIGNAL(newSatAntennas(QList<t_satAntenna>)),
            this, SLOT(slotNewSatAntennas(QList<t_satAntenna>)),conType);

    connect(BNC_CMB, SIGNAL(newOrbCorrections(QList<t_orbCorr>)),
            this, SLOT(slotNewOrbCorrections(QList<t_orbCorr>)),conType);

    connect(BNC_CMB, SIGNAL(newClkCorrections(QList<t_clkCorr>)),
            this, SLOT(slotNewClkCorrections(QList<t_clkCorr>)),conType);

    connect(BNC_CMB, SIGNAL(newCodeBiases(QList<t_satCodeBias>)),
            this, SLOT(slotNewCodeBiases(QList<t_satCodeBias>)),conType);

    connect(BNC_CORE, SIGNAL(providerIDChanged(QString)),
            this, SLOT(slotProviderIDChanged(QString)));
  }
  else {
    _rnxObsFile = 0;
    _rnxNavFile = 0;
    _corrFile   = 0;
    _biasFile   = 0;
    _ionoFile   = 0;
    _speed      = settings.value("PPP/mapSpeedSlider").toInt();
    connect(this, SIGNAL(progressRnxPPP(int)), BNC_CORE, SIGNAL(progressRnxPPP(int)));
    connect(this, SIGNAL(finishedRnxPPP()),    BNC_CORE, SIGNAL(finishedRnxPPP()));
    connect(BNC_CORE, SIGNAL(mapSpeedSliderChanged(int)),
            this, SLOT(slotSetSpeed(int)));
    connect(BNC_CORE, SIGNAL(stopRinexPPP()), this, SLOT(slotSetStopFlag()));
  }

  _stopFlag = false;

  QString roverName(_opt->_roverName.c_str());
  QString country;
  QString monNum = "0";
  QString recNum = "0";
  if (roverName.length() >= 9) {
    monNum = QChar(roverName[4]);
    recNum = QChar(roverName[5]);
  }
  QString intr = "1 day";
  int     sampl  = 0;
  QString distStr;
  int num = 0;
  int statIDlength = roverName.size() -1;
  QString ID4 = roverName.left(4);
  QString ID9 = roverName.left(9);
  ID4 = ID4.toUpper();
  ID4 = ID9.toUpper();
  QListIterator<QString> it(settings.value("mountPoints").toStringList());
  while (it.hasNext()) {
    QString mp = it.next();
    if (mp.indexOf(roverName.left(statIDlength)) != -1) {
      ++num;
    }
    QStringList hlp = mp.split(" ");
    if (hlp.size() < 7)
      continue;
    if (hlp.join(" ").indexOf(roverName, 0) != -1) {
      country = hlp[2];
    }
  }

  if (num > 1) {
    distStr = "." + roverName.right(1);
  }

  if (ID9.size() < 9) {
    ID9 = ID4
        + QString("%1").arg(monNum, 1, 10)
        + QString("%1").arg(recNum, 1, 10)
        +  country;
  }
  QString logFileSkl = settings.value("PPP/logPath").toString();
  int l = logFileSkl.length();
  if (logFileSkl.isEmpty()) {
    _logFile = 0;
  }
  else {
    if (l && logFileSkl[l-1] != QDir::separator() ) {
      logFileSkl += QDir::separator();
    }

    logFileSkl = logFileSkl + ID9 + "${V3PROD}" + distStr + ".ppp";

    _logFile = new bncoutf(logFileSkl, intr, sampl);
  }

  QString nmeaFileSkl = settings.value("PPP/nmeaPath").toString();
  l = nmeaFileSkl.length();
  if (nmeaFileSkl.isEmpty()) {
    _nmeaFile = 0;
  }
  else {
    if (l > 0 && nmeaFileSkl[l-1] != QDir::separator() ) {
      nmeaFileSkl += QDir::separator();
    }
    nmeaFileSkl = nmeaFileSkl + ID9 + "${V3PROD}" + distStr + ".nmea";

    _nmeaFile = new bncoutf(nmeaFileSkl, intr, sampl);
  }
  QString snxtroFileSkl = settings.value("PPP/snxtroPath").toString();
  l = snxtroFileSkl.length();
  if (snxtroFileSkl.isEmpty()) {
    _snxtroFile = 0;
  }
  else {
    if (l > 0 && snxtroFileSkl[l-1] != QDir::separator() ) {
      snxtroFileSkl += QDir::separator();
    }
    snxtroFileSkl = snxtroFileSkl + ID9 + '_' + "${V3PROD}" + distStr + ".TRO";
    sampl = settings.value("PPP/snxtroSampl").toString().split("sec").first().toInt();
    intr  = settings.value("PPP/snxtroIntr").toString();
    _snxtroFile = new bncSinexTro(_opt, snxtroFileSkl, intr, sampl);
  }
}

// Destructor
////////////////////////////////////////////////////////////////////////////
t_pppRun::~t_pppRun() {
  delete _logFile;
  delete _nmeaFile;
  delete _snxtroFile;
  while (!_epoData.empty()) {
    delete _epoData.front();
    _epoData.pop_front();
  }
}

// Ephemeris from a stream selected for PPP (PPP/ephMount)?
////////////////////////////////////////////////////////////////////////////
bool t_pppRun::useEph(const t_eph& eph) const {
  return _ephMounts.isEmpty() || _ephMounts.contains(eph.receptStaID());
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewGPSEph(t_ephGPS eph) {
  if (!useEph(eph)) {
    return;
  }
  QMutexLocker locker(&_mutex);
  _pppClient->putEphemeris(&eph);
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewGlonassEph(t_ephGlo eph) {
  if (!useEph(eph)) {
    return;
  }
  QMutexLocker locker(&_mutex);
  _pppClient->putEphemeris(&eph);
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewGalileoEph(t_ephGal eph) {
  if (!useEph(eph)) {
    return;
  }
  QMutexLocker locker(&_mutex);
  _pppClient->putEphemeris(&eph);
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewBDSEph(t_ephBDS eph) {
  if (!useEph(eph)) {
    return;
  }
  QMutexLocker locker(&_mutex);
  _pppClient->putEphemeris(&eph);
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewObs(QByteArray staID, QList<t_satObs> obsList) {
  QMutexLocker locker(&_mutex);

  //// beg test
  // if (obsList.size() > 0) {
  //   const double sampling = 5.0;
  //   const bncTime& tt = obsList[0]._time;
  //   double mod = remainder(tt.gpssec(), sampling);
  //   if (fabs(mod) > 0.05) {
  //     return;
  //   }
  // }
  //// end test

  if (string(staID.data()) != _opt->_roverName) {
    return;
  }
/*
  int    currentWeek = 0;
  double currentSec  = 0.0;
  currentGPSWeeks(currentWeek, currentSec);
  bncTime currentTime(currentWeek, currentSec);
*/
  // Loop over all observations (possible different epochs)
  // -----------------------------------------------------
  QListIterator<t_satObs> it(obsList);
  while (it.hasNext()) {
    const t_satObs& oldObs = it.next();
    t_satObs*       newObs = new t_satObs(oldObs);

    /* Check regarding current time
    // ----------------------------
    if (OPT->_realTime && BNC_CORE->mode() != t_bncCore::batchPostProcessing ) {
      if ((newObs->_time >= currentTime) ||        // future time stamp
          (currentTime - newObs->_time) > 60.0) {  // very old data sets
        delete newObs;
        continue;
      }
    }
    */
    // Find the corresponding data epoch or create a new one
    // -----------------------------------------------------
    t_epoData* epoch = 0;
    deque<t_epoData*>::const_iterator it;
    for (it = _epoData.begin(); it != _epoData.end(); it++) {
      if (newObs->_time == (*it)->_time) {
        epoch = *it;
        break;
      }
    }
    if (epoch == 0) {
      if (_epoData.empty() || newObs->_time > _epoData.back()->_time) {
        epoch = new t_epoData;
        epoch->_time = newObs->_time;
        _epoData.push_back(epoch);
      }
    }

    // Put data into the epoch
    // -----------------------
    if (epoch != 0) {
      epoch->_satObs.push_back(newObs);
    }
    else {
      delete newObs;
    }
  }

  // Make sure the buffer does not grow beyond any limit
  // ---------------------------------------------------
  const unsigned MAX_EPODATA_SIZE = 120;
  if (_epoData.size() > MAX_EPODATA_SIZE) {
    delete _epoData.front();
    _epoData.pop_front();
  }

  // Process the oldest epochs
  // ------------------------
  while (_epoData.size()) {

    const vector<t_satObs*>& satObs = _epoData.front()->_satObs;

    // No corrections yet, skip the epoch
    // ----------------------------------
    if (_opt->_corrWaitTime && !_lastClkCorrTime.valid()) {
      return;
    }

    // Process the front epoch
    // -----------------------
    if (_opt->_corrWaitTime == 0 ||
        (_epoData.front()->_time - _lastClkCorrTime) < _opt->_corrWaitTime) {

      if (_opt->_corrWaitTime != 0 &&
          (fabs(_epoData.front()->_time - _lastClkCorrTime) > 60.0)) {
        delete _epoData.front();
        _epoData.pop_front();
        continue;
      }

      t_output output;
      _pppClient->processEpoch(satObs, &output);

      if (!output._error) {
        QVector<double> xx(6);
        xx.data()[0] = output._xyzRover[0];
        xx.data()[1] = output._xyzRover[1];
        xx.data()[2] = output._xyzRover[2];
        xx.data()[3] = output._neu[0];
        xx.data()[4] = output._neu[1];
        xx.data()[5] = output._neu[2];
        emit newPosition(staID, output._epoTime, xx);
      }

      delete _epoData.front();
      _epoData.pop_front();

      ostringstream log;
      if (output._error) {
        log << output._log;
      }
      else {
        log.setf(ios::fixed);
        log << string(output._epoTime) << ' ' << staID.data()
            << " X = "  << setprecision(4) << output._xyzRover[0]
            << " Y = "  << setprecision(4) << output._xyzRover[1]
            << " Z = "  << setprecision(4) << output._xyzRover[2]
            << " NEU: " << showpos << setw(8) << setprecision(4) << output._neu[0]
            << " "      << showpos << setw(8) << setprecision(4) << output._neu[1]
            << " "      << showpos << setw(8) << setprecision(4) << output._neu[2]
            << " TRP: " << showpos << setw(8) << setprecision(4) << output._trp0
            << " "      << showpos << setw(8) << setprecision(4) << output._trp;
        if (output._fixRatio > 0.0) {
          log << " fix " <<  int(100*output._fixRatio) << " %";
        }
        else {
          log << " flt ";
        }

      }

      if (_logFile && output._epoTime.valid()) {
          _logFile->write(output._epoTime.gpsw(), output._epoTime.gpssec(),
                        QString(output._log.c_str()));
      }

      if (!output._error) {
        QString rmcStr = nmeaString('R', output);
        QString ggaStr = nmeaString('G', output);
        if (_nmeaFile) {
          _nmeaFile->write(output._epoTime.gpsw(), output._epoTime.gpssec(), rmcStr);
          _nmeaFile->write(output._epoTime.gpsw(), output._epoTime.gpssec(), ggaStr);
        }
        emit newNMEAstr(staID, rmcStr.toLatin1());
        emit newNMEAstr(staID, ggaStr.toLatin1());
        if (_snxtroFile && output._epoTime.valid()) {
          _snxtroFile->write(staID, int(output._epoTime.gpsw()), output._epoTime.gpssec(),
                      output._trp0 + output._trp, output._trpStdev);
        }
      }
      if (output._error) {
        emit newLogRecord(t_logRecord(t_logRecord::warning, "ppp", "ppp_epoch_failed",
                                      staID, QByteArray(log.str().c_str())), true);
      }
      else {
        t_logRecord rec(t_logRecord::info, "ppp", "ppp_position", staID,
                        QByteArray(log.str().c_str()));
        rec.add("epoch", QByteArray((output._epoTime.datestr() + 'T' +
                                     output._epoTime.timestr(3)).c_str()))
           .add("x", output._xyzRover[0], 4)
           .add("y", output._xyzRover[1], 4)
           .add("z", output._xyzRover[2], 4)
           .add("dn", output._neu[0], 4)
           .add("de", output._neu[1], 4)
           .add("du", output._neu[2], 4)
           .add("trp0", output._trp0, 4)
           .add("trp", output._trp, 4)
           .add("trp_stdev", output._trpStdev, 4)
           .add("num_sat", output._numSat)
           .add("hdop", output._hDop, 2)
           .add("fixed", output._fixRatio > 0.0)
           .add("fix_percent", int(100*output._fixRatio));
        emit newLogRecord(rec, true);
      }
    }
    else {
      return;
    }
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewTec(t_vTec vTec) {
  if (vTec._layers.size() == 0) {
    return;
  }

  if (_opt->_realTime) {
    if (_opt->_ionoMount.empty() && _opt->_corrMount.empty()) {
      return;
    }
    if      (!_opt->_ionoMount.empty()) {
      if (_opt->_ionoMount != vTec._staID) {
        return;
      }
    }
    else if (!_opt->_corrMount.empty()) {
      if (_opt->_corrMount != vTec._staID) {
        return;
      }
    }
  }
  _pppClient->putTec(&vTec);
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewMetaData(t_metaData metaData) {
  if (metaData._entries.size() == 0) {
    return;
  }

  if (_opt->_realTime) {
    if (_opt->_corrMount.empty() || _opt->_corrMount != metaData._staID) {
      return;
    }
  }
  _pppClient->putMetaData(&metaData);
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewOrbCorrections(QList<t_orbCorr> orbCorr) {
  if (orbCorr.size() == 0) {
    return;
  }

  if (_opt->_realTime) {
    if (_opt->_corrMount.empty() || _opt->_corrMount != orbCorr[0]._staID) {
      return;
    }
  }
  vector<t_orbCorr*> corrections;
  for (int ii = 0; ii < orbCorr.size(); ii++) {
    corrections.push_back(new t_orbCorr(orbCorr[ii]));
  }

  _pppClient->putOrbCorrections(corrections);

  for (unsigned ii = 0; ii < corrections.size(); ii++) {
    delete corrections[ii];
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewClkCorrections(QList<t_clkCorr> clkCorr) {
  if (clkCorr.size() == 0) {
    return;
  }

  if (_opt->_realTime) {
    if (_opt->_corrMount.empty() || _opt->_corrMount != clkCorr[0]._staID) {
      return;
    }
  }
  vector<t_clkCorr*> corrections;
  for (int ii = 0; ii < clkCorr.size(); ii++) {
    corrections.push_back(new t_clkCorr(clkCorr[ii]));
    _lastClkCorrTime = clkCorr[ii]._time;
  }
  _pppClient->putClkCorrections(corrections);

  for (unsigned ii = 0; ii < corrections.size(); ii++) {
    delete corrections[ii];
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewCodeBiases(QList<t_satCodeBias> codeBiases) {
  if (codeBiases.size() == 0) {
    return;
  }

  if (_opt->_realTime) {
    if (_opt->_corrMount.empty() && _opt->_biasMount.empty()) {
      return;
    }
    if      (!_opt->_biasMount.empty()) {
      if (_opt->_biasMount != codeBiases[0]._staID) {
        return;
      }
    }
    else if (!_opt->_corrMount.empty()) {
      if (_opt->_corrMount != codeBiases[0]._staID) {
        return;
      }
    }
  }
  vector<t_satCodeBias*> biases;
  for (int ii = 0; ii < codeBiases.size(); ii++) {
    biases.push_back(new t_satCodeBias(codeBiases[ii]));
  }

  _pppClient->putCodeBiases(biases);

  for (unsigned ii = 0; ii < biases.size(); ii++) {
    delete biases[ii];
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewSatAntennas(QList<t_satAntenna> satAntennas) {
  if (satAntennas.size() == 0) {
    return;
  }

  // Taken from the same mountpoint as the Metadata (see slotNewMetaData()),
  // since only Metadata of the same provider can reference this antenna
  // information (IOD unique per SSR Provider ID only).
  if (_opt->_realTime) {
    if (_opt->_corrMount.empty() || _opt->_corrMount != satAntennas[0]._staID) {
      return;
    }
  }
  vector<t_satAntenna*> antennas;
  for (int ii = 0; ii < satAntennas.size(); ii++) {
    antennas.push_back(new t_satAntenna(satAntennas[ii]));
  }

  _pppClient->putSatAntennas(antennas);

  for (unsigned ii = 0; ii < antennas.size(); ii++) {
    delete antennas[ii];
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotNewPhaseBiases(QList<t_satPhaseBias> phaseBiases) {
  if (phaseBiases.size() == 0) {
    return;
  }

  if (_opt->_realTime) {
    if (_opt->_corrMount.empty() && _opt->_biasMount.empty()) {
      return;
    }
    if      (!_opt->_biasMount.empty()) {
      if (_opt->_biasMount != phaseBiases[0]._staID) {
        return;
      }
    }
    else if (!_opt->_corrMount.empty()) {
      if (_opt->_corrMount != phaseBiases[0]._staID) {
        return;
      }
    }
  }
  vector<t_satPhaseBias*> biases;
  for (int ii = 0; ii < phaseBiases.size(); ii++) {
    biases.push_back(new t_satPhaseBias(phaseBiases[ii]));
  }

  _pppClient->putPhaseBiases(biases);

  for (unsigned ii = 0; ii < biases.size(); ii++) {
    delete biases[ii];
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::processFiles() {

  try {
    _rnxObsFile = new t_rnxObsFile(QString(_opt->_rinexObs.c_str()), t_rnxObsFile::input);
  }
  catch (...) {
    delete _rnxObsFile; _rnxObsFile = 0;
    emit finishedRnxPPP();
    return;
  }

  try {
    _rnxNavFile = new t_rnxNavFile(QString(_opt->_rinexNav.c_str()), t_rnxNavFile::input);
  }
  catch (...) {
    delete _rnxNavFile; _rnxNavFile = 0;
    emit finishedRnxPPP();
    return;
  }

  if (!_opt->_corrFile.empty()) {
    _corrFile = new t_corrFile(QString(_opt->_corrFile.c_str()));
    connect(_corrFile, SIGNAL(newOrbCorrections(QList<t_orbCorr>)), this, SLOT(slotNewOrbCorrections(QList<t_orbCorr>)));
    connect(_corrFile, SIGNAL(newClkCorrections(QList<t_clkCorr>)), this, SLOT(slotNewClkCorrections(QList<t_clkCorr>)));
    connect(_corrFile, SIGNAL(newMetaData(t_metaData)), this, SLOT(slotNewMetaData(t_metaData)));
    if (_opt->_biasFile.empty()) {
      connect(_corrFile, SIGNAL(newCodeBiases(QList<t_satCodeBias>)), this, SLOT(slotNewCodeBiases(QList<t_satCodeBias>)));
      connect(_corrFile, SIGNAL(newPhaseBiases(QList<t_satPhaseBias>)), this, SLOT(slotNewPhaseBiases(QList<t_satPhaseBias>)));
    }
    connect(_corrFile, SIGNAL(newSatAntennas(QList<t_satAntenna>)), this, SLOT(slotNewSatAntennas(QList<t_satAntenna>)));
    if (_opt->_ionoFile.empty()) {
      connect(_corrFile, SIGNAL(newTec(t_vTec)), this, SLOT(slotNewTec(t_vTec)));
    }
  }

  if (!_opt->_biasFile.empty()) {
    _biasFile = new t_corrFile(QString(_opt->_biasFile.c_str()));
    connect(_biasFile, SIGNAL(newCodeBiases(QList<t_satCodeBias>)), this, SLOT(slotNewCodeBiases(QList<t_satCodeBias>)));
    connect(_biasFile, SIGNAL(newPhaseBiases(QList<t_satPhaseBias>)), this, SLOT(slotNewPhaseBiases(QList<t_satPhaseBias>)));
  }

  if (!_opt->_ionoFile.empty()) {
    _ionoFile = new t_corrFile(QString(_opt->_corrFile.c_str()));
    connect(_ionoFile, SIGNAL(newTec(t_vTec)), this, SLOT(slotNewTec(t_vTec)));
  }

  // Read/Process Observations
  // -------------------------
  int   nEpo = 0;
  const t_rnxObsFile::t_rnxEpo* epo = 0;
  while ( !_stopFlag && (epo = _rnxObsFile->nextEpoch()) != 0 ) {
    ++nEpo;

    if (_speed < 100) {
      double sleepTime = 2.0 / _speed;
      t_pppThread::msleep(int(sleepTime*1.e3));
    }

    // Get Corrections
    // ---------------
    for (t_corrFile* file : {_corrFile, _biasFile, _ionoFile}) {
      if (file) {
        try {
          file->syncRead(epo->tt);
        }
        catch (const char* msg) {
          emit newMessage(QByteArray(msg), true);
          break;
        }
        catch (const string& msg) {
          emit newMessage(QByteArray(msg.c_str()), true);
          break;
        }
        catch (...) {
          emit newMessage("unknown exceptions in file " + file->name().toLatin1(), true);
          break;
        }
      }
    }

    // Get Ephemerides
    // ----------------
    t_eph* eph = 0;
    const QMap<QString, unsigned int>* corrIODs = _corrFile ? &_corrFile->corrIODs() : 0;
    while ( (eph = _rnxNavFile->getNextEph(epo->tt, corrIODs)) != 0 ) {
      _pppClient->putEphemeris(eph);
      delete eph; eph = 0;
    }

    // Create list of observations and start epoch processing
    // ------------------------------------------------------
    QList<t_satObs> obsList;
    for (unsigned iObs = 0; iObs < epo->rnxSat.size(); iObs++) {
      const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iObs];

      t_satObs obs;
      t_rnxObsFile::setObsFromRnx(_rnxObsFile, epo, rnxSat, obs);
      obsList << obs;
    }
    slotNewObs(QByteArray(_opt->_roverName.c_str()), obsList);


    if (nEpo % 10 == 0) {
      emit progressRnxPPP(nEpo);
    }

    QCoreApplication::processEvents();
  }

  emit finishedRnxPPP();

  if (BNC_CORE->mode() != t_bncCore::interactive) {
    qApp->exit(6);
  }
  else {
    BNC_CORE->stopPPP();
  }
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotSetSpeed(int speed) {
  QMutexLocker locker(&_mutex);
  _speed = speed;
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotSetStopFlag() {
  QMutexLocker locker(&_mutex);
  _stopFlag = true;
}

//
////////////////////////////////////////////////////////////////////////////
QString t_pppRun::nmeaString(char strType, const t_output& output) {

  double ell[3];
  xyz2ell(output._xyzRover, ell);
  double phiDeg = ell[0] * 180 / M_PI;
  double lamDeg = ell[1] * 180 / M_PI;

  unsigned year, month, day;
  output._epoTime.civil_date(year, month, day);
  double gps_utc = gnumleap(year, month, day);

  char phiCh = 'N';
  if (phiDeg < 0) {
    phiDeg = -phiDeg;
    phiCh  =  'S';
  }
  char lamCh = 'E';
  if (lamDeg < 0) {
    lamDeg = -lamDeg;
    lamCh  =  'W';
  }

  ostringstream out;
  out.setf(ios::fixed);

  if      (strType == 'R') {
    string datestr = (output._epoTime - gps_utc).datestr(0); // yyyymmdd
    out << "GPRMC,"
        << (output._epoTime - gps_utc).timestr(3,0) << ",A,"
        << setw(2) << setfill('0') << int(phiDeg)
        << setw(7) << setprecision(4) << setfill('0') << fmod(60*phiDeg,60) << ',' << phiCh << ','
        << setw(3) << setfill('0') << int(lamDeg)
        << setw(7) << setprecision(4) << setfill('0') << fmod(60*lamDeg,60) << ',' << lamCh << ",,,"
        << datestr[6] << datestr[7] << datestr[4] << datestr[5] << datestr[2] << datestr[3] << ",,";
  }
  else if (strType == 'G') {
    out << "GPGGA,"
        << (output._epoTime - gps_utc).timestr(2,0) << ','
        << setw(2) << setfill('0') << int(phiDeg)
        << setw(10) << setprecision(7) << setfill('0')
        << fmod(60*phiDeg,60) << ',' << phiCh << ','
        << setw(3) << setfill('0') << int(lamDeg)
        << setw(10) << setprecision(7) << setfill('0')
        << fmod(60*lamDeg,60) << ',' << lamCh
        << ",1," << setw(2) << setfill('0') << output._numSat << ','
        << setw(3) << setprecision(1) << output._hDop << ','
        << setprecision(3) << ell[2] << ",M,0.0,M,,";
  }
  else {
    return "";
  }

  QString nmStr(out.str().c_str());
  unsigned char XOR = 0;
  for (int ii = 0; ii < nmStr.length(); ii++) {
    XOR ^= (unsigned char) nmStr[ii].toLatin1();
  }

  return '$' + nmStr + QString("*%1").arg(XOR, 2, 16, QLatin1Char('0')) + "\n";
}

//
////////////////////////////////////////////////////////////////////////////
void t_pppRun::slotProviderIDChanged(QString mountPoint) {
  QMutexLocker locker(&_mutex);

  if (mountPoint.toStdString() !=_opt->_corrMount) {
    return;
  }

  QString msg = "pppRun " + QString(_opt->_roverName.c_str()) + ": SSR Provider or Service changed: " + mountPoint;
  emit newMessage(msg.toLatin1(), true);

  _pppClient->reset();

  while (!_epoData.empty()) {
    delete _epoData.front();
    _epoData.pop_front();
  }
}