gnss-js 1.15.0 → 1.17.0

package/dist/ubx.d.ts

@@ -1,3 +1,122 @@
+import { E as Ephemeris } from './nav-BAI1a9vK.js';
+
+/**
+ * u-blox UBX framing shared by the measurement and navigation decoders:
+ * 0xB5 0x62 sync, message class and ID, little-endian payload length,
+ * payload, and an 8-bit Fletcher checksum over class..payload.
+ */
+/** One checksum-valid UBX frame located in a byte stream. */
+interface UbxFrame {
+    /** Message class, e.g. 0x02 for RXM. */
+    msgClass: number;
+    /** Message ID within the class, e.g. 0x15 for RXM-RAWX. */
+    msgId: number;
+    /** Payload bytes (a view into the input buffer, not a copy). */
+    payload: Uint8Array;
+    /** Byte offset of the payload within the input buffer. */
+    payloadStart: number;
+}
+/**
+ * Iterate every checksum-valid UBX frame in `data`. Bad checksums are
+ * counted in `stats.badChecksums` and a resync continues at the next
+ * byte; a frame candidate that overruns the buffer ends the scan
+ * (truncated capture).
+ */
+declare function ubxFrames(data: Uint8Array, stats?: {
+    badChecksums: number;
+}): Generator<UbxFrame>;
+
+/**
+ * u-blox UBX navigation-message decoding: GPS/QZSS LNAV broadcast
+ * ephemerides from RXM-SFRBX (class 0x02, id 0x13).
+ *
+ * Ported from RTKLIB demo5 (rtklibexplorer fork, src/rcv/ublox.c:
+ * decode_rxmsfrbx / decode_nav / decode_eph, Copyright (c) 2007-2020
+ * T. Takasu, BSD-2-Clause); the subframe fields themselves are decoded
+ * by `decodeGpsLnavFrame` (the decode_frame port in src/navbits).
+ *
+ * On D30* polarity: the receiver has already checked parity and
+ * resolved the polarity inversion of the data bits before writing the
+ * 30-bit words into the 30 LSBs of each `dwrd`, so — exactly like
+ * RTKLIB's decode_nav (`U4(p)>>6`, "24 x 10 bits w/o parity") — the
+ * conversion to parity-stripped 24-bit words is a plain 6-bit shift
+ * with NO D30* re-inversion. (Verified against the TLM preamble 0x8B
+ * at bits 29..22 of word 1 in every GPS/QZSS LNAV message of the
+ * reference F9P capture.)
+ *
+ * Other constellations in RXM-SFRBX (Galileo I/NAV, BDS D1/D2, GLONASS
+ * strings, SBAS) and GPS/QZSS CNAV are out of scope and skipped
+ * silently. LNAV subframes 4/5 are skipped here except that
+ * `parseUbxIonoUtc` (./iono.ts) decodes the subframe-4 page-18
+ * iono/UTC parameters through the shared `readLnavSubframe` helper.
+ */
+
+interface UbxNavOptions {
+    /**
+     * Full GPS week used to resolve the 10-bit broadcast week. When
+     * omitted, the week is harvested from the first RXM-RAWX message in
+     * the same stream (u-blox raw logs carry both); a stream with
+     * neither yields no ephemerides — the system clock is never used.
+     */
+    refWeek?: number;
+}
+interface UbxNavResult {
+    /** Broadcast ephemerides in stream order, duplicates suppressed. */
+    ephemerides: Ephemeris[];
+    /**
+     * Subframes rejected for inconsistency: an out-of-range subframe ID
+     * in the HOW word, or a complete subframe-1/2/3 set whose
+     * IODC/IODE cross-check failed (mixed issues of data).
+     */
+    badParity: number;
+}
+/**
+ * Decode every GPS/QZSS LNAV ephemeris broadcast in a UBX byte stream
+ * (RXM-SFRBX). Subframes 1–3 are assembled per satellite; a frame is
+ * decoded when subframe 3 arrives with 1 and 2 buffered (RTKLIB
+ * decode_nav), and repeated broadcasts of an unchanged ephemeris are
+ * suppressed by issue of data and clock epoch like `parseNovatelNav`.
+ */
+declare function parseUbxNav(data: Uint8Array, opts?: UbxNavOptions): UbxNavResult;
+
+/**
+ * u-blox UBX ionosphere / UTC decoding: the GPS LNAV subframe 4
+ * page 18 broadcast (Klobuchar coefficients + GPS-UTC parameters)
+ * from RXM-SFRBX.
+ *
+ * Ported from RTKLIB demo5 (rtklibexplorer fork, src/rcvraw.c:
+ * decode_frame_ion / decode_frame_utc, Copyright (c) 2007-2020
+ * T. Takasu, BSD-2-Clause); frame extraction is shared with
+ * `parseUbxNav` through `readLnavSubframe` (./nav.ts).
+ *
+ * Scope is GPS satellites (gnssId 0) only: QZSS also broadcasts
+ * page-18-format iono parameters, but those are the separate
+ * wide-area/Japan-area QZSS coefficient set, not the GPS one a RINEX
+ * `GPSA`/`GPSB` header carries. Scale factors follow IS-GPS-200
+ * §20.3.3.5.1 (alpha_n in s/semicircle^n, beta_n in s/semicircle^n) —
+ * the same semicircle-based units the RINEX nav header prints, so the
+ * output matches `parseNavFile` on a converted header. Like
+ * `parseUbxNav`, the decoder never consults the system clock (nothing
+ * here needs a week reference).
+ */
+interface UbxIonoUtcResult {
+    /**
+     * Iono coefficient sets keyed like `NavHeader.ionoCorrections`:
+     * `GPSA` (alpha_0..3) and `GPSB` (beta_0..3). The page repeats
+     * (nominally every 12.5 min per satellite); the last broadcast in
+     * the stream wins.
+     */
+    ionoCorrections: Record<string, number[]>;
+    /** GPS-UTC ΔtLS from the last page-18 broadcast, if any. */
+    leapSeconds: number | null;
+}
+/**
+ * Decode every GPS LNAV subframe 4 page 18 (data ID 1, SV ID 56) in a
+ * UBX byte stream (RXM-SFRBX) into Klobuchar iono coefficients and the
+ * current leap-second count.
+ */
+declare function parseUbxIonoUtc(data: Uint8Array): UbxIonoUtcResult;
+
 /**
  * u-blox UBX raw-measurement decoding (RXM-RAWX) — the path from a
  * receiver log to RINEX-grade observables.
@@ -12,6 +131,7 @@
  * RINEX attributes chosen to match RTKLIB's u-blox conversion (the de
  * facto reference: GPS L2 CL/CM as 2X, Galileo E1 as 1X, E5b as 7X).
  */
+
 interface UbxMeasurement {
     /** RINEX PRN, e.g. "G04", "R11", "S23". */
     prn: string;
@@ -53,4 +173,4 @@ interface UbxParseResult {
  */
 declare function parseUbxRawx(data: Uint8Array): UbxParseResult;
 
-export { type UbxMeasurement, type UbxParseResult, type UbxRawxEpoch, parseUbxRawx };
+export { type UbxFrame, type UbxIonoUtcResult, type UbxMeasurement, type UbxNavOptions, type UbxNavResult, type UbxParseResult, type UbxRawxEpoch, parseUbxIonoUtc, parseUbxNav, parseUbxRawx, ubxFrames };

package/dist/ubx.js

@@ -1,6 +1,13 @@
 import {
-  parseUbxRawx
-} from "./chunk-5IXK25MX.js";
+  parseUbxIonoUtc,
+  parseUbxNav,
+  parseUbxRawx,
+  ubxFrames
+} from "./chunk-IG5CDNDS.js";
+import "./chunk-REYOYF7O.js";
 export {
-  parseUbxRawx
+  parseUbxIonoUtc,
+  parseUbxNav,
+  parseUbxRawx,
+  ubxFrames
 };

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.15.0",
+  "version": "1.17.0",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,

package/dist/chunk-5IXK25MX.js

@@ -1,116 +0,0 @@
-// src/ubx/index.ts
-var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
-var MS_PER_WEEK = 7 * 864e5;
-var SIGNALS = {
-  0: { 0: ["G", "1C"], 3: ["G", "2X"], 4: ["G", "2X"] },
-  // GPS L1C/A, L2CL, L2CM
-  1: { 0: ["S", "1C"] },
-  // SBAS L1C/A
-  2: { 0: ["E", "1X"], 1: ["E", "1X"], 5: ["E", "7X"], 6: ["E", "7X"] },
-  // E1C/B, E5bI/Q
-  3: { 0: ["C", "2I"], 1: ["C", "2I"], 2: ["C", "7I"], 3: ["C", "7I"] },
-  // B1I D1/D2, B2I D1/D2
-  5: { 0: ["J", "1C"], 4: ["J", "2X"], 5: ["J", "2X"] },
-  // QZSS L1C/A, L2CM/CL
-  6: { 0: ["R", "1C"], 2: ["R", "2C"] }
-  // GLONASS L1OF, L2OF
-};
-function prnFor(gnssId, svId) {
-  const two = (n) => String(n).padStart(2, "0");
-  switch (gnssId) {
-    case 0:
-      return svId >= 1 && svId <= 32 ? `G${two(svId)}` : null;
-    case 1:
-      return svId >= 120 && svId <= 158 ? `S${two(svId - 100)}` : null;
-    case 2:
-      return svId >= 1 && svId <= 36 ? `E${two(svId)}` : null;
-    case 3:
-      return svId >= 1 && svId <= 63 ? `C${two(svId)}` : null;
-    case 5:
-      return svId >= 1 && svId <= 10 ? `J${two(svId)}` : null;
-    case 6:
-      return svId >= 1 && svId <= 32 ? `R${two(svId)}` : null;
-    default:
-      return null;
-  }
-}
-function parseUbxRawx(data) {
-  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
-  const epochs = [];
-  const messageCounts = {};
-  const obsCodes = {};
-  let badChecksums = 0;
-  let i = 0;
-  while (i + 8 <= data.length) {
-    if (data[i] !== 181 || data[i + 1] !== 98) {
-      i++;
-      continue;
-    }
-    const cls = data[i + 2];
-    const id = data[i + 3];
-    const len = data[i + 4] | data[i + 5] << 8;
-    const end = i + 6 + len + 2;
-    if (end > data.length) break;
-    let ckA = 0;
-    let ckB = 0;
-    for (let j = i + 2; j < i + 6 + len; j++) {
-      ckA = ckA + data[j] & 255;
-      ckB = ckB + ckA & 255;
-    }
-    if (ckA !== data[i + 6 + len] || ckB !== data[i + 7 + len]) {
-      badChecksums++;
-      i++;
-      continue;
-    }
-    const key = `${cls.toString(16).padStart(2, "0")}-${id.toString(16).padStart(2, "0")}`;
-    messageCounts[key] = (messageCounts[key] ?? 0) + 1;
-    if (cls === 2 && id === 21 && len >= 16) {
-      const p = i + 6;
-      const rcvTow = view.getFloat64(p, true);
-      const week = view.getUint16(p + 8, true);
-      const leapS = view.getInt8(p + 10);
-      const numMeas = data[p + 11];
-      const recStat = data[p + 12];
-      const leapValid = (recStat & 1) !== 0;
-      const meas = [];
-      for (let k = 0; k < numMeas; k++) {
-        const off = p + 16 + 32 * k;
-        if (off + 32 > p + len) break;
-        const gnssId = data[off + 20];
-        const svId = data[off + 21];
-        const sigId = data[off + 22];
-        const trkStat = data[off + 30];
-        const prn = prnFor(gnssId, svId);
-        const sig = SIGNALS[gnssId]?.[sigId];
-        if (!prn || !sig) continue;
-        const prValid = (trkStat & 1) !== 0;
-        const cpValid = (trkStat & 2) !== 0;
-        const cp = view.getFloat64(off + 8, true);
-        meas.push({
-          prn,
-          code: sig[1],
-          pr: prValid ? view.getFloat64(off, true) : null,
-          cp: cpValid && cp !== 0 ? cp : null,
-          doppler: view.getFloat32(off + 16, true),
-          cn0: data[off + 26],
-          halfCycleAmbiguous: (trkStat & 4) !== 0,
-          lockTimeMs: view.getUint16(off + 24, true)
-        });
-        const sys = sig[0];
-        const codes = obsCodes[sys] ??= [];
-        if (!codes.includes(sig[1])) codes.push(sig[1]);
-      }
-      epochs.push({
-        timeMs: GPS_EPOCH_MS + week * MS_PER_WEEK + Math.round(rcvTow * 1e3),
-        leapS: leapValid ? leapS : null,
-        meas
-      });
-    }
-    i = end;
-  }
-  return { epochs, messageCounts, obsCodes, badChecksums };
-}
-
-export {
-  parseUbxRawx
-};