gnss-js 2.10.1 → 2.12.0

package/README.md

@@ -16,7 +16,8 @@ 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, GLONASS strings, and Galileo HAS (E6 C/NAV) corrections.
+- **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/binex.d.cts

@@ -1,4 +1,4 @@
-import { E as Ephemeris } from './nav-DImXUdbp.cjs';
+import { E as Ephemeris } from './nav-COp0Vmlw.cjs';
 
 /**
  * BINEX (BINary EXchange) framing, shared by the observation and

package/dist/binex.d.ts

@@ -1,4 +1,4 @@
-import { E as Ephemeris } from './nav-DImXUdbp.js';
+import { E as Ephemeris } from './nav-COp0Vmlw.js';
 
 /**
  * BINEX (BINary EXchange) framing, shared by the observation and

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

@@ -1,7 +1,10 @@
+import {
+  bdsCnavAref
+} from "./chunk-NL2ORX75.js";
 import {
   CNAV_A_REF,
   CNAV_OMEGA_DOT_REF
-} from "./chunk-YUW33VXK.js";
+} from "./chunk-HUXQABJI.js";
 import {
   GPS_PI
 } from "./chunk-V2D5KUKH.js";
@@ -661,10 +664,11 @@ var R4_MSG_LINES = {
   // RINEX 4 adds orbit-4 (status flags) vs 3 in RINEX 3
   // EPH type parsed into RinexCnavEphemeris for G/J (skipped otherwise)
   CNAV: 9,
-  // EPH types we skip (different field layout)
+  // BeiDou-3 B-CNAV, parsed into BdsCnavEphemeris
   CNV1: 10,
   CNV2: 10,
   CNV3: 9,
+  // EPH types we skip (different field layout)
   L1NV: 8,
   L1OC: 9,
   L3OC: 9
@@ -873,6 +877,68 @@ function buildCnavEphemeris(sys, prn, date, epochVals, d) {
     tow: num(7, 0)
   };
 }
+var BDT_EPOCH_MS = Date.UTC(2006, 0, 1);
+function buildBdsCnavEphemeris(message, prn, date, epochVals, d) {
+  const num = (r, c) => d[r]?.[c] ?? 0;
+  const opt = (r, c) => d[r]?.[c] ?? null;
+  const cnv3 = message === "CNV3";
+  const t = (date.getTime() - BDT_EPOCH_MS) / 1e3;
+  const weekToc = Math.floor(t / SEC_PER_WEEK);
+  const toc = t - weekToc * SEC_PER_WEEK;
+  const toe = num(2, 0);
+  const week = toe - toc > SEC_PER_WEEK / 2 ? weekToc - 1 : toc - toe > SEC_PER_WEEK / 2 ? weekToc + 1 : weekToc;
+  const sqrtA = num(1, 3);
+  const a = sqrtA * sqrtA;
+  const satType = num(4, 2);
+  return {
+    system: "C",
+    prn,
+    message,
+    week,
+    toe,
+    toeDate: new Date(BDT_EPOCH_MS + (week * SEC_PER_WEEK + toe) * 1e3),
+    toc,
+    tocDate: date,
+    satType,
+    a,
+    deltaA: a - bdsCnavAref(satType),
+    aDot: num(0, 0),
+    deltaN0: num(0, 2),
+    deltaN0Dot: num(4, 1),
+    m0: num(0, 3),
+    e: num(1, 1),
+    omega: num(3, 2),
+    omega0: num(2, 2),
+    i0: num(3, 0),
+    omegaDot: num(3, 3),
+    i0Dot: num(4, 0),
+    cis: num(2, 3),
+    cic: num(2, 1),
+    crs: num(0, 1),
+    crc: num(3, 1),
+    cus: num(1, 2),
+    cuc: num(1, 0),
+    af0: epochVals[0] ?? 0,
+    af1: epochVals[1] ?? 0,
+    af2: epochVals[2] ?? 0,
+    iode: cnv3 ? null : opt(8, 3),
+    iodc: cnv3 ? null : opt(7, 3),
+    tgdB1Cp: cnv3 ? null : opt(6, 2),
+    tgdB2ap: cnv3 ? null : opt(6, 3),
+    iscB1Cd: cnv3 ? null : opt(6, 0),
+    iscB2ad: cnv3 ? null : opt(6, 1),
+    tgdB2bI: cnv3 ? opt(6, 3) : null,
+    top: opt(4, 3),
+    sisaiOe: opt(5, 0),
+    sisaiOcb: opt(5, 1),
+    sisaiOc1: opt(5, 2),
+    sisaiOc2: opt(5, 3),
+    sismai: cnv3 ? opt(6, 0) : opt(7, 0),
+    health: cnv3 ? opt(6, 1) : opt(7, 1),
+    integrityFlags: cnv3 ? opt(6, 2) : opt(7, 2),
+    tow: cnv3 ? num(7, 0) : num(8, 0)
+  };
+}
 function parseNavFile(text) {
   const lines = text.split("\n");
   const header = {
@@ -883,6 +949,7 @@ function parseNavFile(text) {
   };
   const ephemerides = [];
   const cnav = [];
+  const bdsCnav = [];
   let inHeader = true;
   let i = 0;
   while (i < lines.length) {
@@ -951,7 +1018,8 @@ function parseNavFile(text) {
         continue;
       }
       const isCnav = msgType === "CNAV" && (svSys === "G" || svSys === "J");
-      if (!SUPPORTED_EPH_MSGS.has(msgType) && !isCnav) {
+      const isBdsCnav = svSys === "C" && (msgType === "CNV1" || msgType === "CNV2" || msgType === "CNV3");
+      if (!SUPPORTED_EPH_MSGS.has(msgType) && !isCnav && !isBdsCnav) {
         const totalLines = R4_MSG_LINES[msgType] ?? 2;
         i += totalLines + 1;
         continue;
@@ -963,17 +1031,34 @@ function parseNavFile(text) {
       const sys2 = epochLine.charAt(0);
       const parsed = parseNavEpochV3(epochLine);
       const prn2 = `${sys2}${parsed.prn.substring(1)}`;
-      if (isCnav) {
+      if (isCnav || isBdsCnav) {
         const slotLines = [];
         for (let j = 0; j < numDataLines2; j++) {
           i++;
           if (i >= lines.length) break;
           slotLines.push(parseDataLineSlots(lines[i]));
         }
-        if (slotLines.length === numDataLines2 && (sys2 === "G" || sys2 === "J")) {
-          cnav.push(
-            buildCnavEphemeris(sys2, prn2, parsed.date, parsed.values, slotLines)
-          );
+        if (slotLines.length === numDataLines2) {
+          if (isBdsCnav && sys2 === "C")
+            bdsCnav.push(
+              buildBdsCnavEphemeris(
+                msgType,
+                prn2,
+                parsed.date,
+                parsed.values,
+                slotLines
+              )
+            );
+          else if (sys2 === "G" || sys2 === "J")
+            cnav.push(
+              buildCnavEphemeris(
+                sys2,
+                prn2,
+                parsed.date,
+                parsed.values,
+                slotLines
+              )
+            );
         }
         i++;
         continue;
@@ -1051,7 +1136,12 @@ function parseNavFile(text) {
     }
     i++;
   }
-  return cnav.length > 0 ? { header, ephemerides, cnav } : { header, ephemerides };
+  return {
+    header,
+    ephemerides,
+    ...cnav.length > 0 && { cnav },
+    ...bdsCnav.length > 0 && { bdsCnav }
+  };
 }
 
 // src/rinex/nav-writer.ts
@@ -1247,6 +1337,30 @@ function cnavRecord(eph) {
     orbitLine4([eph.tow, resolveWnOp(eph.week, eph.wnOp)])
   ].join("\n");
 }
+function bdsCnavRecord(eph) {
+  const lines = [
+    epochLine4(eph.prn, eph.tocDate, [eph.af0, eph.af1, eph.af2]),
+    orbitLine4([eph.aDot, eph.crs, eph.deltaN0, eph.m0]),
+    orbitLine4([eph.cuc, eph.e, eph.cus, Math.sqrt(eph.a)]),
+    orbitLine4([eph.toe, eph.cic, eph.omega0, eph.cis]),
+    orbitLine4([eph.i0, eph.crc, eph.omega, eph.omegaDot]),
+    orbitLine4([eph.i0Dot, eph.deltaN0Dot, eph.satType, eph.top]),
+    orbitLine4([eph.sisaiOe, eph.sisaiOcb, eph.sisaiOc1, eph.sisaiOc2])
+  ];
+  if (eph.message === "CNV3") {
+    lines.push(
+      orbitLine4([eph.sismai, eph.health, eph.integrityFlags, eph.tgdB2bI]),
+      orbitLine4([eph.tow])
+    );
+  } else {
+    lines.push(
+      orbitLine4([eph.iscB1Cd, eph.iscB2ad, eph.tgdB1Cp, eph.tgdB2ap]),
+      orbitLine4([eph.sismai, eph.health, eph.integrityFlags, eph.iodc]),
+      orbitLine4([eph.tow, null, null, eph.iode])
+    );
+  }
+  return lines.join("\n");
+}
 var KLOBUCHAR_SOURCES = [
   ["GPSA", "GPSB", "G", "LNAV"],
   ["QZSA", "QZSB", "J", "LNAV"],
@@ -1324,6 +1438,15 @@ function writeRinexNav4(nav) {
 ` + cnavRecord(eph)
     });
   }
+  for (const eph of nav.bdsCnav ?? []) {
+    records.push({
+      prn: eph.prn,
+      time: eph.tocDate.getTime(),
+      label: eph.message,
+      body: `> EPH ${padL(eph.prn, 3)} ${eph.message}
+` + bdsCnavRecord(eph)
+    });
+  }
   records.sort(
     (a, b) => a.prn.localeCompare(b.prn) || a.time - b.time || a.label.localeCompare(b.label)
   );
@@ -1494,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() {
@@ -1898,6 +2065,7 @@ export {
   parseBiasSinex,
   biasMetres,
   findSatBias,
+  SatOsbTable,
   createGzipLineSink,
   stationHeaderLines,
   writeRinexObsBlob,

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

@@ -1,19 +1,21 @@
 import {
   BdsAssembler,
-  GalInavAssembler,
   GloStringAssembler,
   bdsSubframeParityOk,
-  galInavPageCrcOk,
   testGloString
-} from "./chunk-3KSDL3CO.js";
+} from "./chunk-NIJ6WFQJ.js";
+import {
+  GalInavAssembler,
+  galInavPageCrcOk
+} from "./chunk-TWGZTRBL.js";
 import {
   decodeSbasGeoNav,
   sbasCrcOk
-} from "./chunk-WQX354KD.js";
+} from "./chunk-H3AKFD25.js";
 import {
   CnavAssembler,
   cnavCrcOk
-} from "./chunk-YUW33VXK.js";
+} from "./chunk-HUXQABJI.js";
 import {
   GpsLnavAssembler,
   getBitS,

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
+};