gnss-js 2.10.1 → 2.12.0

package/dist/rinex.cjs

@@ -22,6 +22,7 @@ var rinex_exports = {};
 __export(rinex_exports, {
   EMPTY_WARNINGS: () => EMPTY_WARNINGS,
   SYSTEM_ORDER: () => SYSTEM_ORDER,
+  SatOsbTable: () => SatOsbTable,
   WarningAccumulator: () => WarningAccumulator,
   biasMetres: () => biasMetres,
   clkBias: () => clkBias,
@@ -59,20 +60,14 @@ module.exports = __toCommonJS(rinex_exports);
 // 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) {
@@ -782,6 +777,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));
@@ -931,7 +935,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] ?? [];
@@ -956,10 +960,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 = [];
@@ -1050,6 +1052,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);
@@ -1067,10 +1070,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;
@@ -1091,11 +1095,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";
@@ -1495,6 +1495,15 @@ var SEC_PER_WEEK2 = 7 * 86400;
 var CNAV_A_REF = 26559710;
 var CNAV_OMEGA_DOT_REF = -26e-10;
 
+// src/navbits/bds-cnav.ts
+var BDT_EPOCH_MS = Date.UTC(2006, 0, 1);
+var SEC_PER_WEEK3 = 7 * 86400;
+var BDS_CNAV_A_REF_MEO = 27906100;
+var BDS_CNAV_A_REF_IGSO_GEO = 42162200;
+function bdsCnavAref(satType) {
+  return satType === 3 ? BDS_CNAV_A_REF_MEO : BDS_CNAV_A_REF_IGSO_GEO;
+}
+
 // src/rinex/nav.ts
 var DATA_LINES = {
   G: 7,
@@ -1520,16 +1529,17 @@ 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
 };
 var GPS_EPOCH_MS3 = Date.UTC(1980, 0, 6);
-var SEC_PER_WEEK3 = 7 * 86400;
+var SEC_PER_WEEK4 = 7 * 86400;
 var QZSS_CNAV_A_REF = 42164200;
 var SUPPORTED_EPH_MSGS = /* @__PURE__ */ new Set([
   "LNAV",
@@ -1674,8 +1684,8 @@ function buildCnavEphemeris(sys, prn, date, epochVals, d) {
   const num = (r, c) => d[r]?.[c] ?? 0;
   const opt = (r, c) => d[r]?.[c] ?? null;
   const t = (date.getTime() - GPS_EPOCH_MS3) / 1e3;
-  const week = Math.floor(t / SEC_PER_WEEK3);
-  const toc = t - week * SEC_PER_WEEK3;
+  const week = Math.floor(t / SEC_PER_WEEK4);
+  const toc = t - week * SEC_PER_WEEK4;
   const sqrtA = num(1, 3);
   const a = sqrtA * sqrtA;
   const aRef = sys === "J" ? QZSS_CNAV_A_REF : CNAV_A_REF;
@@ -1732,6 +1742,68 @@ function buildCnavEphemeris(sys, prn, date, epochVals, d) {
     tow: num(7, 0)
   };
 }
+var BDT_EPOCH_MS2 = 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_MS2) / 1e3;
+  const weekToc = Math.floor(t / SEC_PER_WEEK4);
+  const toc = t - weekToc * SEC_PER_WEEK4;
+  const toe = num(2, 0);
+  const week = toe - toc > SEC_PER_WEEK4 / 2 ? weekToc - 1 : toc - toe > SEC_PER_WEEK4 / 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_MS2 + (week * SEC_PER_WEEK4 + 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 = {
@@ -1742,6 +1814,7 @@ function parseNavFile(text) {
   };
   const ephemerides = [];
   const cnav = [];
+  const bdsCnav = [];
   let inHeader = true;
   let i = 0;
   while (i < lines.length) {
@@ -1810,7 +1883,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;
@@ -1822,17 +1896,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;
@@ -1910,7 +2001,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
@@ -2106,6 +2202,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"],
@@ -2183,6 +2303,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)
   );
@@ -2529,6 +2658,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() {
@@ -2920,6 +3093,7 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
 0 && (module.exports = {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,
+  SatOsbTable,
   WarningAccumulator,
   biasMetres,
   clkBias,

package/dist/rinex.d.cts

@@ -1,8 +1,8 @@
 import { R as RinexHeader } from './parser-JPjjFgeP.cjs';
 export { E as EpochSummary, P as ParseOptions, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.cjs';
-import { a as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris } from './nav-DImXUdbp.cjs';
-export { G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-DImXUdbp.cjs';
-export { C as ClkFile, I as IonexGrid, S as Sp3File, a as Sp3Sample, c as clkBias, p as parseClk, b as parseIonex, d as parseSp3, s as sp3Position } from './clk-DGVBGRzA.cjs';
+import { b as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris, B as BdsCnavEphemeris } from './nav-COp0Vmlw.cjs';
+export { a as BdsCnavMessage, G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-COp0Vmlw.cjs';
+export { B as BiasRecord, a as BiasSinex, C as ClkFile, I as IonexGrid, S as SatOsbTable, b as Sp3File, c as Sp3Sample, d as biasMetres, e as clkBias, f as findSatBias, p as parseBiasSinex, g as parseClk, h as parseIonex, i as parseSp3, s as sp3Position } from './bias-sinex-OGMw9xEZ.cjs';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -155,7 +155,8 @@ declare function writeRinexNav(nav: NavResult): string;
  *
  * Produces `> EPH` data records for the classic ephemeris types
  * (GPS/QZSS/NavIC LNAV, Galileo INAV, BeiDou D1/D2, GLONASS FDMA,
- * SBAS), `> EPH ... CNAV` records for GPS/QZSS CNAV ephemerides, and
+ * SBAS), `> EPH ... CNAV` records for GPS/QZSS CNAV ephemerides, BeiDou-3
+ * `CNV1`/`CNV2`/`CNV3` records for B-CNAV ephemerides, and
  * `> ION` records for the header's Klobuchar/NeQuick-G coefficients,
  * per the RINEX 4.01 record tables (A9-A18, A21, A26, A32, A33).
  */
@@ -165,13 +166,15 @@ interface Nav4Input {
     header: NavHeader;
     ephemerides: Ephemeris[];
     cnav?: RinexCnavEphemeris[];
+    bdsCnav?: BdsCnavEphemeris[];
 }
 /**
  * Write a RINEX 4.01 mixed navigation file.
  *
  * Emits `> EPH` records for every classic ephemeris (message labels
- * per constellation: LNAV, INAV, D1/D2, FDMA, SBAS) and every CNAV
- * ephemeris (`CNAV`), sorted by PRN, epoch, then label; followed by
+ * per constellation: LNAV, INAV, D1/D2, FDMA, SBAS), every GPS/QZSS
+ * CNAV ephemeris (`CNAV`) and every BeiDou-3 B-CNAV ephemeris (`CNV1`/
+ * `CNV2`/`CNV3`), sorted by PRN, epoch, then label; followed by
  * `> ION` records when the header carries Klobuchar (GPSA/GPSB,
  * QZSA/QZSB, BDSA/BDSB, IRNA/IRNB) or NeQuick-G (GAL) coefficients.
  * ION records need a transmit-time epoch, which RINEX 3-style headers
@@ -182,54 +185,6 @@ interface Nav4Input {
  */
 declare function writeRinexNav4(nav: Nav4Input): string;
 
-/**
- * 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;
-
 /**
  * RINEX 3.04 / 4.01 observation file writer.
  *
@@ -320,4 +275,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { BdsCnavEphemeris, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.d.ts

@@ -1,8 +1,8 @@
 import { R as RinexHeader } from './parser-JPjjFgeP.js';
 export { E as EpochSummary, P as ParseOptions, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.js';
-import { a as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris } from './nav-DImXUdbp.js';
-export { G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-DImXUdbp.js';
-export { C as ClkFile, I as IonexGrid, S as Sp3File, a as Sp3Sample, c as clkBias, p as parseClk, b as parseIonex, d as parseSp3, s as sp3Position } from './clk-DGVBGRzA.js';
+import { b as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris, B as BdsCnavEphemeris } from './nav-COp0Vmlw.js';
+export { a as BdsCnavMessage, G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-COp0Vmlw.js';
+export { B as BiasRecord, a as BiasSinex, C as ClkFile, I as IonexGrid, S as SatOsbTable, b as Sp3File, c as Sp3Sample, d as biasMetres, e as clkBias, f as findSatBias, p as parseBiasSinex, g as parseClk, h as parseIonex, i as parseSp3, s as sp3Position } from './bias-sinex-OGMw9xEZ.js';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -155,7 +155,8 @@ declare function writeRinexNav(nav: NavResult): string;
  *
  * Produces `> EPH` data records for the classic ephemeris types
  * (GPS/QZSS/NavIC LNAV, Galileo INAV, BeiDou D1/D2, GLONASS FDMA,
- * SBAS), `> EPH ... CNAV` records for GPS/QZSS CNAV ephemerides, and
+ * SBAS), `> EPH ... CNAV` records for GPS/QZSS CNAV ephemerides, BeiDou-3
+ * `CNV1`/`CNV2`/`CNV3` records for B-CNAV ephemerides, and
  * `> ION` records for the header's Klobuchar/NeQuick-G coefficients,
  * per the RINEX 4.01 record tables (A9-A18, A21, A26, A32, A33).
  */
@@ -165,13 +166,15 @@ interface Nav4Input {
     header: NavHeader;
     ephemerides: Ephemeris[];
     cnav?: RinexCnavEphemeris[];
+    bdsCnav?: BdsCnavEphemeris[];
 }
 /**
  * Write a RINEX 4.01 mixed navigation file.
  *
  * Emits `> EPH` records for every classic ephemeris (message labels
- * per constellation: LNAV, INAV, D1/D2, FDMA, SBAS) and every CNAV
- * ephemeris (`CNAV`), sorted by PRN, epoch, then label; followed by
+ * per constellation: LNAV, INAV, D1/D2, FDMA, SBAS), every GPS/QZSS
+ * CNAV ephemeris (`CNAV`) and every BeiDou-3 B-CNAV ephemeris (`CNV1`/
+ * `CNV2`/`CNV3`), sorted by PRN, epoch, then label; followed by
  * `> ION` records when the header carries Klobuchar (GPSA/GPSB,
  * QZSA/QZSB, BDSA/BDSB, IRNA/IRNB) or NeQuick-G (GAL) coefficients.
  * ION records need a transmit-time epoch, which RINEX 3-style headers
@@ -182,54 +185,6 @@ interface Nav4Input {
  */
 declare function writeRinexNav4(nav: Nav4Input): string;
 
-/**
- * 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;
-
 /**
  * RINEX 3.04 / 4.01 observation file writer.
  *
@@ -320,4 +275,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { BdsCnavEphemeris, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.js

@@ -1,5 +1,6 @@
 import {
   EMPTY_WARNINGS,
+  SatOsbTable,
   WarningAccumulator,
   biasMetres,
   createGzipLineSink,
@@ -21,7 +22,7 @@ import {
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-56HM5B6T.js";
+} from "./chunk-2DJ4XENO.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -30,19 +31,21 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-3U5AX7PY.js";
+} from "./chunk-FEABDV5P.js";
+import "./chunk-NL2ORX75.js";
 import {
   clkBias,
   parseClk,
   parseSp3,
   sp3Position
 } from "./chunk-GPEP67TG.js";
-import "./chunk-YUW33VXK.js";
+import "./chunk-HUXQABJI.js";
 import "./chunk-V2D5KUKH.js";
 import "./chunk-FIEWO4J4.js";
 export {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,
+  SatOsbTable,
   WarningAccumulator,
   biasMetres,
   clkBias,