gnss-js 1.16.0 → 1.17.0

package/dist/ubx.cjs

@@ -20,6 +20,7 @@ var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: tru
 // src/ubx/index.ts
 var ubx_exports = {};
 __export(ubx_exports, {
+  parseUbxIonoUtc: () => parseUbxIonoUtc,
   parseUbxNav: () => parseUbxNav,
   parseUbxRawx: () => parseUbxRawx,
   ubxFrames: () => ubxFrames
@@ -201,6 +202,30 @@ function decodeGpsLnavFrame(subframes, opts = {}) {
 
 // src/ubx/nav.ts
 var PREAMB_CNAV = 139;
+function readLnavSubframe(view, f) {
+  const p = f.payload;
+  if (p.length < 8) return null;
+  const gnssId = p[0];
+  const svId = p[1];
+  let prn;
+  if (gnssId === 0 && svId >= 1 && svId <= 32) {
+    prn = `G${String(svId).padStart(2, "0")}`;
+  } else if (gnssId === 5 && svId >= 1 && svId <= 10) {
+    if (p.length === 44) return null;
+    prn = `J${String(svId).padStart(2, "0")}`;
+  } else {
+    return null;
+  }
+  if (p.length < 8 + 40) return null;
+  const base = f.payloadStart + 8;
+  if (view.getUint32(base, true) >>> 24 === PREAMB_CNAV) return null;
+  const buff = new Uint8Array(30);
+  for (let k = 0; k < 10; k++) {
+    const dwrd = view.getUint32(base + 4 * k, true);
+    setBitU(buff, 24 * k, 24, dwrd >>> 6 & 16777215);
+  }
+  return { prn, gnssId, svId, buff, id: getBitU(buff, 43, 3) };
+}
 function parseUbxNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
@@ -222,28 +247,9 @@ function parseUbxNav(data, opts = {}) {
   const last = /* @__PURE__ */ new Map();
   for (const f of ubxFrames(data)) {
     if (f.msgClass !== 2 || f.msgId !== 19) continue;
-    const p = f.payload;
-    if (p.length < 8) continue;
-    const gnssId = p[0];
-    const svId = p[1];
-    let prn;
-    if (gnssId === 0 && svId >= 1 && svId <= 32) {
-      prn = `G${String(svId).padStart(2, "0")}`;
-    } else if (gnssId === 5 && svId >= 1 && svId <= 10) {
-      if (p.length === 44) continue;
-      prn = `J${String(svId).padStart(2, "0")}`;
-    } else {
-      continue;
-    }
-    if (p.length < 8 + 40) continue;
-    const base = f.payloadStart + 8;
-    if (view.getUint32(base, true) >>> 24 === PREAMB_CNAV) continue;
-    const buff = new Uint8Array(30);
-    for (let k = 0; k < 10; k++) {
-      const dwrd = view.getUint32(base + 4 * k, true);
-      setBitU(buff, 24 * k, 24, dwrd >>> 6 & 16777215);
-    }
-    const id = getBitU(buff, 43, 3);
+    const lnav = readLnavSubframe(view, f);
+    if (!lnav) continue;
+    const { prn, gnssId, svId, buff, id } = lnav;
     if (id < 1 || id > 5) {
       badParity++;
       continue;
@@ -273,6 +279,34 @@ function parseUbxNav(data, opts = {}) {
   return { ephemerides, badParity };
 }
 
+// src/ubx/iono.ts
+function parseUbxIonoUtc(data) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const ionoCorrections = {};
+  let leapSeconds = null;
+  for (const f of ubxFrames(data)) {
+    if (f.msgClass !== 2 || f.msgId !== 19) continue;
+    const sf = readLnavSubframe(view, f);
+    if (!sf || sf.gnssId !== 0 || sf.id !== 4) continue;
+    const b = sf.buff;
+    if (getBitU(b, 48, 2) !== 1 || getBitU(b, 50, 6) !== 56) continue;
+    ionoCorrections["GPSA"] = [
+      getBitS(b, 56, 8) * 2 ** -30,
+      getBitS(b, 64, 8) * 2 ** -27,
+      getBitS(b, 72, 8) * 2 ** -24,
+      getBitS(b, 80, 8) * 2 ** -24
+    ];
+    ionoCorrections["GPSB"] = [
+      getBitS(b, 88, 8) * 2 ** 11,
+      getBitS(b, 96, 8) * 2 ** 14,
+      getBitS(b, 104, 8) * 2 ** 16,
+      getBitS(b, 112, 8) * 2 ** 16
+    ];
+    leapSeconds = getBitS(b, 192, 8);
+  }
+  return { ionoCorrections, leapSeconds };
+}
+
 // src/ubx/index.ts
 var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK = 7 * 864e5;
@@ -367,6 +401,7 @@ function parseUbxRawx(data) {
 }
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
+  parseUbxIonoUtc,
   parseUbxNav,
   parseUbxRawx,
   ubxFrames

package/dist/ubx.d.cts

@@ -46,7 +46,9 @@ declare function ubxFrames(data: Uint8Array, stats?: {
  *
  * 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, as are LNAV subframes 4/5 (iono/UTC/almanac pages).
+ * 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 {
@@ -77,6 +79,44 @@ interface UbxNavResult {
  */
 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.
@@ -133,4 +173,4 @@ interface UbxParseResult {
  */
 declare function parseUbxRawx(data: Uint8Array): UbxParseResult;
 
-export { type UbxFrame, type UbxMeasurement, type UbxNavOptions, type UbxNavResult, type UbxParseResult, type UbxRawxEpoch, parseUbxNav, parseUbxRawx, ubxFrames };
+export { type UbxFrame, type UbxIonoUtcResult, type UbxMeasurement, type UbxNavOptions, type UbxNavResult, type UbxParseResult, type UbxRawxEpoch, parseUbxIonoUtc, parseUbxNav, parseUbxRawx, ubxFrames };

package/dist/ubx.d.ts

@@ -46,7 +46,9 @@ declare function ubxFrames(data: Uint8Array, stats?: {
  *
  * 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, as are LNAV subframes 4/5 (iono/UTC/almanac pages).
+ * 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 {
@@ -77,6 +79,44 @@ interface UbxNavResult {
  */
 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.
@@ -133,4 +173,4 @@ interface UbxParseResult {
  */
 declare function parseUbxRawx(data: Uint8Array): UbxParseResult;
 
-export { type UbxFrame, type UbxMeasurement, type UbxNavOptions, type UbxNavResult, type UbxParseResult, type UbxRawxEpoch, parseUbxNav, parseUbxRawx, ubxFrames };
+export { type UbxFrame, type UbxIonoUtcResult, type UbxMeasurement, type UbxNavOptions, type UbxNavResult, type UbxParseResult, type UbxRawxEpoch, parseUbxIonoUtc, parseUbxNav, parseUbxRawx, ubxFrames };

package/dist/ubx.js

@@ -1,10 +1,12 @@
 import {
+  parseUbxIonoUtc,
   parseUbxNav,
   parseUbxRawx,
   ubxFrames
-} from "./chunk-6CX6EXRK.js";
-import "./chunk-4LVD4DYL.js";
+} from "./chunk-IG5CDNDS.js";
+import "./chunk-REYOYF7O.js";
 export {
+  parseUbxIonoUtc,
   parseUbxNav,
   parseUbxRawx,
   ubxFrames

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.16.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,