Changeset 7943 in ntrip for trunk

Timestamp:

Jun 2, 2016, 10:47:19 AM (8 years ago)

Author:

stuerze

Message:

the approach how to wait for clock corrections in PPP mode, which was well proven in BNC verson 2.11 is now re-implemented

Location:

trunk/BNC/src

Files:

• bnchelp.html (modified) (1 diff)
• pppMain.cpp (modified) (1 diff)
• pppRun.cpp (modified) (2 diffs)
• pppRun.h (modified) (1 diff)
• pppWidgets.cpp (modified) (2 diffs)

trunk/BNC/src/bnchelp.html

@@ -3584 +3584 @@
 <p><h4>2.13.3.7 <a name="pppwaitclockcorr">Wait for Clock Corrections - optional</h4></p>
 <p>
-Specifying 'no' for option 'Wait for clock corr.' means that BNC processes each epoch of data immediately after its arrival using satellite clock corrections available at that time. A non-zero value means that epochs of data are buffered and the processing of each epoch is postponed until satellite clock corrections older than 'Wait for clock corr.' seconds are available. Specifying a value of half the update rate of the clock corrections (e.g. 5 sec) may be appropriate. Note that this causes an additional delay of the PPP solutions in the amount of half of the update rate.
+Specifying 'no' for option 'Wait for clock corr.' means that BNC processes each epoch of data immediately after its arrival using satellite clock corrections available at that time. A non-zero value means that epochs of data are buffered and the processing of each epoch is postponed until satellite clock corrections not older than 'Wait for clock corr.' seconds are available. Specifying a value of half the update rate of the clock corrections (e.g. 5 sec) may be appropriate. Note that this causes an additional delay of the PPP solutions in the amount of half of the update rate.
 </p>
 <p>

trunk/BNC/src/pppMain.cpp

@@ -178 +178 @@
     opt->_sigmaL1      = settings.value("PPP/sigmaL1").toDouble(); if (opt->_sigmaL1 <= 0.0) opt->_sigmaL1 = 0.01;
     opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
+    if (!_realTime || opt->_corrMount.empty()) {
+      opt->_corrWaitTime = 0;
+    }
 
     if      (settings.value("PPP/lcGPS").toString() == "P1") {

trunk/BNC/src/pppRun.cpp

@@ -301 +301 @@
   }
 
+  // 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() && !waitForCorr(_epoData.front()->_time)) {
+  while (_epoData.size()) {
 
     const vector<t_satObs*>& satObs = _epoData.front()->_satObs;
 
-    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;
+    // 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) {
+
+      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 (_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.toAscii());
+        emit newNMEAstr(staID, ggaStr.toAscii());
+        if (_snxtroFile && output._epoTime.valid()) {
+          _snxtroFile->write(staID, int(output._epoTime.gpsw()), output._epoTime.gpssec(),
+                      output._trp0 + output._trp, output._trpStdev);
+        }
+      }
+      emit newMessage(QByteArray(log.str().c_str()), true);
     }
     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 (_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.toAscii());
-      emit newNMEAstr(staID, ggaStr.toAscii());
-      if (_snxtroFile && output._epoTime.valid()) {
-        _snxtroFile->write(staID, int(output._epoTime.gpsw()), output._epoTime.gpssec(),
-                    output._trp0 + output._trp, output._trpStdev);
-      }
-    }
-
-    emit newMessage(QByteArray(log.str().c_str()), true);
+      return;
+    }
   }
 }
@@ -655 +678 @@
   return '$' + nmStr + QString("*%1\n").arg(int(XOR), 0, 16).toUpper();
 }
-
-//
-////////////////////////////////////////////////////////////////////////////
-bool t_pppRun::waitForCorr(const bncTime& epoTime) const {
-
-  if (!_opt->_realTime || _opt->_corrMount.empty()) {
-    return false;
-  }
-  else if (!_lastClkCorrTime.valid()) {
-    return true;
-  }
-  else {
-    double dt = epoTime - _lastClkCorrTime;
-    if (dt > 1.0 && dt < _opt->_corrWaitTime) {
-      return true;
-    }
-    else {
-      return false;
-    }
-  }
-  return false;
-}

trunk/BNC/src/pppRun.h

@@ -62 +62 @@
   };
 
-  bool waitForCorr(const bncTime& epoTime) const;
-
   QMutex                 _mutex;
   const t_pppOptions*    _opt;

trunk/BNC/src/pppWidgets.cpp

@@ -163 +163 @@
   _corrWaitTime->setSingleStep(1);
   _corrWaitTime->setSuffix(" sec");
-  _corrWaitTime->setSpecialValueText("no");
 
   _staTable->setColumnCount(10);
@@ -194 +193 @@
   // WhatsThis, PPP (3)
   // ------------------
-  _corrWaitTime->setWhatsThis(tr("<p>Specifying 'no' means that BNC will not wait for satellite clock corrections and process each epoch of observation data immediately after arrival using the corrections available at that time.</p><p>Specifying a non-zero value (i.e. 5 sec) means that data will be buffered and their processing is postponed till corrections not older than '5 sec' (example) become available.</p>"));
+  _corrWaitTime->setWhatsThis(tr("<p>Zero value means that BNC processes each epoch of data immediately after its arrival using satellite clock corrections available at that time.</p><p> Specifying a non-zero value (i.e. 5 sec) means that the epochs of data are buffered and the processing of each epoch is postponed till the satellite clock corrections not older than '5 sec' (example) become available.<p>"));
   _seedingTime->setWhatsThis(tr("<p>Enter the length of a startup period in seconds for which you want to fix the PPP solutions to known a priori coordinates as introduced through option 'Coordinates file'. Adjust 'Sigma N/E/H' in the PPP Stations table according to the coordinate's precision. Fixing a priori coordinates is done in BNC through setting 'Sigma N/E/H' temporarily to zero.</p><p>This option allows the PPP solution to rapidly converge. It requires that the antenna remains unmoved on the a priori known position throughout the startup period.</p><p>A value of 60 is likely to be an appropriate choice.</p><p>Default is an empty option field, meaning that you don't want BNC to fix PPP solutions during startup to an a priori coordinate.</p>"));