gnss-js 2.11.0 → 2.12.0

package/README.md

@@ -17,6 +17,7 @@ Everything is exposed through subpath exports (`gnss-js/rinex`, `gnss-js/positio
 - **Formats in** — RINEX 2/3/4 obs + nav (plain, Hatanaka `.crx`, `.gz`), SP3, RINEX CLK, IONEX, SINEX bias; RTCM3 (MSM + legacy obs, ephemeris, SSR/IGS-SSR, station); RTCM 2.x; and the receiver binaries SBF, UBX, NovAtel OEM4, Trimble, BINEX.
 - **Formats out** — RINEX 2/3/4 observation + navigation writers and Hatanaka (CRINEX) compression.
 - **Navigation** — Keplerian/GLONASS ephemeris handling plus raw-frame decoding of GPS LNAV/CNAV, Galileo I/NAV + F/NAV, BeiDou D1/D2 and B-CNAV1/2/3 (B1C/B2a/B2b), GLONASS strings, and Galileo HAS (E6 C/NAV) corrections.
+- **Authentication** — Galileo OSNMA: the full chain of trust from the Merkle root through ECDSA-signed TESLA keys to per-satellite MAC tags, flagging navigation data Galileo did not sign.
 - **Orbits** — satellite position/velocity, Doppler, DOP, azimuth/elevation, pass prediction and almanac propagation.
 - **Positioning** — SPP, static and kinematic float PPP, PPP-AR (LAMBDA + wide-/narrow-lane FCB estimation), post-processed RTK, DGNSS, and SBAS/DFMC augmentation.
 - **Quality** — multipath, cycle slips, completeness and slant-TEC ionosphere analysis, with SINEX-DCB calibration.
@@ -99,6 +100,20 @@ const nav = parseSbfNav(bytes); // GPS/GAL/BDS/GLO ephemeris from raw blocks
 const has = parseSbfHas(bytes); // Galileo HAS (E6 C/NAV) SSR corrections
 ```
 
+### Galileo OSNMA (navigation message authentication)
+
+`OsnmaReceiver` takes Galileo E1-B I/NAV pages and verifies them the way the OSNMA SIS ICD v1.1 prescribes: public key against the Merkle root (built in, or broadcast in DSM-PKR), the DSM-KROOT's ECDSA signature, every TESLA chain key back to the root key, and the MAC tags over each satellite's ephemeris/clock (ADKD 0/12) and timing parameters (ADKD 4). A tag that fails means the data differ from what Galileo signed. Validated against the official ICD test vectors and live receivers.
+
+```ts
+import { OsnmaReceiver, verifySbfOsnma } from 'gnss-js/osnma';
+
+const s = verifySbfOsnma(sbfBytes); // whole capture
+s.chains; // verified TESLA chains; s.tags → { ok, failed, ... }
+s.satellites; // per PRN: ced.authenticatedIodNav, timing.lastAuthenticated, failures
+
+const rx = new OsnmaReceiver(); // live: rx.pushPage(svid, wn, tow, page); rx.status()
+```
+
 ### RTCM (RTCM3, RTCM2, SSR)
 
 ```ts

package/dist/analysis.cjs

@@ -1113,20 +1113,14 @@ function applyIonoDcb(iono, satDcb) {
 // src/rinex/crx.ts
 function crxRepair(old, diff) {
   const chars = old.split("");
-  for (let i = 0; i < diff.length; i++) {
+  const n = Math.min(old.length, diff.length);
+  for (let i = 0; i < n; i++) {
     const dc = diff[i];
-    if (dc === " ") {
-    } else if (dc === "&") {
-      chars[i] = " ";
-    } else {
-      chars[i] = dc;
-    }
-  }
-  if (diff.length > chars.length) {
-    for (let i = chars.length; i < diff.length; i++) {
-      chars[i] = diff[i] === "&" ? " " : diff[i];
-    }
+    if (dc === " ") continue;
+    chars[i] = dc === "&" ? " " : dc;
   }
+  for (let i = old.length; i < diff.length; i++)
+    chars[i] = diff[i] === "&" ? " " : diff[i];
   return chars.join("");
 }
 function parseCrxDataLine(line, ntype) {
@@ -1371,6 +1365,15 @@ function readObsValue(line, index) {
   const v = parseFloat(raw);
   return isNaN(v) ? null : v;
 }
+function normalizePrn(id) {
+  let sys = id[0];
+  let num = id.substring(1);
+  if (sys === " " || /\d/.test(sys)) {
+    num = id;
+    sys = "G";
+  }
+  return sys + num.replace(/\s/g, "").padStart(2, "0");
+}
 function parseEpochLine3(line) {
   if (line[0] !== ">") return null;
   const yr = parseInt(line.substring(2, 6));
@@ -1520,7 +1523,7 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
   function processSatLineV3(line) {
     if (!header || line.length < 3) return;
     const sys = line[0];
-    const prn = line.substring(0, 3);
+    const prn = normalizePrn(line.substring(0, 3));
     trackSat(sys, prn);
     const obsLine = line.substring(3);
     const codes = header.obsTypes[sys] ?? [];
@@ -1545,10 +1548,8 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
   let v2ObsAccum = [];
   function processSatLineV2(line) {
     if (!header || v2SatIndex >= v2SatIds.length) return;
-    const satId = v2SatIds[v2SatIndex];
-    let sys = satId[0];
-    if (sys === " " || /\d/.test(sys)) sys = "G";
-    const prn = sys + satId.substring(1).padStart(2, "0");
+    const prn = normalizePrn(v2SatIds[v2SatIndex]);
+    const sys = prn[0];
     if (v2CurrentSatLine === 0) {
       trackSat(sys, prn);
       if (onSatObs) v2ObsAccum = [];
@@ -1639,6 +1640,7 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
         } else {
           epochLine = crxRepair(crxPrevEpochLine, line);
         }
+        epochLine = epochLine.replace(/ +$/, "");
         if (isV3) {
           if (epochLine[0] !== ">") return;
           const info = parseEpochLine3(epochLine);
@@ -1656,10 +1658,11 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
             return;
           }
           resetEpochAccum(info);
-          const satPart = epochLine.substring(41).replace(/\s+$/, "");
           crxSatList = [];
-          for (let i = 0; i + 3 <= satPart.length; i += 3) {
-            crxSatList.push(satPart.substring(i, i + 3));
+          for (let k = 0; k < info.numSats; k++) {
+            crxSatList.push(
+              normalizePrn(epochLine.substring(41 + 3 * k, 44 + 3 * k))
+            );
           }
           crxSatIndex = 0;
           crxPrevEpochLine = epochLine;
@@ -1680,11 +1683,7 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
             return;
           }
           resetEpochAccum(info);
-          crxSatList = (info.satIds ?? []).map((id) => {
-            let sys = id[0];
-            if (sys === " " || /\d/.test(sys)) sys = "G";
-            return sys + id.substring(1).padStart(2, "0");
-          });
+          crxSatList = (info.satIds ?? []).slice(0, info.numSats).map(normalizePrn);
           crxSatIndex = 0;
           crxPrevEpochLine = epochLine;
           crxPhase = "clock";

package/dist/analysis.js

@@ -8,8 +8,8 @@ import {
   computeIonoRate,
   detrendIonoArcs,
   parseSinexBiasDcb
-} from "./chunk-KUSUN7DP.js";
-import "./chunk-3U5AX7PY.js";
+} from "./chunk-MFCPYUEI.js";
+import "./chunk-FEABDV5P.js";
 import "./chunk-FIEWO4J4.js";
 export {
   CompletenessAccumulator,

package/dist/{clk-DGVBGRzA.d.cts → bias-sinex-OGMw9xEZ.d.cts}

@@ -93,4 +93,69 @@ declare function parseClk(text: string): ClkFile;
  */
 declare function clkBias(clk: ClkFile, prn: string, tMs: number): number | null;
 
-export { type ClkFile as C, type IonexGrid as I, type Sp3File as S, type Sp3Sample as a, parseIonex as b, clkBias as c, parseSp3 as d, parseClk as p, sp3Position as s };
+/**
+ * Bias-SINEX (`.BIA`) parser — GNSS code/phase biases published by the IGS
+ * analysis centres. Two record kinds matter here:
+ *
+ *   - **DSB** differential signal bias (OBS1−OBS2), the classic inter-code
+ *     bias used for ionosphere/DCB corrections in SPP.
+ *   - **OSB** observable-specific (absolute) bias for a single OBS, the
+ *     ingredient PPP-AR needs: the *phase* OSBs are the satellite
+ *     fractional-cycle biases that make the ambiguities integer-recoverable.
+ *
+ * Values keep their file unit (`ns` for code/most phase biases, `cyc` for
+ * cycle-based phase biases); `biasMetres` converts ns → m via c.
+ */
+interface BiasRecord {
+    /** 'DSB' | 'OSB' | 'ISB' | … */
+    type: string;
+    /** Satellite PRN (e.g. 'G01', 'E11'); '' for a pure station bias. */
+    prn: string;
+    /** Station (9-char), '' for satellite biases. */
+    station: string;
+    /** First observable code (e.g. 'C1W', 'L1C'). */
+    obs1: string;
+    /** Second observable (DSB only; '' for OSB). */
+    obs2: string;
+    /** Validity start / end (GPS-scale ms); ±Infinity when open. */
+    startMs: number;
+    endMs: number;
+    /** Estimated value in the file unit. */
+    value: number;
+    /** File unit ('ns' | 'cyc' | …). */
+    unit: string;
+}
+interface BiasSinex {
+    records: BiasRecord[];
+    /** All satellite records grouped by PRN (station === ''). */
+    bySat: Map<string, BiasRecord[]>;
+}
+/** Parse a Bias-SINEX file. */
+declare function parseBiasSinex(text: string): BiasSinex;
+/** A satellite bias record's value in metres (ns → m via c; cyc left as-is). */
+declare function biasMetres(r: BiasRecord): number;
+/**
+ * Look up a satellite bias: a DSB by its (obs1,obs2) pair, or an OSB by obs1
+ * alone (pass obs2 = undefined), valid at `tMs` (default: any). Returns the
+ * record or null.
+ */
+declare function findSatBias(bias: BiasSinex, prn: string, obs1: string, obs2?: string, tMs?: number): BiasRecord | null;
+/**
+ * Fast time-aware lookup of satellite OSBs (metres) — the form PPP needs.
+ * Products with integer-recoverable clocks (CNES/GFZ "ABS", CODE) publish
+ * phase OSBs that change every 30 s–15 min, so one PRN/observable can carry
+ * thousands of records; they're sorted once and found by binary search.
+ * Values are as published (observed = true + bias), so the corrected
+ * observable is `observed − osb(...)`.
+ */
+declare class SatOsbTable {
+    private readonly byKey;
+    constructor(bias: BiasSinex | BiasRecord[]);
+    /** Whether any phase (L…) OSB is present — i.e. the product supports AR. */
+    get hasPhase(): boolean;
+    /** OSB (m) of `prn`'s observable `obs` (e.g. 'C1W', 'L2W') at `tMs`;
+     *  null when the product has none valid then. */
+    get(prn: string, obs: string, tMs: number): number | null;
+}
+
+export { type BiasRecord as B, type ClkFile as C, type IonexGrid as I, SatOsbTable as S, type BiasSinex as a, type Sp3File as b, type Sp3Sample as c, biasMetres as d, clkBias as e, findSatBias as f, parseClk as g, parseIonex as h, parseSp3 as i, parseBiasSinex as p, sp3Position as s };

package/dist/{clk-DGVBGRzA.d.ts → bias-sinex-OGMw9xEZ.d.ts}

@@ -93,4 +93,69 @@ declare function parseClk(text: string): ClkFile;
  */
 declare function clkBias(clk: ClkFile, prn: string, tMs: number): number | null;
 
-export { type ClkFile as C, type IonexGrid as I, type Sp3File as S, type Sp3Sample as a, parseIonex as b, clkBias as c, parseSp3 as d, parseClk as p, sp3Position as s };
+/**
+ * Bias-SINEX (`.BIA`) parser — GNSS code/phase biases published by the IGS
+ * analysis centres. Two record kinds matter here:
+ *
+ *   - **DSB** differential signal bias (OBS1−OBS2), the classic inter-code
+ *     bias used for ionosphere/DCB corrections in SPP.
+ *   - **OSB** observable-specific (absolute) bias for a single OBS, the
+ *     ingredient PPP-AR needs: the *phase* OSBs are the satellite
+ *     fractional-cycle biases that make the ambiguities integer-recoverable.
+ *
+ * Values keep their file unit (`ns` for code/most phase biases, `cyc` for
+ * cycle-based phase biases); `biasMetres` converts ns → m via c.
+ */
+interface BiasRecord {
+    /** 'DSB' | 'OSB' | 'ISB' | … */
+    type: string;
+    /** Satellite PRN (e.g. 'G01', 'E11'); '' for a pure station bias. */
+    prn: string;
+    /** Station (9-char), '' for satellite biases. */
+    station: string;
+    /** First observable code (e.g. 'C1W', 'L1C'). */
+    obs1: string;
+    /** Second observable (DSB only; '' for OSB). */
+    obs2: string;
+    /** Validity start / end (GPS-scale ms); ±Infinity when open. */
+    startMs: number;
+    endMs: number;
+    /** Estimated value in the file unit. */
+    value: number;
+    /** File unit ('ns' | 'cyc' | …). */
+    unit: string;
+}
+interface BiasSinex {
+    records: BiasRecord[];
+    /** All satellite records grouped by PRN (station === ''). */
+    bySat: Map<string, BiasRecord[]>;
+}
+/** Parse a Bias-SINEX file. */
+declare function parseBiasSinex(text: string): BiasSinex;
+/** A satellite bias record's value in metres (ns → m via c; cyc left as-is). */
+declare function biasMetres(r: BiasRecord): number;
+/**
+ * Look up a satellite bias: a DSB by its (obs1,obs2) pair, or an OSB by obs1
+ * alone (pass obs2 = undefined), valid at `tMs` (default: any). Returns the
+ * record or null.
+ */
+declare function findSatBias(bias: BiasSinex, prn: string, obs1: string, obs2?: string, tMs?: number): BiasRecord | null;
+/**
+ * Fast time-aware lookup of satellite OSBs (metres) — the form PPP needs.
+ * Products with integer-recoverable clocks (CNES/GFZ "ABS", CODE) publish
+ * phase OSBs that change every 30 s–15 min, so one PRN/observable can carry
+ * thousands of records; they're sorted once and found by binary search.
+ * Values are as published (observed = true + bias), so the corrected
+ * observable is `observed − osb(...)`.
+ */
+declare class SatOsbTable {
+    private readonly byKey;
+    constructor(bias: BiasSinex | BiasRecord[]);
+    /** Whether any phase (L…) OSB is present — i.e. the product supports AR. */
+    get hasPhase(): boolean;
+    /** OSB (m) of `prn`'s observable `obs` (e.g. 'C1W', 'L2W') at `tMs`;
+     *  null when the product has none valid then. */
+    get(prn: string, obs: string, tMs: number): number | null;
+}
+
+export { type BiasRecord as B, type ClkFile as C, type IonexGrid as I, SatOsbTable as S, type BiasSinex as a, type Sp3File as b, type Sp3Sample as c, biasMetres as d, clkBias as e, findSatBias as f, parseClk as g, parseIonex as h, parseSp3 as i, parseBiasSinex as p, sp3Position as s };

package/dist/{chunk-MDO77CCC.js → chunk-2DJ4XENO.js}

@@ -1617,6 +1617,50 @@ function findSatBias(bias, prn, obs1, obs2, tMs) {
   }
   return null;
 }
+var SatOsbTable = class {
+  byKey = /* @__PURE__ */ new Map();
+  constructor(bias) {
+    const recs = Array.isArray(bias) ? bias : bias.records;
+    const groups = /* @__PURE__ */ new Map();
+    for (const r of recs) {
+      if (r.type !== "OSB" || !r.prn || r.station) continue;
+      if (r.unit !== "ns") continue;
+      const k = `${r.prn} ${r.obs1}`;
+      let g = groups.get(k);
+      if (!g) groups.set(k, g = []);
+      g.push(r);
+    }
+    for (const [k, g] of groups) {
+      g.sort((a, b) => a.startMs - b.startMs);
+      this.byKey.set(k, {
+        start: Float64Array.from(g, (r) => r.startMs),
+        end: Float64Array.from(g, (r) => r.endMs),
+        m: Float64Array.from(g, biasMetres)
+      });
+    }
+  }
+  /** Whether any phase (L…) OSB is present — i.e. the product supports AR. */
+  get hasPhase() {
+    for (const k of this.byKey.keys())
+      if (k.split(" ")[1][0] === "L") return true;
+    return false;
+  }
+  /** OSB (m) of `prn`'s observable `obs` (e.g. 'C1W', 'L2W') at `tMs`;
+   *  null when the product has none valid then. */
+  get(prn, obs, tMs) {
+    const e = this.byKey.get(`${prn} ${obs}`) ?? this.byKey.get(`${prn} ${obs.slice(0, 2)}X`);
+    if (!e) return null;
+    let lo = 0;
+    let hi = e.start.length - 1;
+    if (tMs < e.start[0]) return null;
+    while (lo < hi) {
+      const mid = lo + hi + 1 >> 1;
+      if (e.start[mid] <= tMs) lo = mid;
+      else hi = mid - 1;
+    }
+    return tMs <= e.end[lo] ? e.m[lo] : null;
+  }
+};
 
 // src/rinex/obs-writer-common.ts
 function createGzipLineSink() {
@@ -2021,6 +2065,7 @@ export {
   parseBiasSinex,
   biasMetres,
   findSatBias,
+  SatOsbTable,
   createGzipLineSink,
   stationHeaderLines,
   writeRinexObsBlob,

package/dist/{chunk-IQHQHUOB.js → chunk-2KQMCJ5Z.js}

@@ -1,11 +1,13 @@
 import {
   BdsAssembler,
-  GalInavAssembler,
   GloStringAssembler,
   bdsSubframeParityOk,
-  galInavPageCrcOk,
   testGloString
-} from "./chunk-SSWFX3P3.js";
+} from "./chunk-NIJ6WFQJ.js";
+import {
+  GalInavAssembler,
+  galInavPageCrcOk
+} from "./chunk-TWGZTRBL.js";
 import {
   decodeSbasGeoNav,
   sbasCrcOk

package/dist/chunk-E3Y6YKQM.js

@@ -0,0 +1,65 @@
+// src/sbf/frame.ts
+var CRC_TABLE = (() => {
+  const t = new Uint16Array(256);
+  for (let n = 0; n < 256; n++) {
+    let c = n << 8;
+    for (let k = 0; k < 8; k++)
+      c = c & 32768 ? (c << 1 ^ 4129) & 65535 : c << 1 & 65535;
+    t[n] = c;
+  }
+  return t;
+})();
+function crc16(data, start, end) {
+  let crc = 0;
+  for (let i = start; i < end; i++)
+    crc = (crc << 8 ^ CRC_TABLE[(crc >> 8 ^ data[i]) & 255]) & 65535;
+  return crc;
+}
+function scanSbfFrames(data, view, onBlock) {
+  let badCrc = 0;
+  let i = 0;
+  while (i + 8 <= data.length) {
+    if (data[i] !== 36 || data[i + 1] !== 64) {
+      i++;
+      continue;
+    }
+    const len = view.getUint16(i + 6, true);
+    if (len < 8 || len % 4 !== 0) {
+      i++;
+      continue;
+    }
+    if (i + len > data.length) break;
+    if (crc16(data, i + 4, i + len) !== view.getUint16(i + 2, true)) {
+      badCrc++;
+      i++;
+      continue;
+    }
+    onBlock(view.getUint16(i + 4, true) & 8191, i, len);
+    i += len;
+  }
+  return badCrc;
+}
+var two = (n) => String(n).padStart(2, "0");
+function svidToPrn(svid) {
+  if (svid >= 1 && svid <= 37) return svid <= 32 ? `G${two(svid)}` : null;
+  if (svid <= 61) return svid - 37 <= 27 ? `R${two(svid - 37)}` : null;
+  if (svid <= 62) return null;
+  if (svid <= 68) return svid - 38 <= 27 ? `R${two(svid - 38)}` : null;
+  if (svid <= 70) return null;
+  if (svid <= 106) return svid - 70 <= 36 ? `E${two(svid - 70)}` : null;
+  if (svid <= 119) return null;
+  if (svid <= 140) return `S${two(svid - 100)}`;
+  if (svid <= 180) return svid - 140 <= 50 ? `C${two(svid - 140)}` : null;
+  if (svid <= 190) return `J${two(svid - 180)}`;
+  if (svid <= 197) return svid - 190 <= 14 ? `I${two(svid - 190)}` : null;
+  if (svid <= 215) return `S${two(svid - 157)}`;
+  if (svid <= 222) return svid - 208 <= 14 ? `I${two(svid - 208)}` : null;
+  if (svid <= 245) return svid - 182 <= 50 ? `C${two(svid - 182)}` : null;
+  return null;
+}
+
+export {
+  crc16,
+  scanSbfFrames,
+  svidToPrn
+};

package/dist/{chunk-3U5AX7PY.js → chunk-FEABDV5P.js}

@@ -5,20 +5,14 @@ import {
 // src/rinex/crx.ts
 function crxRepair(old, diff) {
   const chars = old.split("");
-  for (let i = 0; i < diff.length; i++) {
+  const n = Math.min(old.length, diff.length);
+  for (let i = 0; i < n; i++) {
     const dc = diff[i];
-    if (dc === " ") {
-    } else if (dc === "&") {
-      chars[i] = " ";
-    } else {
-      chars[i] = dc;
-    }
-  }
-  if (diff.length > chars.length) {
-    for (let i = chars.length; i < diff.length; i++) {
-      chars[i] = diff[i] === "&" ? " " : diff[i];
-    }
+    if (dc === " ") continue;
+    chars[i] = dc === "&" ? " " : dc;
   }
+  for (let i = old.length; i < diff.length; i++)
+    chars[i] = diff[i] === "&" ? " " : diff[i];
   return chars.join("");
 }
 function parseCrxDataLine(line, ntype) {
@@ -266,6 +260,15 @@ function readObsValue(line, index) {
   const v = parseFloat(raw);
   return isNaN(v) ? null : v;
 }
+function normalizePrn(id) {
+  let sys = id[0];
+  let num = id.substring(1);
+  if (sys === " " || /\d/.test(sys)) {
+    num = id;
+    sys = "G";
+  }
+  return sys + num.replace(/\s/g, "").padStart(2, "0");
+}
 function parseEpochLine3(line) {
   if (line[0] !== ">") return null;
   const yr = parseInt(line.substring(2, 6));
@@ -415,7 +418,7 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
   function processSatLineV3(line) {
     if (!header || line.length < 3) return;
     const sys = line[0];
-    const prn = line.substring(0, 3);
+    const prn = normalizePrn(line.substring(0, 3));
     trackSat(sys, prn);
     const obsLine = line.substring(3);
     const codes = header.obsTypes[sys] ?? [];
@@ -440,10 +443,8 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
   let v2ObsAccum = [];
   function processSatLineV2(line) {
     if (!header || v2SatIndex >= v2SatIds.length) return;
-    const satId = v2SatIds[v2SatIndex];
-    let sys = satId[0];
-    if (sys === " " || /\d/.test(sys)) sys = "G";
-    const prn = sys + satId.substring(1).padStart(2, "0");
+    const prn = normalizePrn(v2SatIds[v2SatIndex]);
+    const sys = prn[0];
     if (v2CurrentSatLine === 0) {
       trackSat(sys, prn);
       if (onSatObs) v2ObsAccum = [];
@@ -534,6 +535,7 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
         } else {
           epochLine = crxRepair(crxPrevEpochLine, line);
         }
+        epochLine = epochLine.replace(/ +$/, "");
         if (isV3) {
           if (epochLine[0] !== ">") return;
           const info = parseEpochLine3(epochLine);
@@ -551,10 +553,11 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
             return;
           }
           resetEpochAccum(info);
-          const satPart = epochLine.substring(41).replace(/\s+$/, "");
           crxSatList = [];
-          for (let i = 0; i + 3 <= satPart.length; i += 3) {
-            crxSatList.push(satPart.substring(i, i + 3));
+          for (let k = 0; k < info.numSats; k++) {
+            crxSatList.push(
+              normalizePrn(epochLine.substring(41 + 3 * k, 44 + 3 * k))
+            );
           }
           crxSatIndex = 0;
           crxPrevEpochLine = epochLine;
@@ -575,11 +578,7 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
             return;
           }
           resetEpochAccum(info);
-          crxSatList = (info.satIds ?? []).map((id) => {
-            let sys = id[0];
-            if (sys === " " || /\d/.test(sys)) sys = "G";
-            return sys + id.substring(1).padStart(2, "0");
-          });
+          crxSatList = (info.satIds ?? []).slice(0, info.numSats).map(normalizePrn);
           crxSatIndex = 0;
           crxPrevEpochLine = epochLine;
           crxPhase = "clock";

package/dist/{chunk-RFPR4LKM.js → chunk-HWHD5E6V.js}

@@ -1,13 +1,19 @@
+import {
+  scanSbfFrames,
+  svidToPrn
+} from "./chunk-E3Y6YKQM.js";
 import {
   BdsAssembler,
-  GalFnavAssembler,
-  GalInavAssembler,
   GloStringAssembler,
   bdsSubframeParityOk,
-  galFnavPageCrcOk,
-  galInavPageCrcOk,
   testGloString
-} from "./chunk-SSWFX3P3.js";
+} from "./chunk-NIJ6WFQJ.js";
+import {
+  GalFnavAssembler,
+  GalInavAssembler,
+  galFnavPageCrcOk,
+  galInavPageCrcOk
+} from "./chunk-TWGZTRBL.js";
 import {
   decodeSbasGeoNav,
   sbasCrcOk
@@ -40,66 +46,6 @@ import {
   setBitU
 } from "./chunk-V2D5KUKH.js";
 
-// src/sbf/frame.ts
-var CRC_TABLE = (() => {
-  const t = new Uint16Array(256);
-  for (let n = 0; n < 256; n++) {
-    let c = n << 8;
-    for (let k = 0; k < 8; k++)
-      c = c & 32768 ? (c << 1 ^ 4129) & 65535 : c << 1 & 65535;
-    t[n] = c;
-  }
-  return t;
-})();
-function crc16(data, start, end) {
-  let crc = 0;
-  for (let i = start; i < end; i++)
-    crc = (crc << 8 ^ CRC_TABLE[(crc >> 8 ^ data[i]) & 255]) & 65535;
-  return crc;
-}
-function scanSbfFrames(data, view, onBlock) {
-  let badCrc = 0;
-  let i = 0;
-  while (i + 8 <= data.length) {
-    if (data[i] !== 36 || data[i + 1] !== 64) {
-      i++;
-      continue;
-    }
-    const len = view.getUint16(i + 6, true);
-    if (len < 8 || len % 4 !== 0) {
-      i++;
-      continue;
-    }
-    if (i + len > data.length) break;
-    if (crc16(data, i + 4, i + len) !== view.getUint16(i + 2, true)) {
-      badCrc++;
-      i++;
-      continue;
-    }
-    onBlock(view.getUint16(i + 4, true) & 8191, i, len);
-    i += len;
-  }
-  return badCrc;
-}
-var two = (n) => String(n).padStart(2, "0");
-function svidToPrn(svid) {
-  if (svid >= 1 && svid <= 37) return svid <= 32 ? `G${two(svid)}` : null;
-  if (svid <= 61) return svid - 37 <= 27 ? `R${two(svid - 37)}` : null;
-  if (svid <= 62) return null;
-  if (svid <= 68) return svid - 38 <= 27 ? `R${two(svid - 38)}` : null;
-  if (svid <= 70) return null;
-  if (svid <= 106) return svid - 70 <= 36 ? `E${two(svid - 70)}` : null;
-  if (svid <= 119) return null;
-  if (svid <= 140) return `S${two(svid - 100)}`;
-  if (svid <= 180) return svid - 140 <= 50 ? `C${two(svid - 140)}` : null;
-  if (svid <= 190) return `J${two(svid - 180)}`;
-  if (svid <= 197) return svid - 190 <= 14 ? `I${two(svid - 190)}` : null;
-  if (svid <= 215) return `S${two(svid - 157)}`;
-  if (svid <= 222) return svid - 208 <= 14 ? `I${two(svid - 208)}` : null;
-  if (svid <= 245) return svid - 182 <= 50 ? `C${two(svid - 182)}` : null;
-  return null;
-}
-
 // src/sbf/nav.ts
 var PI = Math.PI;
 var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
@@ -1981,23 +1927,23 @@ var MEAS2_SIG = [
   ["J", "1E"],
   ["J", "5P"]
 ];
-var two2 = (n) => String(n).padStart(2, "0");
+var two = (n) => String(n).padStart(2, "0");
 function m3Prn(navsys, svid) {
   switch (navsys) {
     case 0:
-      return svid < 32 ? `G${two2(svid + 1)}` : null;
+      return svid < 32 ? `G${two(svid + 1)}` : null;
     case 1:
-      return svid < 27 ? `R${two2(svid + 1)}` : null;
+      return svid < 27 ? `R${two(svid + 1)}` : null;
     case 2:
-      return svid < 36 ? `E${two2(svid + 1)}` : null;
+      return svid < 36 ? `E${two(svid + 1)}` : null;
     case 3:
-      return svid < 50 ? `C${two2(svid + 1)}` : null;
+      return svid < 50 ? `C${two(svid + 1)}` : null;
     case 4:
-      return svid <= 38 ? `S${two2(svid + 20)}` : null;
+      return svid <= 38 ? `S${two(svid + 20)}` : null;
     case 5:
-      return svid < 10 ? `J${two2(svid + 1)}` : null;
+      return svid < 10 ? `J${two(svid + 1)}` : null;
     case 6:
-      return svid < 14 ? `I${two2(svid + 1)}` : null;
+      return svid < 14 ? `I${two(svid + 1)}` : null;
     default:
       return null;
   }
@@ -2548,8 +2494,6 @@ function parseSbfMeas(data) {
 }
 
 export {
-  crc16,
-  scanSbfFrames,
   parseSbfNav,
   parseSbfAlmanac,
   parseSbfIonoUtc,

package/dist/{chunk-KUSUN7DP.js → chunk-MFCPYUEI.js}

@@ -1,6 +1,6 @@
 import {
   parseRinexStream
-} from "./chunk-3U5AX7PY.js";
+} from "./chunk-FEABDV5P.js";
 import {
   ARC_GAP_FACTOR,
   BAND_LABELS,