Changeset 11085 in ntrip for trunk/BNC/src

Timestamp:

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

Author:

stuerze

Message:

yaw angle information added to ssrQC

Location:

trunk/BNC/src

Files:

• bnchelp.html (modified) (1 diff)
• ssrQc.cpp (modified) (6 diffs)
• ssrQc.h (modified) (1 diff)

trunk/BNC/src/bnchelp.html

@@ -5132 +5132 @@
       (see section 'Combination'),</li>
     <li>the phase bias signals provided,</li>
+    <li>the availability of the satellite yaw angle with the phase biases: 'not provided', or the number of
+      satellites and the percentage of phase bias records with yaw angle, listing satellites without or only partly
+      with yaw angle. New RTCM-SSR signals the yaw angle with the Yaw Information Indicator; old RTCM-SSR and
+      IGS-SSR always carry the field, a provider without attitude information sends zeros, which is counted as not
+      provided. Without a yaw angle from the stream, a PPP or combination user depends on its own attitude model
+      (see section 'Combination'), which may differ from the provider's attitude during eclipse seasons and noon/midnight
+      turns,</li>
     <li>clock jumps, i.e. clock changes between consecutive epochs (same IOD) beyond the clock rate of more than
       0.5 m,</li>

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()) {

trunk/BNC/src/ssrQc.h

@@ -66 +66 @@
     bool          refBiasOK = false;
     double        biasIF = 0.0;  // ionosphere-free combination of the reference biases [m]
+    unsigned      numPB  = 0;    // phase bias records
+    unsigned      numYaw = 0;    // phase bias records with yaw angle
   };