Changeset 11137 in ntrip for trunk/BNC/src/ssrQc.cpp

Timestamp:

Oct 28, 2026, 9:48:01 PM (2 days ago)

Author:

stuerze

Message:

eph handling improved

File:

• trunk/BNC/src/ssrQc.cpp (modified) (10 diffs)

trunk/BNC/src/ssrQc.cpp

@@ -56 +56 @@
   const double MAX_CLK_JUMP   = 0.5;  // [m] clock change beyond its rate between consecutive epochs
   const double EPH_WARMUP     = 120.0; // [s] no ephemeris checks at the start
+  const double EPH_EARLY_MAX  = 300.0; // [s] IOD received this much later still counted as early
+  const double EPH_IOD_REUSE  = 21600.0; // [s] an outdated data set this far from the epoch: IOD reused
+  const unsigned EPH_KEPT     = 30;    // data sets kept per satellite
   const double CMP_SAMPL      = 30.0;  // [s] sampling of the stream comparison
   const double CMP_DELAY      = 60.0;  // [s] epochs compared when data this much newer arrived
@@ -119 +122 @@
 ////////////////////////////////////////////////////////////////////////////
 t_ssrQc::t_ssrQc(const QString& dirName, double interval, const QStringList& streams,
-                 const QString& reference) : _ephUser(true) {
+                 const QString& reference) : _ephUser(true, EPH_KEPT) {
   _dirName   = dirName;
   _reference = reference.isEmpty() ? QString("CONSENSUS") : reference;
@@ -359 +362 @@
 // Broadcast ephemeris the correction refers to
 ////////////////////////////////////////////////////////////////////////////
-void t_ssrQc::checkEph(t_satStat& stat, const QString& prnInt, unsigned iod, const bncTime& time,
-                       const string& staID) {
+void t_ssrQc::checkEph(t_satStat& stat, const QString& satKey, const QString& prnInt,
+                       unsigned iod, const bncTime& time, const string& staID) {
   // not before all broadcast ephemerides may have been received
   if (time - _staStart[QString::fromStdString(staID)] < EPH_WARMUP) {
     return;
   }
-  t_eph* ephLast = _ephUser.ephLast(prnInt);
-  t_eph* ephPrev = _ephUser.ephPrev(prnInt);
-  t_eph* eph = 0;
-  if      (ephLast && ephLast->IOD() == iod) eph = ephLast;
-  else if (ephPrev && ephPrev->IOD() == iod) eph = ephPrev;
-  if (!ephLast) {
+  if (!_ephUser.ephLast(prnInt)) {
     stat.noEph += 1;
-  }
-  else if (!eph) {
-    stat.iodUnknown += 1;
+    return;
+  }
+  // any data set kept, not only the last two: Galileo has a new IODnav every
+  // 10 minutes, and ACs keep one for up to an hour
+  t_eph* eph = _ephUser.ephWithIOD(prnInt, iod);
+  // an IOD used again later (GLONASS tb every day, GPS IODE) is another data set
+  if (eph && outDatedBcep(eph, time) && fabs(time - eph->TOC()) > EPH_IOD_REUSE) {
+    eph = 0;
+  }
+  if (!eph) {
+    // received later from the ephemeris stream(s)? decided in resolvePendingIod
+    t_pendingIod& pending = _pendingIod[satKey][iod];
+    if (pending.num == 0) {
+      pending.prnInt = prnInt;
+      pending.first  = currentTime();
+    }
+    pending.num += 1;
   }
   else if (eph->checkState() == t_eph::unhealthy || eph->checkState() == t_eph::bad) {
     stat.ephUnhealthy += 1;
+  }
+  else if (outDatedBcep(eph, time)) {
+    stat.ephOutdated += 1;
+  }
+  if (_pendingIod.contains(satKey)) {
+    resolvePendingIod(satKey, false);
+  }
+}
+
+// Corrections referring to an IOD not received before: counted as early if
+// the ephemeris stream(s) delivered the data set later, as unknown if not
+// within EPH_EARLY_MAX (or at the end of the report interval, final)
+////////////////////////////////////////////////////////////////////////////
+void t_ssrQc::resolvePendingIod(const QString& satKey, bool final) {
+  t_satStat& stat = _satStat[satKey];
+  QMap<unsigned, t_pendingIod>& pend = _pendingIod[satKey];
+  QMutableMapIterator<unsigned, t_pendingIod> it(pend);
+  while (it.hasNext()) {
+    it.next();
+    const t_pendingIod& pending = it.value();
+    t_eph* eph = _ephUser.ephWithIOD(pending.prnInt, it.key());
+    if (eph && eph->receptDateTime().isValid()) {
+      bncTime recept(eph->receptDateTime().toString(Qt::ISODate).toStdString());
+      if (recept >= pending.first) {
+        stat.iodEarly += pending.num;
+        stat.maxEarly  = max(stat.maxEarly, recept - pending.first);
+        it.remove();
+        continue;
+      }
+    }
+    if (final || currentTime() - pending.first > EPH_EARLY_MAX) {
+      stat.iodUnknown += pending.num;
+      it.remove();
+    }
+  }
+  if (pend.isEmpty()) {
+    _pendingIod.remove(satKey);
   }
 }
@@ -395 +444 @@
     t_satStat& stat = _satStat[QString::fromStdString(corr._staID) + ' ' + prn];
     stat.numOrb += 1;
-    checkEph(stat, QString::fromStdString(corr._prn.toInternalString()), corr._iod, corr._time, corr._staID);
+    checkEph(stat, QString::fromStdString(corr._staID) + ' ' + prn,
+             QString::fromStdString(corr._prn.toInternalString()), corr._iod, corr._time, corr._staID);
     double radial = fabs(corr._xr[0]);
     double alongCross = max(fabs(corr._xr[1]), fabs(corr._xr[2]));
@@ -417 +467 @@
     t_satStat& stat = _satStat[QString::fromStdString(corr._staID) + ' ' + prn];
     stat.numClk += 1;
-    checkEph(stat, QString::fromStdString(corr._prn.toInternalString()), corr._iod, corr._time, corr._staID);
+    checkEph(stat, QString::fromStdString(corr._staID) + ' ' + prn,
+             QString::fromStdString(corr._prn.toInternalString()), corr._iod, corr._time, corr._staID);
 
     // jump: change between consecutive epochs (same IOD) beyond the clock rate
@@ -837 +888 @@
       }
       const t_orbCorr& orb = (--itO).value();
-      t_eph* eph = 0;
-      t_eph* ephLast = _ephUser.ephLast(prnInt);
-      t_eph* ephPrev = _ephUser.ephPrev(prnInt);
-      if      (ephLast && ephLast->IOD() == orb._iod) eph = ephLast;
-      else if (ephPrev && ephPrev->IOD() == orb._iod) eph = ephPrev;
+      t_eph* eph = _ephUser.ephWithIOD(prnInt, orb._iod);
+      if (eph && outDatedBcep(eph, tt) && fabs(tt - eph->TOC()) > EPH_IOD_REUSE) {
+        eph = 0;
+      }
       if (!eph || clk._iod != orb._iod) {
         continue;
@@ -1363 +1413 @@
   bncTime endTime = final ? _lastTime : _intervalStart + _interval;
 
+  // IODs not received: early or unknown (pending ones younger than
+  // EPH_EARLY_MAX are decided in the next interval)
+  for (const QString& satKey : _pendingIod.keys()) {
+    resolvePendingIod(satKey, final);
+  }
+
   // Healthy and current broadcast ephemerides per system
   QMap<char, QSet<QString> > ephSats;
@@ -1475 +1531 @@
         continue;
       }
+      QMap<QString, unsigned> iodEarly, outdated;
+      double maxEarly = 0.0;
       QMap<QString, unsigned> missOrb, missClk, missCB, missPB, iodUnknown, unhealthy, noEph,
                               jumps, missRef, noYaw, yawPartly;
@@ -1507 +1565 @@
         const t_satStat& stat = itSat.value();
         if (stat.iodUnknown)   iodUnknown[prn] = stat.iodUnknown;
+        if (stat.iodEarly) {
+          iodEarly[prn] = stat.iodEarly;
+          maxEarly = max(maxEarly, stat.maxEarly);
+        }
+        if (stat.ephOutdated)  outdated[prn]   = stat.ephOutdated;
         if (stat.ephUnhealthy) unhealthy[prn]  = stat.ephUnhealthy;
         if (stat.noEph)        noEph[prn]      = stat.noEph;
@@ -1568 +1631 @@
       addSats("phase biases missing (healthy broadcast ephemeris available)", missPB);
       addSats("corrections referring to an IOD of no received broadcast ephemeris", iodUnknown);
+      if (!iodEarly.isEmpty()) {
+        addSats(QString("corrections referring to an IOD received from the ephemeris stream(s) "
+                        "only later (up to %1 s)").arg(maxEarly, 0, 'f', 0), iodEarly);
+      }
+      addSats("corrections referring to a broadcast ephemeris outside its validity", outdated);
       addSats("corrections for unhealthy broadcast ephemeris", unhealthy);
       addSats("corrections for satellites without broadcast ephemeris", noEph);