gnss-js 1.5.0 → 1.5.1

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

@@ -30,6 +30,18 @@ __export(analysis_exports, {
 });
 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 +229,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,12 +712,6 @@ var CompletenessAccumulator = class {
 // src/analysis/ionosphere.ts
 var TEC_FACTOR = 403e15;
 var GF_JUMP_M = 0.15;
-var MIN_ARC_LENGTH2 = 10;
-function median2(values) {
-  const s = [...values].sort((a, b) => a - b);
-  const m = s.length >> 1;
-  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
-}
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
@@ -814,8 +814,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) {
@@ -846,7 +846,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;
@@ -877,7 +877,11 @@ var IonoAccumulator = class {
       });
     }
     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
+    };
   }
 };
 
@@ -885,9 +889,9 @@ var IonoAccumulator = class {
 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 +906,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 +929,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 +992,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 +1001,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

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;

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;

package/dist/analysis.js

@@ -6,7 +6,7 @@ import {
   analyzeQuality,
   applyIonoDcb,
   parseSinexBiasDcb
-} from "./chunk-EJQ6NG5U.js";
+} from "./chunk-GDGHTC3E.js";
 import "./chunk-3U5AX7PY.js";
 import "./chunk-FIEWO4J4.js";
 export {

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-EJQ6NG5U.js → chunk-GDGHTC3E.js}

@@ -13,16 +13,22 @@ import {
   getFreq
 } from "./chunk-FIEWO4J4.js";
 
-// src/analysis/multipath.ts
+// src/analysis/stats-util.ts
 var MIN_ARC_LENGTH = 10;
-var MP_JUMP_M = 1.25;
-var MP_JUMP_WINDOW = 5;
-var MP_EDIT_SIGMA = 3;
 function median(values) {
   const s = [...values].sort((a, b) => a - b);
   const m = s.length >> 1;
   return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
 }
+function percentile(values, p) {
+  const s = [...values].sort((a, b) => a - b);
+  return s[Math.min(s.length - 1, Math.floor(p * s.length))] ?? 0;
+}
+
+// src/analysis/multipath.ts
+var MP_JUMP_M = 1.25;
+var MP_JUMP_WINDOW = 5;
+var MP_EDIT_SIGMA = 3;
 var MultipathAccumulator = class {
   state = /* @__PURE__ */ new Map();
   results = [];
@@ -503,12 +509,6 @@ var CompletenessAccumulator = class {
 // src/analysis/ionosphere.ts
 var TEC_FACTOR = 403e15;
 var GF_JUMP_M = 0.15;
-var MIN_ARC_LENGTH2 = 10;
-function median2(values) {
-  const s = [...values].sort((a, b) => a - b);
-  const m = s.length >> 1;
-  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
-}
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
@@ -611,8 +611,8 @@ var IonoAccumulator = class {
   }
   closeArc(prn, pairKey, ps) {
     const arc = ps.arc;
-    if (arc.times.length >= MIN_ARC_LENGTH2) {
-      const level = median2(arc.l4.map((v, k) => v - arc.p4[k]));
+    if (arc.times.length >= MIN_ARC_LENGTH) {
+      const level = median(arc.l4.map((v, k) => v - arc.p4[k]));
       const toTecu = arc.fi * arc.fi / TEC_FACTOR / (arc.gamma - 1);
       let satArcs = this.closed.get(prn);
       if (!satArcs) {
@@ -643,7 +643,7 @@ var IonoAccumulator = class {
     const series = [];
     let sum = 0;
     let count = 0;
-    let maxStec = 0;
+    let maxStec = -Infinity;
     for (const [prn, satArcs] of this.closed) {
       let bestKey = null;
       let bestLen = 0;
@@ -674,7 +674,11 @@ var IonoAccumulator = class {
       });
     }
     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
+    };
   }
 };
 
@@ -682,9 +686,9 @@ var IonoAccumulator = class {
 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 +703,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 +726,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 +789,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 +798,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

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

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

package/dist/index.cjs

@@ -3906,6 +3906,11 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1111 && msgType <= 1117) return "QZSS";
   if (msgType >= 1121 && msgType <= 1127) return "BeiDou";
   if (msgType >= 1131 && msgType <= 1137) return "NavIC";
+  if (msgType >= 1015 && msgType <= 1017) return "GPS";
+  if (msgType === 1030) return "GPS";
+  if (msgType === 1031) return "GLONASS";
+  if (msgType >= 1057 && msgType <= 1062) return "GPS";
+  if (msgType >= 1063 && msgType <= 1068) return "GLONASS";
   if (msgType === 1230) return "GLONASS";
   return null;
 }
@@ -3945,7 +3950,6 @@ function updateStreamStats(stats, frames, rawBytes) {
     entry.totalBytes += frame.length + 6;
     const msmEpoch = decodeMsmFull(frame);
     if (msmEpoch) {
-      const now2 = Date.now();
       for (const obs of msmEpoch.observations) {
         const signals = [];
         let bestCn0 = 0;
@@ -3960,7 +3964,7 @@ function updateStreamStats(stats, frames, rawBytes) {
             prn: obs.prn,
             system: obs.system,
             cn0: bestCn0,
-            lastSeen: now2,
+            lastSeen: now,
             signals
           });
         }
@@ -4203,20 +4207,10 @@ function invert4x4(m) {
   return result;
 }
 function selectEphemeris(ephemerides, prn, timeMs) {
-  let best = null;
-  let bestDt = Infinity;
-  for (const eph of ephemerides) {
-    if (eph.prn !== prn) continue;
-    const ephTime = eph.tocDate.getTime();
-    const dt = Math.abs(timeMs - ephTime);
-    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
-    if (dt > maxAge) continue;
-    if (dt < bestDt) {
-      bestDt = dt;
-      best = eph;
-    }
-  }
-  return best;
+  return selectBest(
+    ephemerides.filter((e) => e.prn === prn),
+    timeMs
+  );
 }
 function computeSatPosition(eph, timeMs) {
   if (eph.system === "R" || eph.system === "S") {
@@ -4309,8 +4303,9 @@ function selectBest(ephs, timeMs) {
   return best;
 }
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var MS_PER_WEEK2 = 7 * 864e5;
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
-var GAL_EPOCH_MS = GPS_EPOCH_MS;
+var GAL_EPOCH_MS = GPS_EPOCH_MS + 1024 * MS_PER_WEEK2;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -4354,6 +4349,7 @@ function ephInfoToEphemeris(info) {
     return null;
   }
   let epochMs;
+  let week = info.week;
   const tocSec = info.toc ?? info.toe;
   if (sys === "C") {
     epochMs = BDS_EPOCH_MS;
@@ -4361,8 +4357,10 @@ function ephInfoToEphemeris(info) {
     epochMs = GAL_EPOCH_MS;
   } else {
     epochMs = GPS_EPOCH_MS;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS) / MS_PER_WEEK2;
+    week += 1024 * Math.round((weeksNow - week) / 1024);
   }
-  const tocDate = new Date(epochMs + info.week * 7 * 864e5 + tocSec * 1e3);
+  const tocDate = new Date(epochMs + week * MS_PER_WEEK2 + tocSec * 1e3);
   return {
     system: sys,
     prn: info.prn,
@@ -4388,7 +4386,7 @@ function ephInfoToEphemeris(info) {
     omega: info.argPerigee,
     omegaDot: info.omegaDot ?? 0,
     idot: info.idot ?? 0,
-    week: info.week,
+    week,
     svHealth: info.health,
     tgd: 0
   };
@@ -5044,16 +5042,22 @@ async function connectToMountpoint(proxyUrl, info, signal) {
   };
 }
 
-// src/analysis/multipath.ts
+// src/analysis/stats-util.ts
 var MIN_ARC_LENGTH = 10;
-var MP_JUMP_M = 1.25;
-var MP_JUMP_WINDOW = 5;
-var MP_EDIT_SIGMA = 3;
 function median(values) {
   const s = [...values].sort((a, b) => a - b);
   const m = s.length >> 1;
   return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
 }
+function percentile(values, p) {
+  const s = [...values].sort((a, b) => a - b);
+  return s[Math.min(s.length - 1, Math.floor(p * s.length))] ?? 0;
+}
+
+// src/analysis/multipath.ts
+var MP_JUMP_M = 1.25;
+var MP_JUMP_WINDOW = 5;
+var MP_EDIT_SIGMA = 3;
 var MultipathAccumulator = class {
   state = /* @__PURE__ */ new Map();
   results = [];
@@ -5534,12 +5538,6 @@ var CompletenessAccumulator = class {
 // src/analysis/ionosphere.ts
 var TEC_FACTOR = 403e15;
 var GF_JUMP_M = 0.15;
-var MIN_ARC_LENGTH2 = 10;
-function median2(values) {
-  const s = [...values].sort((a, b) => a - b);
-  const m = s.length >> 1;
-  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
-}
 var IonoAccumulator = class {
   state = /* @__PURE__ */ new Map();
   closed = /* @__PURE__ */ new Map();
@@ -5642,8 +5640,8 @@ var IonoAccumulator = class {
   }
   closeArc(prn, pairKey, ps) {
     const arc = ps.arc;
-    if (arc.times.length >= MIN_ARC_LENGTH2) {
-      const level = median2(arc.l4.map((v, k) => v - arc.p4[k]));
+    if (arc.times.length >= MIN_ARC_LENGTH) {
+      const level = median(arc.l4.map((v, k) => v - arc.p4[k]));
       const toTecu = arc.fi * arc.fi / TEC_FACTOR / (arc.gamma - 1);
       let satArcs = this.closed.get(prn);
       if (!satArcs) {
@@ -5674,7 +5672,7 @@ var IonoAccumulator = class {
     const series = [];
     let sum = 0;
     let count = 0;
-    let maxStec = 0;
+    let maxStec = -Infinity;
     for (const [prn, satArcs] of this.closed) {
       let bestKey = null;
       let bestLen = 0;
@@ -5705,7 +5703,11 @@ var IonoAccumulator = class {
       });
     }
     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
+    };
   }
 };
 
@@ -5713,9 +5715,9 @@ var IonoAccumulator = class {
 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 +5732,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 +5755,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 +5818,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 +5827,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

package/dist/index.js

@@ -11,7 +11,7 @@ import {
   analyzeQuality,
   applyIonoDcb,
   parseSinexBiasDcb
-} from "./chunk-EJQ6NG5U.js";
+} from "./chunk-GDGHTC3E.js";
 import {
   frequencyLabel,
   parseAntex
@@ -100,7 +100,7 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-PX7TPPTQ.js";
+} from "./chunk-SOJNWF5T.js";
 import {
   computeAllPositions,
   computeDop,
@@ -113,7 +113,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-PQYBTE42.js";
+} from "./chunk-CD5YSYCG.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -199,7 +199,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-65CQINSB.js";
+} from "./chunk-FGRC3PB3.js";
 import {
   DAYS_MJD2000_MJD,
   DAYS_MJD_JULIAN,

package/dist/orbit.cjs

@@ -372,20 +372,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") {
@@ -478,8 +468,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") {
@@ -523,6 +514,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;
@@ -530,8 +522,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,
@@ -557,7 +551,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/orbit.d.cts

@@ -131,7 +131,7 @@ interface VisibilityPass {
 interface VisibilityResult {
     /** Sample epochs (Unix ms). */
     times: number[];
-    /** Elevation (radians) per PRN per epoch; null when below −0.05 rad / no eph. */
+    /** Elevation (radians) per PRN per epoch; null when no valid ephemeris. */
     elevation: Record<string, (number | null)[]>;
     /** Number of satellites at or above the mask, per epoch. */
     visibleCount: number[];

package/dist/orbit.d.ts

@@ -131,7 +131,7 @@ interface VisibilityPass {
 interface VisibilityResult {
     /** Sample epochs (Unix ms). */
     times: number[];
-    /** Elevation (radians) per PRN per epoch; null when below −0.05 rad / no eph. */
+    /** Elevation (radians) per PRN per epoch; null when no valid ephemeris. */
     elevation: Record<string, (number | null)[]>;
     /** Number of satellites at or above the mask, per epoch. */
     visibleCount: number[];

package/dist/orbit.js

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

package/dist/positioning.cjs

@@ -367,7 +367,9 @@ function computeSatPosition(eph, timeMs) {
   return keplerPosition(eph, tow);
 }
 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 + 1024 * MS_PER_WEEK;
 
 // src/positioning/index.ts
 var GM_GPS2 = 3986005e8;

package/dist/positioning.js

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

package/dist/rtcm3.cjs

@@ -1406,6 +1406,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;
 }
@@ -1445,7 +1450,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;
@@ -1460,7 +1464,7 @@ function updateStreamStats(stats, frames, rawBytes) {
             prn: obs.prn,
             system: obs.system,
             cn0: bestCn0,
-            lastSeen: now2,
+            lastSeen: now,
             signals
           });
         }

package/dist/rtcm3.js

@@ -15,7 +15,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-65CQINSB.js";
+} from "./chunk-FGRC3PB3.js";
 import "./chunk-LEEU5OIO.js";
 import "./chunk-FIEWO4J4.js";
 export {

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.5.0",
+  "version": "1.5.1",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,
@@ -106,10 +106,10 @@
     "positioning"
   ],
   "author": "Miguel González <work@miguel.es>",
-  "license": "AGPL-3.0-only",
+  "license": "AGPL-3.0-only OR LicenseRef-Commercial",
   "repository": {
     "type": "git",
-    "url": "https://github.com/MiguelPuntoEs/gnss-js"
+    "url": "git+https://github.com/MiguelPuntoEs/gnss-js.git"
   },
   "homepage": "https://github.com/MiguelPuntoEs/gnss-js#readme",
   "bugs": {