gnss-js 1.4.1 → 1.5.1

package/README.md

@@ -94,12 +94,16 @@ const itrf = transformFrame(
 ```ts
 import { parseRinexStream, parseNavFile } from 'gnss-js/rinex';
 
-// Streaming observation parser (works with ReadableStream)
-const result = await parseRinexStream(stream, {
-  onObservation: (time, prn, codes, values) => {
-    // process each epoch
-  },
-});
+// Streaming observation parser (positional arguments)
+const result = await parseRinexStream(
+  file, // File (plain, Hatanaka .crx, or .gz)
+  (percent) => {}, // optional progress callback
+  undefined, // optional AbortSignal
+  (time, prn, codes, values) => {
+    // per-satellite observation callback
+  }
+);
+// result.header: RinexHeader
 
 // Navigation file parser
 const nav = parseNavFile(navFileText);
@@ -125,28 +129,36 @@ for (const frame of decoder) {
 import { computeAllPositions, computeDop } from 'gnss-js/orbit';
 
 const positions = computeAllPositions(ephemerides, times, receiverPosition);
-// positions.prns: string[], positions.times: number[], positions.az/el/lat/lon: Float64Array
+// positions.prns: string[], positions.times: number[]
+// positions.positions[prn][epochIdx]: { lat, lon, az, el } | null
 ```
 
 ### Signal analysis
 
 ```ts
-import {
-  MultipathAccumulator,
-  CycleSlipAccumulator,
-  CompletenessAccumulator,
-} from 'gnss-js/analysis';
-
-const mp = new MultipathAccumulator(header);
-const cs = new CycleSlipAccumulator(header);
-await parseRinexStream(stream, {
-  onObservation: (time, prn, codes, values) => {
-    mp.onObservation(time, prn, codes, values);
-    cs.onObservation(time, prn, codes, values);
-  },
-});
-const multipathResult = mp.finalize();
-const cycleSlipResult = cs.finalize();
+import { analyzeQuality } from 'gnss-js/analysis';
+
+// One re-parse pass: multipath (per code band, Anubis-style), cycle
+// slips, completeness, and slant-TEC ionosphere series.
+const q = await analyzeQuality(file, header);
+// q.multipath, q.cycleSlips, q.completeness, q.iono
+
+// The accumulators (MultipathAccumulator, CycleSlipAccumulator,
+// CompletenessAccumulator, IonoAccumulator) are also exported for
+// wiring into a custom parseRinexStream pass.
+```
+
+### Ionosphere and DCBs
+
+```ts
+import { parseSinexBiasDcb, applyIonoDcb } from 'gnss-js/analysis';
+
+// q.iono: slant TEC per satellite from the geometry-free phase
+// levelled to the code — DCB-biased until calibrated:
+const satDcb = parseSinexBiasDcb(sinexBiasText, obsEpochMs); // ESA .BIA / CAS .BSX
+const { result, receiverDcbTecu } = applyIonoDcb(q.iono, satDcb);
+// result.series[n].points: calibrated STEC (TECU); receiver bias
+// estimated from the night-time floor per system/signal pair.
 ```
 
 ### NTRIP client
@@ -173,6 +185,17 @@ const { reader, abort } = await connectToMountpoint(
 );
 ```
 
+### NMEA and ANTEX
+
+```ts
+import { parseNmeaFile, computeStats } from 'gnss-js/nmea';
+import { parseAntex } from 'gnss-js/antex';
+
+const track = parseNmeaFile(nmeaText); // fixes, satellites, per-sentence records
+const stats = computeStats(track.fixes); // position/HDOP/speed statistics
+const antex = parseAntex(antexText); // antenna PCO/PCV calibrations
+```
+
 ### Constants
 
 ```ts

package/dist/analysis.cjs

@@ -24,10 +24,24 @@ __export(analysis_exports, {
   CycleSlipAccumulator: () => CycleSlipAccumulator,
   IonoAccumulator: () => IonoAccumulator,
   MultipathAccumulator: () => MultipathAccumulator,
-  analyzeQuality: () => analyzeQuality
+  analyzeQuality: () => analyzeQuality,
+  applyIonoDcb: () => applyIonoDcb,
+  parseSinexBiasDcb: () => parseSinexBiasDcb
 });
 module.exports = __toCommonJS(analysis_exports);
 
+// src/analysis/stats-util.ts
+var MIN_ARC_LENGTH = 10;
+function median(values) {
+  const s = [...values].sort((a, b) => a - b);
+  const m = s.length >> 1;
+  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
+}
+function percentile(values, p) {
+  const s = [...values].sort((a, b) => a - b);
+  return s[Math.min(s.length - 1, Math.floor(p * s.length))] ?? 0;
+}
+
 // src/constants/gnss.ts
 var SYSTEM_NAMES = {
   G: "GPS",
@@ -215,15 +229,9 @@ function buildObsIndices(header) {
 var ARC_GAP_FACTOR = 5;
 
 // src/analysis/multipath.ts
-var MIN_ARC_LENGTH = 10;
 var MP_JUMP_M = 1.25;
 var MP_JUMP_WINDOW = 5;
 var MP_EDIT_SIGMA = 3;
-function median(values) {
-  const s = [...values].sort((a, b) => a - b);
-  const m = s.length >> 1;
-  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
-}
 var MultipathAccumulator = class {
   state = /* @__PURE__ */ new Map();
   results = [];
@@ -704,12 +712,6 @@ var CompletenessAccumulator = class {
 // src/analysis/ionosphere.ts
 var TEC_FACTOR = 403e15;
 var GF_JUMP_M = 0.15;
-var MIN_ARC_LENGTH2 = 10;
-function median2(values) {
-  const s = [...values].sort((a, b) => a - b);
-  const m = s.length >> 1;
-  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
-}
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
@@ -717,13 +719,15 @@ var IonoAccumulator = class {
   obsIndices;
   gloChannels;
   pairLabel = /* @__PURE__ */ new Map();
+  pairCodes = /* @__PURE__ */ new Map();
+  pairMeta = /* @__PURE__ */ new Map();
   constructor(header) {
     this.interval = header.interval ?? 30;
     this.obsIndices = buildObsIndices(header);
     this.gloChannels = buildGloChannelMap(header.glonassSlots);
   }
   /** Observation callback — wire this into parseRinexStream. */
-  onObservation = (time, prn, _codes, values) => {
+  onObservation = (time, prn, codes, values) => {
     const sys = prn[0];
     const bandMap = this.obsIndices.get(sys) ?? this.obsIndices.get("_v2");
     if (!bandMap) return;
@@ -734,7 +738,12 @@ var IonoAccumulator = class {
       const lVal = values[L];
       const freq = getFreq(this.gloChannels, prn, band);
       if (cVal != null && cVal !== 0 && lVal != null && lVal !== 0 && freq) {
-        bandData.set(band, { C: cVal, L: lVal, f: freq });
+        bandData.set(band, {
+          C: cVal,
+          L: lVal,
+          f: freq,
+          code: codes[C] ?? `C${band}`
+        });
       }
     }
     if (bandData.size < 2) return;
@@ -753,7 +762,9 @@ var IonoAccumulator = class {
         const bLabel = BAND_LABELS[sys]?.[bi] ?? bi;
         const rLabel = BAND_LABELS[sys]?.[bj] ?? bj;
         this.pairLabel.set(`${sys}:${pairKey}`, `${bLabel}-${rLabel}`);
+        this.pairCodes.set(`${sys}:${pairKey}`, [di.code, dj.code]);
       }
+      this.pairMeta.set(`${prn}:${pairKey}`, { fi: di.f, gamma });
       this.push(prn, pairKey, time, l4, p4, di.f, gamma);
       break;
     }
@@ -803,8 +814,8 @@ var IonoAccumulator = class {
   }
   closeArc(prn, pairKey, ps) {
     const arc = ps.arc;
-    if (arc.times.length >= MIN_ARC_LENGTH2) {
-      const level = median2(arc.l4.map((v, k) => v - arc.p4[k]));
+    if (arc.times.length >= MIN_ARC_LENGTH) {
+      const level = median(arc.l4.map((v, k) => v - arc.p4[k]));
       const toTecu = arc.fi * arc.fi / TEC_FACTOR / (arc.gamma - 1);
       let satArcs = this.closed.get(prn);
       if (!satArcs) {
@@ -835,7 +846,7 @@ var IonoAccumulator = class {
     const series = [];
     let sum = 0;
     let count = 0;
-    let maxStec = 0;
+    let maxStec = -Infinity;
     for (const [prn, satArcs] of this.closed) {
       let bestKey = null;
       let bestLen = 0;
@@ -854,18 +865,153 @@ var IonoAccumulator = class {
         count++;
         if (p.stec > maxStec) maxStec = p.stec;
       }
+      const meta = this.pairMeta.get(`${prn}:${bestKey}`);
+      const tecuPerNs = meta ? C_LIGHT * 1e-9 / (meta.gamma - 1) * (meta.fi * meta.fi / TEC_FACTOR) : 0;
       series.push({
         prn,
         system: sys,
         label: this.pairLabel.get(`${sys}:${bestKey}`) ?? bestKey,
+        codes: this.pairCodes.get(`${sys}:${bestKey}`) ?? ["", ""],
+        tecuPerNs,
         points
       });
     }
     series.sort((a, b) => a.prn.localeCompare(b.prn));
-    return { series, maxStec, meanStec: count > 0 ? sum / count : 0 };
+    return {
+      series,
+      maxStec: count > 0 ? maxStec : 0,
+      meanStec: count > 0 ? sum / count : 0
+    };
   }
 };
 
+// src/analysis/dcb.ts
+var PRN_RE = /^[A-Z]\d{2}$/;
+var SVN_RE = /^[A-Z]\d{3}$/;
+var OBS_RE = /^[CL]\d[A-Z]$/;
+function sinexEpochMs(s, openValue) {
+  const [y, d, sec] = s.split(":").map(Number);
+  if (!y) return openValue;
+  return Date.UTC(y, 0, 1) + ((d ?? 1) - 1) * 864e5 + (sec ?? 0) * 1e3;
+}
+function parseSinexBiasDcb(text, epochMs) {
+  const rows = /* @__PURE__ */ new Map();
+  for (const line of text.split("\n")) {
+    if (!line.startsWith(" DSB") && !line.startsWith("DSB")) continue;
+    const t = line.trim().split(/\s+/);
+    if (t.length < 9 || !SVN_RE.test(t[1]) || !PRN_RE.test(t[2]) || !OBS_RE.test(t[3]) || !OBS_RE.test(t[4])) {
+      continue;
+    }
+    const unitIdx = t.indexOf("ns");
+    if (unitIdx < 0 || unitIdx + 1 >= t.length) continue;
+    const value = parseFloat(t[unitIdx + 1]);
+    if (!isFinite(value)) continue;
+    const start = sinexEpochMs(t[5] ?? "", -Infinity);
+    const end = sinexEpochMs(t[6] ?? "", Infinity);
+    const covers = epochMs !== void 0 && epochMs >= start && epochMs < end;
+    const prn = t[2];
+    const pair = `${t[3]}-${t[4]}`;
+    let sat = rows.get(prn);
+    if (!sat) {
+      sat = /* @__PURE__ */ new Map();
+      rows.set(prn, sat);
+    }
+    const prev = sat.get(pair);
+    if (!prev || covers && !prev.covers || covers === prev.covers && end > prev.end) {
+      sat.set(pair, { value, end, covers });
+    }
+  }
+  const out = /* @__PURE__ */ new Map();
+  for (const [prn, sat] of rows) {
+    const m = /* @__PURE__ */ new Map();
+    for (const [pair, row] of sat) m.set(pair, row.value);
+    out.set(prn, m);
+  }
+  return out;
+}
+function applyIonoDcb(iono, satDcb) {
+  const missing = [];
+  let corrected = 0;
+  const dcbFor = (s) => {
+    const sat = satDcb.get(s.prn);
+    if (!sat) return null;
+    const fwd = sat.get(`${s.codes[0]}-${s.codes[1]}`);
+    if (fwd !== void 0) return fwd;
+    const rev = sat.get(`${s.codes[1]}-${s.codes[0]}`);
+    if (rev !== void 0) return -rev;
+    return null;
+  };
+  const series = iono.series.map((s) => {
+    const dcbNs = dcbFor(s);
+    if (dcbNs === null) {
+      missing.push(s.prn);
+      return s;
+    }
+    corrected++;
+    const shift = dcbNs * s.tecuPerNs;
+    return {
+      ...s,
+      points: s.points.map((p) => ({ time: p.time, stec: p.stec + shift }))
+    };
+  });
+  const missingSet = new Set(missing);
+  const receiverDcbTecu = {};
+  const groups = /* @__PURE__ */ new Map();
+  for (const s of series) {
+    const key = `${s.system} ${s.label}`;
+    let g = groups.get(key);
+    if (!g) {
+      g = [];
+      groups.set(key, g);
+    }
+    g.push(s);
+  }
+  const out = [];
+  for (const [key, group] of groups) {
+    const samples = [];
+    for (const s of group) {
+      if (missingSet.has(s.prn)) continue;
+      for (const p of s.points) samples.push(p.stec);
+    }
+    if (samples.length >= 100) {
+      const rx = percentile(samples, 0.01);
+      receiverDcbTecu[key] = rx;
+      for (const s of group) {
+        out.push({
+          ...s,
+          points: s.points.map((p) => ({
+            time: p.time,
+            stec: p.stec - rx
+          }))
+        });
+      }
+    } else {
+      out.push(...group);
+    }
+  }
+  out.sort((a, b) => a.prn.localeCompare(b.prn));
+  let sum = 0;
+  let count = 0;
+  let maxStec = -Infinity;
+  for (const s of out) {
+    for (const p of s.points) {
+      sum += p.stec;
+      count++;
+      if (p.stec > maxStec) maxStec = p.stec;
+    }
+  }
+  return {
+    result: {
+      series: out,
+      maxStec: count > 0 ? maxStec : 0,
+      meanStec: count > 0 ? sum / count : 0
+    },
+    satellitesCorrected: corrected,
+    satellitesMissing: [...new Set(missing)].sort(),
+    receiverDcbTecu
+  };
+}
+
 // src/rinex/crx.ts
 function crxRepair(old, diff) {
   const chars = old.split("");
@@ -1687,5 +1833,7 @@ async function analyzeQuality(file, header, onProgress, signal) {
   CycleSlipAccumulator,
   IonoAccumulator,
   MultipathAccumulator,
-  analyzeQuality
+  analyzeQuality,
+  applyIonoDcb,
+  parseSinexBiasDcb
 });

package/dist/analysis.d.cts

@@ -26,7 +26,7 @@ interface MultipathSeries {
     band: string;
     /** Reference band digit used in dual-frequency combination, e.g. "2". */
     refBand: string;
-    /** Human-readable label, e.g. "MP1 L1-L2". */
+    /** Human-readable label, e.g. "G01 MP L1-L2". */
     label: string;
     /** Time series of multipath values. */
     points: MultipathPoint[];
@@ -34,7 +34,7 @@ interface MultipathSeries {
     rms: number;
 }
 interface MultipathSignalStat {
-    /** Human-readable label, e.g. "MP1 L1-L2 (GPS)". */
+    /** Human-readable label, e.g. "MP L1-L2 (GPS)". */
     label: string;
     /** System letter, e.g. "G". */
     system: string;
@@ -231,6 +231,12 @@ interface IonoSeries {
     system: string;
     /** Band pair used, e.g. "L1-L2". */
     label: string;
+    /** Observation codes of the pair, e.g. ["C1C", "C2W"] — the key for
+     *  matching differential code bias products. */
+    codes: [string, string];
+    /** TECU shift caused by 1 ns of geometry-free code bias (B_i − B_j);
+     *  multiply a product DCB in ns by this to correct the series. */
+    tecuPerNs: number;
     /** Time series of slant TEC values (TECU, DCB-biased). */
     points: IonoPoint[];
 }
@@ -249,9 +255,11 @@ declare class IonoAccumulator {
     private obsIndices;
     private gloChannels;
     private pairLabel;
+    private pairCodes;
+    private pairMeta;
     constructor(header: RinexHeader);
     /** Observation callback — wire this into parseRinexStream. */
-    onObservation: (time: number, prn: string, _codes: string[], values: (number | null)[]) => void;
+    onObservation: (time: number, prn: string, codes: string[], values: (number | null)[]) => void;
     private push;
     /** External cycle-slip notification — close affected arcs. */
     notifySlip(_time: number, prn: string, bands: Set<string>): void;
@@ -260,6 +268,48 @@ declare class IonoAccumulator {
     finalize(): IonoResult;
 }
 
+/**
+ * Differential code bias (DCB) products and their application to the
+ * slant-TEC series from the ionosphere module.
+ *
+ * Satellite DCBs come from SINEX_BIAS files published by the IGS
+ * analysis centres (e.g. CAS `CAS0MGXRAP_*_DCB.BSX`). The levelled
+ * geometry-free combination carries B_i − B_j (satellite) plus the
+ * receiver bias; subtracting the product value per satellite and then
+ * estimating the receiver bias from the night-time floor (the levelled
+ * combination is low-noise and TEC is near zero at night — Hans van
+ * der Marel's "no negative values" criterion) yields calibrated STEC.
+ */
+
+/** PRN → "C1C-C2W" → bias (ns), satellite entries only. */
+type SatDcbMap = Map<string, Map<string, number>>;
+/**
+ * Parse satellite DSB entries from a SINEX_BIAS file (ESA .BIA,
+ * CAS/GFZ .BSX). Station entries (which carry a site name between PRN
+ * and the codes) are skipped. Long-history files publish several
+ * validity windows per satellite/pair — with values differing by many
+ * ns across SVN swaps — so pass the observation epoch to select the
+ * covering window; without it, the latest window wins. Values are
+ * returned in ns as published.
+ */
+declare function parseSinexBiasDcb(text: string, epochMs?: number): SatDcbMap;
+interface IonoDcbResult {
+    result: IonoResult;
+    /** Satellites whose product DCB was found and applied. */
+    satellitesCorrected: number;
+    /** PRNs with no matching product entry (left sat-DCB-biased). */
+    satellitesMissing: string[];
+    /** Estimated receiver bias per "system label" group (TECU). */
+    receiverDcbTecu: Record<string, number>;
+}
+/**
+ * Apply satellite DCBs from a product file to an ionosphere result and
+ * estimate the receiver bias per system/pair group from the night-time
+ * floor (1st percentile → 0). Series without a product entry get the
+ * satellite part left in; they are reported in `satellitesMissing`.
+ */
+declare function applyIonoDcb(iono: IonoResult, satDcb: SatDcbMap): IonoDcbResult;
+
 /**
  * Combined quality analysis — runs cycle slip detection, data completeness,
  * multipath and ionosphere analysis in a single file re-parse pass.
@@ -273,4 +323,4 @@ interface QualityResult {
 }
 declare function analyzeQuality(file: File, header: RinexHeader, onProgress?: (percent: number) => void, signal?: AbortSignal): Promise<QualityResult>;
 
-export { type CompleteSatSignal, type CompleteSignalStats, CompletenessAccumulator, type CompletenessResult, CycleSlipAccumulator, type CycleSlipEvent, type CycleSlipResult, type CycleSlipSignalStats, IonoAccumulator, type IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, analyzeQuality };
+export { type CompleteSatSignal, type CompleteSignalStats, CompletenessAccumulator, type CompletenessResult, CycleSlipAccumulator, type CycleSlipEvent, type CycleSlipResult, type CycleSlipSignalStats, IonoAccumulator, type IonoDcbResult, type IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, type SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb };

package/dist/analysis.d.ts

@@ -26,7 +26,7 @@ interface MultipathSeries {
     band: string;
     /** Reference band digit used in dual-frequency combination, e.g. "2". */
     refBand: string;
-    /** Human-readable label, e.g. "MP1 L1-L2". */
+    /** Human-readable label, e.g. "G01 MP L1-L2". */
     label: string;
     /** Time series of multipath values. */
     points: MultipathPoint[];
@@ -34,7 +34,7 @@ interface MultipathSeries {
     rms: number;
 }
 interface MultipathSignalStat {
-    /** Human-readable label, e.g. "MP1 L1-L2 (GPS)". */
+    /** Human-readable label, e.g. "MP L1-L2 (GPS)". */
     label: string;
     /** System letter, e.g. "G". */
     system: string;
@@ -231,6 +231,12 @@ interface IonoSeries {
     system: string;
     /** Band pair used, e.g. "L1-L2". */
     label: string;
+    /** Observation codes of the pair, e.g. ["C1C", "C2W"] — the key for
+     *  matching differential code bias products. */
+    codes: [string, string];
+    /** TECU shift caused by 1 ns of geometry-free code bias (B_i − B_j);
+     *  multiply a product DCB in ns by this to correct the series. */
+    tecuPerNs: number;
     /** Time series of slant TEC values (TECU, DCB-biased). */
     points: IonoPoint[];
 }
@@ -249,9 +255,11 @@ declare class IonoAccumulator {
     private obsIndices;
     private gloChannels;
     private pairLabel;
+    private pairCodes;
+    private pairMeta;
     constructor(header: RinexHeader);
     /** Observation callback — wire this into parseRinexStream. */
-    onObservation: (time: number, prn: string, _codes: string[], values: (number | null)[]) => void;
+    onObservation: (time: number, prn: string, codes: string[], values: (number | null)[]) => void;
     private push;
     /** External cycle-slip notification — close affected arcs. */
     notifySlip(_time: number, prn: string, bands: Set<string>): void;
@@ -260,6 +268,48 @@ declare class IonoAccumulator {
     finalize(): IonoResult;
 }
 
+/**
+ * Differential code bias (DCB) products and their application to the
+ * slant-TEC series from the ionosphere module.
+ *
+ * Satellite DCBs come from SINEX_BIAS files published by the IGS
+ * analysis centres (e.g. CAS `CAS0MGXRAP_*_DCB.BSX`). The levelled
+ * geometry-free combination carries B_i − B_j (satellite) plus the
+ * receiver bias; subtracting the product value per satellite and then
+ * estimating the receiver bias from the night-time floor (the levelled
+ * combination is low-noise and TEC is near zero at night — Hans van
+ * der Marel's "no negative values" criterion) yields calibrated STEC.
+ */
+
+/** PRN → "C1C-C2W" → bias (ns), satellite entries only. */
+type SatDcbMap = Map<string, Map<string, number>>;
+/**
+ * Parse satellite DSB entries from a SINEX_BIAS file (ESA .BIA,
+ * CAS/GFZ .BSX). Station entries (which carry a site name between PRN
+ * and the codes) are skipped. Long-history files publish several
+ * validity windows per satellite/pair — with values differing by many
+ * ns across SVN swaps — so pass the observation epoch to select the
+ * covering window; without it, the latest window wins. Values are
+ * returned in ns as published.
+ */
+declare function parseSinexBiasDcb(text: string, epochMs?: number): SatDcbMap;
+interface IonoDcbResult {
+    result: IonoResult;
+    /** Satellites whose product DCB was found and applied. */
+    satellitesCorrected: number;
+    /** PRNs with no matching product entry (left sat-DCB-biased). */
+    satellitesMissing: string[];
+    /** Estimated receiver bias per "system label" group (TECU). */
+    receiverDcbTecu: Record<string, number>;
+}
+/**
+ * Apply satellite DCBs from a product file to an ionosphere result and
+ * estimate the receiver bias per system/pair group from the night-time
+ * floor (1st percentile → 0). Series without a product entry get the
+ * satellite part left in; they are reported in `satellitesMissing`.
+ */
+declare function applyIonoDcb(iono: IonoResult, satDcb: SatDcbMap): IonoDcbResult;
+
 /**
  * Combined quality analysis — runs cycle slip detection, data completeness,
  * multipath and ionosphere analysis in a single file re-parse pass.
@@ -273,4 +323,4 @@ interface QualityResult {
 }
 declare function analyzeQuality(file: File, header: RinexHeader, onProgress?: (percent: number) => void, signal?: AbortSignal): Promise<QualityResult>;
 
-export { type CompleteSatSignal, type CompleteSignalStats, CompletenessAccumulator, type CompletenessResult, CycleSlipAccumulator, type CycleSlipEvent, type CycleSlipResult, type CycleSlipSignalStats, IonoAccumulator, type IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, analyzeQuality };
+export { type CompleteSatSignal, type CompleteSignalStats, CompletenessAccumulator, type CompletenessResult, CycleSlipAccumulator, type CycleSlipEvent, type CycleSlipResult, type CycleSlipSignalStats, IonoAccumulator, type IonoDcbResult, type IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, type SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb };

package/dist/analysis.js

@@ -3,8 +3,10 @@ import {
   CycleSlipAccumulator,
   IonoAccumulator,
   MultipathAccumulator,
-  analyzeQuality
-} from "./chunk-UQ7D53R2.js";
+  analyzeQuality,
+  applyIonoDcb,
+  parseSinexBiasDcb
+} from "./chunk-GDGHTC3E.js";
 import "./chunk-3U5AX7PY.js";
 import "./chunk-FIEWO4J4.js";
 export {
@@ -12,5 +14,7 @@ export {
   CycleSlipAccumulator,
   IonoAccumulator,
   MultipathAccumulator,
-  analyzeQuality
+  analyzeQuality,
+  applyIonoDcb,
+  parseSinexBiasDcb
 };

package/dist/{chunk-PQYBTE42.js → chunk-CD5YSYCG.js}

@@ -219,20 +219,10 @@ function invert4x4(m) {
   return result;
 }
 function selectEphemeris(ephemerides, prn, timeMs) {
-  let best = null;
-  let bestDt = Infinity;
-  for (const eph of ephemerides) {
-    if (eph.prn !== prn) continue;
-    const ephTime = eph.tocDate.getTime();
-    const dt = Math.abs(timeMs - ephTime);
-    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
-    if (dt > maxAge) continue;
-    if (dt < bestDt) {
-      bestDt = dt;
-      best = eph;
-    }
-  }
-  return best;
+  return selectBest(
+    ephemerides.filter((e) => e.prn === prn),
+    timeMs
+  );
 }
 function computeSatPosition(eph, timeMs) {
   if (eph.system === "R" || eph.system === "S") {
@@ -325,8 +315,9 @@ function selectBest(ephs, timeMs) {
   return best;
 }
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var MS_PER_WEEK = 7 * 864e5;
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
-var GAL_EPOCH_MS = GPS_EPOCH_MS;
+var GAL_EPOCH_MS = GPS_EPOCH_MS + 1024 * MS_PER_WEEK;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -370,6 +361,7 @@ function ephInfoToEphemeris(info) {
     return null;
   }
   let epochMs;
+  let week = info.week;
   const tocSec = info.toc ?? info.toe;
   if (sys === "C") {
     epochMs = BDS_EPOCH_MS;
@@ -377,8 +369,10 @@ function ephInfoToEphemeris(info) {
     epochMs = GAL_EPOCH_MS;
   } else {
     epochMs = GPS_EPOCH_MS;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS) / MS_PER_WEEK;
+    week += 1024 * Math.round((weeksNow - week) / 1024);
   }
-  const tocDate = new Date(epochMs + info.week * 7 * 864e5 + tocSec * 1e3);
+  const tocDate = new Date(epochMs + week * MS_PER_WEEK + tocSec * 1e3);
   return {
     system: sys,
     prn: info.prn,
@@ -404,7 +398,7 @@ function ephInfoToEphemeris(info) {
     omega: info.argPerigee,
     omegaDot: info.omegaDot ?? 0,
     idot: info.idot ?? 0,
-    week: info.week,
+    week,
     svHealth: info.health,
     tgd: 0
   };

package/dist/{chunk-65CQINSB.js → chunk-FGRC3PB3.js}

@@ -1345,6 +1345,11 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1111 && msgType <= 1117) return "QZSS";
   if (msgType >= 1121 && msgType <= 1127) return "BeiDou";
   if (msgType >= 1131 && msgType <= 1137) return "NavIC";
+  if (msgType >= 1015 && msgType <= 1017) return "GPS";
+  if (msgType === 1030) return "GPS";
+  if (msgType === 1031) return "GLONASS";
+  if (msgType >= 1057 && msgType <= 1062) return "GPS";
+  if (msgType >= 1063 && msgType <= 1068) return "GLONASS";
   if (msgType === 1230) return "GLONASS";
   return null;
 }
@@ -1384,7 +1389,6 @@ function updateStreamStats(stats, frames, rawBytes) {
     entry.totalBytes += frame.length + 6;
     const msmEpoch = decodeMsmFull(frame);
     if (msmEpoch) {
-      const now2 = Date.now();
       for (const obs of msmEpoch.observations) {
         const signals = [];
         let bestCn0 = 0;
@@ -1399,7 +1403,7 @@ function updateStreamStats(stats, frames, rawBytes) {
             prn: obs.prn,
             system: obs.system,
             cn0: bestCn0,
-            lastSeen: now2,
+            lastSeen: now,
             signals
           });
         }