gnss-js 2.4.0 → 2.5.0

package/dist/index.cjs

@@ -169,9 +169,11 @@ __export(src_exports, {
   estimateNetworkFcbs: () => estimateNetworkFcbs,
   estimateWidelaneFcb: () => estimateWidelaneFcb,
   euclidean3D: () => euclidean3D,
+  extractSbfNavicSubframe: () => extractSbfNavicSubframe,
   extractWidelaneArcs: () => extractWidelaneArcs,
   feedGeoBlock: () => feedGeoBlock,
   feedGeoL5Block: () => feedGeoL5Block,
+  feedNavicBlock: () => feedNavicBlock,
   fetchSourcetable: () => fetchSourcetable,
   findSatBias: () => findSatBias,
   fixPppPosition: () => fixPppPosition,
@@ -291,6 +293,7 @@ __export(src_exports, {
   parseSbfIonoUtc: () => parseSbfIonoUtc,
   parseSbfMeas: () => parseSbfMeas,
   parseSbfNav: () => parseSbfNav,
+  parseSbfNavic: () => parseSbfNavic,
   parseSbfPvt: () => parseSbfPvt,
   parseSbfReceiverSetup: () => parseSbfReceiverSetup,
   parseSinexBiasDcb: () => parseSinexBiasDcb,
@@ -5270,6 +5273,70 @@ function decodeQzssEphemeris(payload) {
     fitInterval
   };
 }
+function decodeNavicEphemeris(payload) {
+  if (payload.length < 60) return null;
+  const r = new BitReader(payload);
+  r.skip(12);
+  const svId = r.readU(6);
+  const week = r.readU(10);
+  const af0 = r.readS(22) * 2 ** -31;
+  const af1 = r.readS(16) * 2 ** -43;
+  const af2 = r.readS(8) * 2 ** -55;
+  const ura = r.readU(4);
+  const toc = r.readU(16) * 16;
+  const tgd = r.readS(8) * 2 ** -31;
+  const deltaN = r.readS(22) * 2 ** -41 * PI;
+  const iode = r.readU(8);
+  r.skip(10);
+  const l5Health = r.readU(1);
+  const sHealth = r.readU(1);
+  const cuc = r.readS(15) * 2 ** -28;
+  const cus = r.readS(15) * 2 ** -28;
+  const cic = r.readS(15) * 2 ** -28;
+  const cis = r.readS(15) * 2 ** -28;
+  const crc = r.readS(15) * 2 ** -4;
+  const crs = r.readS(15) * 2 ** -4;
+  const idot = r.readS(14) * 2 ** -43 * PI;
+  const m0 = r.readS(32) * 2 ** -31 * PI;
+  const toe = r.readU(16) * 16;
+  const e = r.readU(32) * 2 ** -33;
+  const sqrtA = r.readU(32) * 2 ** -19;
+  const omega0 = r.readS(32) * 2 ** -31 * PI;
+  const argPerigee = r.readS(32) * 2 ** -31 * PI;
+  const omegaDot = r.readS(22) * 2 ** -41 * PI;
+  const i0 = r.readS(32) * 2 ** -31 * PI;
+  return {
+    prn: `I${String(svId).padStart(2, "0")}`,
+    constellation: "NavIC",
+    health: l5Health << 1 | sHealth,
+    lastReceived: Date.now(),
+    messageType: 1041,
+    week,
+    ura,
+    toc,
+    toe,
+    sqrtA,
+    eccentricity: e,
+    inclination: i0,
+    omega0,
+    omegaDot,
+    argPerigee,
+    meanAnomaly: m0,
+    deltaN,
+    idot,
+    crs,
+    crc,
+    cuc,
+    cus,
+    cic,
+    cis,
+    af0,
+    af1,
+    af2,
+    iode,
+    tgd
+  };
+}
 function decodeEphemeris(frame) {
   try {
     switch (frame.messageType) {
@@ -5277,6 +5344,8 @@ function decodeEphemeris(frame) {
         return decodeGpsEphemeris(frame.payload);
       case 1020:
         return decodeGlonassEphemeris(frame.payload);
+      case 1041:
+        return decodeNavicEphemeris(frame.payload);
       case 1042:
         return decodeBdsEphemeris(frame.payload);
       case 1043:
@@ -8379,6 +8448,173 @@ function parseSbfGpsNav(data) {
   return { ephemerides, badParity, messages };
 }
 
+// src/navbits/navic.ts
+var GPS_EPOCH_MS13 = Date.UTC(1980, 0, 6);
+var SEC_PER_WEEK9 = 604800;
+var DATA = 30;
+function navicSubframeCrcOk(sf) {
+  return crc24q(sf, 262) === getBitU(sf, 262, 24);
+}
+function navicSubframeId(sf) {
+  return getBitU(sf, 27, 2) + 1;
+}
+function decodeNavicSubframes(sf1, sf2, opts = {}) {
+  let i = DATA;
+  getBitU(sf1, i, 10);
+  i += 10;
+  const af0 = getBitS(sf1, i, 22) * 2 ** -31;
+  i += 22;
+  const af1 = getBitS(sf1, i, 16) * 2 ** -43;
+  i += 16;
+  const af2 = getBitS(sf1, i, 8) * 2 ** -55;
+  i += 8;
+  i += 4;
+  const tocSec = getBitU(sf1, i, 16) * 16;
+  i += 16;
+  const tgd = getBitS(sf1, i, 8) * 2 ** -31;
+  i += 8;
+  const deltaN = getBitS(sf1, i, 22) * 2 ** -41 * GPS_PI;
+  i += 22;
+  const iode = getBitU(sf1, i, 8);
+  i += 8;
+  i += 10;
+  const l5Health = getBitU(sf1, i, 1);
+  i += 1;
+  const sHealth = getBitU(sf1, i, 1);
+  i += 1;
+  const cuc = getBitS(sf1, i, 15) * 2 ** -28;
+  i += 15;
+  const cus = getBitS(sf1, i, 15) * 2 ** -28;
+  i += 15;
+  const cic = getBitS(sf1, i, 15) * 2 ** -28;
+  i += 15;
+  const cis = getBitS(sf1, i, 15) * 2 ** -28;
+  i += 15;
+  const crc = getBitS(sf1, i, 15) * 2 ** -4;
+  i += 15;
+  const crs = getBitS(sf1, i, 15) * 2 ** -4;
+  i += 15;
+  const idot = getBitS(sf1, i, 14) * 2 ** -43 * GPS_PI;
+  let j = DATA;
+  const m0 = getBitS(sf2, j, 32) * 2 ** -31 * GPS_PI;
+  j += 32;
+  const toeSec = getBitU(sf2, j, 16) * 16;
+  j += 16;
+  const e = getBitU(sf2, j, 32) * 2 ** -33;
+  j += 32;
+  const sqrtA = getBitU(sf2, j, 32) * 2 ** -19;
+  j += 32;
+  const omega0 = getBitS(sf2, j, 32) * 2 ** -31 * GPS_PI;
+  j += 32;
+  const omega = getBitS(sf2, j, 32) * 2 ** -31 * GPS_PI;
+  j += 32;
+  const omegaDot = getBitS(sf2, j, 22) * 2 ** -41 * GPS_PI;
+  j += 22;
+  const i0 = getBitS(sf2, j, 32) * 2 ** -31 * GPS_PI;
+  if (sqrtA <= 0 || e >= 1) return null;
+  const week = opts.refWeek ?? Math.floor((Date.now() - GPS_EPOCH_MS13) / 1e3 / SEC_PER_WEEK9);
+  const prn = opts.prn ?? "I00";
+  const tocDate = new Date(
+    GPS_EPOCH_MS13 + (week * SEC_PER_WEEK9 + tocSec) * 1e3
+  );
+  return {
+    system: "I",
+    prn,
+    tocDate,
+    toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK9,
+    af0,
+    af1,
+    af2,
+    iode,
+    crs,
+    deltaN,
+    m0,
+    cuc,
+    e,
+    cus,
+    sqrtA,
+    toe: toeSec,
+    cic,
+    omega0,
+    cis,
+    i0,
+    crc,
+    omega,
+    omegaDot,
+    idot,
+    week,
+    svHealth: l5Health << 1 | sHealth,
+    tgd
+  };
+}
+var NavicAssembler = class {
+  sf1 = /* @__PURE__ */ new Map();
+  last = /* @__PURE__ */ new Map();
+  /**
+   * Push one CRC-valid subframe. Subframe 1 opens a master frame (subframe 2
+   * immediately follows it in the NavIC frame order), so an ephemeris is only
+   * emitted when a subframe 2 pairs with the subframe 1 that opened the *same*
+   * frame — never a stale subframe 1 left over across CRC drops.
+   */
+  push(prn, subframe, refWeek) {
+    if (getBitU(subframe, 25, 1)) return null;
+    const id = navicSubframeId(subframe);
+    if (id === 1) {
+      this.sf1.set(prn, subframe);
+      return null;
+    }
+    if (id !== 2) return null;
+    const a = this.sf1.get(prn);
+    if (!a) return null;
+    this.sf1.delete(prn);
+    const eph = decodeNavicSubframes(a, subframe, { prn, refWeek });
+    if (!eph) return null;
+    const key = `${eph.toc}|${eph.iode}`;
+    if (this.last.get(prn) === key) return null;
+    this.last.set(prn, key);
+    return eph;
+  }
+};
+
+// src/sbf/rawnav-navic.ts
+var WNC_DNU2 = 65535;
+function extractSbfNavicSubframe(view, b) {
+  const sf = new Uint8Array(40);
+  for (let k = 0; k < 10; k++) {
+    const w = view.getUint32(b + 20 + 4 * k, true);
+    sf[4 * k] = w >>> 24;
+    sf[4 * k + 1] = w >>> 16 & 255;
+    sf[4 * k + 2] = w >>> 8 & 255;
+    sf[4 * k + 3] = w & 255;
+  }
+  return sf;
+}
+function feedNavicBlock(view, b, assembler) {
+  const prn = svidToPrn(view.getUint8(b + 14));
+  if (!prn || prn[0] !== "I") return {};
+  const wnc = view.getUint16(b + 12, true);
+  if (wnc === WNC_DNU2) return {};
+  const sf = extractSbfNavicSubframe(view, b);
+  if (!navicSubframeCrcOk(sf)) return { badCrc: true };
+  const eph = assembler.push(prn, sf, wnc);
+  return eph ? { eph } : {};
+}
+function parseSbfNavic(data) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const ephemerides = [];
+  const assembler = new NavicAssembler();
+  let badCrc = 0;
+  let messages = 0;
+  scanSbfFrames(data, view, (id, b, len) => {
+    if (id !== 4093 || len < 60) return;
+    messages++;
+    const r = feedNavicBlock(view, b, assembler);
+    if (r.eph) ephemerides.push(r.eph);
+    else if (r.badCrc) badCrc++;
+  });
+  return { ephemerides, badCrc, messages };
+}
+
 // src/sbf/rawnav.ts
 var CNAV_SIGNAL_OF_BLOCK = {
   4018: "L2C",
@@ -8477,10 +8713,10 @@ function isFnavDummy(page) {
 }
 
 // src/sbf/rawnav-bds.ts
-var GPS_EPOCH_MS13 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS14 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK6 = 7 * 86400 * 1e3;
 var TOW_DNU = 4294967295;
-var WNC_DNU2 = 65535;
+var WNC_DNU3 = 65535;
 function feedBdsBlock(view, b, assembler) {
   const prn = svidToPrn(view.getUint8(b + 14));
   if (!prn || prn[0] !== "C") return {};
@@ -8520,7 +8756,7 @@ function feedGloBlock(view, b, assembler) {
   if (!prn || prn[0] !== "R") return {};
   const towMs = view.getUint32(b + 8, true);
   const wnc = view.getUint16(b + 12, true);
-  if (towMs === TOW_DNU || wnc === WNC_DNU2) return {};
+  if (towMs === TOW_DNU || wnc === WNC_DNU3) return {};
   const str2 = new Uint8Array(11);
   for (let k = 0; k < 3; k++) {
     const w = view.getUint32(b + 20 + 4 * k, true);
@@ -8538,7 +8774,7 @@ function feedGloBlock(view, b, assembler) {
   const eph = assembler.push(
     prn,
     str2,
-    new Date(GPS_EPOCH_MS13 + wnc * MS_PER_WEEK6 + towMs),
+    new Date(GPS_EPOCH_MS14 + wnc * MS_PER_WEEK6 + towMs),
     view.getUint8(b + 18) - 8
     // FreqNr, offset by 8
   );
@@ -8713,6 +8949,7 @@ function decodeSbfNavigation(data, opts = {}) {
     bdsRaw: 0,
     sbasRaw: 0,
     sbasL5Raw: 0,
+    navicRaw: 0,
     badFrames: 0
   };
   const dfmcByType = {};
@@ -8722,6 +8959,7 @@ function decodeSbfNavigation(data, opts = {}) {
   const gloAsm = new GloStringAssembler();
   const bdsAsm = new BdsAssembler();
   const lnavAsm = new GpsLnavAssembler();
+  const navicAsm = new NavicAssembler();
   const addEph = (e) => {
     if (!e) return;
     const key = ephKey(e);
@@ -8782,6 +9020,14 @@ function decodeSbfNavigation(data, opts = {}) {
         else if (r.badParity) counts.badFrames++;
         return;
       }
+      case 4093: {
+        if (len < 60) return;
+        counts.navicRaw++;
+        const r = feedNavicBlock(view, b, navicAsm);
+        if (r.eph) addEph(r.eph);
+        else if (r.badCrc) counts.badFrames++;
+        return;
+      }
       case 4018:
       // GPSRawL2C
       case 4019: {
@@ -9402,7 +9648,7 @@ function parseSbfHas(data) {
 
 // src/sbf/index.ts
 var CLIGHT = 299792458;
-var GPS_EPOCH_MS14 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS15 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK7 = 7 * 864e5;
 var M3_SYS = ["G", "R", "E", "C", "S", "J", "I"];
 var M3_PR_BASE = [19e6, 19e6, 22e6, 2e7, 34e6, 34e6, 34e6];
@@ -10056,7 +10302,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_MS14 + wnc * MS_PER_WEEK7 + tow);
+        ensureEpoch(GPS_EPOCH_MS15 + wnc * MS_PER_WEEK7 + tow);
         if (id === 4109) decodeMeas3Ranges(i, tow);
         else if (id === 4110) decodeMeas3CN0(i, len);
         else if (id === 4111) decodeMeas3Doppler(i, len);
@@ -10123,15 +10369,15 @@ var ID_GPSEPHEM = 7;
 var ID_BDSEPHEMERIS = 1696;
 var ID_RAWWAASFRAME = 287;
 var ID_RAWSBASFRAME = 973;
-var GPS_EPOCH_MS15 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS16 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS3 = Date.UTC(2006, 0, 1);
-var SEC_PER_WEEK9 = 7 * 86400;
+var SEC_PER_WEEK10 = 7 * 86400;
 var SEC_PER_DAY2 = 86400;
-var MS_PER_WEEK8 = SEC_PER_WEEK9 * 1e3;
-var gpsMs2 = (week, sec) => GPS_EPOCH_MS15 + week * MS_PER_WEEK8 + sec * 1e3;
-var sowOf3 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK9;
+var MS_PER_WEEK8 = SEC_PER_WEEK10 * 1e3;
+var gpsMs2 = (week, sec) => GPS_EPOCH_MS16 + week * MS_PER_WEEK8 + sec * 1e3;
+var sowOf3 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK10;
 function nearWeekMs(refMs, sow) {
-  const refWeek = Math.floor((refMs - GPS_EPOCH_MS15) / MS_PER_WEEK8);
+  const refWeek = Math.floor((refMs - GPS_EPOCH_MS16) / MS_PER_WEEK8);
   let ms = gpsMs2(refWeek, sow);
   if (ms < refMs - MS_PER_WEEK8 / 2) ms += MS_PER_WEEK8;
   else if (ms > refMs + MS_PER_WEEK8 / 2) ms -= MS_PER_WEEK8;
@@ -10158,9 +10404,9 @@ function decodeGloEphemeris(view, p) {
   tof += Math.floor(tow / SEC_PER_DAY2) * SEC_PER_DAY2;
   if (tof < tow - 43200) tof += SEC_PER_DAY2;
   else if (tof > tow + 43200) tof -= SEC_PER_DAY2;
-  const toeGpsMs = GPS_EPOCH_MS15 + (week * SEC_PER_WEEK9 + tow) * 1e3;
-  const tofGpsMs = GPS_EPOCH_MS15 + (week * SEC_PER_WEEK9 + tof) * 1e3;
-  const tofUtcSec = (getUtcDate(new Date(tofGpsMs)).getTime() - GPS_EPOCH_MS15) / 1e3;
+  const toeGpsMs = GPS_EPOCH_MS16 + (week * SEC_PER_WEEK10 + tow) * 1e3;
+  const tofGpsMs = GPS_EPOCH_MS16 + (week * SEC_PER_WEEK10 + tof) * 1e3;
+  const tofUtcSec = (getUtcDate(new Date(tofGpsMs)).getTime() - GPS_EPOCH_MS16) / 1e3;
   return {
     system: "R",
     prn: `R${String(slot).padStart(2, "0")}`,
@@ -10168,7 +10414,7 @@ function decodeGloEphemeris(view, p) {
     tauN: -view.getFloat64(p + 100, true),
     // RINEX stores −τ_n
     gammaN: view.getFloat64(p + 116, true),
-    messageFrameTime: (tofUtcSec % SEC_PER_WEEK9 + SEC_PER_WEEK9) % SEC_PER_WEEK9,
+    messageFrameTime: (tofUtcSec % SEC_PER_WEEK10 + SEC_PER_WEEK10) % SEC_PER_WEEK10,
     x: view.getFloat64(p + 28, true) / 1e3,
     y: view.getFloat64(p + 36, true) / 1e3,
     z: view.getFloat64(p + 44, true) / 1e3,
@@ -10225,7 +10471,7 @@ function decodeGalEphemeris(view, p, headerMs) {
       omega: view.getFloat64(p + 60, true),
       omegaDot: view.getFloat64(p + 140, true),
       idot: view.getFloat64(p + 124, true),
-      week: Math.floor((toeMs - GPS_EPOCH_MS15) / MS_PER_WEEK8),
+      week: Math.floor((toeMs - GPS_EPOCH_MS16) / MS_PER_WEEK8),
       // 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,
@@ -10412,7 +10658,7 @@ function parseNovatelNav(data, opts = {}) {
 }
 
 // src/novatel/pvt.ts
-var GPS_EPOCH_MS16 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS17 = Date.UTC(1980, 0, 6);
 var BESTPOS = 42;
 function posTypeToMode(t) {
   if (t === 0) return "no-pvt";
@@ -10444,7 +10690,7 @@ function parseNovatelPvt(data) {
     const nrSV = view.getUint8(p + 65);
     const hasFix = posType !== 0 && Number.isFinite(lat) && Number.isFinite(lon);
     records.push({
-      timeMs: GPS_EPOCH_MS16 + f.week * 6048e5 + f.towMs,
+      timeMs: GPS_EPOCH_MS17 + f.week * 6048e5 + f.towMs,
       week: f.week,
       tow: f.towMs / 1e3,
       mode: posTypeToMode(posType),
@@ -10461,7 +10707,7 @@ function parseNovatelPvt(data) {
 }
 
 // src/novatel/index.ts
-var GPS_EPOCH_MS17 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS18 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK9 = 7 * 864e5;
 var HDR = 28;
 var ID_RANGE = 43;
@@ -10631,7 +10877,7 @@ function parseNovatelRange(data) {
           if (!codes.includes(code)) codes.push(code);
         }
         epochs.push({
-          timeMs: GPS_EPOCH_MS17 + week * MS_PER_WEEK9 + towMs,
+          timeMs: GPS_EPOCH_MS18 + week * MS_PER_WEEK9 + towMs,
           meas
         });
       }
@@ -10683,7 +10929,7 @@ function parseNovatelRange(data) {
           if (!codes.includes(code)) codes.push(code);
         }
         epochs.push({
-          timeMs: GPS_EPOCH_MS17 + week * MS_PER_WEEK9 + towMs,
+          timeMs: GPS_EPOCH_MS18 + week * MS_PER_WEEK9 + towMs,
           meas
         });
       }
@@ -10734,15 +10980,15 @@ function* trimbleFrames(data, stats) {
 }
 
 // src/trimble/nav.ts
-var GPS_EPOCH_MS18 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS19 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS4 = Date.UTC(2006, 0, 1);
-var SEC_PER_WEEK10 = 7 * 86400;
-var MS_PER_WEEK10 = SEC_PER_WEEK10 * 1e3;
+var SEC_PER_WEEK11 = 7 * 86400;
+var MS_PER_WEEK10 = SEC_PER_WEEK11 * 1e3;
 var GPS_PI2 = 3.1415926535898;
 var BDS_WEEK_OFFSET = 1356;
 var BDS_LEAP_SEC = 14;
-var gpsMs3 = (week, sec) => GPS_EPOCH_MS18 + week * MS_PER_WEEK10 + sec * 1e3;
-var sowOf4 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK10;
+var gpsMs3 = (week, sec) => GPS_EPOCH_MS19 + week * MS_PER_WEEK10 + sec * 1e3;
+var sowOf4 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK11;
 var SUB_GPS_EPHEMERIS = 1;
 var SUB_ION_UTC = 3;
 var SUB_GLO_EPHEMERIS = 9;
@@ -10894,7 +11140,7 @@ function decodeGloEphemeris2(view, base) {
     // RINEX stores −τn as the clock bias; the record carries τn (ICD sign).
     tauN: -view.getFloat64(base + 113, false),
     gammaN: view.getFloat64(base + 121, false),
-    messageFrameTime: (tocMs - GPS_EPOCH_MS18) / 1e3 % SEC_PER_WEEK10,
+    messageFrameTime: (tocMs - GPS_EPOCH_MS19) / 1e3 % SEC_PER_WEEK11,
     // State vector: metres / m·s⁻¹ / m·s⁻² → km / km·s⁻¹ / km·s⁻².
     x: view.getFloat64(base + 33, false) / 1e3,
     xDot: view.getFloat64(base + 41, false) / 1e3,
@@ -10965,7 +11211,7 @@ function parseTrimbleNav(data) {
 }
 
 // src/trimble/index.ts
-var GPS_EPOCH_MS19 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS20 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK11 = 7 * 864e5;
 var REC_TYPE_27 = 6;
 var SYS = {
@@ -11163,7 +11409,7 @@ function decodeEpoch(epoch, epochs, obsCodes) {
     o += c.length;
   }
   const dv = new DataView(M.buffer, M.byteOffset, M.byteLength);
-  const timeMs = GPS_EPOCH_MS19 + epoch.week * MS_PER_WEEK11 + epoch.towMs;
+  const timeMs = GPS_EPOCH_MS20 + epoch.week * MS_PER_WEEK11 + epoch.towMs;
   const meas = [];
   let q = 0;
   for (let sv = 0; sv < epoch.nsv && q + 8 <= M.length; sv++) {
@@ -11344,14 +11590,14 @@ function* binexRecords(data, view, stats) {
 }
 
 // src/binex/nav.ts
-var GPS_EPOCH_MS20 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS21 = Date.UTC(1980, 0, 6);
 var GLO_DAY0_MS = Date.UTC(1980, 0, 1);
 var BDT_EPOCH_MS5 = Date.UTC(2006, 0, 1);
-var SEC_PER_WEEK11 = 7 * 86400;
-var MS_PER_WEEK12 = SEC_PER_WEEK11 * 1e3;
+var SEC_PER_WEEK12 = 7 * 86400;
+var MS_PER_WEEK12 = SEC_PER_WEEK12 * 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 gpsMs4 = (week, sec) => GPS_EPOCH_MS21 + week * MS_PER_WEEK12 + sec * 1e3;
+var sowOf5 = (dateMs) => (dateMs / 1e3 % SEC_PER_WEEK12 + SEC_PER_WEEK12) % SEC_PER_WEEK12;
 var two9 = (n) => String(n).padStart(2, "0");
 var Reader = class {
   constructor(view, start, le) {
@@ -11638,7 +11884,7 @@ function decodeBinexEph(view, body, len, le) {
 }
 
 // src/binex/index.ts
-var GPS_EPOCH_MS21 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS22 = Date.UTC(1980, 0, 6);
 var CLIGHT2 = 299792458;
 var OBS_CODES = {
   0: [
@@ -11896,7 +12142,7 @@ function decodeObsRecord(data, view, body, len, obsCodes) {
   p += 4;
   const msec = view.getUint16(p);
   p += 2;
-  const timeMs = GPS_EPOCH_MS21 + min * 6e4 + msec;
+  const timeMs = GPS_EPOCH_MS22 + min * 6e4 + msec;
   const flag = data[p++];
   const nsat = (flag & 63) + 1;
   if (flag & 128) p += 3;
@@ -11982,7 +12228,7 @@ function parseBinex(data) {
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
-var GPS_EPOCH_MS22 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS23 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS6 = Date.UTC(2006, 0, 1);
 var WEEK_MS = 6048e5;
 var BDT_MINUS_GPST_MS = 14e3;
@@ -11991,7 +12237,7 @@ function almanacEpochMs(alm) {
   if (alm.system === "C") {
     return BDT_EPOCH_MS6 + alm.weekAlm * WEEK_MS + alm.toaSec * 1e3 + BDT_MINUS_GPST_MS;
   }
-  return GPS_EPOCH_MS22 + alm.weekAlm * WEEK_MS + alm.toaSec * 1e3;
+  return GPS_EPOCH_MS23 + alm.weekAlm * WEEK_MS + alm.toaSec * 1e3;
 }
 function almanacSatPosition(alm, timeMs) {
   const GM = alm.system === "E" ? GM_GAL : alm.system === "C" ? GM_BDS : GM_GPS;
@@ -12500,10 +12746,10 @@ function selectBest(ephs, timeMs) {
   }
   return best;
 }
-var GPS_EPOCH_MS23 = START_GPS_TIME.getTime();
+var GPS_EPOCH_MS24 = START_GPS_TIME.getTime();
 var MS_PER_WEEK13 = 7 * 864e5;
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
-var GAL_EPOCH_MS = GPS_EPOCH_MS23 + 1024 * MS_PER_WEEK13;
+var GAL_EPOCH_MS = GPS_EPOCH_MS24 + 1024 * MS_PER_WEEK13;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -12554,8 +12800,8 @@ function ephInfoToEphemeris(info) {
   } else if (sys === "E") {
     epochMs = GAL_EPOCH_MS;
   } else {
-    epochMs = GPS_EPOCH_MS23;
-    const weeksNow = (info.lastReceived - GPS_EPOCH_MS23) / MS_PER_WEEK13;
+    epochMs = GPS_EPOCH_MS24;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS24) / MS_PER_WEEK13;
     week += 1024 * Math.round((weeksNow - week) / 1024);
   }
   const tocDate = new Date(epochMs + week * MS_PER_WEEK13 + tocSec * 1e3);
@@ -14396,7 +14642,7 @@ var P2_31 = 2 ** -31;
 var P2_39 = 2 ** -39;
 var D2R = Math.PI / 180;
 var R2D = 180 / Math.PI;
-var SEC_PER_WEEK12 = 604800;
+var SEC_PER_WEEK13 = 604800;
 var MAXSBSAGEF = 30;
 var MAXSBSAGEL = 360;
 var MAXBAND = 10;
@@ -14544,7 +14790,7 @@ var igpband2 = [
   [{ coord: 60, arr: x5, bits: 1, bite: 72 }, { coord: 65, arr: x6, bits: 73, bite: 108 }, { coord: 70, arr: x6, bits: 109, bite: 144 }, { coord: 75, arr: x6, bits: 145, bite: 180 }, { coord: 85, arr: x7, bits: 181, bite: 192 }],
   [{ coord: -60, arr: x5, bits: 1, bite: 72 }, { coord: -65, arr: x6, bits: 73, bite: 108 }, { coord: -70, arr: x6, bits: 109, bite: 144 }, { coord: -75, arr: x6, bits: 145, bite: 180 }, { coord: -85, arr: x8, bits: 181, bite: 192 }]
 ];
-var gpsSeconds = (week, tow) => week * SEC_PER_WEEK12 + tow;
+var gpsSeconds = (week, tow) => week * SEC_PER_WEEK13 + tow;
 function maskPrn(i) {
   if (i >= 1 && i <= 37) return `G${two11(i)}`;
   if (i >= 38 && i <= 61) return `R${two11(i - 37)}`;
@@ -19197,9 +19443,11 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   estimateNetworkFcbs,
   estimateWidelaneFcb,
   euclidean3D,
+  extractSbfNavicSubframe,
   extractWidelaneArcs,
   feedGeoBlock,
   feedGeoL5Block,
+  feedNavicBlock,
   fetchSourcetable,
   findSatBias,
   fixPppPosition,
@@ -19319,6 +19567,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   parseSbfIonoUtc,
   parseSbfMeas,
   parseSbfNav,
+  parseSbfNavic,
   parseSbfPvt,
   parseSbfReceiverSetup,
   parseSinexBiasDcb,

package/dist/index.d.cts

@@ -13,7 +13,7 @@ export { RTCM2_MESSAGE_NAMES, Rtcm2Correction, Rtcm2Decoder, Rtcm2Frame, Rtcm2Ob
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.cjs';
 export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, ObsEpoch, ObsSatObs, ObsSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxPvtResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxPvt, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.cjs';
-export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, feedGeoBlock, feedGeoL5Block, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.cjs';
+export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavicResult, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, extractSbfNavicSubframe, feedGeoBlock, feedGeoL5Block, feedNavicBlock, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfNavic, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.cjs';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.cjs';
 export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.cjs';
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';

package/dist/index.d.ts

@@ -13,7 +13,7 @@ export { RTCM2_MESSAGE_NAMES, Rtcm2Correction, Rtcm2Decoder, Rtcm2Frame, Rtcm2Ob
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.js';
 export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, ObsEpoch, ObsSatObs, ObsSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxPvtResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxPvt, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.js';
-export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, feedGeoBlock, feedGeoL5Block, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.js';
+export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavicResult, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, extractSbfNavicSubframe, feedGeoBlock, feedGeoL5Block, feedNavicBlock, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfNavic, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.js';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.js';
 export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.js';
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';

package/dist/index.js

@@ -77,7 +77,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-XPQ4BU6W.js";
+} from "./chunk-7GNRNKBR.js";
 import {
   parseUbxCnav,
   parseUbxIonoUtc,
@@ -92,8 +92,10 @@ import {
   SBAS_L5_PREAMBLE_UW,
   crc16,
   decodeSbfNavigation,
+  extractSbfNavicSubframe,
   feedGeoBlock,
   feedGeoL5Block,
+  feedNavicBlock,
   isSbasL1Preamble,
   isSbasL5Correction,
   parseSbfAlmanac,
@@ -108,12 +110,13 @@ import {
   parseSbfIonoUtc,
   parseSbfMeas,
   parseSbfNav,
+  parseSbfNavic,
   parseSbfPvt,
   parseSbfReceiverSetup,
   sbasL5MessageName,
   sbasL5MessageType,
   scanSbfFrames
-} from "./chunk-WW3DKCDN.js";
+} from "./chunk-74IVYOQV.js";
 import "./chunk-3KSDL3CO.js";
 import {
   OEM4_HLEN,
@@ -542,9 +545,11 @@ export {
   estimateNetworkFcbs,
   estimateWidelaneFcb,
   euclidean3D,
+  extractSbfNavicSubframe,
   extractWidelaneArcs,
   feedGeoBlock,
   feedGeoL5Block,
+  feedNavicBlock,
   fetchSourcetable,
   findSatBias,
   fixPppPosition,
@@ -664,6 +669,7 @@ export {
   parseSbfIonoUtc,
   parseSbfMeas,
   parseSbfNav,
+  parseSbfNavic,
   parseSbfPvt,
   parseSbfReceiverSetup,
   parseSinexBiasDcb,