gnss-js 1.15.0 → 1.17.0

package/dist/index.cjs

@@ -222,11 +222,14 @@ __export(src_exports, {
   parseNovatelRange: () => parseNovatelRange,
   parseRinexStream: () => parseRinexStream,
   parseSbfAlmanac: () => parseSbfAlmanac,
+  parseSbfIonoUtc: () => parseSbfIonoUtc,
   parseSbfMeas: () => parseSbfMeas,
   parseSbfNav: () => parseSbfNav,
   parseSinexBiasDcb: () => parseSinexBiasDcb,
   parseSourcetable: () => parseSourcetable,
   parseSp3: () => parseSp3,
+  parseUbxIonoUtc: () => parseUbxIonoUtc,
+  parseUbxNav: () => parseUbxNav,
   parseUbxRawx: () => parseUbxRawx,
   phiAltBOC: () => phiAltBOC,
   phiBOCc: () => phiBOCc,
@@ -240,6 +243,8 @@ __export(src_exports, {
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   satClockCorrection: () => satClockCorrection,
+  sbfCrc16: () => crc16,
+  scanSbfFrames: () => scanSbfFrames,
   selectEphemeris: () => selectEphemeris,
   setGloFreqNumber: () => setGloFreqNumber,
   setRtcm3DebugHandler: () => setRtcm3DebugHandler,
@@ -249,6 +254,7 @@ __export(src_exports, {
   systemCmp: () => systemCmp,
   systemName: () => systemName,
   transformFrame: () => transformFrame,
+  ubxFrames: () => ubxFrames,
   updateStationMeta: () => updateStationMeta,
   updateStreamStats: () => updateStreamStats,
   verifyChecksum: () => verifyChecksum,
@@ -4598,8 +4604,288 @@ function updateStreamStats(stats, frames, rawBytes) {
   }
 }
 
-// src/ubx/index.ts
+// src/ubx/frame.ts
+function* ubxFrames(data, stats = { badChecksums: 0 }) {
+  let i = 0;
+  while (i + 8 <= data.length) {
+    if (data[i] !== 181 || data[i + 1] !== 98) {
+      i++;
+      continue;
+    }
+    const len = data[i + 4] | data[i + 5] << 8;
+    const end = i + 6 + len + 2;
+    if (end > data.length) break;
+    let ckA = 0;
+    let ckB = 0;
+    for (let j = i + 2; j < i + 6 + len; j++) {
+      ckA = ckA + data[j] & 255;
+      ckB = ckB + ckA & 255;
+    }
+    if (ckA !== data[i + 6 + len] || ckB !== data[i + 7 + len]) {
+      stats.badChecksums++;
+      i++;
+      continue;
+    }
+    yield {
+      msgClass: data[i + 2],
+      msgId: data[i + 3],
+      payload: data.subarray(i + 6, i + 6 + len),
+      payloadStart: i + 6
+    };
+    i = end;
+  }
+}
+
+// src/navbits/index.ts
+var GPS_PI = 3.1415926535898;
 var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
+var SEC_PER_WEEK = 7 * 86400;
+function getBitU(buff, pos, len) {
+  let bits = 0;
+  for (let i = pos; i < pos + len; i++) {
+    bits = bits * 2 + (buff[i >> 3] >> 7 - (i & 7) & 1);
+  }
+  return bits;
+}
+function getBitS(buff, pos, len) {
+  const bits = getBitU(buff, pos, len);
+  if (len <= 0 || bits < 2 ** (len - 1)) return bits;
+  return bits - 2 ** len;
+}
+function setBitU(buff, pos, len, data) {
+  let mask = 2 ** (len - 1);
+  for (let i = pos; i < pos + len; i++, mask /= 2) {
+    if (data >= mask) {
+      buff[i >> 3] |= 1 << 7 - (i & 7);
+      data -= mask;
+    } else {
+      buff[i >> 3] &= ~(1 << 7 - (i & 7));
+    }
+  }
+}
+function decodeGpsLnavFrame(subframes, opts = {}) {
+  if (subframes.length < 90) return null;
+  const b = subframes;
+  let i = 24;
+  const tow1 = getBitU(b, i, 17) * 6;
+  i += 17 + 2;
+  const id1 = getBitU(b, i, 3);
+  i += 3 + 2;
+  const week10 = getBitU(b, i, 10);
+  i += 10;
+  i += 2;
+  i += 4;
+  const svHealth = getBitU(b, i, 6);
+  i += 6;
+  const iodc0 = getBitU(b, i, 2);
+  i += 2;
+  i += 1 + 87;
+  const tgdRaw = getBitS(b, i, 8);
+  i += 8;
+  const iodc1 = getBitU(b, i, 8);
+  i += 8;
+  const tocSec = getBitU(b, i, 16) * 16;
+  i += 16;
+  const af2 = getBitS(b, i, 8) * 2 ** -55;
+  i += 8;
+  const af1 = getBitS(b, i, 16) * 2 ** -43;
+  i += 16;
+  const af0 = getBitS(b, i, 22) * 2 ** -31;
+  i = 240 + 24;
+  i += 17 + 2;
+  const id2 = getBitU(b, i, 3);
+  i += 3 + 2;
+  const iode = getBitU(b, i, 8);
+  i += 8;
+  const crs = getBitS(b, i, 16) * 2 ** -5;
+  i += 16;
+  const deltaN = getBitS(b, i, 16) * 2 ** -43 * GPS_PI;
+  i += 16;
+  const m0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
+  i += 32;
+  const cuc = getBitS(b, i, 16) * 2 ** -29;
+  i += 16;
+  const e = getBitU(b, i, 32) * 2 ** -33;
+  i += 32;
+  const cus = getBitS(b, i, 16) * 2 ** -29;
+  i += 16;
+  const sqrtA = getBitU(b, i, 32) * 2 ** -19;
+  i += 32;
+  const toes = getBitU(b, i, 16) * 16;
+  i = 480 + 24;
+  i += 17 + 2;
+  const id3 = getBitU(b, i, 3);
+  i += 3 + 2;
+  const cic = getBitS(b, i, 16) * 2 ** -29;
+  i += 16;
+  const omega0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
+  i += 32;
+  const cis = getBitS(b, i, 16) * 2 ** -29;
+  i += 16;
+  const i0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
+  i += 32;
+  const crc = getBitS(b, i, 16) * 2 ** -5;
+  i += 16;
+  const omega = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
+  i += 32;
+  const omegaDot = getBitS(b, i, 24) * 2 ** -43 * GPS_PI;
+  i += 24;
+  const iode3 = getBitU(b, i, 8);
+  i += 8;
+  const idot = getBitS(b, i, 14) * 2 ** -43 * GPS_PI;
+  const iodc = iodc0 * 256 + iodc1;
+  if (id1 !== 1 || id2 !== 2 || id3 !== 3) return null;
+  if (iode3 !== iode || iode !== (iodc & 255)) return null;
+  const ref = opts.refWeek ?? Math.floor((Date.now() - GPS_EPOCH_MS) / 1e3 / SEC_PER_WEEK);
+  let week = week10 + 1024 * Math.round((ref - week10) / 1024);
+  if (toes < tow1 - 302400) week++;
+  else if (toes > tow1 + 302400) week--;
+  const prn = opts.prn ?? "G00";
+  const tocDate = new Date(
+    GPS_EPOCH_MS + (week * SEC_PER_WEEK + tocSec) * 1e3
+  );
+  return {
+    system: prn[0] === "J" ? "J" : "G",
+    prn,
+    tocDate,
+    // Same seconds-of-week convention as parseNavFile (rinex/nav.ts).
+    toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK,
+    af0,
+    af1,
+    af2,
+    iode,
+    crs,
+    deltaN,
+    m0,
+    cuc,
+    e,
+    cus,
+    sqrtA,
+    toe: toes,
+    cic,
+    omega0,
+    cis,
+    i0,
+    crc,
+    omega,
+    omegaDot,
+    idot,
+    week,
+    svHealth,
+    tgd: tgdRaw === -128 ? 0 : tgdRaw * 2 ** -31
+    // IS-GPS-200: -128 reserved
+  };
+}
+
+// src/ubx/nav.ts
+var PREAMB_CNAV = 139;
+function readLnavSubframe(view, f) {
+  const p = f.payload;
+  if (p.length < 8) return null;
+  const gnssId = p[0];
+  const svId = p[1];
+  let prn;
+  if (gnssId === 0 && svId >= 1 && svId <= 32) {
+    prn = `G${String(svId).padStart(2, "0")}`;
+  } else if (gnssId === 5 && svId >= 1 && svId <= 10) {
+    if (p.length === 44) return null;
+    prn = `J${String(svId).padStart(2, "0")}`;
+  } else {
+    return null;
+  }
+  if (p.length < 8 + 40) return null;
+  const base = f.payloadStart + 8;
+  if (view.getUint32(base, true) >>> 24 === PREAMB_CNAV) return null;
+  const buff = new Uint8Array(30);
+  for (let k = 0; k < 10; k++) {
+    const dwrd = view.getUint32(base + 4 * k, true);
+    setBitU(buff, 24 * k, 24, dwrd >>> 6 & 16777215);
+  }
+  return { prn, gnssId, svId, buff, id: getBitU(buff, 43, 3) };
+}
+function parseUbxNav(data, opts = {}) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const ephemerides = [];
+  let badParity = 0;
+  let refWeek = opts.refWeek;
+  if (refWeek === void 0) {
+    for (const f of ubxFrames(data)) {
+      if (f.msgClass === 2 && f.msgId === 21 && f.payload.length >= 16) {
+        const week = view.getUint16(f.payloadStart + 8, true);
+        if (week > 0) {
+          refWeek = week;
+          break;
+        }
+      }
+    }
+  }
+  if (refWeek === void 0) return { ephemerides, badParity };
+  const subframes = /* @__PURE__ */ new Map();
+  const last = /* @__PURE__ */ new Map();
+  for (const f of ubxFrames(data)) {
+    if (f.msgClass !== 2 || f.msgId !== 19) continue;
+    const lnav = readLnavSubframe(view, f);
+    if (!lnav) continue;
+    const { prn, gnssId, svId, buff, id } = lnav;
+    if (id < 1 || id > 5) {
+      badParity++;
+      continue;
+    }
+    if (id > 3) continue;
+    const key = gnssId * 256 + svId;
+    let sf = subframes.get(key);
+    if (!sf) {
+      sf = { buf: new Uint8Array(90), have: 0 };
+      subframes.set(key, sf);
+    }
+    sf.buf.set(buff, (id - 1) * 30);
+    sf.have |= 1 << id - 1;
+    if (id !== 3 || sf.have !== 7) continue;
+    const eph = decodeGpsLnavFrame(sf.buf, { prn, refWeek });
+    if (!eph) {
+      badParity++;
+      continue;
+    }
+    const prev = last.get(prn);
+    if (prev && prev.iode === eph.iode && prev.tocDate.getTime() === eph.tocDate.getTime()) {
+      continue;
+    }
+    last.set(prn, eph);
+    ephemerides.push(eph);
+  }
+  return { ephemerides, badParity };
+}
+
+// src/ubx/iono.ts
+function parseUbxIonoUtc(data) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const ionoCorrections = {};
+  let leapSeconds = null;
+  for (const f of ubxFrames(data)) {
+    if (f.msgClass !== 2 || f.msgId !== 19) continue;
+    const sf = readLnavSubframe(view, f);
+    if (!sf || sf.gnssId !== 0 || sf.id !== 4) continue;
+    const b = sf.buff;
+    if (getBitU(b, 48, 2) !== 1 || getBitU(b, 50, 6) !== 56) continue;
+    ionoCorrections["GPSA"] = [
+      getBitS(b, 56, 8) * 2 ** -30,
+      getBitS(b, 64, 8) * 2 ** -27,
+      getBitS(b, 72, 8) * 2 ** -24,
+      getBitS(b, 80, 8) * 2 ** -24
+    ];
+    ionoCorrections["GPSB"] = [
+      getBitS(b, 88, 8) * 2 ** 11,
+      getBitS(b, 96, 8) * 2 ** 14,
+      getBitS(b, 104, 8) * 2 ** 16,
+      getBitS(b, 112, 8) * 2 ** 16
+    ];
+    leapSeconds = getBitS(b, 192, 8);
+  }
+  return { ionoCorrections, leapSeconds };
+}
+
+// src/ubx/index.ts
+var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK2 = 7 * 864e5;
 var SIGNALS = {
   0: { 0: ["G", "1C"], 3: ["G", "2X"], 4: ["G", "2X"] },
@@ -4639,33 +4925,14 @@ function parseUbxRawx(data) {
   const epochs = [];
   const messageCounts = {};
   const obsCodes = {};
-  let badChecksums = 0;
-  let i = 0;
-  while (i + 8 <= data.length) {
-    if (data[i] !== 181 || data[i + 1] !== 98) {
-      i++;
-      continue;
-    }
-    const cls = data[i + 2];
-    const id = data[i + 3];
-    const len = data[i + 4] | data[i + 5] << 8;
-    const end = i + 6 + len + 2;
-    if (end > data.length) break;
-    let ckA = 0;
-    let ckB = 0;
-    for (let j = i + 2; j < i + 6 + len; j++) {
-      ckA = ckA + data[j] & 255;
-      ckB = ckB + ckA & 255;
-    }
-    if (ckA !== data[i + 6 + len] || ckB !== data[i + 7 + len]) {
-      badChecksums++;
-      i++;
-      continue;
-    }
+  const stats = { badChecksums: 0 };
+  for (const frame of ubxFrames(data, stats)) {
+    const { msgClass: cls, msgId: id } = frame;
+    const len = frame.payload.length;
     const key = `${cls.toString(16).padStart(2, "0")}-${id.toString(16).padStart(2, "0")}`;
     messageCounts[key] = (messageCounts[key] ?? 0) + 1;
     if (cls === 2 && id === 21 && len >= 16) {
-      const p = i + 6;
+      const p = frame.payloadStart;
       const rcvTow = view.getFloat64(p, true);
       const week = view.getUint16(p + 8, true);
       const leapS = view.getInt8(p + 10);
@@ -4701,14 +4968,13 @@ function parseUbxRawx(data) {
         if (!codes.includes(sig[1])) codes.push(sig[1]);
       }
       epochs.push({
-        timeMs: GPS_EPOCH_MS + week * MS_PER_WEEK2 + Math.round(rcvTow * 1e3),
+        timeMs: GPS_EPOCH_MS2 + week * MS_PER_WEEK2 + Math.round(rcvTow * 1e3),
         leapS: leapValid ? leapS : null,
         meas
       });
     }
-    i = end;
   }
-  return { epochs, messageCounts, obsCodes, badChecksums };
+  return { epochs, messageCounts, obsCodes, badChecksums: stats.badChecksums };
 }
 
 // src/sbf/frame.ts
@@ -4773,10 +5039,10 @@ function svidToPrn(svid) {
 
 // src/sbf/nav.ts
 var PI2 = Math.PI;
-var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS3 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS = Date.UTC(2006, 0, 1);
-var SEC_PER_WEEK = 7 * 86400;
-var MS_PER_WEEK3 = SEC_PER_WEEK * 1e3;
+var SEC_PER_WEEK2 = 7 * 86400;
+var MS_PER_WEEK3 = SEC_PER_WEEK2 * 1e3;
 var F4_DNU = -2e10;
 var U4_DNU = 4294967295;
 function adjustWeek(refWeek, wn, mod) {
@@ -4785,8 +5051,8 @@ function adjustWeek(refWeek, wn, mod) {
   if (offset < -(mod / 2 - 1)) offset += mod;
   return refWeek - offset;
 }
-var gpsMs = (week, sec) => GPS_EPOCH_MS2 + week * MS_PER_WEEK3 + sec * 1e3;
-var sowOf = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK;
+var gpsMs = (week, sec) => GPS_EPOCH_MS3 + week * MS_PER_WEEK3 + sec * 1e3;
+var sowOf = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK2;
 function decodeGpsQzsNav(view, b, sys) {
   const svid = view.getUint8(b + 14);
   const prn = sys === "G" ? svid : svid - 180;
@@ -4935,7 +5201,7 @@ function decodeGloNav(view, b) {
   const tocDate = new Date(toeGpsMs - leapMs);
   const tofGpsMs = gpsMs(wnc, 0) + view.getUint32(b + 8, true);
   const tofLeapMs = getGpsLeap(new Date(tofGpsMs)) * 1e3;
-  const messageFrameTime = (tofGpsMs - tofLeapMs - GPS_EPOCH_MS2) / 1e3 % SEC_PER_WEEK;
+  const messageFrameTime = (tofGpsMs - tofLeapMs - GPS_EPOCH_MS3) / 1e3 % SEC_PER_WEEK2;
   return {
     system: "R",
     prn,
@@ -5080,9 +5346,42 @@ function parseSbfAlmanac(data) {
   return { almanacs, badCrc };
 }
 
+// src/sbf/iono.ts
+var F4_DNU2 = -2e10;
+function parseSbfIonoUtc(data) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const ionoCorrections = {};
+  let leapSeconds = null;
+  const f4s = (b, n) => {
+    const out = [];
+    for (let k = 0; k < n; k++) {
+      const v = view.getFloat32(b + 4 * k, true);
+      if (v === F4_DNU2) return null;
+      out.push(v);
+    }
+    return out;
+  };
+  scanSbfFrames(data, view, (id, b, len) => {
+    if ((id === 5893 || id === 4120) && len >= 48) {
+      const alpha = f4s(b + 16, 4);
+      const beta = f4s(b + 32, 4);
+      if (!alpha || !beta) return;
+      const sys = id === 5893 ? "GPS" : "BDS";
+      ionoCorrections[`${sys}A`] = alpha;
+      ionoCorrections[`${sys}B`] = beta;
+    } else if (id === 4030 && len >= 28) {
+      const ai = f4s(b + 16, 3);
+      if (ai) ionoCorrections["GAL"] = ai;
+    } else if (id === 5894 && len >= 37) {
+      leapSeconds = view.getInt8(b + 33);
+    }
+  });
+  return { ionoCorrections, leapSeconds };
+}
+
 // src/sbf/index.ts
 var CLIGHT = 299792458;
-var GPS_EPOCH_MS3 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS4 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK4 = 7 * 864e5;
 var M3_SYS = ["G", "R", "E", "C", "S", "J", "I"];
 var M3_PR_BASE = [19e6, 19e6, 22e6, 2e7, 34e6, 34e6, 34e6];
@@ -5736,7 +6035,7 @@ function parseSbfMeas(data) {
       const tow = view.getUint32(i + 8, true);
       const wnc = view.getUint16(i + 12, true);
       if (tow !== 4294967295 && wnc !== 65535) {
-        ensureEpoch(GPS_EPOCH_MS3 + wnc * MS_PER_WEEK4 + tow);
+        ensureEpoch(GPS_EPOCH_MS4 + wnc * MS_PER_WEEK4 + tow);
         if (id === 4109) decodeMeas3Ranges(i, tow);
         else if (id === 4110) decodeMeas3CN0(i, len);
         else if (id === 4111) decodeMeas3Doppler(i, len);
@@ -5792,142 +6091,27 @@ function* oem4Frames(data, view, stats) {
   }
 }
 
-// src/navbits/index.ts
-var GPS_PI = 3.1415926535898;
-var GPS_EPOCH_MS4 = Date.UTC(1980, 0, 6);
-var SEC_PER_WEEK2 = 7 * 86400;
-function getBitU(buff, pos, len) {
-  let bits = 0;
-  for (let i = pos; i < pos + len; i++) {
-    bits = bits * 2 + (buff[i >> 3] >> 7 - (i & 7) & 1);
-  }
-  return bits;
-}
-function getBitS(buff, pos, len) {
-  const bits = getBitU(buff, pos, len);
-  if (len <= 0 || bits < 2 ** (len - 1)) return bits;
-  return bits - 2 ** len;
-}
-function decodeGpsLnavFrame(subframes, opts = {}) {
-  if (subframes.length < 90) return null;
-  const b = subframes;
-  let i = 24;
-  const tow1 = getBitU(b, i, 17) * 6;
-  i += 17 + 2;
-  const id1 = getBitU(b, i, 3);
-  i += 3 + 2;
-  const week10 = getBitU(b, i, 10);
-  i += 10;
-  i += 2;
-  i += 4;
-  const svHealth = getBitU(b, i, 6);
-  i += 6;
-  const iodc0 = getBitU(b, i, 2);
-  i += 2;
-  i += 1 + 87;
-  const tgdRaw = getBitS(b, i, 8);
-  i += 8;
-  const iodc1 = getBitU(b, i, 8);
-  i += 8;
-  const tocSec = getBitU(b, i, 16) * 16;
-  i += 16;
-  const af2 = getBitS(b, i, 8) * 2 ** -55;
-  i += 8;
-  const af1 = getBitS(b, i, 16) * 2 ** -43;
-  i += 16;
-  const af0 = getBitS(b, i, 22) * 2 ** -31;
-  i = 240 + 24;
-  i += 17 + 2;
-  const id2 = getBitU(b, i, 3);
-  i += 3 + 2;
-  const iode = getBitU(b, i, 8);
-  i += 8;
-  const crs = getBitS(b, i, 16) * 2 ** -5;
-  i += 16;
-  const deltaN = getBitS(b, i, 16) * 2 ** -43 * GPS_PI;
-  i += 16;
-  const m0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const cuc = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const e = getBitU(b, i, 32) * 2 ** -33;
-  i += 32;
-  const cus = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const sqrtA = getBitU(b, i, 32) * 2 ** -19;
-  i += 32;
-  const toes = getBitU(b, i, 16) * 16;
-  i = 480 + 24;
-  i += 17 + 2;
-  const id3 = getBitU(b, i, 3);
-  i += 3 + 2;
-  const cic = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const omega0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const cis = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const i0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const crc = getBitS(b, i, 16) * 2 ** -5;
-  i += 16;
-  const omega = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const omegaDot = getBitS(b, i, 24) * 2 ** -43 * GPS_PI;
-  i += 24;
-  const iode3 = getBitU(b, i, 8);
-  i += 8;
-  const idot = getBitS(b, i, 14) * 2 ** -43 * GPS_PI;
-  const iodc = iodc0 * 256 + iodc1;
-  if (id1 !== 1 || id2 !== 2 || id3 !== 3) return null;
-  if (iode3 !== iode || iode !== (iodc & 255)) return null;
-  const ref = opts.refWeek ?? Math.floor((Date.now() - GPS_EPOCH_MS4) / 1e3 / SEC_PER_WEEK2);
-  let week = week10 + 1024 * Math.round((ref - week10) / 1024);
-  if (toes < tow1 - 302400) week++;
-  else if (toes > tow1 + 302400) week--;
-  const prn = opts.prn ?? "G00";
-  const tocDate = new Date(
-    GPS_EPOCH_MS4 + (week * SEC_PER_WEEK2 + tocSec) * 1e3
-  );
-  return {
-    system: prn[0] === "J" ? "J" : "G",
-    prn,
-    tocDate,
-    // Same seconds-of-week convention as parseNavFile (rinex/nav.ts).
-    toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK2,
-    af0,
-    af1,
-    af2,
-    iode,
-    crs,
-    deltaN,
-    m0,
-    cuc,
-    e,
-    cus,
-    sqrtA,
-    toe: toes,
-    cic,
-    omega0,
-    cis,
-    i0,
-    crc,
-    omega,
-    omegaDot,
-    idot,
-    week,
-    svHealth,
-    tgd: tgdRaw === -128 ? 0 : tgdRaw * 2 ** -31
-    // IS-GPS-200: -128 reserved
-  };
-}
-
 // src/novatel/nav.ts
+var ID_IONUTC = 8;
 var ID_RAWEPHEM = 41;
 var ID_GLOEPHEMERIS = 723;
+var ID_GALEPHEMERIS = 1122;
+var ID_QZSSEPHEMERIS = 1336;
+var ID_BDSEPHEMERIS = 1696;
 var GPS_EPOCH_MS5 = Date.UTC(1980, 0, 6);
+var BDT_EPOCH_MS2 = Date.UTC(2006, 0, 1);
 var SEC_PER_WEEK3 = 7 * 86400;
 var SEC_PER_DAY = 86400;
+var MS_PER_WEEK5 = SEC_PER_WEEK3 * 1e3;
+var gpsMs2 = (week, sec) => GPS_EPOCH_MS5 + week * MS_PER_WEEK5 + sec * 1e3;
+var sowOf2 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK3;
+function nearWeekMs(refMs, sow) {
+  const refWeek = Math.floor((refMs - GPS_EPOCH_MS5) / MS_PER_WEEK5);
+  let ms = gpsMs2(refWeek, sow);
+  if (ms < refMs - MS_PER_WEEK5 / 2) ms += MS_PER_WEEK5;
+  else if (ms > refMs + MS_PER_WEEK5 / 2) ms -= MS_PER_WEEK5;
+  return ms;
+}
 function decodeRawEphem(data, view, p) {
   const prn = view.getUint32(p, true);
   if (prn < 1 || prn > 32) return null;
@@ -5973,23 +6157,184 @@ function decodeGloEphemeris(view, p) {
     freqNum
   };
 }
+function decodeGalEphemeris(view, p, headerMs) {
+  const prn = view.getUint32(p, true);
+  if (prn < 1 || prn > 36) return null;
+  const fnav = (view.getUint32(p + 4, true) & 1) !== 0;
+  const svhE1b = view.getUint8(p + 12) & 3;
+  const svhE5a = view.getUint8(p + 13) & 3;
+  const svhE5b = view.getUint8(p + 14) & 3;
+  const dvsE1b = view.getUint8(p + 15) & 1;
+  const dvsE5a = view.getUint8(p + 16) & 1;
+  const dvsE5b = view.getUint8(p + 17) & 1;
+  const iodNav = view.getUint32(p + 20, true);
+  const toes = view.getUint32(p + 24, true);
+  const c = fnav ? p + 148 : p + 176;
+  const tocs = view.getUint32(c, true);
+  const toeMs = nearWeekMs(headerMs, toes);
+  const tocDate = new Date(nearWeekMs(headerMs, tocs));
+  return {
+    fnav,
+    eph: {
+      system: "E",
+      prn: `E${String(prn).padStart(2, "0")}`,
+      toc: sowOf2(tocDate.getTime()),
+      tocDate,
+      af0: view.getFloat64(c + 4, true),
+      af1: view.getFloat64(c + 12, true),
+      af2: view.getFloat64(c + 20, true),
+      iode: iodNav,
+      crs: view.getFloat64(p + 92, true),
+      deltaN: view.getFloat64(p + 36, true),
+      m0: view.getFloat64(p + 44, true),
+      cuc: view.getFloat64(p + 68, true),
+      e: view.getFloat64(p + 52, true),
+      cus: view.getFloat64(p + 76, true),
+      sqrtA: view.getFloat64(p + 28, true),
+      toe: toes,
+      cic: view.getFloat64(p + 100, true),
+      omega0: view.getFloat64(p + 132, true),
+      cis: view.getFloat64(p + 108, true),
+      i0: view.getFloat64(p + 116, true),
+      crc: view.getFloat64(p + 84, true),
+      omega: view.getFloat64(p + 60, true),
+      omegaDot: view.getFloat64(p + 140, true),
+      idot: view.getFloat64(p + 124, true),
+      week: Math.floor((toeMs - GPS_EPOCH_MS5) / MS_PER_WEEK5),
+      // RINEX Galileo SVH bit layout (E1B DVS/HS in bits 0-2, E5a in
+      // 3-5, E5b in 6-8) — same packing as RTKLIB.
+      svHealth: svhE5b << 7 | dvsE5b << 6 | svhE5a << 4 | dvsE5a << 3 | svhE1b << 1 | dvsE1b,
+      tgd: view.getFloat64(p + 204, true)
+      // BGD E5a/E1 (RINEX slot)
+    }
+  };
+}
+function decodeBdsEphemeris2(view, p) {
+  const prn = view.getUint32(p, true);
+  if (prn < 1 || prn > 63) return null;
+  const week = view.getUint32(p + 4, true);
+  const health = view.getUint32(p + 16, true) & 1;
+  const tgd1 = view.getFloat64(p + 20, true);
+  const tocs = view.getUint32(p + 40, true);
+  const toes = view.getUint32(p + 72, true);
+  const tocDate = new Date(BDT_EPOCH_MS2 + week * MS_PER_WEEK5 + tocs * 1e3);
+  return {
+    system: "C",
+    prn: `C${String(prn).padStart(2, "0")}`,
+    toc: sowOf2(tocDate.getTime()),
+    tocDate,
+    af0: view.getFloat64(p + 44, true),
+    af1: view.getFloat64(p + 52, true),
+    af2: view.getFloat64(p + 60, true),
+    iode: view.getUint32(p + 68, true),
+    // AODE
+    crs: view.getFloat64(p + 172, true),
+    deltaN: view.getFloat64(p + 100, true),
+    m0: view.getFloat64(p + 108, true),
+    cuc: view.getFloat64(p + 148, true),
+    e: view.getFloat64(p + 84, true),
+    cus: view.getFloat64(p + 156, true),
+    sqrtA: view.getFloat64(p + 76, true),
+    toe: toes,
+    cic: view.getFloat64(p + 180, true),
+    omega0: view.getFloat64(p + 116, true),
+    cis: view.getFloat64(p + 188, true),
+    i0: view.getFloat64(p + 132, true),
+    crc: view.getFloat64(p + 164, true),
+    omega: view.getFloat64(p + 92, true),
+    omegaDot: view.getFloat64(p + 124, true),
+    idot: view.getFloat64(p + 140, true),
+    week,
+    // RINEX BDS week field is the BDT week
+    svHealth: health,
+    tgd: tgd1
+    // TGD1 (B1) — RINEX slot
+  };
+}
+function decodeQzssEphemeris2(view, p) {
+  const prn = view.getUint32(p, true);
+  if (prn < 193 || prn > 202) return null;
+  const health = view.getUint32(p + 12, true) & 63;
+  const iode1 = view.getUint32(p + 16, true);
+  const week = view.getUint32(p + 24, true);
+  const toes = view.getFloat64(p + 32, true);
+  const toeMs = gpsMs2(week, toes);
+  const tocs = view.getFloat64(p + 164, true);
+  const tocDate = new Date(nearWeekMs(toeMs, tocs));
+  return {
+    system: "J",
+    prn: `J${String(prn - 192).padStart(2, "0")}`,
+    toc: sowOf2(tocDate.getTime()),
+    tocDate,
+    af0: view.getFloat64(p + 180, true),
+    af1: view.getFloat64(p + 188, true),
+    af2: view.getFloat64(p + 196, true),
+    iode: iode1,
+    crs: view.getFloat64(p + 104, true),
+    deltaN: view.getFloat64(p + 48, true),
+    m0: view.getFloat64(p + 56, true),
+    cuc: view.getFloat64(p + 80, true),
+    e: view.getFloat64(p + 64, true),
+    cus: view.getFloat64(p + 88, true),
+    sqrtA: Math.sqrt(view.getFloat64(p + 40, true)),
+    toe: toes,
+    cic: view.getFloat64(p + 112, true),
+    omega0: view.getFloat64(p + 144, true),
+    cis: view.getFloat64(p + 120, true),
+    i0: view.getFloat64(p + 128, true),
+    crc: view.getFloat64(p + 96, true),
+    omega: view.getFloat64(p + 72, true),
+    omegaDot: view.getFloat64(p + 152, true),
+    idot: view.getFloat64(p + 136, true),
+    week,
+    svHealth: health,
+    tgd: view.getFloat64(p + 172, true)
+  };
+}
 function parseNovatelNav(data) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const stats = { badCrc: 0 };
   const ephemerides = [];
-  const lastGps = /* @__PURE__ */ new Map();
+  const ionoCorrections = {};
+  let leapSeconds = null;
+  const lastKep = /* @__PURE__ */ new Map();
   const lastGlo = /* @__PURE__ */ new Map();
+  const pushKepler = (key, eph) => {
+    const prev = lastKep.get(key);
+    if (prev && prev.iode === eph.iode && prev.toe === eph.toe && prev.tocDate.getTime() === eph.tocDate.getTime()) {
+      return;
+    }
+    lastKep.set(key, eph);
+    ephemerides.push(eph);
+  };
   for (const frame of oem4Frames(data, view, stats)) {
     if (!frame.binary) continue;
     if (frame.id === ID_RAWEPHEM && frame.msgLen >= 102) {
       const eph = decodeRawEphem(data, view, frame.payload);
-      if (!eph) continue;
-      const prev = lastGps.get(eph.prn);
-      if (prev && prev.iode === eph.iode && prev.tocDate.getTime() === eph.tocDate.getTime()) {
-        continue;
+      if (eph) pushKepler(eph.prn, eph);
+    } else if (frame.id === ID_GALEPHEMERIS && frame.msgLen >= 220) {
+      if (frame.week === 0) continue;
+      const headerMs = gpsMs2(frame.week, frame.towMs / 1e3);
+      const dec = decodeGalEphemeris(view, frame.payload, headerMs);
+      if (dec)
+        pushKepler(`${dec.eph.prn}:${dec.fnav ? "fnav" : "inav"}`, dec.eph);
+    } else if (frame.id === ID_BDSEPHEMERIS && frame.msgLen >= 196) {
+      const eph = decodeBdsEphemeris2(view, frame.payload);
+      if (eph) pushKepler(eph.prn, eph);
+    } else if (frame.id === ID_QZSSEPHEMERIS && frame.msgLen >= 204) {
+      const eph = decodeQzssEphemeris2(view, frame.payload);
+      if (eph) pushKepler(eph.prn, eph);
+    } else if (frame.id === ID_IONUTC && frame.msgLen >= 108) {
+      const p = frame.payload;
+      const alpha = [];
+      const beta = [];
+      for (let i = 0; i < 4; i++) {
+        alpha.push(view.getFloat64(p + i * 8, true));
+        beta.push(view.getFloat64(p + 32 + i * 8, true));
       }
-      lastGps.set(eph.prn, eph);
-      ephemerides.push(eph);
+      ionoCorrections["GPSA"] = alpha;
+      ionoCorrections["GPSB"] = beta;
+      leapSeconds = view.getInt32(p + 96, true);
     } else if (frame.id === ID_GLOEPHEMERIS && frame.msgLen >= 144) {
       const eph = decodeGloEphemeris(view, frame.payload);
       if (!eph) continue;
@@ -6001,12 +6346,12 @@ function parseNovatelNav(data) {
       ephemerides.push(eph);
     }
   }
-  return { ephemerides, badCrc: stats.badCrc };
+  return { ephemerides, ionoCorrections, leapSeconds, badCrc: stats.badCrc };
 }
 
 // src/novatel/index.ts
 var GPS_EPOCH_MS6 = Date.UTC(1980, 0, 6);
-var MS_PER_WEEK5 = 7 * 864e5;
+var MS_PER_WEEK6 = 7 * 864e5;
 var HDR = 28;
 var ID_RANGE = 43;
 var ID_RANGECMP = 140;
@@ -6175,7 +6520,7 @@ function parseNovatelRange(data) {
           if (!codes.includes(code)) codes.push(code);
         }
         epochs.push({
-          timeMs: GPS_EPOCH_MS6 + week * MS_PER_WEEK5 + towMs,
+          timeMs: GPS_EPOCH_MS6 + week * MS_PER_WEEK6 + towMs,
           meas
         });
       }
@@ -6227,7 +6572,7 @@ function parseNovatelRange(data) {
           if (!codes.includes(code)) codes.push(code);
         }
         epochs.push({
-          timeMs: GPS_EPOCH_MS6 + week * MS_PER_WEEK5 + towMs,
+          timeMs: GPS_EPOCH_MS6 + week * MS_PER_WEEK6 + towMs,
           meas
         });
       }
@@ -6587,9 +6932,9 @@ function selectBest(ephs, timeMs) {
   return best;
 }
 var GPS_EPOCH_MS7 = START_GPS_TIME.getTime();
-var MS_PER_WEEK6 = 7 * 864e5;
+var MS_PER_WEEK7 = 7 * 864e5;
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
-var GAL_EPOCH_MS = GPS_EPOCH_MS7 + 1024 * MS_PER_WEEK6;
+var GAL_EPOCH_MS = GPS_EPOCH_MS7 + 1024 * MS_PER_WEEK7;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -6641,10 +6986,10 @@ function ephInfoToEphemeris(info) {
     epochMs = GAL_EPOCH_MS;
   } else {
     epochMs = GPS_EPOCH_MS7;
-    const weeksNow = (info.lastReceived - GPS_EPOCH_MS7) / MS_PER_WEEK6;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS7) / MS_PER_WEEK7;
     week += 1024 * Math.round((weeksNow - week) / 1024);
   }
-  const tocDate = new Date(epochMs + week * MS_PER_WEEK6 + tocSec * 1e3);
+  const tocDate = new Date(epochMs + week * MS_PER_WEEK7 + tocSec * 1e3);
   return {
     system: sys,
     prn: info.prn,
@@ -9743,11 +10088,14 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   parseNovatelRange,
   parseRinexStream,
   parseSbfAlmanac,
+  parseSbfIonoUtc,
   parseSbfMeas,
   parseSbfNav,
   parseSinexBiasDcb,
   parseSourcetable,
   parseSp3,
+  parseUbxIonoUtc,
+  parseUbxNav,
   parseUbxRawx,
   phiAltBOC,
   phiBOCc,
@@ -9761,6 +10109,8 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
+  sbfCrc16,
+  scanSbfFrames,
   selectEphemeris,
   setGloFreqNumber,
   setRtcm3DebugHandler,
@@ -9770,6 +10120,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   systemCmp,
   systemName,
   transformFrame,
+  ubxFrames,
   updateStationMeta,
   updateStreamStats,
   verifyChecksum,