gnss-js 1.3.1 → 1.4.0

package/dist/analysis.cjs

@@ -22,6 +22,7 @@ var analysis_exports = {};
 __export(analysis_exports, {
   CompletenessAccumulator: () => CompletenessAccumulator,
   CycleSlipAccumulator: () => CycleSlipAccumulator,
+  IonoAccumulator: () => IonoAccumulator,
   MultipathAccumulator: () => MultipathAccumulator,
   analyzeQuality: () => analyzeQuality
 });
@@ -215,6 +216,14 @@ 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 = [];
@@ -271,6 +280,11 @@ var MultipathAccumulator = class {
     const gap = bandState.lastTime > 0 ? (time - bandState.lastTime) / 1e3 : 0;
     if (bandState.lastTime > 0 && gap > this.interval * ARC_GAP_FACTOR) {
       this.closeArc(prn, band, refBand, bandState);
+    } else if (bandState.arc.rawMp.length > 0) {
+      const recent = bandState.arc.rawMp.slice(-MP_JUMP_WINDOW);
+      if (Math.abs(mp - median(recent)) > MP_JUMP_M) {
+        this.closeArc(prn, band, refBand, bandState);
+      }
     }
     bandState.arc.times.push(time);
     bandState.arc.rawMp.push(mp);
@@ -290,13 +304,25 @@ var MultipathAccumulator = class {
   closeArc(prn, band, refBand, state) {
     const arc = state.arc;
     if (arc.times.length >= MIN_ARC_LENGTH) {
-      const mean = arc.rawMp.reduce((a, b) => a + b, 0) / arc.rawMp.length;
+      let keep = arc.times.map((_, i) => i);
+      let mean = 0;
+      for (let pass = 0; pass < 3; pass++) {
+        mean = keep.reduce((s, i) => s + arc.rawMp[i], 0) / keep.length;
+        const sigma = Math.sqrt(
+          keep.reduce((s, i) => s + (arc.rawMp[i] - mean) ** 2, 0) / keep.length
+        );
+        const kept = keep.filter(
+          (i) => Math.abs(arc.rawMp[i] - mean) <= MP_EDIT_SIGMA * sigma
+        );
+        if (kept.length === keep.length || kept.length < MIN_ARC_LENGTH) break;
+        keep = kept;
+      }
       const sys = prn[0];
       const bLabel = BAND_LABELS[sys]?.[band] ?? band;
       const rLabel = BAND_LABELS[sys]?.[refBand] ?? refBand;
       const label = `${prn} MP ${bLabel}-${rLabel}`;
-      const points = arc.times.map((t, i) => ({
-        time: t,
+      const points = keep.map((i) => ({
+        time: arc.times[i],
         mp: arc.rawMp[i] - mean
       }));
       const rms = Math.sqrt(
@@ -350,7 +376,7 @@ var MultipathAccumulator = class {
         refBand,
         rms,
         count,
-        satellites: seriesList.length
+        satellites: new Set(seriesList.map((s) => s.prn)).size
       });
     }
     const sysOrder = "GRECIJS";
@@ -668,6 +694,171 @@ 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();
+  interval;
+  obsIndices;
+  gloChannels;
+  pairLabel = /* @__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) => {
+    const sys = prn[0];
+    const bandMap = this.obsIndices.get(sys) ?? this.obsIndices.get("_v2");
+    if (!bandMap) return;
+    const bandData = /* @__PURE__ */ new Map();
+    for (const [band, { C, L }] of bandMap) {
+      if (C === null) continue;
+      const cVal = values[C];
+      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 });
+      }
+    }
+    if (bandData.size < 2) return;
+    const pairs = DUAL_FREQ_PAIRS[sys] ?? [];
+    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 gamma = di.f * di.f / (dj.f * dj.f);
+      const l4 = di.L * \u03BBi - dj.L * \u03BBj;
+      const p4 = dj.C - di.C;
+      const pairKey = `${bi}-${bj}`;
+      if (!this.pairLabel.has(`${sys}:${pairKey}`)) {
+        const bLabel = BAND_LABELS[sys]?.[bi] ?? bi;
+        const rLabel = BAND_LABELS[sys]?.[bj] ?? bj;
+        this.pairLabel.set(`${sys}:${pairKey}`, `${bLabel}-${rLabel}`);
+      }
+      this.push(prn, pairKey, time, l4, p4, di.f, gamma);
+      break;
+    }
+  };
+  push(prn, pairKey, time, l4, p4, fi, gamma) {
+    if (!isFinite(l4) || !isFinite(p4)) return;
+    let satStates = this.state.get(prn);
+    if (!satStates) {
+      satStates = /* @__PURE__ */ new Map();
+      this.state.set(prn, satStates);
+    }
+    let ps = satStates.get(pairKey);
+    if (!ps) {
+      ps = {
+        arc: { times: [], l4: [], p4: [], fi, gamma },
+        lastTime: 0,
+        lastL4: 0
+      };
+      satStates.set(pairKey, ps);
+    }
+    const gap = ps.lastTime > 0 ? (time - ps.lastTime) / 1e3 : 0;
+    if (ps.lastTime > 0 && gap > this.interval * ARC_GAP_FACTOR) {
+      this.closeArc(prn, pairKey, ps);
+    } else if (ps.arc.times.length > 0 && Math.abs(l4 - ps.lastL4) > GF_JUMP_M) {
+      this.closeArc(prn, pairKey, ps);
+    }
+    if (ps.arc.times.length === 0) {
+      ps.arc.fi = fi;
+      ps.arc.gamma = gamma;
+    }
+    ps.arc.times.push(time);
+    ps.arc.l4.push(l4);
+    ps.arc.p4.push(p4);
+    ps.lastTime = time;
+    ps.lastL4 = l4;
+  }
+  /** External cycle-slip notification — close affected arcs. */
+  notifySlip(_time, prn, bands) {
+    const satStates = this.state.get(prn);
+    if (!satStates) return;
+    for (const [pairKey, ps] of satStates) {
+      const [bi, bj] = pairKey.split("-");
+      if (bands.has(bi) || bands.has(bj)) {
+        this.closeArc(prn, pairKey, ps);
+      }
+    }
+  }
+  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]));
+      const toTecu = arc.fi * arc.fi / TEC_FACTOR / (arc.gamma - 1);
+      let satArcs = this.closed.get(prn);
+      if (!satArcs) {
+        satArcs = /* @__PURE__ */ new Map();
+        this.closed.set(prn, satArcs);
+      }
+      let points = satArcs.get(pairKey);
+      if (!points) {
+        points = [];
+        satArcs.set(pairKey, points);
+      }
+      for (let k = 0; k < arc.times.length; k++) {
+        points.push({
+          time: arc.times[k],
+          stec: (arc.l4[k] - level) * toTecu
+        });
+      }
+    }
+    ps.arc = { times: [], l4: [], p4: [], fi: arc.fi, gamma: arc.gamma };
+  }
+  /** Finalize: close remaining arcs, keep one pair per satellite. */
+  finalize() {
+    for (const [prn, satStates] of this.state) {
+      for (const [pairKey, ps] of satStates) {
+        this.closeArc(prn, pairKey, ps);
+      }
+    }
+    const series = [];
+    let sum = 0;
+    let count = 0;
+    let maxStec = 0;
+    for (const [prn, satArcs] of this.closed) {
+      let bestKey = null;
+      let bestLen = 0;
+      for (const [pairKey, points2] of satArcs) {
+        if (points2.length > bestLen) {
+          bestLen = points2.length;
+          bestKey = pairKey;
+        }
+      }
+      if (!bestKey) continue;
+      const points = satArcs.get(bestKey);
+      points.sort((a, b) => a.time - b.time);
+      const sys = prn[0];
+      for (const p of points) {
+        sum += p.stec;
+        count++;
+        if (p.stec > maxStec) maxStec = p.stec;
+      }
+      series.push({
+        prn,
+        system: sys,
+        label: this.pairLabel.get(`${sys}:${bestKey}`) ?? bestKey,
+        points
+      });
+    }
+    series.sort((a, b) => a.prn.localeCompare(b.prn));
+    return { series, maxStec, meanStec: count > 0 ? sum / count : 0 };
+  }
+};
+
 // src/rinex/crx.ts
 function crxRepair(old, diff) {
   const chars = old.split("");
@@ -1457,8 +1648,10 @@ function computeStats(header, epochs, satellitesSeen) {
 // src/analysis/quality-analysis.ts
 async function analyzeQuality(file, header, onProgress, signal) {
   const mpAccum = new MultipathAccumulator(header);
+  const ionoAccum = new IonoAccumulator(header);
   const csAccum = new CycleSlipAccumulator(header, (time, prn, bands) => {
     mpAccum.notifySlip(time, prn, bands);
+    ionoAccum.notifySlip(time, prn, bands);
   });
   const compAccum = new CompletenessAccumulator(header);
   await parseRinexStream(
@@ -1468,6 +1661,7 @@ async function analyzeQuality(file, header, onProgress, signal) {
     (time, prn, codes, values) => {
       csAccum.onObservation(time, prn, codes, values);
       mpAccum.onObservation(time, prn, codes, values);
+      ionoAccum.onObservation(time, prn, codes, values);
       compAccum.onObservation(time, prn, codes, values);
     },
     true
@@ -1476,13 +1670,15 @@ async function analyzeQuality(file, header, onProgress, signal) {
   return {
     cycleSlips: csAccum.finalize(),
     completeness: compAccum.finalize(),
-    multipath: mpAccum.finalize()
+    multipath: mpAccum.finalize(),
+    iono: ionoAccum.finalize()
   };
 }
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
   CompletenessAccumulator,
   CycleSlipAccumulator,
+  IonoAccumulator,
   MultipathAccumulator,
   analyzeQuality
 });

package/dist/analysis.d.cts

@@ -1,5 +1,5 @@
 import { R as RinexHeader } from './parser-JPjjFgeP.cjs';
-import { O as OnSlipDetected } from './gnss-BT6ulR17.cjs';
+import { i as OnSlipDetected } from './gnss-BtXrG3zH.cjs';
 
 /**
  * Multipath analysis from RINEX dual-frequency observations.
@@ -196,16 +196,81 @@ declare class CompletenessAccumulator {
     finalize(): CompletenessResult;
 }
 
+/**
+ * Slant ionospheric delay from dual-frequency observations.
+ *
+ * Geometry-free phase combination levelled to the geometry-free code,
+ * per continuous arc:
+ *   L4 = λi·Li − λj·Lj = (γ−1)·I_i + B      (phase; B = ambiguity bias)
+ *   P4 = C_j − C_i    = (γ−1)·I_i + DCBs    (code)
+ * with γ = (f_i/f_j)² and I_i the slant ionospheric delay on band i.
+ * The arc levelling constant B̂ = median(L4 − P4) transfers the code's
+ * absolute level onto the low-noise phase; slant TEC follows as
+ *   STEC [TECU] = (L4 − B̂)/(γ−1) · f_i² / 40.3e16.
+ *
+ * The result carries the receiver+satellite differential code biases
+ * (not removed here — several TECU per satellite), so series are
+ * DCB-biased but shape-faithful: diurnal variation, gradients and
+ * disturbances read correctly.
+ *
+ * The geometry-free phase is also a sensitive cycle-slip detector:
+ * arcs split on epoch-to-epoch L4 jumps beyond GF_JUMP_M, in addition
+ * to time gaps and external slip notifications.
+ */
+
+interface IonoPoint {
+    /** Epoch time in milliseconds. */
+    time: number;
+    /** Slant TEC in TECU (DCB-biased). */
+    stec: number;
+}
+interface IonoSeries {
+    /** Satellite PRN, e.g. "G01". */
+    prn: string;
+    /** System letter, e.g. "G". */
+    system: string;
+    /** Band pair used, e.g. "L1-L2". */
+    label: string;
+    /** Time series of slant TEC values (TECU, DCB-biased). */
+    points: IonoPoint[];
+}
+interface IonoResult {
+    /** Per-satellite slant TEC time series. */
+    series: IonoSeries[];
+    /** Highest single STEC sample across all series (TECU). */
+    maxStec: number;
+    /** Mean STEC across all samples (TECU). */
+    meanStec: number;
+}
+declare class IonoAccumulator {
+    private state;
+    private closed;
+    private interval;
+    private obsIndices;
+    private gloChannels;
+    private pairLabel;
+    constructor(header: RinexHeader);
+    /** Observation callback — wire this into parseRinexStream. */
+    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;
+    private closeArc;
+    /** Finalize: close remaining arcs, keep one pair per satellite. */
+    finalize(): IonoResult;
+}
+
 /**
  * Combined quality analysis — runs cycle slip detection, data completeness,
- * and multipath analysis in a single file re-parse pass.
+ * multipath and ionosphere analysis in a single file re-parse pass.
  */
 
 interface QualityResult {
     cycleSlips: CycleSlipResult;
     completeness: CompletenessResult;
     multipath: MultipathResult;
+    iono: IonoResult;
 }
 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, 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 IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, analyzeQuality };

package/dist/analysis.d.ts

@@ -1,5 +1,5 @@
 import { R as RinexHeader } from './parser-JPjjFgeP.js';
-import { O as OnSlipDetected } from './gnss-C-tgoYNa.js';
+import { i as OnSlipDetected } from './gnss-DhnEr0Dy.js';
 
 /**
  * Multipath analysis from RINEX dual-frequency observations.
@@ -196,16 +196,81 @@ declare class CompletenessAccumulator {
     finalize(): CompletenessResult;
 }
 
+/**
+ * Slant ionospheric delay from dual-frequency observations.
+ *
+ * Geometry-free phase combination levelled to the geometry-free code,
+ * per continuous arc:
+ *   L4 = λi·Li − λj·Lj = (γ−1)·I_i + B      (phase; B = ambiguity bias)
+ *   P4 = C_j − C_i    = (γ−1)·I_i + DCBs    (code)
+ * with γ = (f_i/f_j)² and I_i the slant ionospheric delay on band i.
+ * The arc levelling constant B̂ = median(L4 − P4) transfers the code's
+ * absolute level onto the low-noise phase; slant TEC follows as
+ *   STEC [TECU] = (L4 − B̂)/(γ−1) · f_i² / 40.3e16.
+ *
+ * The result carries the receiver+satellite differential code biases
+ * (not removed here — several TECU per satellite), so series are
+ * DCB-biased but shape-faithful: diurnal variation, gradients and
+ * disturbances read correctly.
+ *
+ * The geometry-free phase is also a sensitive cycle-slip detector:
+ * arcs split on epoch-to-epoch L4 jumps beyond GF_JUMP_M, in addition
+ * to time gaps and external slip notifications.
+ */
+
+interface IonoPoint {
+    /** Epoch time in milliseconds. */
+    time: number;
+    /** Slant TEC in TECU (DCB-biased). */
+    stec: number;
+}
+interface IonoSeries {
+    /** Satellite PRN, e.g. "G01". */
+    prn: string;
+    /** System letter, e.g. "G". */
+    system: string;
+    /** Band pair used, e.g. "L1-L2". */
+    label: string;
+    /** Time series of slant TEC values (TECU, DCB-biased). */
+    points: IonoPoint[];
+}
+interface IonoResult {
+    /** Per-satellite slant TEC time series. */
+    series: IonoSeries[];
+    /** Highest single STEC sample across all series (TECU). */
+    maxStec: number;
+    /** Mean STEC across all samples (TECU). */
+    meanStec: number;
+}
+declare class IonoAccumulator {
+    private state;
+    private closed;
+    private interval;
+    private obsIndices;
+    private gloChannels;
+    private pairLabel;
+    constructor(header: RinexHeader);
+    /** Observation callback — wire this into parseRinexStream. */
+    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;
+    private closeArc;
+    /** Finalize: close remaining arcs, keep one pair per satellite. */
+    finalize(): IonoResult;
+}
+
 /**
  * Combined quality analysis — runs cycle slip detection, data completeness,
- * and multipath analysis in a single file re-parse pass.
+ * multipath and ionosphere analysis in a single file re-parse pass.
  */
 
 interface QualityResult {
     cycleSlips: CycleSlipResult;
     completeness: CompletenessResult;
     multipath: MultipathResult;
+    iono: IonoResult;
 }
 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, 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 IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, analyzeQuality };

package/dist/analysis.js

@@ -1,14 +1,16 @@
 import {
   CompletenessAccumulator,
   CycleSlipAccumulator,
+  IonoAccumulator,
   MultipathAccumulator,
   analyzeQuality
-} from "./chunk-G3N4S3DM.js";
-import "./chunk-OZCYOM5D.js";
-import "./chunk-W5WKEV7U.js";
+} from "./chunk-3ZAFWOT6.js";
+import "./chunk-3U5AX7PY.js";
+import "./chunk-FIEWO4J4.js";
 export {
   CompletenessAccumulator,
   CycleSlipAccumulator,
+  IonoAccumulator,
   MultipathAccumulator,
   analyzeQuality
 };

package/dist/{chunk-5SPJH4MG.js → chunk-2K3FCJBX.js}

@@ -4,7 +4,7 @@ import {
   GLO_F1_STEP,
   GLO_F2_BASE,
   GLO_F2_STEP
-} from "./chunk-W5WKEV7U.js";
+} from "./chunk-FIEWO4J4.js";
 
 // src/signals/definitions.ts
 var MHz = 1e6;

package/dist/{chunk-OZCYOM5D.js → chunk-3U5AX7PY.js}

@@ -1,6 +1,6 @@
 import {
   SYSTEM_NAMES
-} from "./chunk-W5WKEV7U.js";
+} from "./chunk-FIEWO4J4.js";
 
 // src/rinex/crx.ts
 function crxRepair(old, diff) {