gnss-js 1.5.0 → 1.6.0

package/README.md

@@ -94,12 +94,16 @@ const itrf = transformFrame(
 ```ts
 import { parseRinexStream, parseNavFile } from 'gnss-js/rinex';
 
-// Streaming observation parser (works with ReadableStream)
-const result = await parseRinexStream(stream, {
-  onObservation: (time, prn, codes, values) => {
-    // process each epoch
-  },
-});
+// Streaming observation parser (positional arguments)
+const result = await parseRinexStream(
+  file, // File (plain, Hatanaka .crx, or .gz)
+  (percent) => {}, // optional progress callback
+  undefined, // optional AbortSignal
+  (time, prn, codes, values) => {
+    // per-satellite observation callback
+  }
+);
+// result.header: RinexHeader
 
 // Navigation file parser
 const nav = parseNavFile(navFileText);
@@ -125,28 +129,36 @@ for (const frame of decoder) {
 import { computeAllPositions, computeDop } from 'gnss-js/orbit';
 
 const positions = computeAllPositions(ephemerides, times, receiverPosition);
-// positions.prns: string[], positions.times: number[], positions.az/el/lat/lon: Float64Array
+// positions.prns: string[], positions.times: number[]
+// positions.positions[prn][epochIdx]: { lat, lon, az, el } | null
 ```
 
 ### Signal analysis
 
 ```ts
-import {
-  MultipathAccumulator,
-  CycleSlipAccumulator,
-  CompletenessAccumulator,
-} from 'gnss-js/analysis';
-
-const mp = new MultipathAccumulator(header);
-const cs = new CycleSlipAccumulator(header);
-await parseRinexStream(stream, {
-  onObservation: (time, prn, codes, values) => {
-    mp.onObservation(time, prn, codes, values);
-    cs.onObservation(time, prn, codes, values);
-  },
-});
-const multipathResult = mp.finalize();
-const cycleSlipResult = cs.finalize();
+import { analyzeQuality } from 'gnss-js/analysis';
+
+// One re-parse pass: multipath (per code band, Anubis-style), cycle
+// slips, completeness, and slant-TEC ionosphere series.
+const q = await analyzeQuality(file, header);
+// q.multipath, q.cycleSlips, q.completeness, q.iono
+
+// The accumulators (MultipathAccumulator, CycleSlipAccumulator,
+// CompletenessAccumulator, IonoAccumulator) are also exported for
+// wiring into a custom parseRinexStream pass.
+```
+
+### Ionosphere and DCBs
+
+```ts
+import { parseSinexBiasDcb, applyIonoDcb } from 'gnss-js/analysis';
+
+// q.iono: slant TEC per satellite from the geometry-free phase
+// levelled to the code — DCB-biased until calibrated:
+const satDcb = parseSinexBiasDcb(sinexBiasText, obsEpochMs); // ESA .BIA / CAS .BSX
+const { result, receiverDcbTecu } = applyIonoDcb(q.iono, satDcb);
+// result.series[n].points: calibrated STEC (TECU); receiver bias
+// estimated from the night-time floor per system/signal pair.
 ```
 
 ### NTRIP client
@@ -173,6 +185,17 @@ const { reader, abort } = await connectToMountpoint(
 );
 ```
 
+### NMEA and ANTEX
+
+```ts
+import { parseNmeaFile, computeStats } from 'gnss-js/nmea';
+import { parseAntex } from 'gnss-js/antex';
+
+const track = parseNmeaFile(nmeaText); // fixes, satellites, per-sentence records
+const stats = computeStats(track.fixes); // position/HDOP/speed statistics
+const antex = parseAntex(antexText); // antenna PCO/PCV calibrations
+```
+
 ### Constants
 
 ```ts

package/dist/analysis.cjs

@@ -26,10 +26,24 @@ __export(analysis_exports, {
   MultipathAccumulator: () => MultipathAccumulator,
   analyzeQuality: () => analyzeQuality,
   applyIonoDcb: () => applyIonoDcb,
+  computeIonoRate: () => computeIonoRate,
+  detrendIonoArcs: () => detrendIonoArcs,
   parseSinexBiasDcb: () => parseSinexBiasDcb
 });
 module.exports = __toCommonJS(analysis_exports);
 
+// src/analysis/stats-util.ts
+var MIN_ARC_LENGTH = 10;
+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/constants/gnss.ts
 var SYSTEM_NAMES = {
   G: "GPS",
@@ -217,15 +231,9 @@ function buildObsIndices(header) {
 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 = [];
@@ -706,15 +714,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;
@@ -814,8 +817,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) {
@@ -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],
@@ -846,7 +860,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;
@@ -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,21 +887,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) {
@@ -902,8 +1002,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]}`;
@@ -925,10 +1025,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;
@@ -992,7 +1088,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;
@@ -1001,7 +1097,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
@@ -1831,5 +1931,7 @@ async function analyzeQuality(file, header, onProgress, signal) {
   MultipathAccumulator,
   analyzeQuality,
   applyIonoDcb,
+  computeIonoRate,
+  detrendIonoArcs,
   parseSinexBiasDcb
 });

package/dist/analysis.d.cts

@@ -26,7 +26,7 @@ interface MultipathSeries {
     band: string;
     /** Reference band digit used in dual-frequency combination, e.g. "2". */
     refBand: string;
-    /** Human-readable label, e.g. "MP1 L1-L2". */
+    /** Human-readable label, e.g. "G01 MP L1-L2". */
     label: string;
     /** Time series of multipath values. */
     points: MultipathPoint[];
@@ -34,7 +34,7 @@ interface MultipathSeries {
     rms: number;
 }
 interface MultipathSignalStat {
-    /** Human-readable label, e.g. "MP1 L1-L2 (GPS)". */
+    /** Human-readable label, e.g. "MP L1-L2 (GPS)". */
     label: string;
     /** System letter, e.g. "G". */
     system: string;
@@ -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

@@ -26,7 +26,7 @@ interface MultipathSeries {
     band: string;
     /** Reference band digit used in dual-frequency combination, e.g. "2". */
     refBand: string;
-    /** Human-readable label, e.g. "MP1 L1-L2". */
+    /** Human-readable label, e.g. "G01 MP L1-L2". */
     label: string;
     /** Time series of multipath values. */
     points: MultipathPoint[];
@@ -34,7 +34,7 @@ interface MultipathSeries {
     rms: number;
 }
 interface MultipathSignalStat {
-    /** Human-readable label, e.g. "MP1 L1-L2 (GPS)". */
+    /** Human-readable label, e.g. "MP L1-L2 (GPS)". */
     label: string;
     /** System letter, e.g. "G". */
     system: string;
@@ -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-EJQ6NG5U.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-PQYBTE42.js → chunk-CD5YSYCG.js}

@@ -219,20 +219,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") {
@@ -325,8 +315,9 @@ function selectBest(ephs, timeMs) {
   return best;
 }
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var MS_PER_WEEK = 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_WEEK;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -370,6 +361,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;
@@ -377,8 +369,10 @@ function ephInfoToEphemeris(info) {
     epochMs = GAL_EPOCH_MS;
   } else {
     epochMs = GPS_EPOCH_MS;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS) / MS_PER_WEEK;
+    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_WEEK + tocSec * 1e3);
   return {
     system: sys,
     prn: info.prn,
@@ -404,7 +398,7 @@ function ephInfoToEphemeris(info) {
     omega: info.argPerigee,
     omegaDot: info.omegaDot ?? 0,
     idot: info.idot ?? 0,
-    week: info.week,
+    week,
     svHealth: info.health,
     tgd: 0
   };

package/dist/{chunk-65CQINSB.js → chunk-FGRC3PB3.js}

@@ -1345,6 +1345,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;
 }
@@ -1384,7 +1389,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;
@@ -1399,7 +1403,7 @@ function updateStreamStats(stats, frames, rawBytes) {
             prn: obs.prn,
             system: obs.system,
             cn0: bestCn0,
-            lastSeen: now2,
+            lastSeen: now,
             signals
           });
         }

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