gnss-js 1.0.1 → 1.1.0

package/dist/constants.js

@@ -1,29 +1,4 @@
 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 {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -45,13 +20,38 @@ import {
   formatUTCTime,
   getFreq,
   gloFreq
-} from "./chunk-LWNTWBHB.js";
+} from "./chunk-W5WKEV7U.js";
 import {
   WGS84_ECCENTRICITY_SQUARED,
   WGS84_FLATTENING,
   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";
 export {
   ARC_GAP_FACTOR,
   BAND_LABELS,

package/dist/{ephemeris-C10stHhM.d.cts → ephemeris-Ohl6NAfw.d.cts}

@@ -7,6 +7,16 @@ interface Rtcm3Frame {
     length: number;
     payload: Uint8Array;
 }
+/** Called when a message decoder swallows an exception. */
+type Rtcm3DebugHandler = (context: string, error: unknown) => void;
+/**
+ * Install a handler to observe decode errors that are otherwise
+ * silently mapped to null results (corrupt frames, truncated payloads).
+ * Pass null to remove.
+ */
+declare function setRtcm3DebugHandler(fn: Rtcm3DebugHandler | null): void;
+/** @internal Report a swallowed decode error to the debug handler. */
+declare function reportDecodeError(context: string, error: unknown): void;
 /** Simple bit reader for parsing RTCM3 payloads at bit level. */
 declare class BitReader {
     private data;
@@ -111,4 +121,4 @@ declare function readString(r: BitReader, nChars: number): string;
 /** Decode any supported ephemeris message. Returns null for non-ephemeris or decode errors. */
 declare function decodeEphemeris(frame: Rtcm3Frame): EphemerisInfo | null;
 
-export { BitReader as B, type EphemerisInfo as E, Rtcm3Decoder as R, type Rtcm3Frame as a, decodeEphemeris as d, readString as r };
+export { BitReader as B, type EphemerisInfo as E, type Rtcm3DebugHandler as R, Rtcm3Decoder as a, type Rtcm3Frame as b, reportDecodeError as c, decodeEphemeris as d, readString as r, setRtcm3DebugHandler as s };

package/dist/{ephemeris-C10stHhM.d.ts → ephemeris-Ohl6NAfw.d.ts}

@@ -7,6 +7,16 @@ interface Rtcm3Frame {
     length: number;
     payload: Uint8Array;
 }
+/** Called when a message decoder swallows an exception. */
+type Rtcm3DebugHandler = (context: string, error: unknown) => void;
+/**
+ * Install a handler to observe decode errors that are otherwise
+ * silently mapped to null results (corrupt frames, truncated payloads).
+ * Pass null to remove.
+ */
+declare function setRtcm3DebugHandler(fn: Rtcm3DebugHandler | null): void;
+/** @internal Report a swallowed decode error to the debug handler. */
+declare function reportDecodeError(context: string, error: unknown): void;
 /** Simple bit reader for parsing RTCM3 payloads at bit level. */
 declare class BitReader {
     private data;
@@ -111,4 +121,4 @@ declare function readString(r: BitReader, nChars: number): string;
 /** Decode any supported ephemeris message. Returns null for non-ephemeris or decode errors. */
 declare function decodeEphemeris(frame: Rtcm3Frame): EphemerisInfo | null;
 
-export { BitReader as B, type EphemerisInfo as E, Rtcm3Decoder as R, type Rtcm3Frame as a, decodeEphemeris as d, readString as r };
+export { BitReader as B, type EphemerisInfo as E, type Rtcm3DebugHandler as R, Rtcm3Decoder as a, type Rtcm3Frame as b, reportDecodeError as c, decodeEphemeris as d, readString as r, setRtcm3DebugHandler as s };

package/dist/index.cjs

@@ -206,10 +206,13 @@ __export(src_exports, {
   phiBPSK: () => phiBPSK,
   rad2deg: () => rad2deg,
   readString: () => readString,
+  reportDecodeError: () => reportDecodeError,
   resetGloFreqCache: () => resetGloFreqCache,
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   selectEphemeris: () => selectEphemeris,
+  setGloFreqNumber: () => setGloFreqNumber,
+  setRtcm3DebugHandler: () => setRtcm3DebugHandler,
   systemCmp: () => systemCmp,
   systemName: () => systemName,
   updateStationMeta: () => updateStationMeta,
@@ -298,7 +301,8 @@ var SYS_SHORT = {
 var C_LIGHT = 299792458;
 var FREQ = {
   G: { "1": 157542e4, "2": 12276e5, "5": 117645e4 },
-  R: { "3": 1202025e3 },
+  // R bands 4/6 are the CDMA L1OC/L2OC center frequencies
+  R: { "3": 1202025e3, "4": 1600995e3, "6": 124806e4 },
   E: {
     "1": 157542e4,
     "5": 117645e4,
@@ -392,11 +396,11 @@ function buildGloChannelMap(slots) {
   return map;
 }
 function gloFreq(gloChannels, prn, band) {
+  if (band === "3" || band === "4" || band === "6") return FREQ.R[band];
   const k = gloChannels[prn];
   if (k === void 0) return void 0;
   if (band === "1") return GLO_F1_BASE + k * GLO_F1_STEP;
   if (band === "2") return GLO_F2_BASE + k * GLO_F2_STEP;
-  if (band === "3") return GLO_F3;
   return void 0;
 }
 function getFreq(gloChannels, prn, band) {
@@ -2556,6 +2560,13 @@ function writeRinexNav(nav) {
 }
 
 // src/rtcm3/decoder.ts
+var debugHandler = null;
+function setRtcm3DebugHandler(fn) {
+  debugHandler = fn;
+}
+function reportDecodeError(context, error) {
+  if (debugHandler) debugHandler(context, error);
+}
 var BitReader = class {
   data;
   bitPos = 0;
@@ -3090,31 +3101,32 @@ function decodeEphemeris(frame) {
       default:
         return null;
     }
-  } catch {
+  } catch (err) {
+    reportDecodeError("ephemeris", err);
     return null;
   }
 }
 
 // src/rtcm3/msm.ts
-var GPS_L1 = 157542e4;
-var GPS_L2 = 12276e5;
-var GPS_L5 = 117645e4;
-var GLO_L1a = 1600995e3;
-var GLO_L2a = 124806e4;
-var GLO_L3 = 1202025e3;
-var GAL_E1 = 157542e4;
-var GAL_E5a = 117645e4;
-var GAL_E5b = 120714e4;
-var GAL_E5 = 1191795e3;
-var GAL_E6 = 127875e4;
-var BDS_B1 = 1561098e3;
-var BDS_B3 = 126852e4;
-var BDS_B2 = 120714e4;
-var BDS_B1C = 157542e4;
-var BDS_B2a = 117645e4;
-var BDS_B2b = 120714e4;
-var QZSS_L6 = 127875e4;
-var NAVIC_S = 2492028e3;
+var GPS_L1 = FREQ.G["1"];
+var GPS_L2 = FREQ.G["2"];
+var GPS_L5 = FREQ.G["5"];
+var GLO_L1a = FREQ.R["4"];
+var GLO_L2a = FREQ.R["6"];
+var GLO_L3 = FREQ.R["3"];
+var GAL_E1 = FREQ.E["1"];
+var GAL_E5a = FREQ.E["5"];
+var GAL_E5b = FREQ.E["7"];
+var GAL_E5 = FREQ.E["8"];
+var GAL_E6 = FREQ.E["6"];
+var BDS_B1 = FREQ.C["2"];
+var BDS_B3 = FREQ.C["6"];
+var BDS_B2 = FREQ.C["7"];
+var BDS_B1C = FREQ.C["1"];
+var BDS_B2a = FREQ.C["5"];
+var BDS_B2b = FREQ.C["7"];
+var QZSS_L6 = FREQ.J["6"];
+var NAVIC_S = FREQ.I["9"];
 var EMPTY = { code: "", freq: 0 };
 function sd(code, freq) {
   return { code, freq };
@@ -3366,40 +3378,6 @@ function signalTable(sys) {
       return GPS_SIGNALS;
   }
 }
-var BitReader2 = class {
-  data;
-  pos;
-  constructor(data, startBit = 0) {
-    this.data = data;
-    this.pos = startBit;
-  }
-  /** Read unsigned value of numBits bits (max 53) */
-  u(numBits) {
-    let val = 0;
-    for (let i = 0; i < numBits; i++) {
-      const byteIdx = this.pos + i >>> 3;
-      const bitIdx = 7 - (this.pos + i & 7);
-      const bit = byteIdx < this.data.length ? this.data[byteIdx] >>> bitIdx & 1 : 0;
-      val = val * 2 + bit;
-    }
-    this.pos += numBits;
-    return val;
-  }
-  /** Read signed value (two's complement) */
-  s(numBits) {
-    const raw = this.u(numBits);
-    const half = 2 ** (numBits - 1);
-    return raw >= half ? raw - 2 ** numBits : raw;
-  }
-  /** Skip bits */
-  skip(n) {
-    this.pos += n;
-  }
-  /** Current bit position */
-  get position() {
-    return this.pos;
-  }
-};
 var gloFreqNum = new Int8Array(64).fill(-128);
 function gloWavelength(satIdx, freqMarker) {
   const k = gloFreqNum[satIdx];
@@ -3434,28 +3412,9 @@ function lockTimeSec(msmType, lti) {
   if (variant <= 5) {
     return LTI_TABLE_4BIT[lti] ?? 0;
   }
-  if (lti < 64) return lti;
-  if (lti < 96) return (lti - 64) * 2 + 64;
-  if (lti < 128) return (lti - 96) * 4 + 128;
-  if (lti < 160) return (lti - 128) * 8 + 256;
-  if (lti < 192) return (lti - 160) * 16 + 512;
-  if (lti < 224) return (lti - 192) * 32 + 1024;
-  if (lti < 256) return (lti - 224) * 64 + 2048;
-  if (lti < 288) return (lti - 256) * 128 + 4096;
-  if (lti < 320) return (lti - 288) * 256 + 8192;
-  if (lti < 352) return (lti - 320) * 512 + 16384;
-  if (lti < 384) return (lti - 352) * 1024 + 32768;
-  if (lti < 416) return (lti - 384) * 2048 + 65536;
-  if (lti < 448) return (lti - 416) * 4096 + 131072;
-  if (lti < 480) return (lti - 448) * 8192 + 262144;
-  if (lti < 512) return (lti - 480) * 16384 + 524288;
-  if (lti < 544) return (lti - 512) * 32768 + 1048576;
-  if (lti < 576) return (lti - 544) * 65536 + 2097152;
-  if (lti < 608) return (lti - 576) * 131072 + 4194304;
-  if (lti < 640) return (lti - 608) * 262144 + 8388608;
-  if (lti < 672) return (lti - 640) * 524288 + 16777216;
-  if (lti < 704) return (lti - 672) * 1048576 + 33554432;
-  return (lti - 704) * 2097152 + 67108864;
+  if (lti < 64) return lti / 1e3;
+  const n = Math.min(20, Math.floor((lti - 64) / 32));
+  return ((lti - 64 - 32 * n) * 2 ** (n + 1) + 2 ** (n + 6)) / 1e3;
 }
 function msmSystem(type) {
   if (type >= 1071 && type <= 1077) return "G";
@@ -3481,10 +3440,10 @@ function decodeMsmFull(frame) {
   const payload = frame.payload;
   if (!payload || payload.length < 10) return null;
   try {
-    const bits = new BitReader2(payload);
+    const bits = new BitReader(payload);
     bits.skip(12);
     bits.skip(12);
-    const epochMs = bits.u(30);
+    const epochMs = bits.readU(30);
     bits.skip(1);
     bits.skip(3);
     bits.skip(7);
@@ -3492,8 +3451,8 @@ function decodeMsmFull(frame) {
     bits.skip(2);
     bits.skip(1);
     bits.skip(3);
-    const satMaskHi = bits.u(32) >>> 0;
-    const satMaskLo = bits.u(32) >>> 0;
+    const satMaskHi = bits.readU(32) >>> 0;
+    const satMaskLo = bits.readU(32) >>> 0;
     const satIndices = [];
     for (let i = 0; i < 32; i++) {
       if (satMaskHi & 1 << 31 - i) satIndices.push(i + 1);
@@ -3503,7 +3462,7 @@ function decodeMsmFull(frame) {
     }
     const numSat = satIndices.length;
     if (numSat === 0) return null;
-    const sigMask = bits.u(32) >>> 0;
+    const sigMask = bits.readU(32) >>> 0;
     const sigIndices = [];
     for (let i = 0; i < 32; i++) {
       if (sigMask & 1 << 31 - i) sigIndices.push(i);
@@ -3513,7 +3472,7 @@ function decodeMsmFull(frame) {
     const numCells = numSat * numSig;
     const cellMask = [];
     for (let i = 0; i < numCells; i++) {
-      cellMask.push(bits.u(1) === 1);
+      cellMask.push(bits.readU(1) === 1);
     }
     let activeCells = 0;
     for (const c of cellMask) if (c) activeCells++;
@@ -3523,13 +3482,13 @@ function decodeMsmFull(frame) {
     const extsat = new Int8Array(numSat);
     const rdop = new Float64Array(numSat);
     if (variant === 4 || variant === 6) {
-      for (let j = 0; j < numSat; j++) rrint[j] = bits.u(8);
-      for (let j = 0; j < numSat; j++) rrmod[j] = bits.u(10) / 1024;
+      for (let j = 0; j < numSat; j++) rrint[j] = bits.readU(8);
+      for (let j = 0; j < numSat; j++) rrmod[j] = bits.readU(10) / 1024;
     } else {
-      for (let j = 0; j < numSat; j++) rrint[j] = bits.u(8);
-      for (let j = 0; j < numSat; j++) extsat[j] = bits.u(4);
-      for (let j = 0; j < numSat; j++) rrmod[j] = bits.u(10) / 1024;
-      for (let j = 0; j < numSat; j++) rdop[j] = bits.s(14) * 1;
+      for (let j = 0; j < numSat; j++) rrint[j] = bits.readU(8);
+      for (let j = 0; j < numSat; j++) extsat[j] = bits.readU(4);
+      for (let j = 0; j < numSat; j++) rrmod[j] = bits.readU(10) / 1024;
+      for (let j = 0; j < numSat; j++) rdop[j] = bits.readS(14) * 1;
     }
     if (sys === "R") {
       for (let j = 0; j < numSat; j++) {
@@ -3549,31 +3508,31 @@ function decodeMsmFull(frame) {
     cp.fill(-1e30);
     dop.fill(-1e30);
     if (variant === 4) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(15) / 2 ** 24;
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(22) / 2 ** 29;
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(4);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(6);
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(6);
     } else if (variant === 5) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(15) / 2 ** 24;
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(22) / 2 ** 29;
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(4);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(6);
-      for (let i = 0; i < activeCells; i++) dop[i] = bits.s(15) * 1e-4;
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(6);
+      for (let i = 0; i < activeCells; i++) dop[i] = bits.readS(15) * 1e-4;
     } else if (variant === 6) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(20) / 2 ** 29;
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(24) / 2 ** 31;
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(10);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(10) / 16;
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(20) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(24) / 2 ** 31;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(10);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(10) / 16;
     } else {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(20) / 2 ** 29;
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(24) / 2 ** 31;
-      for (let i = 0; i < activeCells; i++) ll[i] = bits.u(10);
-      for (let i = 0; i < activeCells; i++) hc[i] = bits.u(1);
-      for (let i = 0; i < activeCells; i++) cnr[i] = bits.u(10) / 16;
-      for (let i = 0; i < activeCells; i++) dop[i] = bits.s(15) * 1e-4;
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(20) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(24) / 2 ** 31;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(10);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+      for (let i = 0; i < activeCells; i++) cnr[i] = bits.readU(10) / 16;
+      for (let i = 0; i < activeCells; i++) dop[i] = bits.readS(15) * 1e-4;
     }
     const sigTable = signalTable(sys);
     const observations = [];
@@ -3607,11 +3566,11 @@ function decodeMsmFull(frame) {
           wavelength
         };
         const psrVal = psr[cellIdx];
-        if (psrVal > -1 / (1 << 10)) {
+        if (psrVal > -(2 ** -10)) {
           signal.pseudorange = psrVal * C_LIGHT / 1e3 + roughRange_m;
         }
         const cpVal = cp[cellIdx];
-        if (cpVal > -1 / (1 << 8) && wavelength > 0) {
+        if (cpVal > -(2 ** -8) && wavelength > 0) {
           signal.phase = cpVal * C_LIGHT / 1e3 / wavelength + roughRange_m / wavelength;
         }
         if (variant === 5 || variant === 7) {
@@ -3639,12 +3598,13 @@ function decodeMsmFull(frame) {
       system: sys,
       observations
     };
-  } catch {
+  } catch (err) {
+    reportDecodeError("msm", err);
     return null;
   }
 }
-var GPS_EPOCH = Date.UTC(1980, 0, 6);
-var MS_PER_WEEK = 6048e5;
+var GPS_EPOCH = START_GPS_TIME.getTime();
+var MS_PER_WEEK = MILLISECONDS_IN_WEEK;
 function msmEpochToDate(sys, epochMs, refTime = /* @__PURE__ */ new Date()) {
   const refMs = refTime.getTime();
   if (sys === "R") {
@@ -3659,11 +3619,12 @@ function msmEpochToDate(sys, epochMs, refTime = /* @__PURE__ */ new Date()) {
     return new Date(t2);
   }
   if (sys === "C") {
-    const BDS_EPOCH = Date.UTC(2006, 0, 1) - 14e3;
+    const BDS_EPOCH = START_BDS_TIME.getTime();
     const weeksSinceEpoch2 = Math.floor((refMs - BDS_EPOCH) / MS_PER_WEEK);
     let t2 = BDS_EPOCH + weeksSinceEpoch2 * MS_PER_WEEK + epochMs;
     if (t2 - refMs > MS_PER_WEEK / 2) t2 -= MS_PER_WEEK;
     else if (refMs - t2 > MS_PER_WEEK / 2) t2 += MS_PER_WEEK;
+    t2 -= 18e3;
     return new Date(t2);
   }
   const weeksSinceEpoch = Math.floor((refMs - GPS_EPOCH) / MS_PER_WEEK);
@@ -3673,6 +3634,11 @@ function msmEpochToDate(sys, epochMs, refTime = /* @__PURE__ */ new Date()) {
   t -= 18e3;
   return new Date(t);
 }
+function setGloFreqNumber(slot, k) {
+  if (slot >= 1 && slot <= 64 && k >= -7 && k <= 13) {
+    gloFreqNum[slot - 1] = k;
+  }
+}
 function resetGloFreqCache() {
   gloFreqNum.fill(-128);
 }
@@ -3806,7 +3772,8 @@ function updateStationMeta(meta, frame) {
       default:
         return false;
     }
-  } catch {
+  } catch (err) {
+    reportDecodeError("station", err);
     return false;
   }
 }
@@ -3994,6 +3961,9 @@ function updateStreamStats(stats, frames, rawBytes) {
     const eph = decodeEphemeris(frame);
     if (eph) {
       stats.ephemerides.set(eph.prn, eph);
+      if (eph.freqChannel !== void 0 && eph.prn.startsWith("R")) {
+        setGloFreqNumber(parseInt(eph.prn.slice(1), 10), eph.freqChannel);
+      }
     }
     updateStationMeta(stats.stationMeta, frame);
   }
@@ -4013,6 +3983,8 @@ 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;
@@ -4050,7 +4022,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);
@@ -4231,9 +4203,10 @@ function computeSatPosition(eph, timeMs) {
   if (eph.system === "R" || eph.system === "S") {
     return glonassPosition(eph, timeMs / 1e3);
   }
-  const GPS_EPOCH2 = Date.UTC(1980, 0, 6);
+  const GPS_EPOCH2 = START_GPS_TIME.getTime();
   const gpsSeconds = (timeMs - GPS_EPOCH2) / 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) {
@@ -4315,8 +4288,8 @@ 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);
@@ -4357,6 +4330,9 @@ function ephInfoToEphemeris(info) {
   }
   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") {
@@ -5487,32 +5463,33 @@ function parseNmeaFile(content) {
 }
 
 // src/signals/definitions.ts
-var FREQ_GPS_L1 = 1575.42;
-var FREQ_GPS_L2 = 1227.6;
-var FREQ_GPS_L5 = 1176.45;
-var FREQ_GLO_L1 = 1602;
-var FREQ_GLO_L2 = 1246;
-var DELTA_GLO_L1 = 0.5625;
-var DELTA_GLO_L2 = 0.4375;
-var FREQ_GLO_L1OC = 1600.995;
-var FREQ_GLO_L2OC = 1248.06;
-var FREQ_GLO_L3OC = 1202.025;
-var FREQ_GAL_E1 = 1575.42;
-var FREQ_GAL_E5 = 1191.795;
-var FREQ_GAL_E5a = 1176.45;
-var FREQ_GAL_E5b = 1207.14;
-var FREQ_GAL_E6 = 1278.75;
-var FREQ_BDS_B1I = 1561.098;
-var FREQ_BDS_B2I = 1207.14;
-var FREQ_BDS_B3I = 1268.52;
-var FREQ_BDS_B1A = 1575.42;
-var FREQ_BDS_B1C = 1575.42;
-var FREQ_BDS_B3A = 1268.52;
-var FREQ_QZS_L1 = 1575.42;
-var FREQ_QZS_L2 = 1227.6;
-var FREQ_QZS_L5 = 1176.45;
-var FREQ_QZS_L6 = 1278.75;
-var FREQ_NAVIC_L5 = 1176.45;
+var MHz = 1e6;
+var FREQ_GPS_L1 = FREQ.G["1"] / MHz;
+var FREQ_GPS_L2 = FREQ.G["2"] / MHz;
+var FREQ_GPS_L5 = FREQ.G["5"] / MHz;
+var FREQ_GLO_L1 = GLO_F1_BASE / MHz;
+var FREQ_GLO_L2 = GLO_F2_BASE / MHz;
+var DELTA_GLO_L1 = GLO_F1_STEP / MHz;
+var DELTA_GLO_L2 = GLO_F2_STEP / MHz;
+var FREQ_GLO_L1OC = FREQ.R["4"] / MHz;
+var FREQ_GLO_L2OC = FREQ.R["6"] / MHz;
+var FREQ_GLO_L3OC = FREQ.R["3"] / MHz;
+var FREQ_GAL_E1 = FREQ.E["1"] / MHz;
+var FREQ_GAL_E5 = FREQ.E["8"] / MHz;
+var FREQ_GAL_E5a = FREQ.E["5"] / MHz;
+var FREQ_GAL_E5b = FREQ.E["7"] / MHz;
+var FREQ_GAL_E6 = FREQ.E["6"] / MHz;
+var FREQ_BDS_B1I = FREQ.C["2"] / MHz;
+var FREQ_BDS_B2I = FREQ.C["7"] / MHz;
+var FREQ_BDS_B3I = FREQ.C["6"] / MHz;
+var FREQ_BDS_B1A = FREQ.C["1"] / MHz;
+var FREQ_BDS_B1C = FREQ.C["1"] / MHz;
+var FREQ_BDS_B3A = FREQ.C["6"] / MHz;
+var FREQ_QZS_L1 = FREQ.J["1"] / MHz;
+var FREQ_QZS_L2 = FREQ.J["2"] / MHz;
+var FREQ_QZS_L5 = FREQ.J["5"] / MHz;
+var FREQ_QZS_L6 = FREQ.J["6"] / MHz;
+var FREQ_NAVIC_L5 = FREQ.I["5"] / MHz;
 var GPS_COLOR = "#4ade80";
 var GPS_M_COLOR = "#ef4444";
 var GPS_P_COLOR = "#dbc51f";
@@ -6438,10 +6415,13 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   phiBPSK,
   rad2deg,
   readString,
+  reportDecodeError,
   resetGloFreqCache,
   rhumbLine,
   rtcm3Constellation,
   selectEphemeris,
+  setGloFreqNumber,
+  setRtcm3DebugHandler,
   systemCmp,
   systemName,
   updateStationMeta,

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-C10stHhM.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-C10stHhM.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';