gnss-js 1.0.0 → 1.1.0

package/dist/{chunk-SDRRAJT5.js → chunk-PGHDJFQK.js}

@@ -1,6 +1,10 @@
 import {
   ecefToGeodetic
 } from "./chunk-37QNKGTC.js";
+import {
+  START_BDS_TIME,
+  START_GPS_TIME
+} from "./chunk-LEEU5OIO.js";
 
 // src/orbit/index.ts
 var GM_GPS = 3986005e8;
@@ -8,6 +12,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;
@@ -33,10 +39,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);
@@ -48,7 +51,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);
@@ -199,7 +202,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;
@@ -226,9 +232,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) {
@@ -310,16 +317,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;
@@ -344,7 +357,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/{chunk-IS4UUDBV.js → chunk-W4YMQKWH.js}

@@ -39,10 +39,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)");
@@ -51,7 +59,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) {
@@ -118,7 +130,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",
@@ -208,7 +225,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"));
 }
@@ -415,9 +440,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++;
@@ -453,9 +494,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++;

package/dist/{chunk-354IRDOG.js → chunk-W5WKEV7U.js}

@@ -20,9 +20,22 @@ var SYS_SHORT = {
 var C_LIGHT = 299792458;
 var FREQ = {
   G: { "1": 157542e4, "2": 12276e5, "5": 117645e4 },
-  R: { "3": 1202025e3 },
-  E: { "1": 157542e4, "5": 117645e4, "6": 127875e4, "7": 120714e4, "8": 1191795e3 },
-  C: { "1": 157542e4, "2": 1561098e3, "5": 117645e4, "6": 126852e4, "7": 120714e4 },
+  // R bands 4/6 are the CDMA L1OC/L2OC center frequencies
+  R: { "3": 1202025e3, "4": 1600995e3, "6": 124806e4 },
+  E: {
+    "1": 157542e4,
+    "5": 117645e4,
+    "6": 127875e4,
+    "7": 120714e4,
+    "8": 1191795e3
+  },
+  C: {
+    "1": 157542e4,
+    "2": 1561098e3,
+    "5": 117645e4,
+    "6": 126852e4,
+    "7": 120714e4
+  },
   J: { "1": 157542e4, "2": 12276e5, "5": 117645e4, "6": 127875e4 },
   I: { "5": 117645e4, "9": 2492028e3 },
   S: { "1": 157542e4, "5": 117645e4 }
@@ -37,11 +50,28 @@ var BAND_LABELS = {
   S: { "1": "L1", "5": "L5" }
 };
 var DUAL_FREQ_PAIRS = {
-  G: [["1", "2"], ["1", "5"]],
-  R: [["1", "2"], ["1", "3"]],
-  E: [["1", "5"], ["1", "7"], ["1", "6"]],
-  C: [["2", "7"], ["2", "6"], ["1", "5"]],
-  J: [["1", "2"], ["1", "5"]],
+  G: [
+    ["1", "2"],
+    ["1", "5"]
+  ],
+  R: [
+    ["1", "2"],
+    ["1", "3"]
+  ],
+  E: [
+    ["1", "5"],
+    ["1", "7"],
+    ["1", "6"]
+  ],
+  C: [
+    ["2", "7"],
+    ["2", "6"],
+    ["1", "5"]
+  ],
+  J: [
+    ["1", "2"],
+    ["1", "5"]
+  ],
   I: [["5", "9"]],
   S: [["1", "5"]]
 };
@@ -85,11 +115,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) {

package/dist/constants.cjs

@@ -148,9 +148,22 @@ var SYS_SHORT = {
 var C_LIGHT = 299792458;
 var FREQ = {
   G: { "1": 157542e4, "2": 12276e5, "5": 117645e4 },
-  R: { "3": 1202025e3 },
-  E: { "1": 157542e4, "5": 117645e4, "6": 127875e4, "7": 120714e4, "8": 1191795e3 },
-  C: { "1": 157542e4, "2": 1561098e3, "5": 117645e4, "6": 126852e4, "7": 120714e4 },
+  // R bands 4/6 are the CDMA L1OC/L2OC center frequencies
+  R: { "3": 1202025e3, "4": 1600995e3, "6": 124806e4 },
+  E: {
+    "1": 157542e4,
+    "5": 117645e4,
+    "6": 127875e4,
+    "7": 120714e4,
+    "8": 1191795e3
+  },
+  C: {
+    "1": 157542e4,
+    "2": 1561098e3,
+    "5": 117645e4,
+    "6": 126852e4,
+    "7": 120714e4
+  },
   J: { "1": 157542e4, "2": 12276e5, "5": 117645e4, "6": 127875e4 },
   I: { "5": 117645e4, "9": 2492028e3 },
   S: { "1": 157542e4, "5": 117645e4 }
@@ -165,11 +178,28 @@ var BAND_LABELS = {
   S: { "1": "L1", "5": "L5" }
 };
 var DUAL_FREQ_PAIRS = {
-  G: [["1", "2"], ["1", "5"]],
-  R: [["1", "2"], ["1", "3"]],
-  E: [["1", "5"], ["1", "7"], ["1", "6"]],
-  C: [["2", "7"], ["2", "6"], ["1", "5"]],
-  J: [["1", "2"], ["1", "5"]],
+  G: [
+    ["1", "2"],
+    ["1", "5"]
+  ],
+  R: [
+    ["1", "2"],
+    ["1", "3"]
+  ],
+  E: [
+    ["1", "5"],
+    ["1", "7"],
+    ["1", "6"]
+  ],
+  C: [
+    ["2", "7"],
+    ["2", "6"],
+    ["1", "5"]
+  ],
+  J: [
+    ["1", "2"],
+    ["1", "5"]
+  ],
   I: [["5", "9"]],
   S: [["1", "5"]]
 };
@@ -213,11 +243,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) {

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-354IRDOG.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/coordinates.cjs

@@ -156,7 +156,10 @@ function vincenty(lat1, lon1, lat2, lon2) {
     const cosLat2 = Math.cos(lat2), sinLat2 = Math.sin(lat2);
     const dLon = lon2 - lon1;
     sigma = Math.acos(
-      Math.min(1, Math.max(-1, sinLat1 * sinLat2 + cosLat1 * cosLat2 * Math.cos(dLon)))
+      Math.min(
+        1,
+        Math.max(-1, sinLat1 * sinLat2 + cosLat1 * cosLat2 * Math.cos(dLon))
+      )
     );
     const distance2 = WGS84_SEMI_MAJOR_AXIS * sigma;
     const ib = Math.atan2(

package/dist/coordinates.js

@@ -10,7 +10,7 @@ import {
   rad2deg,
   rhumbLine,
   vincenty
-} from "./chunk-K7WZQFBV.js";
+} from "./chunk-MIIM4LDY.js";
 import {
   clampUnit,
   ecefToGeodetic,

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

@@ -7,12 +7,22 @@ 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;
     private bitPos;
     constructor(data: Uint8Array);
-    /** Read `n` bits as unsigned integer (max 32). */
+    /** Read `n` bits as unsigned integer (max 53). */
     readU(n: number): number;
     /** Read `n` bits as signed (two's complement) integer. */
     readS(n: number): number;
@@ -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-BUWzfmBy.d.ts → ephemeris-Ohl6NAfw.d.ts}

@@ -7,12 +7,22 @@ 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;
     private bitPos;
     constructor(data: Uint8Array);
-    /** Read `n` bits as unsigned integer (max 32). */
+    /** Read `n` bits as unsigned integer (max 53). */
     readU(n: number): number;
     /** Read `n` bits as signed (two's complement) integer. */
     readS(n: number): number;
@@ -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 };