gnss-js 1.29.0 → 1.31.0

package/dist/rinex.cjs

@@ -23,14 +23,19 @@ __export(rinex_exports, {
   EMPTY_WARNINGS: () => EMPTY_WARNINGS,
   SYSTEM_ORDER: () => SYSTEM_ORDER,
   WarningAccumulator: () => WarningAccumulator,
+  biasMetres: () => biasMetres,
+  clkBias: () => clkBias,
   createGzipLineSink: () => createGzipLineSink,
   crxDecompress: () => crxDecompress,
   crxRepair: () => crxRepair,
+  findSatBias: () => findSatBias,
   fmtD: () => fmtD,
   fmtF: () => fmtF,
   hdrLine: () => hdrLine,
   padL: () => padL,
   padR: () => padR,
+  parseBiasSinex: () => parseBiasSinex,
+  parseClk: () => parseClk,
   parseCrxDataLine: () => parseCrxDataLine,
   parseIonex: () => parseIonex,
   parseNavFile: () => parseNavFile,
@@ -595,6 +600,7 @@ var SYSTEM_NAMES = {
   I: "NavIC",
   S: "SBAS"
 };
+var C_LIGHT = 299792458;
 
 // src/rinex/parser.ts
 var CHUNK_SIZE = 2 * 1024 * 1024;
@@ -2388,6 +2394,141 @@ function sp3Position(sp3, prn, tMs, order = 9) {
   return { x, y, z, clk };
 }
 
+// src/rinex/clk.ts
+var GPS_EPOCH_MS4 = Date.UTC(1980, 0, 6);
+function parseClk(text) {
+  const sats = {};
+  let inData = false;
+  for (const line of text.split("\n")) {
+    if (!inData) {
+      if (line.includes("END OF HEADER")) inData = true;
+      continue;
+    }
+    if (line.length < 3 || line[0] !== "A" || line[1] !== "S") continue;
+    const p = line.slice(2).trim().split(/\s+/);
+    if (p.length < 9) continue;
+    const prn = p[0];
+    const y = +p[1];
+    const mo = +p[2];
+    const d = +p[3];
+    const h = +p[4];
+    const mi = +p[5];
+    const s = +p[6];
+    const bias = +p[8];
+    if (!Number.isFinite(bias)) continue;
+    const t = Date.UTC(y, mo - 1, d, h, mi, Math.floor(s)) + s % 1 * 1e3;
+    let e = sats[prn];
+    if (!e) {
+      e = { t: [], c: [] };
+      sats[prn] = e;
+    }
+    e.t.push(t);
+    e.c.push(bias);
+  }
+  let intervalSec = 30;
+  for (const e of Object.values(sats)) {
+    if (e.t.length >= 2) {
+      intervalSec = Math.round((e.t[1] - e.t[0]) / 1e3);
+      break;
+    }
+  }
+  return { sats, intervalSec };
+}
+function clkBias(clk, prn, tMs) {
+  const e = clk.sats[prn];
+  if (!e || e.t.length === 0) return null;
+  const { t, c } = e;
+  let lo = 0;
+  let hi = t.length - 1;
+  if (tMs < t[0] || tMs > t[hi]) return null;
+  while (lo < hi) {
+    const mid = lo + hi + 1 >> 1;
+    if (t[mid] <= tMs) lo = mid;
+    else hi = mid - 1;
+  }
+  if (t[lo] === tMs || lo === t.length - 1) return c[lo];
+  const t0 = t[lo];
+  const t1 = t[lo + 1];
+  if (t1 - t0 > clk.intervalSec * 1e3 * 4) return null;
+  const f = (tMs - t0) / (t1 - t0);
+  return c[lo] + f * (c[lo + 1] - c[lo]);
+}
+
+// src/rinex/bias-sinex.ts
+function parseEpoch(s, openTo) {
+  const m = /^(\d{4}):(\d{3}):(\d{5})$/.exec(s.trim());
+  if (!m) return openTo;
+  const y = +m[1];
+  if (y === 0) return openTo;
+  const doy = +m[2];
+  const sod = +m[3];
+  return Date.UTC(y, 0, 1) + (doy - 1) * 864e5 + sod * 1e3;
+}
+function parseBiasSinex(text) {
+  const records = [];
+  let inBlock = false;
+  for (const line of text.split("\n")) {
+    if (line.startsWith("+BIAS/SOLUTION")) {
+      inBlock = true;
+      continue;
+    }
+    if (line.startsWith("-BIAS/SOLUTION")) {
+      inBlock = false;
+      continue;
+    }
+    if (!inBlock || line.length < 65 || line.startsWith("*")) continue;
+    const type = line.slice(1, 4).trim();
+    if (type !== "DSB" && type !== "OSB" && type !== "ISB") continue;
+    const prn = line.slice(11, 14).trim();
+    const station = line.slice(15, 24).trim();
+    const obs1 = line.slice(25, 29).trim();
+    const obs2 = line.slice(30, 34).trim();
+    const startMs = parseEpoch(line.slice(35, 49), -Infinity);
+    const endMs = parseEpoch(line.slice(50, 64), Infinity);
+    const tail = line.slice(65).trim().split(/\s+/);
+    if (tail.length < 2) continue;
+    const unit = tail[0];
+    const value = Number(tail[1]);
+    if (!Number.isFinite(value)) continue;
+    records.push({
+      type,
+      prn,
+      station,
+      obs1,
+      obs2,
+      startMs,
+      endMs,
+      value,
+      unit
+    });
+  }
+  const bySat = /* @__PURE__ */ new Map();
+  for (const r of records) {
+    if (!r.prn || r.station) continue;
+    let a = bySat.get(r.prn);
+    if (!a) {
+      a = [];
+      bySat.set(r.prn, a);
+    }
+    a.push(r);
+  }
+  return { records, bySat };
+}
+function biasMetres(r) {
+  return r.unit === "ns" ? r.value * 1e-9 * C_LIGHT : r.value;
+}
+function findSatBias(bias, prn, obs1, obs2, tMs) {
+  const recs = bias.bySat.get(prn);
+  if (!recs) return null;
+  for (const r of recs) {
+    if (r.obs1 !== obs1) continue;
+    if (obs2 === void 0 ? r.obs2 !== "" : r.obs2 !== obs2) continue;
+    if (tMs != null && (tMs < r.startMs || tMs > r.endMs)) continue;
+    return r;
+  }
+  return null;
+}
+
 // src/rinex/obs-writer-common.ts
 function createGzipLineSink() {
   const encoder = new TextEncoder();
@@ -2766,14 +2907,19 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,
   WarningAccumulator,
+  biasMetres,
+  clkBias,
   createGzipLineSink,
   crxDecompress,
   crxRepair,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
+  parseClk,
   parseCrxDataLine,
   parseIonex,
   parseNavFile,

package/dist/rinex.d.cts

@@ -2,7 +2,7 @@ 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 { I as IonexGrid, S as Sp3File, a as Sp3Sample, p as parseIonex, b as parseSp3, s as sp3Position } from './sp3-BAq8l39b.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';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -182,6 +182,54 @@ 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.
  *
@@ -263,4 +311,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { 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, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+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, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.d.ts

@@ -2,7 +2,7 @@ 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 { I as IonexGrid, S as Sp3File, a as Sp3Sample, p as parseIonex, b as parseSp3, s as sp3Position } from './sp3-BAq8l39b.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';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -182,6 +182,54 @@ 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.
  *
@@ -263,4 +311,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { 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, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+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, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.js

@@ -1,12 +1,15 @@
 import {
   EMPTY_WARNINGS,
   WarningAccumulator,
+  biasMetres,
   createGzipLineSink,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
   parseIonex,
   parseNavFile,
   stationHeaderLines,
@@ -17,7 +20,7 @@ import {
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-6GTS5MVI.js";
+} from "./chunk-5RNDQQ6J.js";
 import "./chunk-726LQBGM.js";
 import "./chunk-IQ5OIMQD.js";
 import {
@@ -30,22 +33,29 @@ import {
   systemName
 } from "./chunk-3U5AX7PY.js";
 import {
+  clkBias,
+  parseClk,
   parseSp3,
   sp3Position
-} from "./chunk-4KMWWCUR.js";
+} from "./chunk-GPEP67TG.js";
 import "./chunk-FIEWO4J4.js";
 export {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,
   WarningAccumulator,
+  biasMetres,
+  clkBias,
   createGzipLineSink,
   crxDecompress,
   crxRepair,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
+  parseClk,
   parseCrxDataLine,
   parseIonex,
   parseNavFile,

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.29.0",
+  "version": "1.31.0",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,