gnss-js 1.5.1 → 1.6.1

package/dist/analysis.cjs

@@ -26,6 +26,8 @@ __export(analysis_exports, {
   MultipathAccumulator: () => MultipathAccumulator,
   analyzeQuality: () => analyzeQuality,
   applyIonoDcb: () => applyIonoDcb,
+  computeIonoRate: () => computeIonoRate,
+  detrendIonoArcs: () => detrendIonoArcs,
   parseSinexBiasDcb: () => parseSinexBiasDcb
 });
 module.exports = __toCommonJS(analysis_exports);
@@ -715,6 +717,7 @@ var GF_JUMP_M = 0.15;
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
+  closedArcStarts = /* @__PURE__ */ new Map();
   interval;
   obsIndices;
   gloChannels;
@@ -827,6 +830,17 @@ var IonoAccumulator = class {
         points = [];
         satArcs.set(pairKey, points);
       }
+      let satStarts = this.closedArcStarts.get(prn);
+      if (!satStarts) {
+        satStarts = /* @__PURE__ */ new Map();
+        this.closedArcStarts.set(prn, satStarts);
+      }
+      let starts = satStarts.get(pairKey);
+      if (!starts) {
+        starts = [];
+        satStarts.set(pairKey, starts);
+      }
+      starts.push(points.length);
       for (let k = 0; k < arc.times.length; k++) {
         points.push({
           time: arc.times[k],
@@ -858,7 +872,7 @@ var IonoAccumulator = class {
       }
       if (!bestKey) continue;
       const points = satArcs.get(bestKey);
-      points.sort((a, b) => a.time - b.time);
+      const arcStarts = this.closedArcStarts.get(prn)?.get(bestKey) ?? [0];
       const sys = prn[0];
       for (const p of points) {
         sum += p.stec;
@@ -873,7 +887,8 @@ var IonoAccumulator = class {
         label: this.pairLabel.get(`${sys}:${bestKey}`) ?? bestKey,
         codes: this.pairCodes.get(`${sys}:${bestKey}`) ?? ["", ""],
         tecuPerNs,
-        points
+        points,
+        arcStarts
       });
     }
     series.sort((a, b) => a.prn.localeCompare(b.prn));
@@ -884,6 +899,87 @@ var IonoAccumulator = class {
     };
   }
 };
+function computeIonoRate(result, intervalSec, order = 1) {
+  const out = [];
+  for (const s of result.series) {
+    const pts = [];
+    const bounds = [...s.arcStarts, s.points.length];
+    for (let a = 0; a < bounds.length - 1; a++) {
+      const lo = bounds[a];
+      const hi = bounds[a + 1];
+      if (order === 2) {
+        for (let i = lo + 2; i < hi; i++) {
+          const dt1 = s.points[i - 1].time - s.points[i - 2].time;
+          const dt2 = s.points[i].time - s.points[i - 1].time;
+          if (Math.abs(dt1 - dt2) > 0.1 * Math.max(dt1, dt2)) continue;
+          pts.push({
+            time: s.points[i].time,
+            value: s.points[i].stec - 2 * s.points[i - 1].stec + s.points[i - 2].stec
+          });
+        }
+        continue;
+      }
+      if (intervalSec === void 0) {
+        for (let i = lo + 1; i < hi; i++) {
+          const dtMin = (s.points[i].time - s.points[i - 1].time) / 6e4;
+          if (dtMin <= 0) continue;
+          pts.push({
+            time: s.points[i].time,
+            value: (s.points[i].stec - s.points[i - 1].stec) / dtMin
+          });
+        }
+      } else {
+        const targetMs = intervalSec * 1e3;
+        let j = lo;
+        for (let i = lo; i < hi; i++) {
+          if (j <= i) j = i + 1;
+          while (j < hi && s.points[j].time - s.points[i].time < targetMs)
+            j++;
+          if (j >= hi) break;
+          const dtMs = s.points[j].time - s.points[i].time;
+          if (dtMs > 1.5 * targetMs) continue;
+          pts.push({
+            time: s.points[j].time,
+            value: (s.points[j].stec - s.points[i].stec) / dtMs * 6e4
+          });
+        }
+      }
+    }
+    if (pts.length > 0) out.push({ prn: s.prn, system: s.system, points: pts });
+  }
+  return out;
+}
+function detrendIonoArcs(result) {
+  let sum = 0;
+  let count = 0;
+  let maxStec = -Infinity;
+  const series = result.series.map((s) => {
+    const bounds = [...s.arcStarts, s.points.length];
+    const points = new Array(s.points.length);
+    for (let a = 0; a < bounds.length - 1; a++) {
+      const lo = bounds[a];
+      const hi = bounds[a + 1];
+      const anchor = s.points[lo]?.stec ?? 0;
+      for (let i = lo; i < hi; i++) {
+        points[i] = {
+          time: s.points[i].time,
+          stec: s.points[i].stec - anchor
+        };
+      }
+    }
+    for (const p of points) {
+      sum += p.stec;
+      count++;
+      if (p.stec > maxStec) maxStec = p.stec;
+    }
+    return { ...s, points };
+  });
+  return {
+    series,
+    maxStec: count > 0 ? maxStec : 0,
+    meanStec: count > 0 ? sum / count : 0
+  };
+}
 
 // src/analysis/dcb.ts
 var PRN_RE = /^[A-Z]\d{2}$/;
@@ -1835,5 +1931,7 @@ async function analyzeQuality(file, header, onProgress, signal) {
   MultipathAccumulator,
   analyzeQuality,
   applyIonoDcb,
+  computeIonoRate,
+  detrendIonoArcs,
   parseSinexBiasDcb
 });

package/dist/analysis.d.cts

@@ -239,6 +239,12 @@ interface IonoSeries {
     tecuPerNs: number;
     /** Time series of slant TEC values (TECU, DCB-biased). */
     points: IonoPoint[];
+    /**
+     * Indices into `points` where a new continuous arc begins (gap,
+     * cycle slip, or geometry-free jump). Differencing and per-arc
+     * detrending must not cross these boundaries.
+     */
+    arcStarts: number[];
 }
 interface IonoResult {
     /** Per-satellite slant TEC time series. */
@@ -251,6 +257,7 @@ interface IonoResult {
 declare class IonoAccumulator {
     private state;
     private closed;
+    private closedArcStarts;
     private interval;
     private obsIndices;
     private gloChannels;
@@ -267,6 +274,39 @@ declare class IonoAccumulator {
     /** Finalize: close remaining arcs, keep one pair per satellite. */
     finalize(): IonoResult;
 }
+interface IonoRatePoint {
+    /** Epoch time (ms) — the later epoch of the differenced pair. */
+    time: number;
+    /** ROT in TECU/min (order 1) or undivided ΔΔSTEC in TECU (order 2). */
+    value: number;
+}
+interface IonoRateSeries {
+    prn: string;
+    system: string;
+    points: IonoRatePoint[];
+}
+/**
+ * Sequential time differences of the slant TEC series — the biases
+ * (ambiguities, DCBs) cancel, leaving ionospheric rate of TEC, phase
+ * noise, and scintillation. Differences never cross arc boundaries,
+ * so cycle slips do not appear as outliers (arcs split there).
+ *
+ * @param intervalSec Standardized differencing baseline in seconds.
+ *   Omit for native-rate sequential differences; set e.g. 60 for a
+ *   sample-rate-independent picture (all pairs ~1 min apart are used).
+ * @param order 1 (default) = first difference in TECU/min;
+ *   2 = second undivided difference in TECU (native rate only,
+ *   gradients removed, noise and scintillation amplified).
+ */
+declare function computeIonoRate(result: IonoResult, intervalSec?: number, order?: 1 | 2): IonoRateSeries[];
+/**
+ * Remove the per-arc bias from a slant TEC result — each arc is
+ * shifted so its first observation reads zero, leaving only the TEC
+ * *variation* along the arc (Hans van der Marel's suggestion; using
+ * the first observation keeps the shape untouched at the small risk
+ * of anchoring on an outlier).
+ */
+declare function detrendIonoArcs(result: IonoResult): IonoResult;
 
 /**
  * Differential code bias (DCB) products and their application to the
@@ -323,4 +363,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 IonoDcbResult, type IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, type SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb };
+export { type CompleteSatSignal, type CompleteSignalStats, CompletenessAccumulator, type CompletenessResult, CycleSlipAccumulator, type CycleSlipEvent, type CycleSlipResult, type CycleSlipSignalStats, IonoAccumulator, type IonoDcbResult, type IonoPoint, type IonoRatePoint, type IonoRateSeries, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, type SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb };

package/dist/analysis.d.ts

@@ -239,6 +239,12 @@ interface IonoSeries {
     tecuPerNs: number;
     /** Time series of slant TEC values (TECU, DCB-biased). */
     points: IonoPoint[];
+    /**
+     * Indices into `points` where a new continuous arc begins (gap,
+     * cycle slip, or geometry-free jump). Differencing and per-arc
+     * detrending must not cross these boundaries.
+     */
+    arcStarts: number[];
 }
 interface IonoResult {
     /** Per-satellite slant TEC time series. */
@@ -251,6 +257,7 @@ interface IonoResult {
 declare class IonoAccumulator {
     private state;
     private closed;
+    private closedArcStarts;
     private interval;
     private obsIndices;
     private gloChannels;
@@ -267,6 +274,39 @@ declare class IonoAccumulator {
     /** Finalize: close remaining arcs, keep one pair per satellite. */
     finalize(): IonoResult;
 }
+interface IonoRatePoint {
+    /** Epoch time (ms) — the later epoch of the differenced pair. */
+    time: number;
+    /** ROT in TECU/min (order 1) or undivided ΔΔSTEC in TECU (order 2). */
+    value: number;
+}
+interface IonoRateSeries {
+    prn: string;
+    system: string;
+    points: IonoRatePoint[];
+}
+/**
+ * Sequential time differences of the slant TEC series — the biases
+ * (ambiguities, DCBs) cancel, leaving ionospheric rate of TEC, phase
+ * noise, and scintillation. Differences never cross arc boundaries,
+ * so cycle slips do not appear as outliers (arcs split there).
+ *
+ * @param intervalSec Standardized differencing baseline in seconds.
+ *   Omit for native-rate sequential differences; set e.g. 60 for a
+ *   sample-rate-independent picture (all pairs ~1 min apart are used).
+ * @param order 1 (default) = first difference in TECU/min;
+ *   2 = second undivided difference in TECU (native rate only,
+ *   gradients removed, noise and scintillation amplified).
+ */
+declare function computeIonoRate(result: IonoResult, intervalSec?: number, order?: 1 | 2): IonoRateSeries[];
+/**
+ * Remove the per-arc bias from a slant TEC result — each arc is
+ * shifted so its first observation reads zero, leaving only the TEC
+ * *variation* along the arc (Hans van der Marel's suggestion; using
+ * the first observation keeps the shape untouched at the small risk
+ * of anchoring on an outlier).
+ */
+declare function detrendIonoArcs(result: IonoResult): IonoResult;
 
 /**
  * Differential code bias (DCB) products and their application to the
@@ -323,4 +363,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 IonoDcbResult, type IonoPoint, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, type SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb };
+export { type CompleteSatSignal, type CompleteSignalStats, CompletenessAccumulator, type CompletenessResult, CycleSlipAccumulator, type CycleSlipEvent, type CycleSlipResult, type CycleSlipSignalStats, IonoAccumulator, type IonoDcbResult, type IonoPoint, type IonoRatePoint, type IonoRateSeries, type IonoResult, type IonoSeries, MultipathAccumulator, type MultipathPoint, type MultipathResult, type MultipathSeries, type MultipathSignalStat, type QualityResult, type SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb };

package/dist/analysis.js

@@ -5,8 +5,10 @@ import {
   MultipathAccumulator,
   analyzeQuality,
   applyIonoDcb,
+  computeIonoRate,
+  detrendIonoArcs,
   parseSinexBiasDcb
-} from "./chunk-GDGHTC3E.js";
+} from "./chunk-Z5UZUNCB.js";
 import "./chunk-3U5AX7PY.js";
 import "./chunk-FIEWO4J4.js";
 export {
@@ -16,5 +18,7 @@ export {
   MultipathAccumulator,
   analyzeQuality,
   applyIonoDcb,
+  computeIonoRate,
+  detrendIonoArcs,
   parseSinexBiasDcb
 };

package/dist/{chunk-SOJNWF5T.js → chunk-4PXLRDK6.js}

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

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

@@ -98,7 +98,7 @@ function gloDerivatives(state, acc) {
   const r3 = r2 * r;
   const r5 = r2 * r3;
   const mu_r3 = GM_GLO / r3;
-  const j2_term = 1.5 * J2_GLO * AE_GLO * AE_GLO / r5;
+  const j2_term = 1.5 * J2_GLO * GM_GLO * AE_GLO * AE_GLO / r5;
   const dvx = -mu_r3 * x + j2_term * x * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * x + 2 * OMEGA_E * vy + ax;
   const dvy = -mu_r3 * y + j2_term * y * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * y - 2 * OMEGA_E * vx + ay;
   const dvz = -mu_r3 * z + j2_term * z * (5 * z * z / r2 - 3) + az;

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

@@ -512,6 +512,7 @@ var GF_JUMP_M = 0.15;
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
+  closedArcStarts = /* @__PURE__ */ new Map();
   interval;
   obsIndices;
   gloChannels;
@@ -624,6 +625,17 @@ var IonoAccumulator = class {
         points = [];
         satArcs.set(pairKey, points);
       }
+      let satStarts = this.closedArcStarts.get(prn);
+      if (!satStarts) {
+        satStarts = /* @__PURE__ */ new Map();
+        this.closedArcStarts.set(prn, satStarts);
+      }
+      let starts = satStarts.get(pairKey);
+      if (!starts) {
+        starts = [];
+        satStarts.set(pairKey, starts);
+      }
+      starts.push(points.length);
       for (let k = 0; k < arc.times.length; k++) {
         points.push({
           time: arc.times[k],
@@ -655,7 +667,7 @@ var IonoAccumulator = class {
       }
       if (!bestKey) continue;
       const points = satArcs.get(bestKey);
-      points.sort((a, b) => a.time - b.time);
+      const arcStarts = this.closedArcStarts.get(prn)?.get(bestKey) ?? [0];
       const sys = prn[0];
       for (const p of points) {
         sum += p.stec;
@@ -670,7 +682,8 @@ var IonoAccumulator = class {
         label: this.pairLabel.get(`${sys}:${bestKey}`) ?? bestKey,
         codes: this.pairCodes.get(`${sys}:${bestKey}`) ?? ["", ""],
         tecuPerNs,
-        points
+        points,
+        arcStarts
       });
     }
     series.sort((a, b) => a.prn.localeCompare(b.prn));
@@ -681,6 +694,87 @@ var IonoAccumulator = class {
     };
   }
 };
+function computeIonoRate(result, intervalSec, order = 1) {
+  const out = [];
+  for (const s of result.series) {
+    const pts = [];
+    const bounds = [...s.arcStarts, s.points.length];
+    for (let a = 0; a < bounds.length - 1; a++) {
+      const lo = bounds[a];
+      const hi = bounds[a + 1];
+      if (order === 2) {
+        for (let i = lo + 2; i < hi; i++) {
+          const dt1 = s.points[i - 1].time - s.points[i - 2].time;
+          const dt2 = s.points[i].time - s.points[i - 1].time;
+          if (Math.abs(dt1 - dt2) > 0.1 * Math.max(dt1, dt2)) continue;
+          pts.push({
+            time: s.points[i].time,
+            value: s.points[i].stec - 2 * s.points[i - 1].stec + s.points[i - 2].stec
+          });
+        }
+        continue;
+      }
+      if (intervalSec === void 0) {
+        for (let i = lo + 1; i < hi; i++) {
+          const dtMin = (s.points[i].time - s.points[i - 1].time) / 6e4;
+          if (dtMin <= 0) continue;
+          pts.push({
+            time: s.points[i].time,
+            value: (s.points[i].stec - s.points[i - 1].stec) / dtMin
+          });
+        }
+      } else {
+        const targetMs = intervalSec * 1e3;
+        let j = lo;
+        for (let i = lo; i < hi; i++) {
+          if (j <= i) j = i + 1;
+          while (j < hi && s.points[j].time - s.points[i].time < targetMs)
+            j++;
+          if (j >= hi) break;
+          const dtMs = s.points[j].time - s.points[i].time;
+          if (dtMs > 1.5 * targetMs) continue;
+          pts.push({
+            time: s.points[j].time,
+            value: (s.points[j].stec - s.points[i].stec) / dtMs * 6e4
+          });
+        }
+      }
+    }
+    if (pts.length > 0) out.push({ prn: s.prn, system: s.system, points: pts });
+  }
+  return out;
+}
+function detrendIonoArcs(result) {
+  let sum = 0;
+  let count = 0;
+  let maxStec = -Infinity;
+  const series = result.series.map((s) => {
+    const bounds = [...s.arcStarts, s.points.length];
+    const points = new Array(s.points.length);
+    for (let a = 0; a < bounds.length - 1; a++) {
+      const lo = bounds[a];
+      const hi = bounds[a + 1];
+      const anchor = s.points[lo]?.stec ?? 0;
+      for (let i = lo; i < hi; i++) {
+        points[i] = {
+          time: s.points[i].time,
+          stec: s.points[i].stec - anchor
+        };
+      }
+    }
+    for (const p of points) {
+      sum += p.stec;
+      count++;
+      if (p.stec > maxStec) maxStec = p.stec;
+    }
+    return { ...s, points };
+  });
+  return {
+    series,
+    maxStec: count > 0 ? maxStec : 0,
+    meanStec: count > 0 ? sum / count : 0
+  };
+}
 
 // src/analysis/dcb.ts
 var PRN_RE = /^[A-Z]\d{2}$/;
@@ -844,6 +938,8 @@ export {
   CycleSlipAccumulator,
   CompletenessAccumulator,
   IonoAccumulator,
+  computeIonoRate,
+  detrendIonoArcs,
   parseSinexBiasDcb,
   applyIonoDcb,
   analyzeQuality

package/dist/index.cjs

@@ -111,6 +111,7 @@ __export(src_exports, {
   clampUnit: () => clampUnit,
   computeAllPositions: () => computeAllPositions,
   computeDop: () => computeDop,
+  computeIonoRate: () => computeIonoRate,
   computeLiveSkyPositions: () => computeLiveSkyPositions,
   computePsdDb: () => computePsdDb,
   computeSatPosition: () => computeSatPosition,
@@ -126,6 +127,7 @@ __export(src_exports, {
   decodeMsmFull: () => decodeMsmFull,
   deg2dms: () => deg2dms,
   deg2rad: () => deg2rad,
+  detrendIonoArcs: () => detrendIonoArcs,
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
   ephInfoToEphemeris: () => ephInfoToEphemeris,
@@ -4086,7 +4088,7 @@ function gloDerivatives(state, acc) {
   const r3 = r2 * r;
   const r5 = r2 * r3;
   const mu_r3 = GM_GLO / r3;
-  const j2_term = 1.5 * J2_GLO * AE_GLO * AE_GLO / r5;
+  const j2_term = 1.5 * J2_GLO * GM_GLO * AE_GLO * AE_GLO / r5;
   const dvx = -mu_r3 * x + j2_term * x * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * x + 2 * OMEGA_E * vy + ax;
   const dvy = -mu_r3 * y + j2_term * y * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * y - 2 * OMEGA_E * vx + ay;
   const dvz = -mu_r3 * z + j2_term * z * (5 * z * z / r2 - 3) + az;
@@ -5541,6 +5543,7 @@ var GF_JUMP_M = 0.15;
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
+  closedArcStarts = /* @__PURE__ */ new Map();
   interval;
   obsIndices;
   gloChannels;
@@ -5653,6 +5656,17 @@ var IonoAccumulator = class {
         points = [];
         satArcs.set(pairKey, points);
       }
+      let satStarts = this.closedArcStarts.get(prn);
+      if (!satStarts) {
+        satStarts = /* @__PURE__ */ new Map();
+        this.closedArcStarts.set(prn, satStarts);
+      }
+      let starts = satStarts.get(pairKey);
+      if (!starts) {
+        starts = [];
+        satStarts.set(pairKey, starts);
+      }
+      starts.push(points.length);
       for (let k = 0; k < arc.times.length; k++) {
         points.push({
           time: arc.times[k],
@@ -5684,7 +5698,7 @@ var IonoAccumulator = class {
       }
       if (!bestKey) continue;
       const points = satArcs.get(bestKey);
-      points.sort((a, b) => a.time - b.time);
+      const arcStarts = this.closedArcStarts.get(prn)?.get(bestKey) ?? [0];
       const sys = prn[0];
       for (const p of points) {
         sum += p.stec;
@@ -5699,7 +5713,8 @@ var IonoAccumulator = class {
         label: this.pairLabel.get(`${sys}:${bestKey}`) ?? bestKey,
         codes: this.pairCodes.get(`${sys}:${bestKey}`) ?? ["", ""],
         tecuPerNs,
-        points
+        points,
+        arcStarts
       });
     }
     series.sort((a, b) => a.prn.localeCompare(b.prn));
@@ -5710,6 +5725,87 @@ var IonoAccumulator = class {
     };
   }
 };
+function computeIonoRate(result, intervalSec, order = 1) {
+  const out = [];
+  for (const s of result.series) {
+    const pts = [];
+    const bounds = [...s.arcStarts, s.points.length];
+    for (let a = 0; a < bounds.length - 1; a++) {
+      const lo = bounds[a];
+      const hi = bounds[a + 1];
+      if (order === 2) {
+        for (let i = lo + 2; i < hi; i++) {
+          const dt1 = s.points[i - 1].time - s.points[i - 2].time;
+          const dt2 = s.points[i].time - s.points[i - 1].time;
+          if (Math.abs(dt1 - dt2) > 0.1 * Math.max(dt1, dt2)) continue;
+          pts.push({
+            time: s.points[i].time,
+            value: s.points[i].stec - 2 * s.points[i - 1].stec + s.points[i - 2].stec
+          });
+        }
+        continue;
+      }
+      if (intervalSec === void 0) {
+        for (let i = lo + 1; i < hi; i++) {
+          const dtMin = (s.points[i].time - s.points[i - 1].time) / 6e4;
+          if (dtMin <= 0) continue;
+          pts.push({
+            time: s.points[i].time,
+            value: (s.points[i].stec - s.points[i - 1].stec) / dtMin
+          });
+        }
+      } else {
+        const targetMs = intervalSec * 1e3;
+        let j = lo;
+        for (let i = lo; i < hi; i++) {
+          if (j <= i) j = i + 1;
+          while (j < hi && s.points[j].time - s.points[i].time < targetMs)
+            j++;
+          if (j >= hi) break;
+          const dtMs = s.points[j].time - s.points[i].time;
+          if (dtMs > 1.5 * targetMs) continue;
+          pts.push({
+            time: s.points[j].time,
+            value: (s.points[j].stec - s.points[i].stec) / dtMs * 6e4
+          });
+        }
+      }
+    }
+    if (pts.length > 0) out.push({ prn: s.prn, system: s.system, points: pts });
+  }
+  return out;
+}
+function detrendIonoArcs(result) {
+  let sum = 0;
+  let count = 0;
+  let maxStec = -Infinity;
+  const series = result.series.map((s) => {
+    const bounds = [...s.arcStarts, s.points.length];
+    const points = new Array(s.points.length);
+    for (let a = 0; a < bounds.length - 1; a++) {
+      const lo = bounds[a];
+      const hi = bounds[a + 1];
+      const anchor = s.points[lo]?.stec ?? 0;
+      for (let i = lo; i < hi; i++) {
+        points[i] = {
+          time: s.points[i].time,
+          stec: s.points[i].stec - anchor
+        };
+      }
+    }
+    for (const p of points) {
+      sum += p.stec;
+      count++;
+      if (p.stec > maxStec) maxStec = p.stec;
+    }
+    return { ...s, points };
+  });
+  return {
+    series,
+    maxStec: count > 0 ? maxStec : 0,
+    meanStec: count > 0 ? sum / count : 0
+  };
+}
 
 // src/analysis/dcb.ts
 var PRN_RE = /^[A-Z]\d{2}$/;
@@ -7150,6 +7246,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   clampUnit,
   computeAllPositions,
   computeDop,
+  computeIonoRate,
   computeLiveSkyPositions,
   computePsdDb,
   computeSatPosition,
@@ -7165,6 +7262,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   decodeMsmFull,
   deg2dms,
   deg2rad,
+  detrendIonoArcs,
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,

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, IonoDcbResult, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb } from './analysis.cjs';
+export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, 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, IonoDcbResult, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, parseSinexBiasDcb } from './analysis.js';
+export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, 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

@@ -10,8 +10,10 @@ import {
   MultipathAccumulator,
   analyzeQuality,
   applyIonoDcb,
+  computeIonoRate,
+  detrendIonoArcs,
   parseSinexBiasDcb
-} from "./chunk-GDGHTC3E.js";
+} from "./chunk-Z5UZUNCB.js";
 import {
   frequencyLabel,
   parseAntex
@@ -100,7 +102,7 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-SOJNWF5T.js";
+} from "./chunk-4PXLRDK6.js";
 import {
   computeAllPositions,
   computeDop,
@@ -113,7 +115,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-CD5YSYCG.js";
+} from "./chunk-JG2HDJGK.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -340,6 +342,7 @@ export {
   clampUnit,
   computeAllPositions,
   computeDop,
+  computeIonoRate,
   computeLiveSkyPositions,
   computePsdDb,
   computeSatPosition,
@@ -355,6 +358,7 @@ export {
   decodeMsmFull,
   deg2dms,
   deg2rad,
+  detrendIonoArcs,
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,

package/dist/orbit.cjs

@@ -251,7 +251,7 @@ function gloDerivatives(state, acc) {
   const r3 = r2 * r;
   const r5 = r2 * r3;
   const mu_r3 = GM_GLO / r3;
-  const j2_term = 1.5 * J2_GLO * AE_GLO * AE_GLO / r5;
+  const j2_term = 1.5 * J2_GLO * GM_GLO * AE_GLO * AE_GLO / r5;
   const dvx = -mu_r3 * x + j2_term * x * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * x + 2 * OMEGA_E * vy + ax;
   const dvy = -mu_r3 * y + j2_term * y * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * y - 2 * OMEGA_E * vx + ay;
   const dvz = -mu_r3 * z + j2_term * z * (5 * z * z / r2 - 3) + az;

package/dist/orbit.js

@@ -10,7 +10,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-CD5YSYCG.js";
+} from "./chunk-JG2HDJGK.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,

package/dist/positioning.cjs

@@ -235,7 +235,7 @@ function gloDerivatives(state, acc) {
   const r3 = r2 * r;
   const r5 = r2 * r3;
   const mu_r3 = GM_GLO / r3;
-  const j2_term = 1.5 * J2_GLO * AE_GLO * AE_GLO / r5;
+  const j2_term = 1.5 * J2_GLO * GM_GLO * AE_GLO * AE_GLO / r5;
   const dvx = -mu_r3 * x + j2_term * x * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * x + 2 * OMEGA_E * vy + ax;
   const dvy = -mu_r3 * y + j2_term * y * (5 * z * z / r2 - 1) + OMEGA_E * OMEGA_E * y - 2 * OMEGA_E * vx + ay;
   const dvz = -mu_r3 * z + j2_term * z * (5 * z * z / r2 - 3) + az;

package/dist/positioning.js

@@ -2,8 +2,8 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-SOJNWF5T.js";
-import "./chunk-CD5YSYCG.js";
+} from "./chunk-4PXLRDK6.js";
+import "./chunk-JG2HDJGK.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.5.1",
+  "version": "1.6.1",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,
@@ -118,14 +118,14 @@
   "devDependencies": {
     "@eslint/js": "^10.0.1",
     "@types/node": "^20.11.17",
-    "eslint": "^10.1.0",
+    "eslint": "^10.7.0",
     "husky": "^9.1.7",
     "lint-staged": "^16.4.0",
-    "prettier": "^3.5.3",
+    "prettier": "^3.9.5",
     "tsup": "^8.0.2",
     "typescript": "^5.8.3",
-    "typescript-eslint": "^8.57.2",
-    "vitest": "^4.0.0"
+    "typescript-eslint": "^8.64.0",
+    "vitest": "^4.1.10"
   },
   "lint-staged": {
     "*.{ts,js}": [