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

Timestamp:

Oct 14, 2026, 11:47:57 AM (3 days ago)

Author:

stuerze

Message:

yaw angle information added to ssrQC

File:

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

trunk/BNC/src/ssrQc.cpp

@@ -182 +182 @@
            << "Checks: epochs and update intervals, gaps (interval > 2 x max(declared, median)), latency,"
            << " completeness w.r.t. healthy broadcast ephemerides," << endl
-           << "IOD matching, code and phase bias signals, code biases of the reference signals, clock jumps > "
+           << "IOD matching, code and phase bias signals, yaw angle, code biases of the reference signals, clock jumps > "
            << MAX_CLK_JUMP << " m, largest orbit corrections" << endl;
     }
@@ -416 +416 @@
       _phaseSignals[staSys].insert(QString::fromStdString(frqBias._rnxType2ch));
     }
+    // Yaw angle: new RTCM-SSR has an indicator, old RTCM-SSR and IGS-SSR always
+    // carry the field, providers without attitude information send zeros
+    t_satStat& stat = _satStat[QString::fromStdString(bias._staID) + ' ' +
+                               QString::fromStdString(bias._prn.toString())];
+    stat.numPB++;
+    bool yawProvided = (bias._ssrFormat == ssrRtcmNew) ? (bias._satYawInfoInd != 0)
+                                                       : (bias._yaw != 0.0 || bias._yawRate != 0.0);
+    if (yawProvided) {
+      stat.numYaw++;
+    }
   }
 }
@@ -707 +717 @@
     if (final) {
       endTime = _staLast.value(staID, _lastTime);
+      if (endTime < _intervalStart) { // stream without data in the last interval
+        endTime = _lastTime;
+      }
     }
     out << endl << "=== SSR QC " << staID.toStdString() << ' '
@@ -746 +759 @@
     for (char sys : systems) {
       QMap<QString, unsigned> missOrb, missClk, missCB, missPB, iodUnknown, unhealthy, noEph,
-                              jumps, missRef;
+                              jumps, missRef, noYaw, yawPartly;
+      unsigned numPB = 0, numYaw = 0, yawSats = 0;
       double maxRadial = 0.0, maxAlongCross = 0.0;
       QString maxRadialPrn, maxAlongCrossPrn;
@@ -792 +806 @@
         if (stat.refBiasOK) {
           biasIF.push_back(stat.biasIF);
+        }
+        if (stat.numPB) {
+          numPB  += stat.numPB;
+          numYaw += stat.numYaw;
+          if (stat.numYaw) {
+            yawSats++;
+            if (stat.numYaw < stat.numPB) {
+              yawPartly[prn] = stat.numPB - stat.numYaw;
+            }
+          }
+          else {
+            noYaw[prn] = 0;
+          }
         }
       }
@@ -843 +870 @@
         lines << "phase bias signals: " + sigs.join(" ");
       }
+      if (numPB) {
+        if (numYaw == 0) {
+          lines << "yaw angle: not provided";
+        }
+        else {
+          lines << QString("yaw angle: provided for %1 of %2 satellites (%3 % of the phase bias records)")
+                   .arg(yawSats).arg(yawSats + noYaw.size()).arg(100.0 * numYaw / numPB, 0, 'f', 1);
+          addSats("phase biases without yaw angle", noYaw);
+          addSats("phase biases partly without yaw angle (records)", yawPartly);
+        }
+      }
       addSats(QString("clock jumps > %1 m").arg(MAX_CLK_JUMP), jumps);
       if (!maxRadialPrn.isEmpty()) {