gnss-js 1.4.0 → 1.5.0

package/dist/analysis.cjs

@@ -24,7 +24,9 @@ __export(analysis_exports, {
   CycleSlipAccumulator: () => CycleSlipAccumulator,
   IonoAccumulator: () => IonoAccumulator,
   MultipathAccumulator: () => MultipathAccumulator,
-  analyzeQuality: () => analyzeQuality
+  analyzeQuality: () => analyzeQuality,
+  applyIonoDcb: () => applyIonoDcb,
+  parseSinexBiasDcb: () => parseSinexBiasDcb
 });
 module.exports = __toCommonJS(analysis_exports);
 
@@ -251,17 +253,24 @@ var MultipathAccumulator = class {
       }
     }
     if (bandData.size < 2) return;
+    const mpOf = (a, b) => {
+      const \u03BBa = C_LIGHT / a.f;
+      const \u03BBb = C_LIGHT / b.f;
+      const \u03B1 = a.f * a.f / (b.f * b.f);
+      const coeff = 2 / (\u03B1 - 1);
+      return a.C - (1 + coeff) * a.L * \u03BBa + coeff * b.L * \u03BBb;
+    };
     const pairs = DUAL_FREQ_PAIRS[sys] ?? [];
+    let primaryDone = false;
     for (const [bi, bj] of pairs) {
       const di = bandData.get(bi);
       const dj = bandData.get(bj);
       if (!di || !dj) continue;
-      const \u03BBi = C_LIGHT / di.f;
-      const \u03BBj = C_LIGHT / dj.f;
-      const \u03B1 = di.f * di.f / (dj.f * dj.f);
-      const coeff = 2 / (\u03B1 - 1);
-      const mp = di.C - (1 + coeff) * di.L * \u03BBi + coeff * dj.L * \u03BBj;
-      this.pushMp(prn, bi, bj, time, mp);
+      if (!primaryDone) {
+        this.pushMp(prn, bi, bj, time, mpOf(di, dj));
+        primaryDone = true;
+      }
+      this.pushMp(prn, bj, bi, time, mpOf(dj, di));
     }
   };
   pushMp(prn, band, refBand, time, mp) {
@@ -710,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;
@@ -727,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;
@@ -746,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;
     }
@@ -847,10 +865,14 @@ 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
       });
     }
@@ -859,6 +881,133 @@ var IonoAccumulator = class {
   }
 };
 
+// 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) {
+  const [y, d, sec] = s.split(":").map(Number);
+  if (!y) return Infinity;
+  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] ?? "");
+    const end = sinexEpochMs(t[6] ?? "");
+    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 percentile(values, p) {
+  const s = [...values].sort((a, b) => a - b);
+  return s[Math.min(s.length - 1, Math.floor(p * s.length))] ?? 0;
+}
+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 = 0;
+  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, 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("");
@@ -1680,5 +1829,7 @@ async function analyzeQuality(file, header, onProgress, signal) {
   CycleSlipAccumulator,
   IonoAccumulator,
   MultipathAccumulator,
-  analyzeQuality
+  analyzeQuality,
+  applyIonoDcb,
+  parseSinexBiasDcb
 });

package/dist/analysis.d.cts

@@ -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

@@ -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-3ZAFWOT6.js";
+  analyzeQuality,
+  applyIonoDcb,
+  parseSinexBiasDcb
+} from "./chunk-EJQ6NG5U.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-3ZAFWOT6.js → chunk-EJQ6NG5U.js}

@@ -50,17 +50,24 @@ var MultipathAccumulator = class {
       }
     }
     if (bandData.size < 2) return;
+    const mpOf = (a, b) => {
+      const \u03BBa = C_LIGHT / a.f;
+      const \u03BBb = C_LIGHT / b.f;
+      const \u03B1 = a.f * a.f / (b.f * b.f);
+      const coeff = 2 / (\u03B1 - 1);
+      return a.C - (1 + coeff) * a.L * \u03BBa + coeff * b.L * \u03BBb;
+    };
     const pairs = DUAL_FREQ_PAIRS[sys] ?? [];
+    let primaryDone = false;
     for (const [bi, bj] of pairs) {
       const di = bandData.get(bi);
       const dj = bandData.get(bj);
       if (!di || !dj) continue;
-      const \u03BBi = C_LIGHT / di.f;
-      const \u03BBj = C_LIGHT / dj.f;
-      const \u03B1 = di.f * di.f / (dj.f * dj.f);
-      const coeff = 2 / (\u03B1 - 1);
-      const mp = di.C - (1 + coeff) * di.L * \u03BBi + coeff * dj.L * \u03BBj;
-      this.pushMp(prn, bi, bj, time, mp);
+      if (!primaryDone) {
+        this.pushMp(prn, bi, bj, time, mpOf(di, dj));
+        primaryDone = true;
+      }
+      this.pushMp(prn, bj, bi, time, mpOf(dj, di));
     }
   };
   pushMp(prn, band, refBand, time, mp) {
@@ -509,13 +516,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;
@@ -526,7 +535,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;
@@ -545,7 +559,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;
     }
@@ -646,10 +662,14 @@ 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
       });
     }
@@ -658,6 +678,133 @@ var IonoAccumulator = class {
   }
 };
 
+// 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) {
+  const [y, d, sec] = s.split(":").map(Number);
+  if (!y) return Infinity;
+  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] ?? "");
+    const end = sinexEpochMs(t[6] ?? "");
+    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 percentile(values, p) {
+  const s = [...values].sort((a, b) => a - b);
+  return s[Math.min(s.length - 1, Math.floor(p * s.length))] ?? 0;
+}
+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 = 0;
+  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, meanStec: count > 0 ? sum / count : 0 },
+    satellitesCorrected: corrected,
+    satellitesMissing: [...new Set(missing)].sort(),
+    receiverDcbTecu
+  };
+}
+
 // src/analysis/quality-analysis.ts
 async function analyzeQuality(file, header, onProgress, signal) {
   const mpAccum = new MultipathAccumulator(header);
@@ -693,5 +840,7 @@ export {
   CycleSlipAccumulator,
   CompletenessAccumulator,
   IonoAccumulator,
+  parseSinexBiasDcb,
+  applyIonoDcb,
   analyzeQuality
 };

package/dist/index.cjs

@@ -105,6 +105,7 @@ __export(src_exports, {
   WarningAccumulator: () => WarningAccumulator,
   analyzeQuality: () => analyzeQuality,
   applyHelmert: () => applyHelmert,
+  applyIonoDcb: () => applyIonoDcb,
   buildGloChannelMap: () => buildGloChannelMap,
   buildObsIndices: () => buildObsIndices,
   clampUnit: () => clampUnit,
@@ -207,6 +208,7 @@ __export(src_exports, {
   parseNavFile: () => parseNavFile,
   parseNmeaFile: () => parseNmeaFile,
   parseRinexStream: () => parseRinexStream,
+  parseSinexBiasDcb: () => parseSinexBiasDcb,
   parseSourcetable: () => parseSourcetable,
   phiAltBOC: () => phiAltBOC,
   phiBOCc: () => phiBOCc,
@@ -5079,17 +5081,24 @@ var MultipathAccumulator = class {
       }
     }
     if (bandData.size < 2) return;
+    const mpOf = (a, b) => {
+      const \u03BBa = C_LIGHT / a.f;
+      const \u03BBb = C_LIGHT / b.f;
+      const \u03B1 = a.f * a.f / (b.f * b.f);
+      const coeff = 2 / (\u03B1 - 1);
+      return a.C - (1 + coeff) * a.L * \u03BBa + coeff * b.L * \u03BBb;
+    };
     const pairs = DUAL_FREQ_PAIRS[sys] ?? [];
+    let primaryDone = false;
     for (const [bi, bj] of pairs) {
       const di = bandData.get(bi);
       const dj = bandData.get(bj);
       if (!di || !dj) continue;
-      const \u03BBi = C_LIGHT / di.f;
-      const \u03BBj = C_LIGHT / dj.f;
-      const \u03B1 = di.f * di.f / (dj.f * dj.f);
-      const coeff = 2 / (\u03B1 - 1);
-      const mp = di.C - (1 + coeff) * di.L * \u03BBi + coeff * dj.L * \u03BBj;
-      this.pushMp(prn, bi, bj, time, mp);
+      if (!primaryDone) {
+        this.pushMp(prn, bi, bj, time, mpOf(di, dj));
+        primaryDone = true;
+      }
+      this.pushMp(prn, bj, bi, time, mpOf(dj, di));
     }
   };
   pushMp(prn, band, refBand, time, mp) {
@@ -5538,13 +5547,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;
@@ -5555,7 +5566,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;
@@ -5574,7 +5590,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;
     }
@@ -5675,10 +5693,14 @@ 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
       });
     }
@@ -5687,6 +5709,133 @@ var IonoAccumulator = class {
   }
 };
 
+// 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) {
+  const [y, d, sec] = s.split(":").map(Number);
+  if (!y) return Infinity;
+  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] ?? "");
+    const end = sinexEpochMs(t[6] ?? "");
+    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 percentile(values, p) {
+  const s = [...values].sort((a, b) => a - b);
+  return s[Math.min(s.length - 1, Math.floor(p * s.length))] ?? 0;
+}
+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 = 0;
+  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, meanStec: count > 0 ? sum / count : 0 },
+    satellitesCorrected: corrected,
+    satellitesMissing: [...new Set(missing)].sort(),
+    receiverDcbTecu
+  };
+}
+
 // src/analysis/quality-analysis.ts
 async function analyzeQuality(file, header, onProgress, signal) {
   const mpAccum = new MultipathAccumulator(header);
@@ -6993,6 +7142,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   WarningAccumulator,
   analyzeQuality,
   applyHelmert,
+  applyIonoDcb,
   buildGloChannelMap,
   buildObsIndices,
   clampUnit,
@@ -7095,6 +7245,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   parseNavFile,
   parseNmeaFile,
   parseRinexStream,
+  parseSinexBiasDcb,
   parseSourcetable,
   phiAltBOC,
   phiBOCc,

package/dist/index.d.cts

@@ -13,7 +13,7 @@ export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPositi
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
 export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } from './positioning.cjs';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';
-export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.cjs';
+export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';
 export { NmeaFix, NmeaStats, NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum } from './nmea.cjs';
 export { BAND_PRESETS, BandPreset, CONSTELLATIONS, ConstellationRow, DELTA_GLO_L1, DELTA_GLO_L2, FREQ_BDS_B1A, FREQ_BDS_B1C, FREQ_BDS_B1I, FREQ_BDS_B2I, FREQ_BDS_B3A, FREQ_BDS_B3I, FREQ_GAL_E1, FREQ_GAL_E5, FREQ_GAL_E5a, FREQ_GAL_E5b, FREQ_GAL_E6, FREQ_GLO_L1, FREQ_GLO_L1OC, FREQ_GLO_L2, FREQ_GLO_L2OC, FREQ_GLO_L3OC, FREQ_GPS_L1, FREQ_GPS_L2, FREQ_GPS_L5, FREQ_NAVIC_L5, FREQ_QZS_L1, FREQ_QZS_L2, FREQ_QZS_L5, FREQ_QZS_L6, Modulation, SignalDef, computePsdDb, phiAltBOC, phiBOCc, phiBOCs, phiBPSK } from './signals.cjs';

package/dist/index.d.ts

@@ -13,7 +13,7 @@ export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPositi
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
 export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } from './positioning.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';
-export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.js';
+export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';
 export { NmeaFix, NmeaStats, NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum } from './nmea.js';
 export { BAND_PRESETS, BandPreset, CONSTELLATIONS, ConstellationRow, DELTA_GLO_L1, DELTA_GLO_L2, FREQ_BDS_B1A, FREQ_BDS_B1C, FREQ_BDS_B1I, FREQ_BDS_B2I, FREQ_BDS_B3A, FREQ_BDS_B3I, FREQ_GAL_E1, FREQ_GAL_E5, FREQ_GAL_E5a, FREQ_GAL_E5b, FREQ_GAL_E6, FREQ_GLO_L1, FREQ_GLO_L1OC, FREQ_GLO_L2, FREQ_GLO_L2OC, FREQ_GLO_L3OC, FREQ_GPS_L1, FREQ_GPS_L2, FREQ_GPS_L5, FREQ_NAVIC_L5, FREQ_QZS_L1, FREQ_QZS_L2, FREQ_QZS_L5, FREQ_QZS_L6, Modulation, SignalDef, computePsdDb, phiAltBOC, phiBOCc, phiBOCs, phiBPSK } from './signals.js';

package/dist/index.js

@@ -8,8 +8,10 @@ import {
   CycleSlipAccumulator,
   IonoAccumulator,
   MultipathAccumulator,
-  analyzeQuality
-} from "./chunk-3ZAFWOT6.js";
+  analyzeQuality,
+  applyIonoDcb,
+  parseSinexBiasDcb
+} from "./chunk-EJQ6NG5U.js";
 import {
   frequencyLabel,
   parseAntex
@@ -332,6 +334,7 @@ export {
   WarningAccumulator,
   analyzeQuality,
   applyHelmert,
+  applyIonoDcb,
   buildGloChannelMap,
   buildObsIndices,
   clampUnit,
@@ -434,6 +437,7 @@ export {
   parseNavFile,
   parseNmeaFile,
   parseRinexStream,
+  parseSinexBiasDcb,
   parseSourcetable,
   phiAltBOC,
   phiBOCc,

package/package.json

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