gnss-js 1.63.0 → 1.65.0

package/dist/index.cjs

@@ -144,9 +144,12 @@ __export(src_exports, {
   dateToEpoch: () => dateToEpoch,
   decodeBinexEph: () => decodeBinexEph,
   decodeEphemeris: () => decodeEphemeris,
+  decodeIgsSsr: () => decodeIgsSsr,
   decodeLegacyObs: () => decodeLegacyObs,
   decodeMsmFull: () => decodeMsmFull,
   decodeSbfNavigation: () => decodeSbfNavigation,
+  decodeSsr: () => decodeSsr,
+  decodeSystemParams: () => decodeSystemParams,
   deg2dms: () => deg2dms,
   deg2rad: () => deg2rad,
   detrendIonoArcs: () => detrendIonoArcs,
@@ -239,6 +242,7 @@ __export(src_exports, {
   ionoFree: () => ionoFree,
   isSbasL1Preamble: () => isSbasL1Preamble,
   isSbasL5Correction: () => isSbasL5Correction,
+  isSsrMessage: () => isSsrMessage,
   keplerPosition: () => keplerPosition,
   klobucharDelay: () => klobucharDelay,
   lambdaReduction: () => lambdaReduction,
@@ -5020,6 +5024,192 @@ function decodeEphemeris(frame) {
   }
 }
 
+// src/rtcm3/igs-ssr.ts
+var UPDATE_INTERVAL_S = [
+  1,
+  2,
+  5,
+  10,
+  15,
+  30,
+  60,
+  120,
+  240,
+  300,
+  600,
+  900,
+  1800,
+  3600,
+  7200,
+  10800
+];
+var SYSTEM_BY_BASE = {
+  20: "G",
+  40: "R",
+  60: "E",
+  80: "J",
+  100: "C",
+  120: "S"
+};
+var KIND_BY_OFFSET = {
+  1: "orbit",
+  2: "clock",
+  3: "combined",
+  4: "highRateClock",
+  5: "codeBias",
+  6: "phaseBias",
+  7: "ura"
+};
+var two = (n) => n.toString().padStart(2, "0");
+function prnFor(system, id) {
+  if (system === "J") return `J${two(id)}`;
+  if (system === "S") return `S${two(19 + id)}`;
+  return `${system}${two(id)}`;
+}
+function ssrUraMm(v) {
+  if (v === 0) return null;
+  const cls = v >> 3 & 7;
+  const val = v & 7;
+  return 3 ** cls * (1 + val / 4) - 1;
+}
+function decodeIgsSsr(frame) {
+  if (frame.messageType !== 4076) return null;
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  const version = r.readU(3);
+  const im = r.readU(8);
+  if (im === 201) return decodeIono(r, version, im);
+  const base = Math.floor(im / 20) * 20;
+  const system = SYSTEM_BY_BASE[base];
+  const kind = KIND_BY_OFFSET[im - base];
+  if (!system || !kind) return null;
+  const epochS = r.readU(20);
+  const updateIntervalS = UPDATE_INTERVAL_S[r.readU(4)];
+  const multipleMessage = r.readU(1) === 1;
+  const iodSsr = r.readU(4);
+  const providerId = r.readU(16);
+  const solutionId = r.readU(4);
+  const hasDatum = kind === "orbit" || kind === "combined";
+  const referenceDatum = hasDatum ? r.readU(1) : void 0;
+  let dispersiveConsistency;
+  let mwConsistency;
+  if (kind === "phaseBias") {
+    dispersiveConsistency = r.readU(1) === 1;
+    mwConsistency = r.readU(1) === 1;
+  }
+  const nsat = r.readU(6);
+  const readOrbit = (s) => {
+    s.iode = r.readU(8);
+    s.deltaRadial = r.readS(22) * 1e-4;
+    s.deltaAlongTrack = r.readS(20) * 4e-4;
+    s.deltaCrossTrack = r.readS(20) * 4e-4;
+    s.dotRadial = r.readS(21) * 1e-6;
+    s.dotAlongTrack = r.readS(19) * 4e-6;
+    s.dotCrossTrack = r.readS(19) * 4e-6;
+  };
+  const readClock2 = (s) => {
+    s.c0 = r.readS(22) * 1e-4;
+    s.c1 = r.readS(21) * 1e-6;
+    s.c2 = r.readS(27) * 2e-8;
+  };
+  const satellites = [];
+  for (let i = 0; i < nsat; i++) {
+    const id = r.readU(6);
+    const sat = { prn: prnFor(system, id) };
+    if (kind === "orbit" || kind === "combined") readOrbit(sat);
+    if (kind === "clock" || kind === "combined") readClock2(sat);
+    else if (kind === "ura")
+      sat.uraMm = ssrUraMm(r.readU(6));
+    else if (kind === "highRateClock")
+      sat.highRateClock = r.readS(22) * 1e-4;
+    else if (kind === "codeBias") {
+      const n = r.readU(5);
+      const biases = [];
+      for (let b = 0; b < n; b++) {
+        const signal = r.readU(5);
+        biases.push({
+          signal,
+          signalName: `#${signal}`,
+          biasM: r.readS(14) * 0.01
+          // IDF025 0.01 m
+        });
+      }
+      sat.codeBiases = biases;
+    } else if (kind === "phaseBias") {
+      const n = r.readU(5);
+      r.readU(9);
+      r.readS(8);
+      const biases = [];
+      for (let b = 0; b < n; b++) {
+        const signal = r.readU(5);
+        const integer = r.readU(1) === 1;
+        const wideLaneGroup = r.readU(2);
+        const discontinuity = r.readU(4);
+        const biasM = r.readS(20) * 1e-4;
+        biases.push({ signal, integer, wideLaneGroup, discontinuity, biasM });
+      }
+      sat.phaseBiases = biases;
+    }
+    satellites.push(sat);
+  }
+  return {
+    messageType: 4076,
+    version,
+    igsMessageNumber: im,
+    system,
+    kind,
+    epochS,
+    updateIntervalS,
+    multipleMessage,
+    referenceDatum,
+    iodSsr,
+    providerId,
+    solutionId,
+    dispersiveConsistency,
+    mwConsistency,
+    satellites
+  };
+}
+function decodeIono(r, version, im) {
+  const epochS = r.readU(20);
+  const updateIntervalS = UPDATE_INTERVAL_S[r.readU(4)];
+  const multipleMessage = r.readU(1) === 1;
+  const iodSsr = r.readU(4);
+  const providerId = r.readU(16);
+  const solutionId = r.readU(4);
+  const qualityTecu = r.readU(9) * 0.05;
+  const nLayers = r.readU(2) + 1;
+  const layers = [];
+  for (let l = 0; l < nLayers; l++) {
+    const heightKm = r.readU(8) * 10;
+    const degree = r.readU(4) + 1;
+    const order = r.readU(4) + 1;
+    const N = degree;
+    const M = order;
+    const nCoef = (N + 1) * (N + 1) - (N - M) * (N - M + 1);
+    const nCos = nCoef;
+    let nSin = 0;
+    for (let m = 1; m <= M; m++) nSin += N - m + 1;
+    for (let i = 0; i < nCos + nSin; i++) r.readS(16);
+    layers.push({ heightKm, degree, order, coefficientCount: nCos + nSin });
+  }
+  return {
+    messageType: 4076,
+    version,
+    igsMessageNumber: im,
+    system: "I",
+    kind: "ionosphere",
+    epochS,
+    updateIntervalS,
+    multipleMessage,
+    iodSsr,
+    providerId,
+    solutionId,
+    satellites: [],
+    iono: { qualityTecu, layers }
+  };
+}
+
 // src/rtcm3/legacy-obs.ts
 var PRUNIT_GPS = 299792.458;
 var PRUNIT_GLO = 599584.916;
@@ -5036,7 +5226,7 @@ function legacyLockSec(i) {
   if (i < 127) return 32 * i - 3096;
   return 937;
 }
-var two = (n) => n.toString().padStart(2, "0");
+var two2 = (n) => n.toString().padStart(2, "0");
 function decodeLegacyObs(frame) {
   const type = frame.messageType;
   const gps = type >= 1001 && type <= 1004;
@@ -5122,7 +5312,7 @@ function decodeLegacyObs(frame) {
       signals.push(sig);
     }
     observations.push({
-      prn: `${system}${two(satId)}`,
+      prn: `${system}${two2(satId)}`,
       system,
       signals
     });
@@ -5664,6 +5854,168 @@ function resetGloFreqCache() {
   gloFreqNum.fill(-128);
 }
 
+// src/rtcm3/ssr.ts
+var UPDATE_INTERVAL_S2 = [
+  1,
+  2,
+  5,
+  10,
+  15,
+  30,
+  60,
+  120,
+  240,
+  300,
+  600,
+  900,
+  1800,
+  3600,
+  7200,
+  10800
+];
+var GPS_SIGNAL = {
+  0: "L1 C/A",
+  1: "L1 P",
+  2: "L1 Z-tracking",
+  5: "L2 C/A",
+  6: "L2 semi-codeless",
+  7: "L2C(M)",
+  8: "L2C(L)",
+  9: "L2C(M+L)",
+  10: "L2 P",
+  11: "L2 Z-tracking",
+  14: "L5 I",
+  15: "L5 Q"
+};
+var GLO_SIGNAL = {
+  0: "G1 C/A",
+  1: "G1 P",
+  2: "G2 C/A",
+  3: "G2 P"
+};
+var SSR_KIND = {
+  1057: { system: "G", kind: "orbit" },
+  1058: { system: "G", kind: "clock" },
+  1059: { system: "G", kind: "codeBias" },
+  1060: { system: "G", kind: "combined" },
+  1061: { system: "G", kind: "ura" },
+  1062: { system: "G", kind: "highRateClock" },
+  1063: { system: "R", kind: "orbit" },
+  1064: { system: "R", kind: "clock" },
+  1065: { system: "R", kind: "codeBias" },
+  1066: { system: "R", kind: "combined" },
+  1067: { system: "R", kind: "ura" },
+  1068: { system: "R", kind: "highRateClock" }
+};
+function ssrUraMm2(v) {
+  if (v === 0) return null;
+  const cls = v >> 3 & 7;
+  const val = v & 7;
+  return 3 ** cls * (1 + val / 4) - 1;
+}
+var two3 = (n) => n.toString().padStart(2, "0");
+function decodeSsr(frame) {
+  const spec = SSR_KIND[frame.messageType];
+  if (!spec) return null;
+  const { system, kind } = spec;
+  const gps = system === "G";
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  const epochS = r.readU(gps ? 20 : 17);
+  const updateIntervalS = UPDATE_INTERVAL_S2[r.readU(4)];
+  const multipleMessage = r.readU(1) === 1;
+  const hasDatum = kind === "orbit" || kind === "combined";
+  const referenceDatum = hasDatum ? r.readU(1) : void 0;
+  const iodSsr = r.readU(4);
+  const providerId = r.readU(16);
+  const solutionId = r.readU(4);
+  const nsat = r.readU(6);
+  const sigName = gps ? GPS_SIGNAL : GLO_SIGNAL;
+  const satellites = [];
+  const readOrbit = (sat) => {
+    sat.deltaRadial = r.readS(22) * 1e-4;
+    sat.deltaAlongTrack = r.readS(20) * 4e-4;
+    sat.deltaCrossTrack = r.readS(20) * 4e-4;
+    sat.dotRadial = r.readS(21) * 1e-6;
+    sat.dotAlongTrack = r.readS(19) * 4e-6;
+    sat.dotCrossTrack = r.readS(19) * 4e-6;
+  };
+  const readClock2 = (sat) => {
+    sat.c0 = r.readS(22) * 1e-4;
+    sat.c1 = r.readS(21) * 1e-6;
+    sat.c2 = r.readS(27) * 2e-8;
+  };
+  for (let i = 0; i < nsat; i++) {
+    const satId = gps ? r.readU(6) : r.readU(5);
+    const sat = { prn: `${system}${two3(satId)}` };
+    if (kind === "orbit" || kind === "combined") {
+      sat.iode = r.readU(8);
+      readOrbit(sat);
+    }
+    if (kind === "combined") readClock2(sat);
+    else if (kind === "clock") readClock2(sat);
+    else if (kind === "ura")
+      sat.uraMm = ssrUraMm2(r.readU(6));
+    else if (kind === "highRateClock")
+      sat.highRateClock = r.readS(22) * 1e-4;
+    else if (kind === "codeBias") {
+      const ncb = r.readU(5);
+      const biases = [];
+      for (let c = 0; c < ncb; c++) {
+        const signal = r.readU(5);
+        const biasM = r.readS(14) * 0.01;
+        biases.push({
+          signal,
+          signalName: sigName[signal] ?? `#${signal}`,
+          biasM
+        });
+      }
+      sat.codeBiases = biases;
+    }
+    satellites.push(sat);
+  }
+  return {
+    messageType: frame.messageType,
+    system,
+    kind,
+    epochS,
+    updateIntervalS,
+    multipleMessage,
+    referenceDatum,
+    iodSsr,
+    providerId,
+    solutionId,
+    satellites
+  };
+}
+function isSsrMessage(messageType) {
+  return messageType in SSR_KIND;
+}
+
+// src/rtcm3/system-params.ts
+function decodeSystemParams(frame) {
+  if (frame.messageType !== 1013) return null;
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  const referenceStationId = r.readU(12);
+  const mjd = r.readU(16);
+  const secondsOfDay = r.readU(17);
+  const nm = r.readU(5);
+  const leapSeconds = r.readU(8);
+  const messages = [];
+  for (let i = 0; i < nm; i++) {
+    messages.push({
+      messageId: r.readU(12),
+      // DF055
+      sync: r.readU(1) === 1,
+      // DF056
+      intervalS: r.readU(16) * 0.1
+      // DF057 (0.1 s)
+    });
+  }
+  return { referenceStationId, mjd, secondsOfDay, leapSeconds, messages };
+}
+
 // src/rtcm3/station.ts
 function createStationMeta() {
   return {
@@ -5677,7 +6029,8 @@ function createStationMeta() {
     receiverType: null,
     receiverFirmware: null,
     receiverSerial: null,
-    description: null
+    description: null,
+    physicalRefStationId: null
   };
 }
 function decodeStationARP(payload, msgType, meta) {
@@ -5700,6 +6053,18 @@ function decodeStationARP(payload, msgType, meta) {
     meta.antennaHeight = r.readU(16) * 1e-4;
   }
 }
+function decodeStation1032(payload, meta) {
+  if (payload.length < 19) return;
+  const r = new BitReader(payload);
+  r.skip(12);
+  meta.stationId = r.readU(12);
+  meta.physicalRefStationId = r.readU(12);
+  meta.itrf = r.readU(6);
+  const x = r.readS(38) * 1e-4;
+  const y = r.readS(38) * 1e-4;
+  const z = r.readS(38) * 1e-4;
+  if (x !== 0 || y !== 0 || z !== 0) meta.position = [x, y, z];
+}
 function decodeAntennaDescriptor(payload, msgType, meta) {
   if (payload.length < 5) return;
   const r = new BitReader(payload);
@@ -5780,6 +6145,9 @@ function updateStationMeta(meta, frame) {
       case 1006:
         decodeStationARP(frame.payload, frame.messageType, meta);
         return true;
+      case 1032:
+        decodeStation1032(frame.payload, meta);
+        return true;
       case 1007:
       case 1008:
         decodeAntennaDescriptor(frame.payload, frame.messageType, meta);
@@ -6460,7 +6828,7 @@ function weekOfToe(week, tow, toe) {
   if (toe - tow < -HALF_WEEK3) return week + 1;
   return week;
 }
-var two2 = (n) => String(n).padStart(2, "0");
+var two4 = (n) => String(n).padStart(2, "0");
 function decodeGalInavWords(words) {
   if (words.length < 96) return null;
   const b = words;
@@ -6548,7 +6916,7 @@ function decodeGalInavWords(words) {
   const tocDate = gpsDate2(week, tocs);
   return {
     system: "E",
-    prn: `E${two2(svid)}`,
+    prn: `E${two4(svid)}`,
     toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK5,
     tocDate,
     af0,
@@ -6656,7 +7024,7 @@ function decodeGalFnavPages(pages) {
   const tocDate = gpsDate2(week, tocs);
   return {
     system: "E",
-    prn: `E${two2(svid)}`,
+    prn: `E${two4(svid)}`,
     toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK5,
     tocDate,
     af0,
@@ -6712,7 +7080,7 @@ var GalInavAssembler = class {
     sat.words[type * 16 + 15] = getBitU(page, 124, 8);
     if (type !== 5) return null;
     const eph = decodeGalInavWords(sat.words);
-    if (!eph || eph.prn !== `E${two2(prn)}`) return null;
+    if (!eph || eph.prn !== `E${two4(prn)}`) return null;
     const key = dedupKey(eph);
     if (key === sat.lastKey) return null;
     sat.lastKey = key;
@@ -6739,7 +7107,7 @@ var GalFnavAssembler = class {
     sat.words.set(page.subarray(0, 31), (type - 1) * 31);
     if (type !== 4) return null;
     const eph = decodeGalFnavPages(sat.words);
-    if (!eph || eph.prn !== `E${two2(prn)}`) return null;
+    if (!eph || eph.prn !== `E${two4(prn)}`) return null;
     const key = dedupKey(eph);
     if (key === sat.lastKey) return null;
     sat.lastKey = key;
@@ -6913,7 +7281,7 @@ var GloStringAssembler = class {
 // src/navbits/sbas.ts
 var GPS_EPOCH_MS7 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK3 = 7 * 86400 * 1e3;
-var two3 = (n) => String(n).padStart(2, "0");
+var two5 = (n) => String(n).padStart(2, "0");
 var SBAS_MT_GEO_NAV = 9;
 function sbasMessageType(msg) {
   return getBitU(msg, 8, 6);
@@ -6932,7 +7300,7 @@ function decodeSbasGeoNav(msg, prn, week, tow) {
   const ura = getBitU(msg, 35, 4);
   return {
     system: "S",
-    prn: `S${two3(prn - 100)}`,
+    prn: `S${two5(prn - 100)}`,
     tocDate: new Date(GPS_EPOCH_MS7 + week * MS_PER_WEEK3 + t0Sec * 1e3),
     tauN: getBitS(msg, 206, 12) * 2 ** -31,
     // a_Gf0, clock offset (s)
@@ -6956,7 +7324,7 @@ function decodeSbasGeoNav(msg, prn, week, tow) {
 // src/ubx/rawnav.ts
 var GPS_EPOCH_MS8 = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK7 = 7 * 86400;
-var two4 = (n) => String(n).padStart(2, "0");
+var two6 = (n) => String(n).padStart(2, "0");
 function parseUbxRawNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
@@ -7032,7 +7400,7 @@ function parseUbxRawNav(data, opts = {}) {
         badParity++;
         continue;
       }
-      const eph = bds.push(`C${two4(svId)}`, sf);
+      const eph = bds.push(`C${two6(svId)}`, sf);
       if (eph) ephemerides.push(eph);
     } else if (gnssId === 6) {
       if (p.length < 8 + 16) continue;
@@ -7051,7 +7419,7 @@ function parseUbxRawNav(data, opts = {}) {
         continue;
       }
       if (!refDate) continue;
-      const eph = glo.push(`R${two4(svId)}`, str2, refDate, p[3] - 7);
+      const eph = glo.push(`R${two6(svId)}`, str2, refDate, p[3] - 7);
       if (eph) ephemerides.push(eph);
     } else if (gnssId === 1) {
       if (p.length < 8 + 32) continue;
@@ -7152,21 +7520,21 @@ var SIGNALS = {
   6: { 0: ["R", "1C"], 2: ["R", "2C"] }
   // GLONASS L1OF, L2OF
 };
-function prnFor(gnssId, svId) {
-  const two10 = (n) => String(n).padStart(2, "0");
+function prnFor2(gnssId, svId) {
+  const two12 = (n) => String(n).padStart(2, "0");
   switch (gnssId) {
     case 0:
-      return svId >= 1 && svId <= 32 ? `G${two10(svId)}` : null;
+      return svId >= 1 && svId <= 32 ? `G${two12(svId)}` : null;
     case 1:
-      return svId >= 120 && svId <= 158 ? `S${two10(svId - 100)}` : null;
+      return svId >= 120 && svId <= 158 ? `S${two12(svId - 100)}` : null;
     case 2:
-      return svId >= 1 && svId <= 36 ? `E${two10(svId)}` : null;
+      return svId >= 1 && svId <= 36 ? `E${two12(svId)}` : null;
     case 3:
-      return svId >= 1 && svId <= 63 ? `C${two10(svId)}` : null;
+      return svId >= 1 && svId <= 63 ? `C${two12(svId)}` : null;
     case 5:
-      return svId >= 1 && svId <= 10 ? `J${two10(svId)}` : null;
+      return svId >= 1 && svId <= 10 ? `J${two12(svId)}` : null;
     case 6:
-      return svId >= 1 && svId <= 32 ? `R${two10(svId)}` : null;
+      return svId >= 1 && svId <= 32 ? `R${two12(svId)}` : null;
     default:
       return null;
   }
@@ -7198,7 +7566,7 @@ function parseUbxRawx(data) {
         const svId = data[off + 21];
         const sigId = data[off + 22];
         const trkStat = data[off + 30];
-        const prn = prnFor(gnssId, svId);
+        const prn = prnFor2(gnssId, svId);
         const sig = SIGNALS[gnssId]?.[sigId];
         if (!prn || !sig) continue;
         const prValid = (trkStat & 1) !== 0;
@@ -7269,22 +7637,22 @@ function scanSbfFrames(data, view, onBlock) {
   }
   return badCrc;
 }
-var two5 = (n) => String(n).padStart(2, "0");
+var two7 = (n) => String(n).padStart(2, "0");
 function svidToPrn(svid) {
-  if (svid >= 1 && svid <= 37) return svid <= 32 ? `G${two5(svid)}` : null;
-  if (svid <= 61) return svid - 37 <= 27 ? `R${two5(svid - 37)}` : null;
+  if (svid >= 1 && svid <= 37) return svid <= 32 ? `G${two7(svid)}` : null;
+  if (svid <= 61) return svid - 37 <= 27 ? `R${two7(svid - 37)}` : null;
   if (svid <= 62) return null;
-  if (svid <= 68) return svid - 38 <= 27 ? `R${two5(svid - 38)}` : null;
+  if (svid <= 68) return svid - 38 <= 27 ? `R${two7(svid - 38)}` : null;
   if (svid <= 70) return null;
-  if (svid <= 106) return svid - 70 <= 36 ? `E${two5(svid - 70)}` : null;
+  if (svid <= 106) return svid - 70 <= 36 ? `E${two7(svid - 70)}` : null;
   if (svid <= 119) return null;
-  if (svid <= 140) return `S${two5(svid - 100)}`;
-  if (svid <= 180) return svid - 140 <= 50 ? `C${two5(svid - 140)}` : null;
-  if (svid <= 190) return `J${two5(svid - 180)}`;
-  if (svid <= 197) return svid - 190 <= 14 ? `I${two5(svid - 190)}` : null;
-  if (svid <= 215) return `S${two5(svid - 157)}`;
-  if (svid <= 222) return svid - 208 <= 14 ? `I${two5(svid - 208)}` : null;
-  if (svid <= 245) return svid - 182 <= 50 ? `C${two5(svid - 182)}` : null;
+  if (svid <= 140) return `S${two7(svid - 100)}`;
+  if (svid <= 180) return svid - 140 <= 50 ? `C${two7(svid - 140)}` : null;
+  if (svid <= 190) return `J${two7(svid - 180)}`;
+  if (svid <= 197) return svid - 190 <= 14 ? `I${two7(svid - 190)}` : null;
+  if (svid <= 215) return `S${two7(svid - 157)}`;
+  if (svid <= 222) return svid - 208 <= 14 ? `I${two7(svid - 208)}` : null;
+  if (svid <= 245) return svid - 182 <= 50 ? `C${two7(svid - 182)}` : null;
   return null;
 }
 
@@ -8896,23 +9264,23 @@ var MEAS2_SIG = [
   ["J", "1E"],
   ["J", "5P"]
 ];
-var two6 = (n) => String(n).padStart(2, "0");
+var two8 = (n) => String(n).padStart(2, "0");
 function m3Prn(navsys, svid) {
   switch (navsys) {
     case 0:
-      return svid < 32 ? `G${two6(svid + 1)}` : null;
+      return svid < 32 ? `G${two8(svid + 1)}` : null;
     case 1:
-      return svid < 27 ? `R${two6(svid + 1)}` : null;
+      return svid < 27 ? `R${two8(svid + 1)}` : null;
     case 2:
-      return svid < 36 ? `E${two6(svid + 1)}` : null;
+      return svid < 36 ? `E${two8(svid + 1)}` : null;
     case 3:
-      return svid < 50 ? `C${two6(svid + 1)}` : null;
+      return svid < 50 ? `C${two8(svid + 1)}` : null;
     case 4:
-      return svid <= 38 ? `S${two6(svid + 20)}` : null;
+      return svid <= 38 ? `S${two8(svid + 20)}` : null;
     case 5:
-      return svid < 10 ? `J${two6(svid + 1)}` : null;
+      return svid < 10 ? `J${two8(svid + 1)}` : null;
     case 6:
-      return svid < 14 ? `I${two6(svid + 1)}` : null;
+      return svid < 14 ? `I${two8(svid + 1)}` : null;
     default:
       return null;
   }
@@ -10746,7 +11114,7 @@ var MS_PER_WEEK12 = SEC_PER_WEEK11 * 1e3;
 var SC2RAD = GPS_PI;
 var gpsMs4 = (week, sec) => GPS_EPOCH_MS20 + week * MS_PER_WEEK12 + sec * 1e3;
 var sowOf5 = (dateMs) => (dateMs / 1e3 % SEC_PER_WEEK11 + SEC_PER_WEEK11) % SEC_PER_WEEK11;
-var two7 = (n) => String(n).padStart(2, "0");
+var two9 = (n) => String(n).padStart(2, "0");
 var Reader = class {
   constructor(view, start, le) {
     this.view = view;
@@ -10841,7 +11209,7 @@ function decodeGpsQzss(view, p, le, sys) {
   const tocDate = new Date(gpsMs4(week, toes));
   return {
     system: sys,
-    prn: `${sys}${two7(prn)}`,
+    prn: `${sys}${two9(prn)}`,
     toc: sowOf5(tocDate.getTime()),
     tocDate,
     af0,
@@ -10874,7 +11242,7 @@ function decodeGalileo(view, p, le) {
   const tocDate = new Date(gpsMs4(week, toes));
   return {
     system: "E",
-    prn: `E${two7(prn)}`,
+    prn: `E${two9(prn)}`,
     toc: sowOf5(tocDate.getTime()),
     tocDate,
     af0,
@@ -10912,7 +11280,7 @@ function decodeBeidou(view, p, le) {
   const tocDate = new Date(BDT_EPOCH_MS5 + week * MS_PER_WEEK12 + toes * 1e3);
   return {
     system: "C",
-    prn: `C${two7(prn)}`,
+    prn: `C${two9(prn)}`,
     toc: sowOf5(tocDate.getTime()),
     tocDate,
     af0,
@@ -10952,7 +11320,7 @@ function decodeGlonass(view, p, le) {
   );
   return {
     system: "R",
-    prn: `R${two7(slot)}`,
+    prn: `R${two9(slot)}`,
     tocDate,
     tauN,
     gammaN,
@@ -10991,7 +11359,7 @@ function decodeSbas(view, p, le) {
   if (prn < 120 || prn > 158) return null;
   return {
     system: "S",
-    prn: `S${two7(prn - 100)}`,
+    prn: `S${two9(prn - 100)}`,
     tocDate: new Date(gpsMs4(week, toe)),
     tauN: af0,
     gammaN: af1,
@@ -11168,22 +11536,22 @@ var OBS_CODES = {
   ]
   // QZSS
 };
-var two8 = (n) => String(n).padStart(2, "0");
+var two10 = (n) => String(n).padStart(2, "0");
 function obsPrn(sys, prn) {
   switch (sys) {
     case 0:
-      return prn >= 1 && prn <= 32 ? `G${two8(prn)}` : null;
+      return prn >= 1 && prn <= 32 ? `G${two10(prn)}` : null;
     case 1:
-      return prn >= 1 && prn <= 27 ? `R${two8(prn)}` : null;
+      return prn >= 1 && prn <= 27 ? `R${two10(prn)}` : null;
     case 2:
-      return prn >= 120 && prn <= 158 ? `S${two8(prn - 100)}` : null;
+      return prn >= 120 && prn <= 158 ? `S${two10(prn - 100)}` : null;
     case 3:
-      return prn >= 1 && prn <= 36 ? `E${two8(prn)}` : null;
+      return prn >= 1 && prn <= 36 ? `E${two10(prn)}` : null;
     case 4:
-      return prn >= 1 && prn <= 63 ? `C${two8(prn)}` : null;
+      return prn >= 1 && prn <= 63 ? `C${two10(prn)}` : null;
     case 5:
-      if (prn >= 193 && prn <= 202) return `J${two8(prn - 192)}`;
-      return prn >= 1 && prn <= 10 ? `J${two8(prn)}` : null;
+      if (prn >= 193 && prn <= 202) return `J${two10(prn - 192)}`;
+      return prn >= 1 && prn <= 10 ? `J${two10(prn)}` : null;
     default:
       return null;
   }
@@ -13536,7 +13904,7 @@ var SEC_PER_WEEK12 = 604800;
 var MAXSBSAGEF = 30;
 var MAXSBSAGEL = 360;
 var MAXBAND = 10;
-var two9 = (n) => String(n).padStart(2, "0");
+var two11 = (n) => String(n).padStart(2, "0");
 var VAR_FCORR = [
   0.052,
   0.0924,
@@ -13682,10 +14050,10 @@ var igpband2 = [
 ];
 var gpsSeconds = (week, tow) => week * SEC_PER_WEEK12 + tow;
 function maskPrn(i) {
-  if (i >= 1 && i <= 37) return `G${two9(i)}`;
-  if (i >= 38 && i <= 61) return `R${two9(i - 37)}`;
-  if (i >= 120 && i <= 138) return `S${two9(i - 100)}`;
-  if (i >= 193 && i <= 202) return `J${two9(i - 192)}`;
+  if (i >= 1 && i <= 37) return `G${two11(i)}`;
+  if (i >= 38 && i <= 61) return `R${two11(i - 37)}`;
+  if (i >= 120 && i <= 138) return `S${two11(i - 100)}`;
+  if (i >= 193 && i <= 202) return `J${two11(i - 192)}`;
   return null;
 }
 function ionppp(latRad, lonRad, azRad, elRad, re, hion) {
@@ -13724,7 +14092,7 @@ var SbasProcessor = class {
     if (week === 0) return -1;
     const type = sbasMessageType(msg);
     if (geoPrn != null && geoPrn >= 120 && geoPrn <= 158)
-      this.geoPrns.add(`S${two9(geoPrn - 100)}`);
+      this.geoPrns.add(`S${two11(geoPrn - 100)}`);
     switch (type) {
       case 1:
         return this.decodeMask(msg) ? type : -1;
@@ -18273,9 +18641,12 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   dateToEpoch,
   decodeBinexEph,
   decodeEphemeris,
+  decodeIgsSsr,
   decodeLegacyObs,
   decodeMsmFull,
   decodeSbfNavigation,
+  decodeSsr,
+  decodeSystemParams,
   deg2dms,
   deg2rad,
   detrendIonoArcs,
@@ -18368,6 +18739,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ionoFree,
   isSbasL1Preamble,
   isSbasL5Correction,
+  isSsrMessage,
   keplerPosition,
   klobucharDelay,
   lambdaReduction,