gnss-js 1.0.0 → 1.0.1

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

@@ -21,8 +21,20 @@ 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 },
+  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 +49,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"]]
 };

package/dist/{chunk-K7WZQFBV.js → chunk-MIIM4LDY.js}

@@ -63,7 +63,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/{chunk-Y3R57B5P.js → chunk-SO3POWWR.js}

@@ -1,6 +1,6 @@
 import {
   parseRinexStream
-} from "./chunk-CB6EOOLA.js";
+} from "./chunk-HBLU2EJ4.js";
 import {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -11,7 +11,7 @@ import {
   buildGloChannelMap,
   buildObsIndices,
   getFreq
-} from "./chunk-354IRDOG.js";
+} from "./chunk-LWNTWBHB.js";
 
 // src/analysis/multipath.ts
 var MIN_ARC_LENGTH = 10;
@@ -99,8 +99,18 @@ var MultipathAccumulator = class {
         time: t,
         mp: arc.rawMp[i] - mean
       }));
-      const rms = Math.sqrt(points.reduce((s, p) => s + p.mp * p.mp, 0) / points.length);
-      this.results.push({ prn, system: sys, band, refBand, label, points, rms });
+      const rms = Math.sqrt(
+        points.reduce((s, p) => s + p.mp * p.mp, 0) / points.length
+      );
+      this.results.push({
+        prn,
+        system: sys,
+        band,
+        refBand,
+        label,
+        points,
+        rms
+      });
     }
     state.arc = { times: [], rawMp: [] };
   }
@@ -254,10 +264,20 @@ var CycleSlipAccumulator = class {
           const sigKey = `${sys}:${signal}`;
           this.countEpoch(sigKey);
           const stddev = mwSmoothStddev(prevSmooth);
-          const threshold = Math.max(MW_MIN_THRESHOLD, MW_SIGMA_FACTOR * stddev);
+          const threshold = Math.max(
+            MW_MIN_THRESHOLD,
+            MW_SIGMA_FACTOR * stddev
+          );
           const deviation = Math.abs(mw - prevSmooth.mean);
           if (deviation > threshold) {
-            this.recordSlip(sigKey, time, prn, signal, deviation, /* @__PURE__ */ new Set([b1, b2]));
+            this.recordSlip(
+              sigKey,
+              time,
+              prn,
+              signal,
+              deviation,
+              /* @__PURE__ */ new Set([b1, b2])
+            );
             mwSmooth.set(pairKey, mwSmoothInit(mw));
           } else {
             mwSmooth.set(pairKey, mwSmoothUpdate(prevSmooth, mw));
@@ -312,7 +332,13 @@ var CycleSlipAccumulator = class {
         }
       }
     }
-    this.prev.set(prn, { time, mwSmooth, gf: gfMap, phase: phaseM, code: codeM });
+    this.prev.set(prn, {
+      time,
+      mwSmooth,
+      gf: gfMap,
+      phase: phaseM,
+      code: codeM
+    });
   };
   finalize() {
     const satSlipCounts = {};
@@ -320,7 +346,10 @@ var CycleSlipAccumulator = class {
       satSlipCounts[ev.prn] = (satSlipCounts[ev.prn] ?? 0) + 1;
     }
     const signalStats = [];
-    const allKeys = /* @__PURE__ */ new Set([...this.signalEpochs.keys(), ...this.signalSlips.keys()]);
+    const allKeys = /* @__PURE__ */ new Set([
+      ...this.signalEpochs.keys(),
+      ...this.signalSlips.keys()
+    ]);
     for (const sigKey of allKeys) {
       const [sys, ...rest] = sigKey.split(":");
       const signal = rest.join(":");

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-PRSZIWKM.js → chunk-WP2JFDLA.js}

@@ -4,7 +4,7 @@ import {
   GLO_F1_STEP,
   GLO_F2_BASE,
   GLO_F2_STEP
-} from "./chunk-354IRDOG.js";
+} from "./chunk-LWNTWBHB.js";
 
 // src/rtcm3/decoder.ts
 var BitReader = class {
@@ -13,21 +13,18 @@ var BitReader = class {
   constructor(data) {
     this.data = data;
   }
-  /** Read `n` bits as unsigned integer (max 32). */
+  /** Read `n` bits as unsigned integer (max 53). */
   readU(n) {
     if (n === 0) return 0;
     let val = 0;
     for (let i = 0; i < n; i++) {
       const byteIdx = this.bitPos + i >> 3;
       const bitIdx = 7 - (this.bitPos + i & 7);
-      if (byteIdx < this.data.length) {
-        val = val << 1 | this.data[byteIdx] >> bitIdx & 1;
-      } else {
-        val <<= 1;
-      }
+      const bit = byteIdx < this.data.length ? this.data[byteIdx] >> bitIdx & 1 : 0;
+      val = val * 2 + bit;
     }
     this.bitPos += n;
-    return val >>> 0;
+    return val;
   }
   /** Read `n` bits as signed (two's complement) integer. */
   readS(n) {
@@ -38,9 +35,9 @@ var BitReader = class {
   /** Read `n` bits as sign-magnitude integer (MSB = sign, rest = magnitude). */
   readSM(n) {
     const val = this.readU(n);
-    const sign = val >> n - 1;
-    const mag = val & (1 << n - 1) - 1;
-    return sign ? -mag : mag;
+    const half = 2 ** (n - 1);
+    const mag = val % half;
+    return val >= half ? -mag : mag;
   }
   /** Skip `n` bits. */
   skip(n) {
@@ -827,17 +824,14 @@ var BitReader2 = class {
     this.data = data;
     this.pos = startBit;
   }
-  /** Read unsigned value of numBits bits */
+  /** 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);
-      if (byteIdx < this.data.length) {
-        val = val << 1 | this.data[byteIdx] >>> bitIdx & 1;
-      } else {
-        val <<= 1;
-      }
+      const bit = byteIdx < this.data.length ? this.data[byteIdx] >>> bitIdx & 1 : 0;
+      val = val * 2 + bit;
     }
     this.pos += numBits;
     return val;
@@ -845,8 +839,8 @@ var BitReader2 = class {
   /** Read signed value (two's complement) */
   s(numBits) {
     const raw = this.u(numBits);
-    const sign = 1 << numBits - 1;
-    return (raw & sign - 1) - (raw & sign);
+    const half = 2 ** (numBits - 1);
+    return raw >= half ? raw - 2 ** numBits : raw;
   }
   /** Skip bits */
   skip(n) {
@@ -981,7 +975,6 @@ function decodeMsmFull(frame) {
     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++) extsat[j] = bits.u(4);
       for (let j = 0; j < numSat; j++) rrmod[j] = bits.u(10) / 1024;
     } else {
       for (let j = 0; j < numSat; j++) rrint[j] = bits.u(8);
@@ -1007,27 +1000,27 @@ function decodeMsmFull(frame) {
     cp.fill(-1e30);
     dop.fill(-1e30);
     if (variant === 4) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(15) / (1 << 24);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(22) / (1 << 29);
+      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);
     } else if (variant === 5) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(15) / (1 << 24);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(22) / (1 << 29);
+      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;
     } else if (variant === 6) {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(20) / (1 << 29);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(24) / (1 << 31);
+      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;
     } else {
-      for (let i = 0; i < activeCells; i++) psr[i] = bits.s(20) / (1 << 29);
-      for (let i = 0; i < activeCells; i++) cp[i] = bits.s(24) / (1 << 31);
+      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;

package/dist/constants.cjs

@@ -149,8 +149,20 @@ 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 },
+  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 +177,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"]]
 };

package/dist/constants.js

@@ -45,7 +45,7 @@ import {
   formatUTCTime,
   getFreq,
   gloFreq
-} from "./chunk-354IRDOG.js";
+} from "./chunk-LWNTWBHB.js";
 import {
   WGS84_ECCENTRICITY_SQUARED,
   WGS84_FLATTENING,

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-C10stHhM.d.cts}

@@ -12,7 +12,7 @@ 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;

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

@@ -12,7 +12,7 @@ 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;