Changeset 11061 in ntrip for trunk/BNC/src/combination

Timestamp:

Oct 1, 2026, 6:34:50 PM (4 weeks ago)

Author:

stuerze

Message:

some more fixes

Location:

trunk/BNC/src/combination

Files:

• bnccomb.cpp (modified) (4 diffs)
• bnccomb.h (modified) (1 diff)

trunk/BNC/src/combination/bnccomb.cpp

@@ -899 +899 @@
           _newCorr->_satCodeBiasIF += codeCoeff * codeBiasesRefSig[frqType];
         }
+        if (codeBiasesRefSig.size() == 2) {
+          _acBiasIF[acName][prnStr.mid(0,3)] = _newCorr->_satCodeBiasIF;
+        }
+        else {
+          _acBiasIF[acName].remove(prnStr.mid(0,3));
+        }
         _newCorr->_satCodeBias._bias.clear();
       }
@@ -1079 +1085 @@
       out.flush();
     }
+    if (_method == filter) {
+      checkBiasConsistency(epoTime, sys, out);
+    }
     printResults(epoTime, out, masterCorr);
     dumpResults(epoTime, masterCorr);
@@ -1144 +1153 @@
         if (pp->type == cmbParam::offACSat            &&
             pp->AC   == corrs(sys)[maxResIndex-1]->_acName &&
-            pp->prn  == corrs(sys)[maxResIndex-1]->_prn.mid(0,3)) {
+            pp->prn  == corrs(sys)[maxResIndex-1]->_prn) {
           QQ_sav.Row(iPar)    = 0.0;
           QQ_sav.Column(iPar) = 0.0;
@@ -1174 +1183 @@
 
   return success;
+}
+
+// Check the consistency of the AC code biases with the AC clocks
+////////////////////////////////////////////////////////////////////////////
+// An AC's clocks refer to its own clock signals; the code biases of the
+// reference signals refer them to the reference signals (B_IF, ionosphere-
+// free combination of these biases, added to the AC clocks). If the code
+// biases do not belong to the clocks (e.g. taken from another product), the
+// error B_IF is absorbed by the AC's satellite offsets. The satellite offsets
+// of each AC are therefore regressed against its B_IF every 5 minutes (from
+// 5 minutes after the start); a warning is issued (at most once per hour per
+// AC and system) if B_IF varies noticeably between the satellites and the
+// offsets clearly follow it.
+////////////////////////////////////////////////////////////////////////////
+void bncComb::checkBiasConsistency(bncTime epoTime, char sys, QTextStream& out) {
+
+  const double CHECK_INT   = 300.0;   // [s] check interval, first check after start
+  const unsigned MIN_SAT   = 5;
+  const double MIN_SPREAD  = 0.03;    // [m] minimum std. deviation of B_IF
+  const double MIN_SLOPE   = 0.3;     // 0.5 expected for 2 ACs, -> 1 for more ACs
+  const double MIN_CORR    = 0.7;
+
+  if (!_biasCheckStart[sys].valid()) {
+    _biasCheckStart[sys] = epoTime;
+    _biasCheckLast[sys]  = epoTime;
+    return;
+  }
+  if (epoTime - _biasCheckLast[sys] < CHECK_INT) {
+    return;
+  }
+  _biasCheckLast[sys] = epoTime;
+
+  QListIterator<cmbAC*> itAc(_ACs);
+  while (itAc.hasNext()) {
+    const QString& acName = itAc.next()->name;
+    if (!_acBiasIF.contains(acName)) {
+      continue;
+    }
+    const QMap<QString, double>& biasIF = _acBiasIF[acName];
+
+    // satellites of this AC in the current epoch with B_IF available
+    QVector<double> xx, yy;
+    for (int ii = 0; ii < corrs(sys).size(); ii++) {
+      const cmbCorr* corr = corrs(sys)[ii];
+      if (corr->_acName != acName || !biasIF.contains(corr->_prn.mid(0,3))) {
+        continue;
+      }
+      for (int iPar = 1; iPar <= _params[sys].size(); iPar++) {
+        const cmbParam* pp = _params[sys][iPar-1];
+        if (pp->type == cmbParam::offACSat && pp->AC == acName && pp->prn == corr->_prn) {
+          xx.push_back(biasIF[corr->_prn.mid(0,3)]);
+          yy.push_back(pp->xx);
+          break;
+        }
+      }
+    }
+    unsigned nn = xx.size();
+    if (nn < MIN_SAT) {
+      continue;
+    }
+    double mx = 0.0, my = 0.0;
+    for (unsigned ii = 0; ii < nn; ii++) {
+      mx += xx[ii] / nn;
+      my += yy[ii] / nn;
+    }
+    double sxx = 0.0, syy = 0.0, sxy = 0.0;
+    for (unsigned ii = 0; ii < nn; ii++) {
+      sxx += (xx[ii] - mx) * (xx[ii] - mx);
+      syy += (yy[ii] - my) * (yy[ii] - my);
+      sxy += (xx[ii] - mx) * (yy[ii] - my);
+    }
+    double spread = sqrt(sxx / nn);
+    if (spread < MIN_SPREAD || syy <= 0.0) {
+      continue;
+    }
+    double slope = sxy / sxx;
+    double corrCoeff = sxy / sqrt(sxx * syy);
+    if (slope < MIN_SLOPE || corrCoeff < MIN_CORR) {
+      continue;
+    }
+
+    QString refSig;
+    for (unsigned ii = 1; ii < cmbRefSig::cIF; ii++) {
+      t_frequency::type frqType = cmbRefSig::toFreq(sys, static_cast<cmbRefSig::type>(ii));
+      refSig += QString("%1C%2%3").arg(ii > 1 ? "/" : "").arg(t_frequency::toString(frqType)[1])
+                                  .arg(cmbRefSig::toAttrib(sys, static_cast<cmbRefSig::type>(ii)));
+    }
+    QString msg = QString("%1 %2 code biases (%3) seem inconsistent with its clocks: "
+                          "satellite offsets follow the bias combination B_IF "
+                          "(%4 satellites, B_IF spread %5 m, slope %6, correlation %7)")
+                  .arg(acName).arg(sys).arg(refSig).arg(nn)
+                  .arg(spread, 0, 'f', 3).arg(slope, 0, 'f', 2).arg(corrCoeff, 0, 'f', 2);
+    out << epoTime.datestr().c_str() << " " << epoTime.timestr().c_str()
+        << " " << msg << "\n";
+    QString key = acName + sys;
+    bncTime& warned = _biasCheckWarned[key];
+    if (!warned.valid() || epoTime - warned >= 3600.0) {
+      warned = epoTime;
+      emit newMessage("bncComb: " + msg.toLatin1(), true);
+    }
+  }
 }
 

trunk/BNC/src/combination/bnccomb.h

@@ -270 +270 @@
   QMap<QString, bncTime>                     _refBiasWarned; // last warning per AC and system
   QMap<char, bncTime>                        _singleACWarned; // last 'fewer than 2 ACs' message per system
+  QMap<QString, QMap<QString, double> >      _acBiasIF;       // latest B_IF of the reference signals per AC and satellite
+  QMap<char, bncTime>                        _biasCheckStart; // first combined epoch per system
+  QMap<char, bncTime>                        _biasCheckLast;  // last bias/clock consistency check per system
+  QMap<QString, bncTime>                     _biasCheckWarned; // last inconsistency warning per AC and system
+  void checkBiasConsistency(bncTime epoTime, char sys, QTextStream& out);
   bool allACsBeyond(const bncTime& epoTime) const;
   bncTime                                    _resTime;