gnss-js 1.4.1 → 1.5.1

package/dist/{chunk-UQ7D53R2.js → chunk-GDGHTC3E.js}

@@ -13,16 +13,22 @@ import {
   getFreq
 } from "./chunk-FIEWO4J4.js";
 
-// src/analysis/multipath.ts
+// src/analysis/stats-util.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;
 }
+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/analysis/multipath.ts
+var MP_JUMP_M = 1.25;
+var MP_JUMP_WINDOW = 5;
+var MP_EDIT_SIGMA = 3;
 var MultipathAccumulator = class {
   state = /* @__PURE__ */ new Map();
   results = [];
@@ -503,12 +509,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();
@@ -516,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;
@@ -533,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;
@@ -552,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;
     }
@@ -602,8 +611,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) {
@@ -634,7 +643,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;
@@ -653,18 +662,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/analysis/quality-analysis.ts
 async function analyzeQuality(file, header, onProgress, signal) {
   const mpAccum = new MultipathAccumulator(header);
@@ -700,5 +844,7 @@ export {
   CycleSlipAccumulator,
   CompletenessAccumulator,
   IonoAccumulator,
+  parseSinexBiasDcb,
+  applyIonoDcb,
   analyzeQuality
 };

package/dist/{chunk-PX7TPPTQ.js → chunk-SOJNWF5T.js}

@@ -2,7 +2,7 @@ import {
   computeDop,
   computeSatPosition,
   ecefToAzEl
-} from "./chunk-PQYBTE42.js";
+} from "./chunk-CD5YSYCG.js";
 import {
   C_LIGHT,
   OMEGA_E

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,
@@ -3904,6 +3906,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;
 }
@@ -3943,7 +3950,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;
@@ -3958,7 +3964,7 @@ function updateStreamStats(stats, frames, rawBytes) {
             prn: obs.prn,
             system: obs.system,
             cn0: bestCn0,
-            lastSeen: now2,
+            lastSeen: now,
             signals
           });
         }
@@ -4201,20 +4207,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") {
@@ -4307,8 +4303,9 @@ function selectBest(ephs, timeMs) {
   return best;
 }
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var MS_PER_WEEK2 = 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_WEEK2;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -4352,6 +4349,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;
@@ -4359,8 +4357,10 @@ function ephInfoToEphemeris(info) {
     epochMs = GAL_EPOCH_MS;
   } else {
     epochMs = GPS_EPOCH_MS;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS) / MS_PER_WEEK2;
+    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_WEEK2 + tocSec * 1e3);
   return {
     system: sys,
     prn: info.prn,
@@ -4386,7 +4386,7 @@ function ephInfoToEphemeris(info) {
     omega: info.argPerigee,
     omegaDot: info.omegaDot ?? 0,
     idot: info.idot ?? 0,
-    week: info.week,
+    week,
     svHealth: info.health,
     tgd: 0
   };
@@ -5042,16 +5042,22 @@ async function connectToMountpoint(proxyUrl, info, signal) {
   };
 }
 
-// src/analysis/multipath.ts
+// src/analysis/stats-util.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;
 }
+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/analysis/multipath.ts
+var MP_JUMP_M = 1.25;
+var MP_JUMP_WINDOW = 5;
+var MP_EDIT_SIGMA = 3;
 var MultipathAccumulator = class {
   state = /* @__PURE__ */ new Map();
   results = [];
@@ -5532,12 +5538,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();
@@ -5545,13 +5545,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;
@@ -5562,7 +5564,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;
@@ -5581,7 +5588,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;
     }
@@ -5631,8 +5640,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) {
@@ -5663,7 +5672,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;
@@ -5682,18 +5691,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/analysis/quality-analysis.ts
 async function analyzeQuality(file, header, onProgress, signal) {
   const mpAccum = new MultipathAccumulator(header);
@@ -7000,6 +7144,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   WarningAccumulator,
   analyzeQuality,
   applyHelmert,
+  applyIonoDcb,
   buildGloChannelMap,
   buildObsIndices,
   clampUnit,
@@ -7102,6 +7247,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';