Changeset 11087 in ntrip for trunk/BNC/src/combination/bncbiassnx.cpp

Timestamp:

Oct 16, 2026, 10:19:51 AM (29 hours ago)

Author:

stuerze

Message:

Bias check is added in ssrQC, refernce signals are configurable now

File:

• trunk/BNC/src/combination/bncbiassnx.cpp (modified) (3 diffs)

trunk/BNC/src/combination/bncbiassnx.cpp

@@ -40 +40 @@
 
 #include "bncbiassnx.h"
+#include "ssrRefSig.h"
 
 #include <iostream>
@@ -48 +49 @@
 // Constructor
 ////////////////////////////////////////////////////////////////////////////
-bncBiasSnx::bncBiasSnx() {
+bncBiasSnx::bncBiasSnx(bool allGnss) {
   bncSettings settings;
+  if (allGnss) {
+    for (char sys : QString("GRECJSI").toLatin1()) {
+      _useGnss[sys] = true;
+    }
+    return;
+  }
   _useGnss['G'] = (Qt::CheckState(settings.value("cmbGps").toInt()) == Qt::Checked) ? true : false;
   _useGnss['R'] = (Qt::CheckState(settings.value("cmbGlo").toInt()) == Qt::Checked) ? true : false;
@@ -307 +314 @@
 }
 
-//
+// DCBs (OSB_A - OSB_B) used to derive the code biases from the default
+// reference signals (see t_ssrRefSig), in the order of evaluation
+////////////////////////////////////////////////////////////////////////////
+static const QVector<QPair<QString, QString> >& defaultDcbs(char sys) {
+  static QMap<char, QVector<QPair<QString, QString> > > dcbs;
+  static QMutex mutex;
+  QMutexLocker locker(&mutex);
+  if (dcbs.isEmpty()) {
+    auto add = [](QVector<QPair<QString, QString> >& vv, const char* list) {
+      for (const QString& dcb : QString(list).split(' ', Qt::SkipEmptyParts)) {
+        vv.append(qMakePair(dcb.section('-', 0, 0), dcb.section('-', 1, 1)));
+      }
+    };
+    add(dcbs['G'], "1C-1W 1C-5Q 1C-5X 1C-2W 2C-2W 2W-2S 2W-2L 2W-2X");
+    add(dcbs['R'], "1C-1P 2C-2P");
+    add(dcbs['E'], "1C-6C 1C-7Q 1C-8Q");
+    add(dcbs['C'], "2I-7I 1X-6I 1X-7Z 1X-8X 1X-5X 1P-6I 1P-5P 1D-6I 1D-5D");
+    add(dcbs['J'], "1C-1X 1X-2X 1C-5X");
+  }
+  static const QVector<QPair<QString, QString> > none;
+  return dcbs.contains(sys) ? dcbs[sys] : none;
+}
+
+// Satellite code biases (OSB) from the DCBs, referred to the reference
+// signals of the clocks: aIF1 * OSB_ref1 + aIF2 * OSB_ref2 = 0
 ////////////////////////////////////////////////////////////////////////////
 void bncBiasSnx::determineSsrSatCodeBiases(QString prn, double aIF1, double aIF2,
                                            t_satCodeBias& satCodeBias) {
-  std::string rnxType2ch;
-  double value;
-  if (_snxSatCodeBiasMap.contains(prn)) {
-    char sys = prn.at(0).toLatin1();
-    QMap<QString, double>& cb = _snxSatCodeBiasMap[prn]->_biasMap;
-    switch (sys) {
-      case 'G':
-        if  (cb.contains("C1W-C2W")) {
-          double dcb_1W2W = cb.find("C1W-C2W").value();
-          rnxType2ch = "1W";
-          value =  aIF2 * dcb_1W2W;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          rnxType2ch = "2W";
-          value = -aIF1 * dcb_1W2W;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          if   (cb.contains("C1C-C1W")) {
-            double dcb_1C1W = cb.find("C1C-C1W").value();
-            rnxType2ch = "1C";
-            value = dcb_1C1W + aIF2 * dcb_1W2W;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            if  (cb.contains("C1C-C5Q")) {
-              double dcb_1C5Q = cb.find("C1C-C5Q").value();
-              rnxType2ch = "5Q";
-              value = dcb_1C1W + aIF2 * dcb_1W2W - dcb_1C5Q;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-            if  (cb.contains("C1C-C5X")) {
-              double dcb_1C5X = cb.find("C1C-C5X").value();
-              rnxType2ch = "5X";
-              value = dcb_1C1W + aIF2 * dcb_1W2W - dcb_1C5X;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-          }
-          else {
-            if  (cb.contains("C1C-C2W")) {
-              double dcb_1C2W = cb.find("C1C-C2W").value();
-              rnxType2ch = "1C";
-              value = dcb_1C2W - aIF1 * dcb_1W2W;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-          }
-          if   (cb.contains("C2C-C2W")) {
-            double dcb_2C2W = cb.find("C2C-C2W").value();
-            rnxType2ch = "2C";
-            value = dcb_2C2W - aIF1 * dcb_1W2W;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if   (cb.contains("C2W-C2S")) {
-            double dcb_2W2S = cb.find("C2W-C2S").value();
-            rnxType2ch = "2S";
-            value = -aIF1 * dcb_1W2W - dcb_2W2S;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if   (cb.contains("C2W-C2L")) {
-            double dcb_2W2L = cb.find("C2W-C2L").value();
-            rnxType2ch = "2L";
-            value = -aIF1 * dcb_1W2W - dcb_2W2L;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if   (cb.contains("C2W-C2X")) {
-            double dcb_2W2X = cb.find("C2W-C2X").value();
-            rnxType2ch = "2X";
-            value = -aIF1 * dcb_1W2W - dcb_2W2X;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-        }
-        break;
-      case 'R':
-        if  (cb.contains("C1P-C2P")) {
-          double dcb_1P2P = cb.find("C1P-C2P").value();
-          rnxType2ch = "1P";
-          value =  aIF2 * dcb_1P2P;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          rnxType2ch = "2P";
-          value = -aIF1 * dcb_1P2P;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          if   (cb.contains("C1C-C1P")) {
-            double dcb_1C1P = cb.find("C1C-C1P").value();
-            rnxType2ch = "1C";
-            value = dcb_1C1P + aIF2 * dcb_1P2P;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if   (cb.contains("C2C-C2P")) {
-            double dcb_2C2P = cb.find("C2C-C2P").value();
-            rnxType2ch = "2C";
-            value = dcb_2C2P - aIF1 * dcb_1P2P;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-        }
-        break;
-      case 'E':
-        if  (cb.contains("C1C-C5Q")) {
-          double dcb_1C5Q = cb.find("C1C-C5Q").value();
-          rnxType2ch = "5Q";
-          value = -aIF1 * dcb_1C5Q;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          rnxType2ch = "1C";
-          value =  aIF2 * dcb_1C5Q;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          if  (cb.contains("C1C-C6C")) {
-            double dcb_1C6C = cb.find("C1C-C6C").value();
-            rnxType2ch = "6C";
-            value = aIF2 * dcb_1C5Q - dcb_1C6C;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if  (cb.contains("C1C-C7Q")) {
-            double dcb_1C7Q = cb.find("C1C-C7Q").value();
-            rnxType2ch = "7Q";
-            value = aIF2 * dcb_1C5Q - dcb_1C7Q;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if  (cb.contains("C1C-C8Q")) {
-            double dcb_1C8Q = cb.find("C1C-C8Q").value();
-            rnxType2ch = "8Q";
-            value = aIF2 * dcb_1C5Q - dcb_1C8Q;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-        }
-        break;
-      case 'C':
-        if  (cb.contains("C2I-C6I")) {
-          double dcb_2I6I = cb.find("C2I-C6I").value();
-          rnxType2ch = "2I";
-          value =  aIF2 * dcb_2I6I;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          rnxType2ch = "6I";
-          value = -aIF1 * dcb_2I6I;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          if (cb.contains("C2I-C7I")) {
-            double dcb_2I7I = cb.find("C2I-C7I").value();
-            rnxType2ch = "7I";
-            value = aIF2 * dcb_2I6I - dcb_2I7I;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-          if (cb.contains("C1X-C6I")) {
-            double dcb_1X6I = cb.find("C1X-C6I").value();
-            rnxType2ch = "1X";
-            value = dcb_1X6I - aIF1 * dcb_2I6I;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            if (cb.contains("C1X-C7Z")) {
-              double dcb_1X7Z = cb.find("C1X-C7Z").value();
-              rnxType2ch = "7Z";
-              value = dcb_1X6I - aIF1 * dcb_2I6I - dcb_1X7Z;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-            if (cb.contains("C1X-C8X")) {
-              double dcb_1X8X = cb.find("C1X-C8X").value();
-              rnxType2ch = "8X";
-              value = dcb_1X6I - aIF1 * dcb_2I6I - dcb_1X8X;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-            if (cb.contains("C1X-C5X")) {
-              double dcb_1X5X = cb.find("C1X-C5X").value();
-              rnxType2ch = "5X";
-              value = dcb_1X6I - aIF1 * dcb_2I6I - dcb_1X5X;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-          }
-          if (cb.contains("C1P-C6I")) {
-            double dcb_1P6I = cb.find("C1P-C6I").value();
-            rnxType2ch = "1P";
-            value = dcb_1P6I - aIF1 * dcb_2I6I;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            if (cb.contains("C1P-C5P")) {
-              double dcb_1P5P = cb.find("C1P-C5P").value();
-              rnxType2ch = "5P";
-              value = dcb_1P6I - aIF1 * dcb_2I6I - dcb_1P5P;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-          }
-          if (cb.contains("C1D-C6I")) {
-            double dcb_1D6I = cb.find("C1D-C6I").value();
-            rnxType2ch = "1D";
-            value = dcb_1D6I - aIF1 * dcb_2I6I;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            if (cb.contains("C1D-C5D")) {
-              double dcb_1D5D = cb.find("C1D-C5D").value();
-              rnxType2ch = "5D";
-              value = dcb_1D6I - aIF1 * dcb_2I6I - dcb_1D5D;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-          }
-        }
-        break;
-      case 'J':
-        if (cb.contains("C1C-C2L")) {
-          double dcb_1C2L = cb.find("C1C-C2L").value();
-          rnxType2ch = "1C";
-          value =  aIF2 * dcb_1C2L;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          rnxType2ch = "2L";
-          value = -aIF1 * dcb_1C2L;
-          setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          if (cb.contains("C1C-C1X")) {
-            double dcb_1C1X = cb.find("C1C-C1X").value();
-            rnxType2ch = "1X";
-            value =  aIF2 * dcb_1C2L - dcb_1C1X;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            if (cb.contains("C1X-C2X")) {
-              double dcb_1X2X = cb.find("C1X-C2X").value();
-              rnxType2ch = "2X";
-              value =  aIF2 * dcb_1C2L - dcb_1C1X - dcb_1X2X;
-              setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-            }
-          }
-          if (cb.contains("C1C-C5X")) {
-            double dcb_1C5X = cb.find("C1C-C5X").value();
-            rnxType2ch = "5X";
-            value =  aIF2 * dcb_1C2L - dcb_1C5X;
-            setSsrSatCodeBias(rnxType2ch, value, satCodeBias);
-          }
-        }
-        break;
-      default:
-        break;
-    }
-  }
-  return;
+  if (!_snxSatCodeBiasMap.contains(prn)) {
+    return;
+  }
+  char sys = prn.at(0).toLatin1();
+  t_ssrRefSig::t_pair pair;
+  if (!t_ssrRefSig::get(sys, pair)) {
+    return;
+  }
+  const QString ref1 = QString::fromStdString(pair.sig1());
+  const QString ref2 = QString::fromStdString(pair.sig2());
+  const QMap<QString, double>& cb = _snxSatCodeBiasMap[prn]->_biasMap;
+  auto dcb = [&cb](const QString& sig1, const QString& sig2, double& value) {
+    QString key = "C" + sig1 + "-C" + sig2;
+    if (cb.contains(key)) {
+      value = cb[key];
+      return true;
+    }
+    return false;
+  };
+
+  QMap<QString, double> osb;
+  QStringList           order;
+
+  // Default reference signals: the DCB of the reference signals and the
+  // DCBs of the IGS conventions, evaluated in one pass as before
+  if (t_ssrRefSig::isDefault(sys)) {
+    const QVector<QPair<QString, QString> >& dcbs = defaultDcbs(sys);
+    double dcbRef;
+    if (dcbs.isEmpty() || !dcb(ref1, ref2, dcbRef)) {
+      return;
+    }
+    osb[ref1] =  aIF2 * dcbRef;
+    osb[ref2] = -aIF1 * dcbRef;
+    if (sys == 'E') {   // order of the output as before
+      order << ref2 << ref1;
+    }
+    else {
+      order << ref1 << ref2;
+    }
+    for (const QPair<QString, QString>& dd : dcbs) {
+      double value;
+      if (!dcb(dd.first, dd.second, value)) {
+        continue;
+      }
+      if (osb.contains(dd.first) && !osb.contains(dd.second)) {
+        osb[dd.second] = osb[dd.first] - value;
+        order << dd.second;
+      }
+      else if (osb.contains(dd.second) && !osb.contains(dd.first)) {
+        osb[dd.first] = value + osb[dd.second];
+        order << dd.first;
+      }
+    }
+  }
+
+  // Other reference signals: all DCBs of the satellite, the DCB of the
+  // reference signals may also be derived via other signals
+  else {
+    QMap<QString, double> rel;   // OSB - OSB_ref1
+    QStringList           relOrder;
+    rel[ref1] = 0.0;
+    relOrder << ref1;
+    bool changed = true;
+    while (changed) {
+      changed = false;
+      QMapIterator<QString, double> it(cb);
+      while (it.hasNext()) {
+        it.next();
+        QStringList hlp = it.key().split('-');
+        if (hlp.size() != 2 || hlp[0].size() != 3 || hlp[1].size() != 3) {
+          continue;
+        }
+        QString sigA = hlp[0].mid(1);
+        QString sigB = hlp[1].mid(1);
+        if (rel.contains(sigA) && !rel.contains(sigB)) {
+          rel[sigB] = rel[sigA] - it.value();
+          relOrder << sigB;
+          changed = true;
+        }
+        else if (rel.contains(sigB) && !rel.contains(sigA)) {
+          rel[sigA] = rel[sigB] + it.value();
+          relOrder << sigA;
+          changed = true;
+        }
+      }
+    }
+    if (!rel.contains(ref2)) {
+      return;
+    }
+    double osbRef1 = -aIF2 * rel[ref2];
+    order << ref1 << ref2;
+    for (const QString& sig : relOrder) {
+      osb[sig] = osbRef1 + rel[sig];
+      if (sig != ref1 && sig != ref2) {
+        order << sig;
+      }
+    }
+  }
+
+  for (const QString& sig : order) {
+    setSsrSatCodeBias(sig.toStdString(), osb[sig], satCodeBias);
+  }
 }