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

Timestamp:

Oct 22, 2026, 6:20:09 PM (2 weeks ago)

Author:

stuerze

Message:

SSR QC improved

File:

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

trunk/BNC/src/ssrQc.cpp

@@ -1279 +1279 @@
       QString maxRadialPrn, maxAlongCrossPrn;
       QVector<double> biasIF;
+      double  maxBiasIF = 0.0;
+      QString biasIFPrn;
       const QSet<QString>& ephs = ephSats[sys];
       for (int type = ORB; type <= PBIAS; type++) {
@@ -1321 +1323 @@
         if (stat.refBiasOK) {
           biasIF.push_back(stat.biasIF);
+          if (biasIFPrn.isEmpty() || fabs(stat.biasIF) > maxBiasIF) {
+            maxBiasIF  = fabs(stat.biasIF);
+            biasIFPrn  = prn;
+          }
         }
         if (stat.numPB) {
@@ -1373 +1379 @@
           double var = 0.0;
           for (double bb : biasIF) var += (bb - mean) * (bb - mean) / biasIF.size();
-          double maxAbs = max(fabs(*std::min_element(biasIF.begin(), biasIF.end())),
-                              fabs(*std::max_element(biasIF.begin(), biasIF.end())));
+          // the median shows whether the clocks refer to other signals (all
+          // satellites) or single satellites have inconsistent biases
+          QVector<double> absIF;
+          for (double bb : biasIF) absIF << fabs(bb);
           lines << QString("ionosphere-free combination B_IF of the reference signal biases: "
-                           "std. dev. %1 m, max |B_IF| %2 m (%3 satellites)")
-                   .arg(sqrt(var), 0, 'f', 3).arg(maxAbs, 0, 'f', 3).arg(biasIF.size());
+                           "std. dev. %1 m, median |B_IF| %2 m, max |B_IF| %3 m (%4) (%5 satellites)")
+                   .arg(sqrt(var), 0, 'f', 3).arg(median(absIF), 0, 'f', 3).arg(maxBiasIF, 0, 'f', 3)
+                   .arg(biasIFPrn).arg(biasIF.size());
         }
       }