gnss-js 1.0.0 → 1.1.0

package/dist/index.d.cts

@@ -7,8 +7,8 @@ export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer,
 export { CrxField, DiffState, EMPTY_WARNINGS, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, writeRinexNav } from './rinex.cjs';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.cjs';
 export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, p as parseNavFile } from './nav-BAI1a9vK.cjs';
-export { B as BitReader, E as EphemerisInfo, R as Rtcm3Decoder, a as Rtcm3Frame, d as decodeEphemeris, r as readString } from './ephemeris-BUWzfmBy.cjs';
-export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
+export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.cjs';
+export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
 export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.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, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.cjs';

package/dist/index.d.ts

@@ -7,8 +7,8 @@ export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer,
 export { CrxField, DiffState, EMPTY_WARNINGS, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, writeRinexNav } from './rinex.js';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.js';
 export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, p as parseNavFile } from './nav-BAI1a9vK.js';
-export { B as BitReader, E as EphemerisInfo, R as Rtcm3Decoder, a as Rtcm3Frame, d as decodeEphemeris, r as readString } from './ephemeris-BUWzfmBy.js';
-export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, updateStationMeta, updateStreamStats } from './rtcm3.js';
+export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.js';
+export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
 export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.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, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.js';

package/dist/index.js

@@ -3,7 +3,7 @@ import {
   CycleSlipAccumulator,
   MultipathAccumulator,
   analyzeQuality
-} from "./chunk-Y3R57B5P.js";
+} from "./chunk-G3N4S3DM.js";
 import {
   frequencyLabel,
   parseAntex
@@ -48,7 +48,7 @@ import {
   phiBOCc,
   phiBOCs,
   phiBPSK
-} from "./chunk-HKN3PUGN.js";
+} from "./chunk-5SPJH4MG.js";
 import {
   Scale,
   getBdsTime,
@@ -114,33 +114,8 @@ import {
   rad2deg,
   rhumbLine,
   vincenty
-} from "./chunk-K7WZQFBV.js";
+} from "./chunk-MIIM4LDY.js";
 import "./chunk-5S5MONFK.js";
-import {
-  DAYS_MJD2000_MJD,
-  DAYS_MJD_JULIAN,
-  MILLISECONDS_GPS_TAI,
-  MILLISECONDS_IN_DAY,
-  MILLISECONDS_IN_HOUR,
-  MILLISECONDS_IN_MINUTE,
-  MILLISECONDS_IN_SECOND,
-  MILLISECONDS_IN_WEEK,
-  MILLISECONDS_TT_TAI,
-  RINEX_CODES,
-  SECONDS_IN_DAY,
-  SECONDS_IN_HOUR,
-  SECONDS_IN_MINUTE,
-  SECONDS_IN_WEEK,
-  START_BDS_TIME,
-  START_GAL_TIME,
-  START_GLO_LEAP,
-  START_GPS_TIME,
-  START_JULIAN_TAI,
-  START_MJD_UNIX_SECONDS,
-  START_NTP_TIME,
-  START_TAI_TIME,
-  START_UNIX_TIME
-} from "./chunk-LEEU5OIO.js";
 import {
   EMPTY_WARNINGS,
   WarningAccumulator,
@@ -150,7 +125,7 @@ import {
   padR,
   parseNavFile,
   writeRinexNav
-} from "./chunk-IS4UUDBV.js";
+} from "./chunk-W4YMQKWH.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -159,7 +134,7 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-CB6EOOLA.js";
+} from "./chunk-OZCYOM5D.js";
 import {
   BitReader,
   RTCM3_MESSAGE_NAMES,
@@ -170,11 +145,14 @@ import {
   decodeMsmFull,
   msmEpochToDate,
   readString,
+  reportDecodeError,
   resetGloFreqCache,
   rtcm3Constellation,
+  setGloFreqNumber,
+  setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-PRSZIWKM.js";
+} from "./chunk-JDO3LEPC.js";
 import {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -196,7 +174,7 @@ import {
   formatUTCTime,
   getFreq,
   gloFreq
-} from "./chunk-354IRDOG.js";
+} from "./chunk-W5WKEV7U.js";
 import {
   computeAllPositions,
   computeDop,
@@ -207,7 +185,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-SDRRAJT5.js";
+} from "./chunk-PGHDJFQK.js";
 import {
   clampUnit,
   ecefToGeodetic,
@@ -221,11 +199,36 @@ import {
   WGS84_SEMI_MAJOR_AXIS,
   WGS84_SEMI_MINOR_AXIS
 } from "./chunk-6FAL6P4G.js";
+import {
+  DAYS_MJD2000_MJD,
+  DAYS_MJD_JULIAN,
+  MILLISECONDS_GPS_TAI,
+  MILLISECONDS_IN_DAY,
+  MILLISECONDS_IN_HOUR,
+  MILLISECONDS_IN_MINUTE,
+  MILLISECONDS_IN_SECOND,
+  MILLISECONDS_IN_WEEK,
+  MILLISECONDS_TT_TAI,
+  RINEX_CODES,
+  SECONDS_IN_DAY,
+  SECONDS_IN_HOUR,
+  SECONDS_IN_MINUTE,
+  SECONDS_IN_WEEK,
+  START_BDS_TIME,
+  START_GAL_TIME,
+  START_GLO_LEAP,
+  START_GPS_TIME,
+  START_JULIAN_TAI,
+  START_MJD_UNIX_SECONDS,
+  START_NTP_TIME,
+  START_TAI_TIME,
+  START_UNIX_TIME
+} from "./chunk-LEEU5OIO.js";
 import {
   connectToMountpoint,
   fetchSourcetable,
   parseSourcetable
-} from "./chunk-WR7LCB52.js";
+} from "./chunk-BJHTBYNG.js";
 export {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -413,10 +416,13 @@ export {
   phiBPSK,
   rad2deg,
   readString,
+  reportDecodeError,
   resetGloFreqCache,
   rhumbLine,
   rtcm3Constellation,
   selectEphemeris,
+  setGloFreqNumber,
+  setRtcm3DebugHandler,
   systemCmp,
   systemName,
   updateStationMeta,

package/dist/ntrip.cjs

@@ -124,14 +124,18 @@ async function ntripFetch(proxyUrl, path, headers, signal) {
   } catch (err) {
     if (err instanceof Error && err.name === "AbortError") throw err;
     const message = err instanceof Error ? err.message : "The service may be temporarily unavailable.";
-    throw new Error(`Could not reach the NTRIP proxy: ${message}`, { cause: err });
+    throw new Error(`Could not reach the NTRIP proxy: ${message}`, {
+      cause: err
+    });
   }
 }
 async function fetchSourcetable(proxyUrl, info, signal) {
   const headers = ntripHeaders(info);
   const res = await ntripFetch(proxyUrl, "/", headers, signal);
   if (res.status === 401) {
-    throw new Error("Authentication required. Please provide valid credentials.");
+    throw new Error(
+      "Authentication required. Please provide valid credentials."
+    );
   }
   if (!res.ok) {
     throw new Error(`Caster returned ${res.status} ${res.statusText}`);
@@ -143,7 +147,12 @@ async function connectToMountpoint(proxyUrl, info, signal) {
   const headers = ntripHeaders(info);
   const controller = new AbortController();
   const combinedSignal = signal ? AbortSignal.any([signal, controller.signal]) : controller.signal;
-  const res = await ntripFetch(proxyUrl, `/${info.mountpoint}`, headers, combinedSignal);
+  const res = await ntripFetch(
+    proxyUrl,
+    `/${info.mountpoint}`,
+    headers,
+    combinedSignal
+  );
   if (res.status === 401) {
     throw new Error("Authentication required for this mountpoint.");
   }
@@ -154,7 +163,9 @@ async function connectToMountpoint(proxyUrl, info, signal) {
     throw new Error(`Caster returned ${res.status} ${res.statusText}`);
   }
   if (!res.body) {
-    throw new Error("No response body \u2014 streaming not supported by this environment.");
+    throw new Error(
+      "No response body \u2014 streaming not supported by this environment."
+    );
   }
   return {
     reader: res.body.getReader(),

package/dist/ntrip.js

@@ -2,7 +2,7 @@ import {
   connectToMountpoint,
   fetchSourcetable,
   parseSourcetable
-} from "./chunk-WR7LCB52.js";
+} from "./chunk-BJHTBYNG.js";
 export {
   connectToMountpoint,
   fetchSourcetable,

package/dist/orbit.cjs

@@ -71,12 +71,18 @@ function ecefToGeodetic(x, y, z) {
   return [lat, lon, alt];
 }
 
+// src/constants/time.ts
+var START_GPS_TIME = /* @__PURE__ */ new Date("1980-01-06T00:00:00Z");
+var START_BDS_TIME = /* @__PURE__ */ new Date("2006-01-01T00:00:14Z");
+
 // src/orbit/index.ts
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
 var OMEGA_E = 72921151467e-15;
+var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
+var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
 var J2_GLO = 108263e-8;
 var TWO_PI = 2 * Math.PI;
@@ -102,10 +108,7 @@ function keplerPosition(eph, t) {
   }
   const sinE = Math.sin(Ek);
   const cosE = Math.cos(Ek);
-  const vk = Math.atan2(
-    Math.sqrt(1 - eph.e * eph.e) * sinE,
-    cosE - eph.e
-  );
+  const vk = Math.atan2(Math.sqrt(1 - eph.e * eph.e) * sinE, cosE - eph.e);
   const phik = vk + eph.omega;
   const sin2phi = Math.sin(2 * phik);
   const cos2phi = Math.cos(2 * phik);
@@ -117,7 +120,7 @@ function keplerPosition(eph, t) {
   const ik = eph.i0 + dik + eph.idot * tk;
   const xp = rk * Math.cos(uk);
   const yp = rk * Math.sin(uk);
-  const isBdsGeo = eph.system === "C" && Math.abs(eph.i0) < 0.1 && a > 4e7;
+  const isBdsGeo = eph.system === "C" && Math.abs(eph.i0) < BDS_GEO_MAX_INCLINATION_RAD && a > BDS_GEO_MIN_SEMIMAJOR_AXIS_M;
   if (isBdsGeo) {
     const omegak2 = eph.omega0 + eph.omegaDot * tk - OMEGA_E * eph.toe;
     const cosO2 = Math.cos(omegak2);
@@ -268,7 +271,10 @@ function invert4x4(m) {
   let det = a[0] * inv[0] + a[1] * inv[4] + a[2] * inv[8] + a[3] * inv[12];
   if (Math.abs(det) < 1e-20) return null;
   det = 1 / det;
-  const result = Array.from({ length: 4 }, () => new Array(4).fill(0));
+  const result = Array.from(
+    { length: 4 },
+    () => new Array(4).fill(0)
+  );
   for (let i = 0; i < 4; i++) {
     for (let j = 0; j < 4; j++) {
       result[i][j] = inv[i * 4 + j] * det;
@@ -295,9 +301,10 @@ function computeSatPosition(eph, timeMs) {
   if (eph.system === "R" || eph.system === "S") {
     return glonassPosition(eph, timeMs / 1e3);
   }
-  const GPS_EPOCH = Date.UTC(1980, 0, 6);
+  const GPS_EPOCH = START_GPS_TIME.getTime();
   const gpsSeconds = (timeMs - GPS_EPOCH) / 1e3;
-  const tow = gpsSeconds % (7 * 86400);
+  let tow = gpsSeconds % (7 * 86400);
+  if (eph.system === "C") tow = (tow - 14 + 604800) % 604800;
   return keplerPosition(eph, tow);
 }
 function navTimesFromEph(ephs) {
@@ -379,16 +386,22 @@ function selectBest(ephs, timeMs) {
   }
   return best;
 }
-var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
-var BDS_EPOCH_MS = Date.UTC(2006, 0, 1);
+var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var BDS_EPOCH_MS = START_BDS_TIME.getTime();
 var GAL_EPOCH_MS = GPS_EPOCH_MS;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
-    if (info.x === void 0 || info.y === void 0 || info.z === void 0) return null;
-    if (info.vx === void 0 || info.vy === void 0 || info.vz === void 0) return null;
+    if (info.x === void 0 || info.y === void 0 || info.z === void 0)
+      return null;
+    if (info.vx === void 0 || info.vy === void 0 || info.vz === void 0)
+      return null;
     const now = new Date(info.lastReceived);
-    const utcMidnight = Date.UTC(now.getUTCFullYear(), now.getUTCMonth(), now.getUTCDate());
+    const utcMidnight = Date.UTC(
+      now.getUTCFullYear(),
+      now.getUTCMonth(),
+      now.getUTCDate()
+    );
     const tbSec = (info.tb ?? 0) * 900;
     const moscowOffset = 3 * 3600;
     const tocMs = utcMidnight + (tbSec - moscowOffset) * 1e3;
@@ -413,7 +426,11 @@ function ephInfoToEphemeris(info) {
       freqNum: info.freqChannel ?? 0
     };
   }
-  if (info.sqrtA === void 0 || info.eccentricity === void 0 || info.inclination === void 0 || info.omega0 === void 0 || info.argPerigee === void 0 || info.meanAnomaly === void 0 || info.toe === void 0 || info.week === void 0) return null;
+  if (info.sqrtA === void 0 || info.eccentricity === void 0 || info.inclination === void 0 || info.omega0 === void 0 || info.argPerigee === void 0 || info.meanAnomaly === void 0 || info.toe === void 0 || info.week === void 0)
+    return null;
+  if (info.eccentricity < 0 || info.eccentricity >= 1 || info.sqrtA < 3e3 || info.sqrtA > 9e3 || info.toe < 0 || info.toe > 604800) {
+    return null;
+  }
   let epochMs;
   const tocSec = info.toc ?? info.toe;
   if (sys === "C") {

package/dist/orbit.d.cts

@@ -1,5 +1,5 @@
 import { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris } from './nav-BAI1a9vK.cjs';
-import { E as EphemerisInfo } from './ephemeris-BUWzfmBy.cjs';
+import { E as EphemerisInfo } from './ephemeris-Ohl6NAfw.cjs';
 export { e as ecefToGeodetic, g as geodeticToEcef } from './ecef-CF0uAysr.cjs';
 
 /**

package/dist/orbit.d.ts

@@ -1,5 +1,5 @@
 import { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris } from './nav-BAI1a9vK.js';
-import { E as EphemerisInfo } from './ephemeris-BUWzfmBy.js';
+import { E as EphemerisInfo } from './ephemeris-Ohl6NAfw.js';
 export { e as ecefToGeodetic, g as geodeticToEcef } from './ecef-CF0uAysr.js';
 
 /**

package/dist/orbit.js

@@ -8,12 +8,13 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-SDRRAJT5.js";
+} from "./chunk-PGHDJFQK.js";
 import {
   ecefToGeodetic,
   geodeticToEcef
 } from "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
+import "./chunk-LEEU5OIO.js";
 export {
   computeAllPositions,
   computeDop,

package/dist/rinex.cjs

@@ -78,10 +78,20 @@ function parseCrxDataLine(line, ntype) {
     if (token.length >= 3 && token[1] === "&") {
       const arcOrder = parseInt(token[0]);
       const value = parseInt(token.substring(2));
-      fields.push({ empty: false, init: true, arcOrder, value: isNaN(value) ? 0 : value });
+      fields.push({
+        empty: false,
+        init: true,
+        arcOrder,
+        value: isNaN(value) ? 0 : value
+      });
     } else {
       const value = parseInt(token);
-      fields.push({ empty: false, init: false, arcOrder: -1, value: isNaN(value) ? 0 : value });
+      fields.push({
+        empty: false,
+        init: false,
+        arcOrder: -1,
+        value: isNaN(value) ? 0 : value
+      });
     }
   }
   const flagStr = pos < line.length ? line.substring(pos) : "";
@@ -91,7 +101,10 @@ function crxDecompress(prev, field) {
   if (field.init) {
     const a = new Array(Math.max(field.arcOrder + 1, 1)).fill(0);
     a[0] = field.value;
-    return { state: { accum: a, order: 0, arcOrder: field.arcOrder }, result: field.value };
+    return {
+      state: { accum: a, order: 0, arcOrder: field.arcOrder },
+      result: field.value
+    };
   }
   if (!prev) {
     return { state: { accum: [0], order: 0, arcOrder: 0 }, result: 0 };
@@ -104,7 +117,10 @@ function crxDecompress(prev, field) {
   for (let k = order; k >= 1; k--) {
     accum[k - 1] += accum[k];
   }
-  return { state: { accum, order, arcOrder: prev.arcOrder }, result: accum[0] };
+  return {
+    state: { accum, order, arcOrder: prev.arcOrder },
+    result: accum[0]
+  };
 }
 
 // src/rinex/format.ts
@@ -137,7 +153,9 @@ var CHUNK_SIZE = 2 * 1024 * 1024;
 var yieldToMain = () => new Promise((r) => setTimeout(r, 0));
 var noYield = () => Promise.resolve();
 var SYSTEM_ORDER = ["G", "R", "E", "C", "J", "I", "S"];
-var SYSTEM_RANK = Object.fromEntries(SYSTEM_ORDER.map((s, i) => [s, i]));
+var SYSTEM_RANK = Object.fromEntries(
+  SYSTEM_ORDER.map((s, i) => [s, i])
+);
 function systemCmp(a, b) {
   const ra = SYSTEM_RANK[a.charAt(0)] ?? 99;
   const rb = SYSTEM_RANK[b.charAt(0)] ?? 99;
@@ -323,7 +341,11 @@ function parseEpochLine3(line) {
   const wholeSec = Math.floor(sc);
   const ms = Math.round((sc - wholeSec) * 1e3);
   const time = Date.UTC(yr, mo - 1, dy, hr, mn, wholeSec, ms);
-  return { time, flag: isNaN(flag) ? 0 : flag, numSats: isNaN(numSats) ? 0 : numSats };
+  return {
+    time,
+    flag: isNaN(flag) ? 0 : flag,
+    numSats: isNaN(numSats) ? 0 : numSats
+  };
 }
 function parseEpochLine2(line) {
   const yr2 = parseInt(line.substring(1, 3));
@@ -345,7 +367,12 @@ function parseEpochLine2(line) {
     const id = satPart.substring(i, i + 3).trim();
     if (id) satIds.push(id);
   }
-  return { time, flag: isNaN(flag) ? 0 : flag, numSats: isNaN(numSats) ? 0 : numSats, satIds };
+  return {
+    time,
+    flag: isNaN(flag) ? 0 : flag,
+    numSats: isNaN(numSats) ? 0 : numSats,
+    satIds
+  };
 }
 async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode) {
   const skipEpochs = !!workerMode;
@@ -398,16 +425,19 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
     if (!skipEpochs) {
       const snrBySys = {};
       for (const [sys, vals] of Object.entries(snrPerSystemAccum)) {
-        if (vals.length > 0) snrBySys[sys] = vals.reduce((a, b) => a + b, 0) / vals.length;
+        if (vals.length > 0)
+          snrBySys[sys] = vals.reduce((a, b) => a + b, 0) / vals.length;
       }
       const meanSnr = snrValues.length > 0 ? snrValues.reduce((a, b) => a + b, 0) / snrValues.length : null;
       const snrPerSat = {};
       for (const [prn, vals] of Object.entries(snrPerSatAccum)) {
-        if (vals.length > 0) snrPerSat[prn] = vals.reduce((a, b) => a + b, 0) / vals.length;
+        if (vals.length > 0)
+          snrPerSat[prn] = vals.reduce((a, b) => a + b, 0) / vals.length;
       }
       const snrPerSatBand = {};
       for (const [key, vals] of Object.entries(snrPerSatBandAccum)) {
-        if (vals.length > 0) snrPerSatBand[key] = vals.reduce((a, b) => a + b, 0) / vals.length;
+        if (vals.length > 0)
+          snrPerSatBand[key] = vals.reduce((a, b) => a + b, 0) / vals.length;
       }
       epochs.push({
         time: epochInfo.time,
@@ -452,7 +482,8 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
     const codes = header.obsTypes[sys] ?? [];
     if (onSatObs && epochInfo) {
       const values = new Array(codes.length);
-      for (let i = 0; i < codes.length; i++) values[i] = readObsValue(obsLine, i);
+      for (let i = 0; i < codes.length; i++)
+        values[i] = readObsValue(obsLine, i);
       if (!skipEpochs) {
         for (const { idx, band } of getSnrIndices(sys)) {
           const val = values[idx] ?? null;
@@ -689,7 +720,10 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
           v2SnrBandInfo = snrIndicesWithBandV2(header.obsTypes);
           v2SatIndex = 0;
           v2CurrentSatLine = 0;
-          const continuationLines = Math.max(0, Math.ceil((info.numSats - 12) / 12));
+          const continuationLines = Math.max(
+            0,
+            Math.ceil((info.numSats - 12) / 12)
+          );
           if (continuationLines > 0) {
             v2ContinuationSatsRemaining = continuationLines;
           } else {
@@ -729,7 +763,9 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
   const magic = new Uint8Array(await file.slice(0, 2).arrayBuffer());
   const isGz = magic[0] === 31 && magic[1] === 139;
   if (magic[0] === 31 && magic[1] === 157) {
-    throw new Error("Unix compress (.Z) files are not supported. Please decompress first (e.g. uncompress or gzip -d).");
+    throw new Error(
+      "Unix compress (.Z) files are not supported. Please decompress first (e.g. uncompress or gzip -d)."
+    );
   }
   if (isGz) {
     let bytesRead = 0;
@@ -745,7 +781,9 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
       if (done) break;
       bytesRead += value.byteLength;
       processChunkText(decoder.decode(value, { stream: true }), false);
-      onProgress?.(Math.min(99, Math.round(bytesRead / (file.size * 4) * 100)));
+      onProgress?.(
+        Math.min(99, Math.round(bytesRead / (file.size * 4) * 100))
+      );
       await yield_();
     }
     processChunkText(decoder.decode(), true);
@@ -755,7 +793,10 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
       const end = Math.min(offset + CHUNK_SIZE, file.size);
       const slice = file.slice(offset, end);
       const arrayBuf = await slice.arrayBuffer();
-      processChunkText(decoder.decode(arrayBuf, { stream: end < file.size }), false);
+      processChunkText(
+        decoder.decode(arrayBuf, { stream: end < file.size }),
+        false
+      );
       onProgress?.(Math.min(99, Math.round(end / file.size * 100)));
       await yield_();
     }
@@ -763,19 +804,27 @@ async function parseRinexStream(file, onProgress, signal, onSatObs, workerMode)
       processChunkText("", true);
     }
   }
-  if (epochInfo && satLinesRemaining === 0 && crxPhase === "epoch") finishEpoch();
-  if (!header) throw new Error("No valid RINEX header found (missing END OF HEADER).");
+  if (epochInfo && satLinesRemaining === 0 && crxPhase === "epoch")
+    finishEpoch();
+  if (!header)
+    throw new Error("No valid RINEX header found (missing END OF HEADER).");
   onProgress?.(100);
-  return { header, epochs, stats: computeStats(header, epochs, satellitesSeen) };
+  return {
+    header,
+    epochs,
+    stats: computeStats(header, epochs, satellitesSeen)
+  };
 }
 function computeStats(header, epochs, satellitesSeen) {
   const n = epochs.length;
   const startTime = n > 0 ? new Date(epochs[0].time) : header.timeOfFirstObs;
   const endTime = n > 0 ? new Date(epochs[n - 1].time) : header.timeOfLastObs;
   let duration = null;
-  if (startTime && endTime) duration = (endTime.getTime() - startTime.getTime()) / 1e3;
+  if (startTime && endTime)
+    duration = (endTime.getTime() - startTime.getTime()) / 1e3;
   let interval = header.interval;
-  if (interval === null && n >= 2) interval = (epochs[1].time - epochs[0].time) / 1e3;
+  if (interval === null && n >= 2)
+    interval = (epochs[1].time - epochs[0].time) / 1e3;
   const uniqueSatsPerSystem = {};
   let totalUnique = 0;
   const systems = [];
@@ -830,10 +879,18 @@ var WarningAccumulator = class {
       );
     }
     if (!header.antType) {
-      this.add("MISSING_ANT_TYPE", "info", "Missing antenna type (ANT # / TYPE)");
+      this.add(
+        "MISSING_ANT_TYPE",
+        "info",
+        "Missing antenna type (ANT # / TYPE)"
+      );
     }
     if (!header.receiverType) {
-      this.add("MISSING_REC_TYPE", "info", "Missing receiver type (REC # / TYPE / VERS)");
+      this.add(
+        "MISSING_REC_TYPE",
+        "info",
+        "Missing receiver type (REC # / TYPE / VERS)"
+      );
     }
     if (!header.markerName) {
       this.add("MISSING_MARKER", "info", "Missing marker name (MARKER NAME)");
@@ -842,7 +899,11 @@ var WarningAccumulator = class {
       this.add("MISSING_INTERVAL", "info", "No observation interval in header");
     }
     if (Object.keys(header.obsTypes).length === 0) {
-      this.add("NO_OBS_TYPES", "error", "No observation types defined in header");
+      this.add(
+        "NO_OBS_TYPES",
+        "error",
+        "No observation types defined in header"
+      );
     }
   }
   setInterval(intervalMs) {
@@ -909,7 +970,12 @@ var WarningAccumulator = class {
     };
     const range = ranges[sys];
     if (!range) {
-      this.add("UNKNOWN_SYSTEM", "info", "Unknown satellite system letter", prn);
+      this.add(
+        "UNKNOWN_SYSTEM",
+        "info",
+        "Unknown satellite system letter",
+        prn
+      );
     } else if (num < range[0] || num > range[1]) {
       this.add(
         "PRN_OUT_OF_RANGE",
@@ -999,7 +1065,15 @@ var R4_MSG_LINES = {
   EOP: 3,
   ION: 3
 };
-var SUPPORTED_EPH_MSGS = /* @__PURE__ */ new Set(["LNAV", "INAV", "FNAV", "D1", "D2", "SBAS", "FDMA"]);
+var SUPPORTED_EPH_MSGS = /* @__PURE__ */ new Set([
+  "LNAV",
+  "INAV",
+  "FNAV",
+  "D1",
+  "D2",
+  "SBAS",
+  "FDMA"
+]);
 function parseFloat19(s) {
   return parseFloat(s.trim().replace(/[dD]/g, "E"));
 }
@@ -1206,9 +1280,25 @@ function parseNavFile(text) {
       }
       if (dataLines2.length === numDataLines2) {
         if (sys2 === "R" || sys2 === "S") {
-          ephemerides.push(buildStateVectorEphemeris(sys2, prn2, parsed.date, parsed.values, dataLines2));
+          ephemerides.push(
+            buildStateVectorEphemeris(
+              sys2,
+              prn2,
+              parsed.date,
+              parsed.values,
+              dataLines2
+            )
+          );
         } else if (sys2 === "G" || sys2 === "E" || sys2 === "C" || sys2 === "J" || sys2 === "I") {
-          ephemerides.push(buildKeplerEphemeris(sys2, prn2, parsed.date, parsed.values, dataLines2));
+          ephemerides.push(
+            buildKeplerEphemeris(
+              sys2,
+              prn2,
+              parsed.date,
+              parsed.values,
+              dataLines2
+            )
+          );
         }
       }
       i++;
@@ -1244,9 +1334,13 @@ function parseNavFile(text) {
     }
     if (dataLines.length === numDataLines) {
       if (sys === "R" || sys === "S") {
-        ephemerides.push(buildStateVectorEphemeris(sys, prn, date, epochVals, dataLines));
+        ephemerides.push(
+          buildStateVectorEphemeris(sys, prn, date, epochVals, dataLines)
+        );
       } else if (sys === "G" || sys === "E" || sys === "C" || sys === "J" || sys === "I") {
-        ephemerides.push(buildKeplerEphemeris(sys, prn, date, epochVals, dataLines));
+        ephemerides.push(
+          buildKeplerEphemeris(sys, prn, date, epochVals, dataLines)
+        );
       }
     }
     i++;