gnss-js 1.62.0 → 1.64.0

package/dist/index.cjs

@@ -144,8 +144,11 @@ __export(src_exports, {
   dateToEpoch: () => dateToEpoch,
   decodeBinexEph: () => decodeBinexEph,
   decodeEphemeris: () => decodeEphemeris,
+  decodeIgsSsr: () => decodeIgsSsr,
+  decodeLegacyObs: () => decodeLegacyObs,
   decodeMsmFull: () => decodeMsmFull,
   decodeSbfNavigation: () => decodeSbfNavigation,
+  decodeSsr: () => decodeSsr,
   deg2dms: () => deg2dms,
   deg2rad: () => deg2rad,
   detrendIonoArcs: () => detrendIonoArcs,
@@ -238,6 +241,7 @@ __export(src_exports, {
   ionoFree: () => ionoFree,
   isSbasL1Preamble: () => isSbasL1Preamble,
   isSbasL5Correction: () => isSbasL5Correction,
+  isSsrMessage: () => isSsrMessage,
   keplerPosition: () => keplerPosition,
   klobucharDelay: () => klobucharDelay,
   lambdaReduction: () => lambdaReduction,
@@ -5019,6 +5023,302 @@ 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;
+var LAM_GPS_L1 = C_LIGHT / 157542e4;
+var LAM_GPS_L2 = C_LIGHT / 12276e5;
+var PPR_INVALID = -(2 ** 19);
+var PR21_INVALID = -(2 ** 13);
+function legacyLockSec(i) {
+  if (i < 24) return i;
+  if (i < 48) return 2 * i - 24;
+  if (i < 72) return 4 * i - 120;
+  if (i < 96) return 8 * i - 408;
+  if (i < 120) return 16 * i - 1176;
+  if (i < 127) return 32 * i - 3096;
+  return 937;
+}
+var two2 = (n) => n.toString().padStart(2, "0");
+function decodeLegacyObs(frame) {
+  const type = frame.messageType;
+  const gps = type >= 1001 && type <= 1004;
+  const glo = type >= 1009 && type <= 1012;
+  if (!gps && !glo) return null;
+  const ext = type === 1002 || type === 1004 || type === 1010 || type === 1012;
+  const dual = type === 1003 || type === 1004 || type === 1011 || type === 1012;
+  const system = gps ? "G" : "R";
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  r.readU(12);
+  let epochMs;
+  let nsat;
+  if (gps) {
+    epochMs = r.readU(30);
+    r.readU(1);
+    nsat = r.readU(5);
+    r.readU(1);
+    r.readU(3);
+  } else {
+    const tk = r.readU(27);
+    r.readU(1);
+    nsat = r.readU(5);
+    r.readU(1);
+    r.readU(3);
+    const gloDow = new Date(Date.now() + 3 * 36e5).getUTCDay();
+    epochMs = (gloDow << 27 | tk) >>> 0;
+  }
+  const observations = [];
+  for (let s = 0; s < nsat; s++) {
+    const satId = r.readU(6);
+    const code1 = r.readU(1);
+    const freqCh = glo ? r.readU(5) : 0;
+    const pr1raw = r.readU(gps ? 24 : 25);
+    const ppr1 = r.readS(20);
+    const lock1 = r.readU(7);
+    let amb = 0;
+    let cnr1 = 0;
+    if (ext) {
+      amb = r.readU(gps ? 8 : 7);
+      cnr1 = r.readU(8);
+    }
+    let code2 = 0;
+    let pr21 = 0;
+    let ppr2 = 0;
+    let lock2 = 0;
+    let cnr2 = 0;
+    if (dual) {
+      code2 = r.readU(2);
+      pr21 = r.readS(14);
+      ppr2 = r.readS(20);
+      lock2 = r.readU(7);
+      if (ext) cnr2 = r.readU(8);
+    }
+    if (gps && (satId < 1 || satId > 32)) continue;
+    if (glo && (satId < 1 || satId > 24)) continue;
+    const prunit = gps ? PRUNIT_GPS : PRUNIT_GLO;
+    const pr1 = pr1raw * 0.02 + amb * prunit;
+    const k = freqCh - 7;
+    const lamL1 = gps ? LAM_GPS_L1 : C_LIGHT / (GLO_F1_BASE + k * GLO_F1_STEP);
+    const lamL2 = gps ? LAM_GPS_L2 : C_LIGHT / (GLO_F2_BASE + k * GLO_F2_STEP);
+    const signals = [];
+    {
+      const sig = {
+        rinexCode: gps ? code1 ? "1W" : "1C" : code1 ? "1P" : "1C",
+        pseudorange: pr1,
+        wavelength: lamL1,
+        lockTime: legacyLockSec(lock1)
+      };
+      if (ppr1 !== PPR_INVALID) sig.phase = (pr1 + ppr1 * 5e-4) / lamL1;
+      if (ext && cnr1 > 0) sig.cn0 = cnr1 * 0.25;
+      signals.push(sig);
+    }
+    if (dual) {
+      const sig = {
+        rinexCode: gps ? code2 === 0 ? "2X" : "2W" : code2 === 0 ? "2C" : "2P",
+        wavelength: lamL2,
+        lockTime: legacyLockSec(lock2)
+      };
+      if (pr21 !== PR21_INVALID) sig.pseudorange = pr1 + pr21 * 0.02;
+      if (ppr2 !== PPR_INVALID) sig.phase = (pr1 + ppr2 * 5e-4) / lamL2;
+      if (ext && cnr2 > 0) sig.cn0 = cnr2 * 0.25;
+      signals.push(sig);
+    }
+    observations.push({
+      prn: `${system}${two2(satId)}`,
+      system,
+      signals
+    });
+  }
+  return { messageType: type, epochMs, system, observations };
+}
+
 // src/rtcm3/msm.ts
 var GPS_L1 = FREQ.G["1"];
 var GPS_L2 = FREQ.G["2"];
@@ -5553,6 +5853,144 @@ 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/station.ts
 function createStationMeta() {
   return {
@@ -6349,7 +6787,7 @@ function weekOfToe(week, tow, toe) {
   if (toe - tow < -HALF_WEEK3) return week + 1;
   return week;
 }
-var two = (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;
@@ -6437,7 +6875,7 @@ function decodeGalInavWords(words) {
   const tocDate = gpsDate2(week, tocs);
   return {
     system: "E",
-    prn: `E${two(svid)}`,
+    prn: `E${two4(svid)}`,
     toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK5,
     tocDate,
     af0,
@@ -6545,7 +6983,7 @@ function decodeGalFnavPages(pages) {
   const tocDate = gpsDate2(week, tocs);
   return {
     system: "E",
-    prn: `E${two(svid)}`,
+    prn: `E${two4(svid)}`,
     toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK5,
     tocDate,
     af0,
@@ -6601,7 +7039,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${two(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;
@@ -6628,7 +7066,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${two(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;
@@ -6802,7 +7240,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 two2 = (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);
@@ -6821,7 +7259,7 @@ function decodeSbasGeoNav(msg, prn, week, tow) {
   const ura = getBitU(msg, 35, 4);
   return {
     system: "S",
-    prn: `S${two2(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)
@@ -6845,7 +7283,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 two3 = (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 = [];
@@ -6921,7 +7359,7 @@ function parseUbxRawNav(data, opts = {}) {
         badParity++;
         continue;
       }
-      const eph = bds.push(`C${two3(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;
@@ -6940,7 +7378,7 @@ function parseUbxRawNav(data, opts = {}) {
         continue;
       }
       if (!refDate) continue;
-      const eph = glo.push(`R${two3(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;
@@ -7041,21 +7479,21 @@ var SIGNALS = {
   6: { 0: ["R", "1C"], 2: ["R", "2C"] }
   // GLONASS L1OF, L2OF
 };
-function prnFor(gnssId, svId) {
-  const two9 = (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${two9(svId)}` : null;
+      return svId >= 1 && svId <= 32 ? `G${two12(svId)}` : null;
     case 1:
-      return svId >= 120 && svId <= 158 ? `S${two9(svId - 100)}` : null;
+      return svId >= 120 && svId <= 158 ? `S${two12(svId - 100)}` : null;
     case 2:
-      return svId >= 1 && svId <= 36 ? `E${two9(svId)}` : null;
+      return svId >= 1 && svId <= 36 ? `E${two12(svId)}` : null;
     case 3:
-      return svId >= 1 && svId <= 63 ? `C${two9(svId)}` : null;
+      return svId >= 1 && svId <= 63 ? `C${two12(svId)}` : null;
     case 5:
-      return svId >= 1 && svId <= 10 ? `J${two9(svId)}` : null;
+      return svId >= 1 && svId <= 10 ? `J${two12(svId)}` : null;
     case 6:
-      return svId >= 1 && svId <= 32 ? `R${two9(svId)}` : null;
+      return svId >= 1 && svId <= 32 ? `R${two12(svId)}` : null;
     default:
       return null;
   }
@@ -7087,7 +7525,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;
@@ -7158,22 +7596,22 @@ function scanSbfFrames(data, view, onBlock) {
   }
   return badCrc;
 }
-var two4 = (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${two4(svid)}` : null;
-  if (svid <= 61) return svid - 37 <= 27 ? `R${two4(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${two4(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${two4(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${two4(svid - 100)}`;
-  if (svid <= 180) return svid - 140 <= 50 ? `C${two4(svid - 140)}` : null;
-  if (svid <= 190) return `J${two4(svid - 180)}`;
-  if (svid <= 197) return svid - 190 <= 14 ? `I${two4(svid - 190)}` : null;
-  if (svid <= 215) return `S${two4(svid - 157)}`;
-  if (svid <= 222) return svid - 208 <= 14 ? `I${two4(svid - 208)}` : null;
-  if (svid <= 245) return svid - 182 <= 50 ? `C${two4(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;
 }
 
@@ -8785,23 +9223,23 @@ var MEAS2_SIG = [
   ["J", "1E"],
   ["J", "5P"]
 ];
-var two5 = (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${two5(svid + 1)}` : null;
+      return svid < 32 ? `G${two8(svid + 1)}` : null;
     case 1:
-      return svid < 27 ? `R${two5(svid + 1)}` : null;
+      return svid < 27 ? `R${two8(svid + 1)}` : null;
     case 2:
-      return svid < 36 ? `E${two5(svid + 1)}` : null;
+      return svid < 36 ? `E${two8(svid + 1)}` : null;
     case 3:
-      return svid < 50 ? `C${two5(svid + 1)}` : null;
+      return svid < 50 ? `C${two8(svid + 1)}` : null;
     case 4:
-      return svid <= 38 ? `S${two5(svid + 20)}` : null;
+      return svid <= 38 ? `S${two8(svid + 20)}` : null;
     case 5:
-      return svid < 10 ? `J${two5(svid + 1)}` : null;
+      return svid < 10 ? `J${two8(svid + 1)}` : null;
     case 6:
-      return svid < 14 ? `I${two5(svid + 1)}` : null;
+      return svid < 14 ? `I${two8(svid + 1)}` : null;
     default:
       return null;
   }
@@ -10635,7 +11073,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 two6 = (n) => String(n).padStart(2, "0");
+var two9 = (n) => String(n).padStart(2, "0");
 var Reader = class {
   constructor(view, start, le) {
     this.view = view;
@@ -10730,7 +11168,7 @@ function decodeGpsQzss(view, p, le, sys) {
   const tocDate = new Date(gpsMs4(week, toes));
   return {
     system: sys,
-    prn: `${sys}${two6(prn)}`,
+    prn: `${sys}${two9(prn)}`,
     toc: sowOf5(tocDate.getTime()),
     tocDate,
     af0,
@@ -10763,7 +11201,7 @@ function decodeGalileo(view, p, le) {
   const tocDate = new Date(gpsMs4(week, toes));
   return {
     system: "E",
-    prn: `E${two6(prn)}`,
+    prn: `E${two9(prn)}`,
     toc: sowOf5(tocDate.getTime()),
     tocDate,
     af0,
@@ -10801,7 +11239,7 @@ function decodeBeidou(view, p, le) {
   const tocDate = new Date(BDT_EPOCH_MS5 + week * MS_PER_WEEK12 + toes * 1e3);
   return {
     system: "C",
-    prn: `C${two6(prn)}`,
+    prn: `C${two9(prn)}`,
     toc: sowOf5(tocDate.getTime()),
     tocDate,
     af0,
@@ -10841,7 +11279,7 @@ function decodeGlonass(view, p, le) {
   );
   return {
     system: "R",
-    prn: `R${two6(slot)}`,
+    prn: `R${two9(slot)}`,
     tocDate,
     tauN,
     gammaN,
@@ -10880,7 +11318,7 @@ function decodeSbas(view, p, le) {
   if (prn < 120 || prn > 158) return null;
   return {
     system: "S",
-    prn: `S${two6(prn - 100)}`,
+    prn: `S${two9(prn - 100)}`,
     tocDate: new Date(gpsMs4(week, toe)),
     tauN: af0,
     gammaN: af1,
@@ -11057,22 +11495,22 @@ var OBS_CODES = {
   ]
   // QZSS
 };
-var two7 = (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${two7(prn)}` : null;
+      return prn >= 1 && prn <= 32 ? `G${two10(prn)}` : null;
     case 1:
-      return prn >= 1 && prn <= 27 ? `R${two7(prn)}` : null;
+      return prn >= 1 && prn <= 27 ? `R${two10(prn)}` : null;
     case 2:
-      return prn >= 120 && prn <= 158 ? `S${two7(prn - 100)}` : null;
+      return prn >= 120 && prn <= 158 ? `S${two10(prn - 100)}` : null;
     case 3:
-      return prn >= 1 && prn <= 36 ? `E${two7(prn)}` : null;
+      return prn >= 1 && prn <= 36 ? `E${two10(prn)}` : null;
     case 4:
-      return prn >= 1 && prn <= 63 ? `C${two7(prn)}` : null;
+      return prn >= 1 && prn <= 63 ? `C${two10(prn)}` : null;
     case 5:
-      if (prn >= 193 && prn <= 202) return `J${two7(prn - 192)}`;
-      return prn >= 1 && prn <= 10 ? `J${two7(prn)}` : null;
+      if (prn >= 193 && prn <= 202) return `J${two10(prn - 192)}`;
+      return prn >= 1 && prn <= 10 ? `J${two10(prn)}` : null;
     default:
       return null;
   }
@@ -13425,7 +13863,7 @@ var SEC_PER_WEEK12 = 604800;
 var MAXSBSAGEF = 30;
 var MAXSBSAGEL = 360;
 var MAXBAND = 10;
-var two8 = (n) => String(n).padStart(2, "0");
+var two11 = (n) => String(n).padStart(2, "0");
 var VAR_FCORR = [
   0.052,
   0.0924,
@@ -13571,10 +14009,10 @@ var igpband2 = [
 ];
 var gpsSeconds = (week, tow) => week * SEC_PER_WEEK12 + tow;
 function maskPrn(i) {
-  if (i >= 1 && i <= 37) return `G${two8(i)}`;
-  if (i >= 38 && i <= 61) return `R${two8(i - 37)}`;
-  if (i >= 120 && i <= 138) return `S${two8(i - 100)}`;
-  if (i >= 193 && i <= 202) return `J${two8(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) {
@@ -13613,7 +14051,7 @@ var SbasProcessor = class {
     if (week === 0) return -1;
     const type = sbasMessageType(msg);
     if (geoPrn != null && geoPrn >= 120 && geoPrn <= 158)
-      this.geoPrns.add(`S${two8(geoPrn - 100)}`);
+      this.geoPrns.add(`S${two11(geoPrn - 100)}`);
     switch (type) {
       case 1:
         return this.decodeMask(msg) ? type : -1;
@@ -18162,8 +18600,11 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   dateToEpoch,
   decodeBinexEph,
   decodeEphemeris,
+  decodeIgsSsr,
+  decodeLegacyObs,
   decodeMsmFull,
   decodeSbfNavigation,
+  decodeSsr,
   deg2dms,
   deg2rad,
   detrendIonoArcs,
@@ -18256,6 +18697,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ionoFree,
   isSbasL1Preamble,
   isSbasL5Correction,
+  isSsrMessage,
   keplerPosition,
   klobucharDelay,
   lambdaReduction,