gnss-js 2.11.0 → 2.13.0

package/dist/{chunk-RFPR4LKM.js → chunk-HWHD5E6V.js}

@@ -1,13 +1,19 @@
+import {
+  scanSbfFrames,
+  svidToPrn
+} from "./chunk-E3Y6YKQM.js";
 import {
   BdsAssembler,
-  GalFnavAssembler,
-  GalInavAssembler,
   GloStringAssembler,
   bdsSubframeParityOk,
-  galFnavPageCrcOk,
-  galInavPageCrcOk,
   testGloString
-} from "./chunk-SSWFX3P3.js";
+} from "./chunk-NIJ6WFQJ.js";
+import {
+  GalFnavAssembler,
+  GalInavAssembler,
+  galFnavPageCrcOk,
+  galInavPageCrcOk
+} from "./chunk-TWGZTRBL.js";
 import {
   decodeSbasGeoNav,
   sbasCrcOk
@@ -40,66 +46,6 @@ import {
   setBitU
 } from "./chunk-V2D5KUKH.js";
 
-// src/sbf/frame.ts
-var CRC_TABLE = (() => {
-  const t = new Uint16Array(256);
-  for (let n = 0; n < 256; n++) {
-    let c = n << 8;
-    for (let k = 0; k < 8; k++)
-      c = c & 32768 ? (c << 1 ^ 4129) & 65535 : c << 1 & 65535;
-    t[n] = c;
-  }
-  return t;
-})();
-function crc16(data, start, end) {
-  let crc = 0;
-  for (let i = start; i < end; i++)
-    crc = (crc << 8 ^ CRC_TABLE[(crc >> 8 ^ data[i]) & 255]) & 65535;
-  return crc;
-}
-function scanSbfFrames(data, view, onBlock) {
-  let badCrc = 0;
-  let i = 0;
-  while (i + 8 <= data.length) {
-    if (data[i] !== 36 || data[i + 1] !== 64) {
-      i++;
-      continue;
-    }
-    const len = view.getUint16(i + 6, true);
-    if (len < 8 || len % 4 !== 0) {
-      i++;
-      continue;
-    }
-    if (i + len > data.length) break;
-    if (crc16(data, i + 4, i + len) !== view.getUint16(i + 2, true)) {
-      badCrc++;
-      i++;
-      continue;
-    }
-    onBlock(view.getUint16(i + 4, true) & 8191, i, len);
-    i += len;
-  }
-  return badCrc;
-}
-var two = (n) => String(n).padStart(2, "0");
-function svidToPrn(svid) {
-  if (svid >= 1 && svid <= 37) return svid <= 32 ? `G${two(svid)}` : null;
-  if (svid <= 61) return svid - 37 <= 27 ? `R${two(svid - 37)}` : null;
-  if (svid <= 62) return null;
-  if (svid <= 68) return svid - 38 <= 27 ? `R${two(svid - 38)}` : null;
-  if (svid <= 70) return null;
-  if (svid <= 106) return svid - 70 <= 36 ? `E${two(svid - 70)}` : null;
-  if (svid <= 119) return null;
-  if (svid <= 140) return `S${two(svid - 100)}`;
-  if (svid <= 180) return svid - 140 <= 50 ? `C${two(svid - 140)}` : null;
-  if (svid <= 190) return `J${two(svid - 180)}`;
-  if (svid <= 197) return svid - 190 <= 14 ? `I${two(svid - 190)}` : null;
-  if (svid <= 215) return `S${two(svid - 157)}`;
-  if (svid <= 222) return svid - 208 <= 14 ? `I${two(svid - 208)}` : null;
-  if (svid <= 245) return svid - 182 <= 50 ? `C${two(svid - 182)}` : null;
-  return null;
-}
-
 // src/sbf/nav.ts
 var PI = Math.PI;
 var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
@@ -1981,23 +1927,23 @@ var MEAS2_SIG = [
   ["J", "1E"],
   ["J", "5P"]
 ];
-var two2 = (n) => String(n).padStart(2, "0");
+var two = (n) => String(n).padStart(2, "0");
 function m3Prn(navsys, svid) {
   switch (navsys) {
     case 0:
-      return svid < 32 ? `G${two2(svid + 1)}` : null;
+      return svid < 32 ? `G${two(svid + 1)}` : null;
     case 1:
-      return svid < 27 ? `R${two2(svid + 1)}` : null;
+      return svid < 27 ? `R${two(svid + 1)}` : null;
     case 2:
-      return svid < 36 ? `E${two2(svid + 1)}` : null;
+      return svid < 36 ? `E${two(svid + 1)}` : null;
     case 3:
-      return svid < 50 ? `C${two2(svid + 1)}` : null;
+      return svid < 50 ? `C${two(svid + 1)}` : null;
     case 4:
-      return svid <= 38 ? `S${two2(svid + 20)}` : null;
+      return svid <= 38 ? `S${two(svid + 20)}` : null;
     case 5:
-      return svid < 10 ? `J${two2(svid + 1)}` : null;
+      return svid < 10 ? `J${two(svid + 1)}` : null;
     case 6:
-      return svid < 14 ? `I${two2(svid + 1)}` : null;
+      return svid < 14 ? `I${two(svid + 1)}` : null;
     default:
       return null;
   }
@@ -2548,8 +2494,6 @@ function parseSbfMeas(data) {
 }
 
 export {
-  crc16,
-  scanSbfFrames,
   parseSbfNav,
   parseSbfAlmanac,
   parseSbfIonoUtc,

package/dist/{chunk-KUSUN7DP.js → chunk-MFCPYUEI.js}

@@ -1,6 +1,6 @@
 import {
   parseRinexStream
-} from "./chunk-3U5AX7PY.js";
+} from "./chunk-FEABDV5P.js";
 import {
   ARC_GAP_FACTOR,
   BAND_LABELS,

package/dist/{chunk-SSWFX3P3.js → chunk-NIJ6WFQJ.js}

@@ -1,9 +1,6 @@
 import {
   getUtcDate
 } from "./chunk-HVXYFUCB.js";
-import {
-  crc24q
-} from "./chunk-HUXQABJI.js";
 import {
   GPS_PI,
   getBitS,
@@ -288,317 +285,9 @@ var BdsAssembler = class {
   }
 };
 
-// src/navbits/gal.ts
+// src/navbits/glo.ts
 var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK2 = 7 * 86400;
-var HALF_WEEK2 = 302400;
-var GST_GPS_WEEK_OFFSET = 1024;
-function galInavPageCrcOk(page) {
-  if (page.length < 28) return false;
-  return crc24q(page, 196) === getBitU(page, 196, 24);
-}
-function galFnavPageCrcOk(page) {
-  if (page.length < 30) return false;
-  return crc24q(page, 214) === getBitU(page, 214, 24);
-}
-function gpsDate(week, sec) {
-  return new Date(GPS_EPOCH_MS + (week * SEC_PER_WEEK2 + sec) * 1e3);
-}
-function weekOfToe(week, tow, toe) {
-  if (toe - tow > HALF_WEEK2) return week - 1;
-  if (toe - tow < -HALF_WEEK2) return week + 1;
-  return week;
-}
-var two = (n) => String(n).padStart(2, "0");
-function decodeGalInavWords(words) {
-  if (words.length < 96) return null;
-  const b = words;
-  let i = 128;
-  const type1 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav1 = getBitU(b, i, 10);
-  i += 10;
-  const toes = getBitU(b, i, 14) * 60;
-  i += 14;
-  const m0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const e = getBitU(b, i, 32) * 2 ** -33;
-  i += 32;
-  const sqrtA = getBitU(b, i, 32) * 2 ** -19;
-  i = 128 * 2;
-  const type2 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav2 = getBitU(b, i, 10);
-  i += 10;
-  const omega0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const i0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const omega = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const idot = getBitS(b, i, 14) * 2 ** -43 * GPS_PI;
-  i = 128 * 3;
-  const type3 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav3 = getBitU(b, i, 10);
-  i += 10;
-  const omegaDot = getBitS(b, i, 24) * 2 ** -43 * GPS_PI;
-  i += 24;
-  const deltaN = getBitS(b, i, 16) * 2 ** -43 * GPS_PI;
-  i += 16;
-  const cuc = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const cus = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const crc = getBitS(b, i, 16) * 2 ** -5;
-  i += 16;
-  const crs = getBitS(b, i, 16) * 2 ** -5;
-  i = 128 * 4;
-  const type4 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav4 = getBitU(b, i, 10);
-  i += 10;
-  const svid = getBitU(b, i, 6);
-  i += 6;
-  const cic = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const cis = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const tocs = getBitU(b, i, 14) * 60;
-  i += 14;
-  const af0 = getBitS(b, i, 31) * 2 ** -34;
-  i += 31;
-  const af1 = getBitS(b, i, 21) * 2 ** -46;
-  i += 21;
-  const af2 = getBitS(b, i, 6) * 2 ** -59;
-  i = 128 * 5;
-  const type5 = getBitU(b, i, 6);
-  i += 6 + 11 + 11 + 14 + 5;
-  const bgdE5a = getBitS(b, i, 10) * 2 ** -32;
-  i += 10;
-  i += 10;
-  const e5bHs = getBitU(b, i, 2);
-  i += 2;
-  const e1bHs = getBitU(b, i, 2);
-  i += 2;
-  const e5bDvs = getBitU(b, i, 1);
-  i += 1;
-  const e1bDvs = getBitU(b, i, 1);
-  i += 1;
-  const gstWeek = getBitU(b, i, 12);
-  i += 12;
-  const tow = getBitU(b, i, 20);
-  if (type1 !== 1 || type2 !== 2 || type3 !== 3 || type4 !== 4 || type5 !== 5)
-    return null;
-  if (iodNav1 !== iodNav2 || iodNav1 !== iodNav3 || iodNav1 !== iodNav4)
-    return null;
-  if (svid < 1 || svid > 36) return null;
-  const week = weekOfToe(gstWeek, tow, toes) + GST_GPS_WEEK_OFFSET;
-  const tocDate = gpsDate(week, tocs);
-  return {
-    system: "E",
-    prn: `E${two(svid)}`,
-    toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK2,
-    tocDate,
-    af0,
-    af1,
-    af2,
-    iode: iodNav1,
-    crs,
-    deltaN,
-    m0,
-    cuc,
-    e,
-    cus,
-    sqrtA,
-    toe: toes,
-    cic,
-    omega0,
-    cis,
-    i0,
-    crc,
-    omega,
-    omegaDot,
-    idot,
-    week,
-    svHealth: e5bHs << 7 | e5bDvs << 6 | e1bHs << 1 | e1bDvs,
-    tgd: bgdE5a
-  };
-}
-function decodeGalFnavPages(pages) {
-  if (pages.length < 124) return null;
-  const b = pages;
-  let i = 0;
-  const type1 = getBitU(b, i, 6);
-  i += 6;
-  const svid = getBitU(b, i, 6);
-  i += 6;
-  const iodNav1 = getBitU(b, i, 10);
-  i += 10;
-  const tocs = getBitU(b, i, 14) * 60;
-  i += 14;
-  const af0 = getBitS(b, i, 31) * 2 ** -34;
-  i += 31;
-  const af1 = getBitS(b, i, 21) * 2 ** -46;
-  i += 21;
-  const af2 = getBitS(b, i, 6) * 2 ** -59;
-  i += 6;
-  i += 8 + 11 + 11 + 14 + 5;
-  const bgdE5a = getBitS(b, i, 10) * 2 ** -32;
-  i += 10;
-  const e5aHs = getBitU(b, i, 2);
-  i += 2;
-  const gstWeek = getBitU(b, i, 12);
-  i += 12;
-  const tow = getBitU(b, i, 20);
-  i += 20;
-  const e5aDvs = getBitU(b, i, 1);
-  i = 31 * 8;
-  const type2 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav2 = getBitU(b, i, 10);
-  i += 10;
-  const m0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const omegaDot = getBitS(b, i, 24) * 2 ** -43 * GPS_PI;
-  i += 24;
-  const e = getBitU(b, i, 32) * 2 ** -33;
-  i += 32;
-  const sqrtA = getBitU(b, i, 32) * 2 ** -19;
-  i += 32;
-  const omega0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const idot = getBitS(b, i, 14) * 2 ** -43 * GPS_PI;
-  i = 62 * 8;
-  const type3 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav3 = getBitU(b, i, 10);
-  i += 10;
-  const i0 = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const omega = getBitS(b, i, 32) * 2 ** -31 * GPS_PI;
-  i += 32;
-  const deltaN = getBitS(b, i, 16) * 2 ** -43 * GPS_PI;
-  i += 16;
-  const cuc = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const cus = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const crc = getBitS(b, i, 16) * 2 ** -5;
-  i += 16;
-  const crs = getBitS(b, i, 16) * 2 ** -5;
-  i += 16;
-  const toes = getBitU(b, i, 14) * 60;
-  i = 93 * 8;
-  const type4 = getBitU(b, i, 6);
-  i += 6;
-  const iodNav4 = getBitU(b, i, 10);
-  i += 10;
-  const cic = getBitS(b, i, 16) * 2 ** -29;
-  i += 16;
-  const cis = getBitS(b, i, 16) * 2 ** -29;
-  if (type1 !== 1 || type2 !== 2 || type3 !== 3 || type4 !== 4) return null;
-  if (iodNav1 !== iodNav2 || iodNav1 !== iodNav3 || iodNav1 !== iodNav4)
-    return null;
-  if (svid < 1 || svid > 36) return null;
-  const week = weekOfToe(gstWeek, tow, toes) + GST_GPS_WEEK_OFFSET;
-  const tocDate = gpsDate(week, tocs);
-  return {
-    system: "E",
-    prn: `E${two(svid)}`,
-    toc: tocDate.getTime() / 1e3 % SEC_PER_WEEK2,
-    tocDate,
-    af0,
-    af1,
-    af2,
-    iode: iodNav1,
-    crs,
-    deltaN,
-    m0,
-    cuc,
-    e,
-    cus,
-    sqrtA,
-    toe: toes,
-    cic,
-    omega0,
-    cis,
-    i0,
-    crc,
-    omega,
-    omegaDot,
-    idot,
-    week,
-    svHealth: e5aHs << 4 | e5aDvs << 3,
-    tgd: bgdE5a
-  };
-}
-var dedupKey = (eph) => `${eph.iode}:${eph.toe}:${eph.tocDate.getTime()}`;
-var GalInavAssembler = class {
-  sats = /* @__PURE__ */ new Map();
-  /**
-   * Push one I/NAV page pair for satellite `prn` (1-36): the 234-bit
-   * even+odd concatenation, even part first with its 6 tail bits
-   * removed, bit 0 = the even part's Even/Odd bit (≥ 17 bytes).
-   * Returns the newly completed ephemeris, or null. Pages that are
-   * not an even/odd nominal pair, alert pages and word types outside
-   * 1-5 are ignored.
-   */
-  push(prn, page) {
-    if (prn < 1 || prn > 36 || page.length < 17) return null;
-    if (getBitU(page, 0, 1) !== 0 || getBitU(page, 114, 1) !== 1) return null;
-    if (getBitU(page, 1, 1) === 1 || getBitU(page, 115, 1) === 1) return null;
-    const type = getBitU(page, 2, 6);
-    if (type < 1 || type > 5) return null;
-    let sat = this.sats.get(prn);
-    if (!sat) {
-      sat = { words: new Uint8Array(96) };
-      this.sats.set(prn, sat);
-    }
-    for (let k = 0; k < 14; k++)
-      sat.words[type * 16 + k] = getBitU(page, 2 + 8 * k, 8);
-    sat.words[type * 16 + 14] = getBitU(page, 116, 8);
-    sat.words[type * 16 + 15] = getBitU(page, 124, 8);
-    if (type !== 5) return null;
-    const eph = decodeGalInavWords(sat.words);
-    if (!eph || eph.prn !== `E${two(prn)}`) return null;
-    const key = dedupKey(eph);
-    if (key === sat.lastKey) return null;
-    sat.lastKey = key;
-    return eph;
-  }
-};
-var GalFnavAssembler = class {
-  sats = /* @__PURE__ */ new Map();
-  /**
-   * Push one F/NAV page for satellite `prn` (1-36): 244 bits with the
-   * sync field stripped, bit 0 = first page-type bit (≥ 31 bytes).
-   * Returns the newly completed ephemeris, or null. Dummy pages
-   * (type 63) and page types outside 1-4 are ignored.
-   */
-  push(prn, page) {
-    if (prn < 1 || prn > 36 || page.length < 31) return null;
-    const type = getBitU(page, 0, 6);
-    if (type < 1 || type > 4) return null;
-    let sat = this.sats.get(prn);
-    if (!sat) {
-      sat = { words: new Uint8Array(124) };
-      this.sats.set(prn, sat);
-    }
-    sat.words.set(page.subarray(0, 31), (type - 1) * 31);
-    if (type !== 4) return null;
-    const eph = decodeGalFnavPages(sat.words);
-    if (!eph || eph.prn !== `E${two(prn)}`) return null;
-    const key = dedupKey(eph);
-    if (key === sat.lastKey) return null;
-    sat.lastKey = key;
-    return eph;
-  }
-};
-
-// src/navbits/glo.ts
-var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
-var SEC_PER_WEEK3 = 7 * 86400;
 var SEC_PER_DAY = 86400;
 var GLO_STRING_BYTES = 10;
 function getBitG(b, pos, len) {
@@ -691,9 +380,9 @@ function decodeGloStrings(strings, refDate, opts = {}) {
   const slot = getBitU(b, i, 5);
   if (frn1 !== 1 || frn2 !== 2 || frn3 !== 3 || frn4 !== 4) return null;
   if (slot < 1 || slot > 27) return null;
-  const utcSec = (getUtcDate(refDate).getTime() - GPS_EPOCH_MS2) / 1e3;
-  const week = Math.floor(utcSec / SEC_PER_WEEK3);
-  let tow = utcSec - week * SEC_PER_WEEK3;
+  const utcSec = (getUtcDate(refDate).getTime() - GPS_EPOCH_MS) / 1e3;
+  const week = Math.floor(utcSec / SEC_PER_WEEK2);
+  let tow = utcSec - week * SEC_PER_WEEK2;
   const tod = tow % SEC_PER_DAY;
   tow -= tod;
   let tof = tkH * 3600 + tkM * 60 + tkS - 10800;
@@ -706,12 +395,12 @@ function decodeGloStrings(strings, refDate, opts = {}) {
     system: "R",
     prn: `R${String(slot).padStart(2, "0")}`,
     // RINEX GLONASS epochs are UTC: build the UTC toe Date directly.
-    tocDate: new Date(GPS_EPOCH_MS2 + (week * SEC_PER_WEEK3 + tow + toe) * 1e3),
+    tocDate: new Date(GPS_EPOCH_MS + (week * SEC_PER_WEEK2 + tow + toe) * 1e3),
     tauN: -tauN,
     // RINEX stores −τn; the SIS carries τn (ICD sign)
     gammaN,
     // v3 message frame time: seconds of the UTC week (RTKLIB tof).
-    messageFrameTime: ((tow + tof) % SEC_PER_WEEK3 + SEC_PER_WEEK3) % SEC_PER_WEEK3,
+    messageFrameTime: ((tow + tof) % SEC_PER_WEEK2 + SEC_PER_WEEK2) % SEC_PER_WEEK2,
     x,
     xDot,
     xAcc,
@@ -762,10 +451,6 @@ var GloStringAssembler = class {
 export {
   bdsSubframeParityOk,
   BdsAssembler,
-  galInavPageCrcOk,
-  galFnavPageCrcOk,
-  GalInavAssembler,
-  GalFnavAssembler,
   testGloString,
   GloStringAssembler
 };