Changeset 11194 in ntrip for trunk/BNC/src/RTCM3/RTCM3coDecoder.cpp

Timestamp:

Nov 10, 2026, 2:21:34 PM (83 minutes ago)

Author:

stuerze

Message:

IGS-SSR format checked and fixed

File:

• trunk/BNC/src/RTCM3/RTCM3coDecoder.cpp (modified) (7 diffs)

trunk/BNC/src/RTCM3/RTCM3coDecoder.cpp

@@ -375 +375 @@
       continue;
     }
+    num = satNumber(sys, num);
 
     // Orbit correction
@@ -429 +430 @@
 
       t_clkCorr clkCorr;
-      clkCorr._prn.set(sys, _clkOrb.Sat[ii].ID, flag);
+      clkCorr._prn.set(sys, num, flag);
       checkPrnRange(clkCorr._prn);
       clkCorr._staID      = _staID.toStdString();
@@ -454 +455 @@
          _clkOrb.messageType == _ssrCorr->COTYPE_SBASHR ||
          _clkOrb.messageType == _ssrCorr->COTYPE_BDSHR) {
-      t_prn prn(sys, _clkOrb.Sat[ii].ID, flag);
+      t_prn prn(sys, num, flag);
       if (_lastClkCorrections.contains(prn)) {
         t_clkCorr clkCorr;
         clkCorr            = _lastClkCorrections[prn];
+        // The polynomial keeps the reference time t0 of its clock message
+        // (IGS SSR v1.00 4.2, RTCM 10403.3): re-reference it to the t0 of
+        // the High Rate Clock message, then add the High Rate correction.
+        double t0Clk = (clkCorr._updateInt ? 0.5 * ssrUpdateInt[clkCorr._updateInt] : 0.0);
+        double t0HR  = (_clkOrb.UpdateInterval ? 0.5 * ssrUpdateInt[_clkOrb.UpdateInterval] : 0.0);
+        double dt    = (_lastTime + t0HR) - (clkCorr._time + t0Clk);
+        clkCorr._dClk      += clkCorr._dotDClk * dt + clkCorr._dotDotDClk * dt * dt;
+        clkCorr._dotDClk   += 2.0 * clkCorr._dotDotDClk * dt;
         clkCorr._time      = _lastTime;
         clkCorr._updateInt = _clkOrb.UpdateInterval;
@@ -521 +530 @@
       continue;
     }
+    num = satNumber(sys, num);
     t_satCodeBias satCodeBias;
     satCodeBias._prn.set(sys, num, flag);
@@ -614 +624 @@
       continue;
     }
+    num = satNumber(sys, num);
     t_satPhaseBias satPhaseBias;
     satPhaseBias._prn.set(sys, num, flag);
@@ -992 +1003 @@
         corrIsOutOfRange = true;
       }
-      if (coSat.Clock.DeltaA2 < -1.3421772 ||
-          coSat.Clock.DeltaA2 > +1.3421772) {
+      if (coSat.Clock.DeltaA2 < -1.34217726 ||
+          coSat.Clock.DeltaA2 > +1.34217726) {
         ssrParStr = "Clock::DeltaA2";
-        ssrParValue = QString::number(coSat.Clock.DeltaA2, 'f', 7);
+        ssrParValue = QString::number(coSat.Clock.DeltaA2, 'f', 8);
         corrIsOutOfRange = true;
       }
@@ -1368 +1379 @@
 }
 
+// Satellite number from the GNSS Satellite ID on the wire. IGS-SSR (IGS SSR
+// v1.00 IDF011): ID 0 means 64 for GPS, GLONASS and BDS; SBAS ID 1 is PRN 120
+// (RINEX S20). RTCM-SSR IDs are used as they are.
+////////////////////////////////////////////////////////////////////////////
+int RTCM3coDecoder::satNumber(char sys, int id) const {
+  if (_type == IGSssr) {
+    if (id == 0 && (sys == 'G' || sys == 'R' || sys == 'C')) {
+      return 64;
+    }
+    if (sys == 'S') {
+      return id + 19;
+    }
+  }
+  return id;
+}
+
 // Satellite numbers on the wire can exceed BNC's supported range (e.g. 6-bit
 // IDs up to 63, while t_prn supports G01-G32). They are passed on unchanged,