gnss-js 1.5.0 → 1.6.0

package/dist/{chunk-EJQ6NG5U.js → chunk-Z5UZUNCB.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,15 +509,10 @@ 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();
+  closedArcStarts = /* @__PURE__ */ new Map();
   interval;
   obsIndices;
   gloChannels;
@@ -611,8 +612,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) {
@@ -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],
@@ -643,7 +655,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;
@@ -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,21 +682,107 @@ 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));
-    return { series, maxStec, meanStec: count > 0 ? sum / count : 0 };
+    return {
+      series,
+      maxStec: count > 0 ? maxStec : 0,
+      meanStec: count > 0 ? sum / count : 0
+    };
   }
 };
+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}$/;
 var SVN_RE = /^[A-Z]\d{3}$/;
 var OBS_RE = /^[CL]\d[A-Z]$/;
-function sinexEpochMs(s) {
+function sinexEpochMs(s, openValue) {
   const [y, d, sec] = s.split(":").map(Number);
-  if (!y) return Infinity;
+  if (!y) return openValue;
   return Date.UTC(y, 0, 1) + ((d ?? 1) - 1) * 864e5 + (sec ?? 0) * 1e3;
 }
 function parseSinexBiasDcb(text, epochMs) {
@@ -699,8 +797,8 @@ function parseSinexBiasDcb(text, epochMs) {
     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 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]}`;
@@ -722,10 +820,6 @@ function parseSinexBiasDcb(text, epochMs) {
   }
   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;
@@ -789,7 +883,7 @@ function applyIonoDcb(iono, satDcb) {
   out.sort((a, b) => a.prn.localeCompare(b.prn));
   let sum = 0;
   let count = 0;
-  let maxStec = 0;
+  let maxStec = -Infinity;
   for (const s of out) {
     for (const p of s.points) {
       sum += p.stec;
@@ -798,7 +892,11 @@ function applyIonoDcb(iono, satDcb) {
     }
   }
   return {
-    result: { series: out, maxStec, meanStec: count > 0 ? sum / count : 0 },
+    result: {
+      series: out,
+      maxStec: count > 0 ? maxStec : 0,
+      meanStec: count > 0 ? sum / count : 0
+    },
     satellitesCorrected: corrected,
     satellitesMissing: [...new Set(missing)].sort(),
     receiverDcbTecu
@@ -840,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,
@@ -3906,6 +3908,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;
 }
@@ -3945,7 +3952,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;
@@ -3960,7 +3966,7 @@ function updateStreamStats(stats, frames, rawBytes) {
             prn: obs.prn,
             system: obs.system,
             cn0: bestCn0,
-            lastSeen: now2,
+            lastSeen: now,
             signals
           });
         }
@@ -4203,20 +4209,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") {
@@ -4309,8 +4305,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") {
@@ -4354,6 +4351,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;
@@ -4361,8 +4359,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,
@@ -4388,7 +4388,7 @@ function ephInfoToEphemeris(info) {
     omega: info.argPerigee,
     omegaDot: info.omegaDot ?? 0,
     idot: info.idot ?? 0,
-    week: info.week,
+    week,
     svHealth: info.health,
     tgd: 0
   };
@@ -5044,16 +5044,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 = [];
@@ -5534,15 +5540,10 @@ 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();
+  closedArcStarts = /* @__PURE__ */ new Map();
   interval;
   obsIndices;
   gloChannels;
@@ -5642,8 +5643,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) {
@@ -5655,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],
@@ -5674,7 +5686,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;
@@ -5686,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;
@@ -5701,21 +5713,107 @@ 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));
-    return { series, maxStec, meanStec: count > 0 ? sum / count : 0 };
+    return {
+      series,
+      maxStec: count > 0 ? maxStec : 0,
+      meanStec: count > 0 ? sum / count : 0
+    };
   }
 };
+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}$/;
 var SVN_RE = /^[A-Z]\d{3}$/;
 var OBS_RE = /^[CL]\d[A-Z]$/;
-function sinexEpochMs(s) {
+function sinexEpochMs(s, openValue) {
   const [y, d, sec] = s.split(":").map(Number);
-  if (!y) return Infinity;
+  if (!y) return openValue;
   return Date.UTC(y, 0, 1) + ((d ?? 1) - 1) * 864e5 + (sec ?? 0) * 1e3;
 }
 function parseSinexBiasDcb(text, epochMs) {
@@ -5730,8 +5828,8 @@ function parseSinexBiasDcb(text, epochMs) {
     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 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]}`;
@@ -5753,10 +5851,6 @@ function parseSinexBiasDcb(text, epochMs) {
   }
   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;
@@ -5820,7 +5914,7 @@ function applyIonoDcb(iono, satDcb) {
   out.sort((a, b) => a.prn.localeCompare(b.prn));
   let sum = 0;
   let count = 0;
-  let maxStec = 0;
+  let maxStec = -Infinity;
   for (const s of out) {
     for (const p of s.points) {
       sum += p.stec;
@@ -5829,7 +5923,11 @@ function applyIonoDcb(iono, satDcb) {
     }
   }
   return {
-    result: { series: out, maxStec, meanStec: count > 0 ? sum / count : 0 },
+    result: {
+      series: out,
+      maxStec: count > 0 ? maxStec : 0,
+      meanStec: count > 0 ? sum / count : 0
+    },
     satellitesCorrected: corrected,
     satellitesMissing: [...new Set(missing)].sort(),
     receiverDcbTecu
@@ -7148,6 +7246,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   clampUnit,
   computeAllPositions,
   computeDop,
+  computeIonoRate,
   computeLiveSkyPositions,
   computePsdDb,
   computeSatPosition,
@@ -7163,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-EJQ6NG5U.js";
+} from "./chunk-Z5UZUNCB.js";
 import {
   frequencyLabel,
   parseAntex
@@ -100,7 +102,7 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-PX7TPPTQ.js";
+} from "./chunk-SOJNWF5T.js";
 import {
   computeAllPositions,
   computeDop,
@@ -113,7 +115,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-PQYBTE42.js";
+} from "./chunk-CD5YSYCG.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -199,7 +201,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-65CQINSB.js";
+} from "./chunk-FGRC3PB3.js";
 import {
   DAYS_MJD2000_MJD,
   DAYS_MJD_JULIAN,
@@ -340,6 +342,7 @@ export {
   clampUnit,
   computeAllPositions,
   computeDop,
+  computeIonoRate,
   computeLiveSkyPositions,
   computePsdDb,
   computeSatPosition,
@@ -355,6 +358,7 @@ export {
   decodeMsmFull,
   deg2dms,
   deg2rad,
+  detrendIonoArcs,
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,