gnss-js 1.25.0 → 1.27.0

package/dist/rinex.d.cts

@@ -49,6 +49,32 @@ declare function crxDecompress(prev: DiffState | null, field: CrxField): {
     result: number;
 };
 
+/**
+ * Compact RINEX (Hatanaka / CRX) *compression* — the inverse of crx.ts.
+ *
+ * This is a faithful TypeScript port of RNXCMP's `rnx2crx.c` (v4.2.0,
+ * Y. Hatanaka / GSI Japan), operating on RINEX text in and CRINEX text
+ * out. It reproduces the reference compressor byte-for-byte in the data
+ * section (the CRINEX header's PROG/DATE line carries a timestamp, so it
+ * is the only line that legitimately differs from a given run of the C
+ * tool). Validated against the `rnx2crx`/`crx2rnx` oracles on real
+ * multi-GNSS RINEX 3.03/3.04 observation files.
+ *
+ * Supports RINEX 2.x (CRINEX 1.0) and RINEX 3.x / 4.0x (CRINEX 3.0).
+ * The RINEX 4.02 pico-second epoch record (CRINEX 3.1) is handled too.
+ */
+interface WriteCrxOptions {
+    /** CRINEX PROG / DATE program label (≤20 chars). */
+    prog?: string;
+    /** CRINEX PROG / DATE timestamp (as `dd-MMM-yy HH:mm`). Defaults to now (UTC). */
+    date?: string;
+}
+/**
+ * Compress a full RINEX observation file (text) to Compact RINEX (text).
+ * The inverse of the crx.ts decoder. Throws on structurally invalid input.
+ */
+declare function writeCrx(rinexText: string, opts?: WriteCrxOptions): string;
+
 /**
  * Shared RINEX formatting helpers.
  *
@@ -284,4 +310,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.d.ts

@@ -49,6 +49,32 @@ declare function crxDecompress(prev: DiffState | null, field: CrxField): {
     result: number;
 };
 
+/**
+ * Compact RINEX (Hatanaka / CRX) *compression* — the inverse of crx.ts.
+ *
+ * This is a faithful TypeScript port of RNXCMP's `rnx2crx.c` (v4.2.0,
+ * Y. Hatanaka / GSI Japan), operating on RINEX text in and CRINEX text
+ * out. It reproduces the reference compressor byte-for-byte in the data
+ * section (the CRINEX header's PROG/DATE line carries a timestamp, so it
+ * is the only line that legitimately differs from a given run of the C
+ * tool). Validated against the `rnx2crx`/`crx2rnx` oracles on real
+ * multi-GNSS RINEX 3.03/3.04 observation files.
+ *
+ * Supports RINEX 2.x (CRINEX 1.0) and RINEX 3.x / 4.0x (CRINEX 3.0).
+ * The RINEX 4.02 pico-second epoch record (CRINEX 3.1) is handled too.
+ */
+interface WriteCrxOptions {
+    /** CRINEX PROG / DATE program label (≤20 chars). */
+    prog?: string;
+    /** CRINEX PROG / DATE timestamp (as `dd-MMM-yy HH:mm`). Defaults to now (UTC). */
+    date?: string;
+}
+/**
+ * Compress a full RINEX observation file (text) to Compact RINEX (text).
+ * The inverse of the crx.ts decoder. Throws on structurally invalid input.
+ */
+declare function writeCrx(rinexText: string, opts?: WriteCrxOptions): string;
+
 /**
  * Shared RINEX formatting helpers.
  *
@@ -284,4 +310,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.js

@@ -12,14 +12,16 @@ import {
   parseSp3,
   sp3Position,
   stationHeaderLines,
+  writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-NOQAAUGY.js";
-import "./chunk-WZDFZ76K.js";
+} from "./chunk-QVNNHACQ.js";
+import "./chunk-726LQBGM.js";
+import "./chunk-IQ5OIMQD.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -51,6 +53,7 @@ export {
   stationHeaderLines,
   systemCmp,
   systemName,
+  writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,

package/dist/sbf.js

@@ -10,9 +10,10 @@ import {
   parseSbfMeas,
   parseSbfNav,
   scanSbfFrames
-} from "./chunk-PZNVFCKJ.js";
-import "./chunk-VTX3VCL4.js";
-import "./chunk-WZDFZ76K.js";
+} from "./chunk-AFEDNWXI.js";
+import "./chunk-5AUK42AL.js";
+import "./chunk-726LQBGM.js";
+import "./chunk-IQ5OIMQD.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";
 export {

package/dist/trimble.cjs

@@ -0,0 +1,452 @@
+"use strict";
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// src/trimble/index.ts
+var trimble_exports = {};
+__export(trimble_exports, {
+  parseTrimble: () => parseTrimble,
+  parseTrimbleNav: () => parseTrimbleNav,
+  trimbleFrames: () => trimbleFrames
+});
+module.exports = __toCommonJS(trimble_exports);
+
+// src/trimble/frame.ts
+var STX = 2;
+var ETX = 3;
+var RAWDATA = 87;
+var RETSVDATA = 85;
+function* trimbleFrames(data, stats) {
+  let i = 0;
+  while (i + 6 <= data.length) {
+    if (data[i] !== STX) {
+      i++;
+      continue;
+    }
+    const len = data[i + 3];
+    const total = 4 + len + 2;
+    if (i + total > data.length) {
+      i++;
+      continue;
+    }
+    if (data[i + total - 1] !== ETX) {
+      i++;
+      continue;
+    }
+    let cs = 0;
+    for (let k = 1; k <= 3 + len; k++) cs = cs + data[i + k] & 255;
+    if (cs !== data[i + 4 + len]) {
+      stats.badChecksum++;
+      i++;
+      continue;
+    }
+    yield {
+      start: i,
+      status: data[i + 1],
+      type: data[i + 2],
+      len,
+      payload: i + 4
+    };
+    i += total;
+  }
+}
+
+// src/trimble/nav.ts
+var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
+var SEC_PER_WEEK = 7 * 86400;
+var MS_PER_WEEK = SEC_PER_WEEK * 1e3;
+var GPS_PI = 3.1415926535898;
+var gpsMs = (week, sec) => GPS_EPOCH_MS + week * MS_PER_WEEK + sec * 1e3;
+var sowOf = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK;
+var SUB_GPS_EPHEMERIS = 1;
+var SUB_ION_UTC = 3;
+function decodeGpsEphemeris(data, view, base) {
+  const prn = data[base + 5];
+  if (prn < 1 || prn > 32) return null;
+  const week = view.getUint16(base + 6, false);
+  const iode = data[base + 11];
+  const toc = view.getInt32(base + 16, false);
+  const toe = view.getUint32(base + 20, false);
+  const tgd = view.getFloat64(base + 24, false);
+  const af2 = view.getFloat64(base + 32, false);
+  const af1 = view.getFloat64(base + 40, false);
+  const af0 = view.getFloat64(base + 48, false);
+  const crs = view.getFloat64(base + 56, false);
+  const deltaN = view.getFloat64(base + 64, false) * GPS_PI;
+  const m0 = view.getFloat64(base + 72, false) * GPS_PI;
+  const cuc = view.getFloat64(base + 80, false) * GPS_PI;
+  const e = view.getFloat64(base + 88, false);
+  const cus = view.getFloat64(base + 96, false) * GPS_PI;
+  const sqrtA = view.getFloat64(base + 104, false);
+  const cic = view.getFloat64(base + 112, false) * GPS_PI;
+  const omega0 = view.getFloat64(base + 120, false) * GPS_PI;
+  const cis = view.getFloat64(base + 128, false) * GPS_PI;
+  const i0 = view.getFloat64(base + 136, false) * GPS_PI;
+  const crc = view.getFloat64(base + 144, false);
+  const omega = view.getFloat64(base + 152, false) * GPS_PI;
+  const omegaDot = view.getFloat64(base + 160, false) * GPS_PI;
+  const idot = view.getFloat64(base + 168, false) * GPS_PI;
+  const flags = view.getUint32(base + 176, false);
+  return {
+    system: "G",
+    prn: `G${String(prn).padStart(2, "0")}`,
+    toc: sowOf(gpsMs(week, toc)),
+    tocDate: new Date(gpsMs(week, toc)),
+    af0,
+    af1,
+    af2,
+    iode,
+    crs,
+    deltaN,
+    m0,
+    cuc,
+    e,
+    cus,
+    sqrtA,
+    toe,
+    cic,
+    omega0,
+    cis,
+    i0,
+    crc,
+    omega,
+    omegaDot,
+    idot,
+    week,
+    svHealth: flags >>> 4 & 127,
+    tgd
+  };
+}
+function parseTrimbleNav(data) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const stats = { badChecksum: 0 };
+  const ephemerides = [];
+  const ionoCorrections = {};
+  const retsvCounts = {};
+  let leapSeconds = null;
+  const lastIode = /* @__PURE__ */ new Map();
+  for (const f of trimbleFrames(data, stats)) {
+    if (f.type !== RETSVDATA) continue;
+    const sub = data[f.payload];
+    retsvCounts[sub] = (retsvCounts[sub] ?? 0) + 1;
+    if (sub === SUB_GPS_EPHEMERIS && f.len >= 178) {
+      const eph = decodeGpsEphemeris(data, view, f.start);
+      if (!eph) continue;
+      if (lastIode.get(eph.prn) === eph.iode) continue;
+      lastIode.set(eph.prn, eph.iode);
+      ephemerides.push(eph);
+    } else if (sub === SUB_ION_UTC && f.len >= 125) {
+      const p = f.start;
+      const alpha = [];
+      const beta = [];
+      for (let i = 0; i < 4; i++) {
+        alpha.push(view.getFloat64(p + 6 + i * 8, false));
+        beta.push(view.getFloat64(p + 38 + i * 8, false));
+      }
+      ionoCorrections["GPSA"] = alpha;
+      ionoCorrections["GPSB"] = beta;
+      leapSeconds = Math.trunc(view.getFloat64(p + 94, false));
+    }
+  }
+  return {
+    ephemerides,
+    ionoCorrections,
+    leapSeconds,
+    retsvCounts,
+    badChecksum: stats.badChecksum
+  };
+}
+
+// src/trimble/index.ts
+var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
+var MS_PER_WEEK2 = 7 * 864e5;
+var REC_TYPE_27 = 6;
+var SYS = {
+  0: "G",
+  // GPS
+  1: "S",
+  // SBAS
+  2: "R",
+  // GLONASS
+  3: "E",
+  // Galileo
+  4: "J",
+  // QZSS
+  5: "C",
+  // BeiDou (pre-ICD numbering)
+  7: "C",
+  // BeiDou (ICD numbering)
+  9: "I",
+  // NavIC / IRNSS
+  10: "C"
+  // BeiDou
+};
+var CODES = {
+  "G0-0": "1C",
+  "G0-20": "1L",
+  "G1-2": "2W",
+  "G1-5": "2L",
+  "G2-6": "5I",
+  "G2-7": "5Q",
+  "G2-8": "5X",
+  "R0-0": "1C",
+  "R0-1": "1P",
+  "R1-0": "2C",
+  "R1-1": "2P",
+  "E0-23": "1C",
+  "E2-11": "5Q",
+  "E3-11": "7Q",
+  "E4-14": "8Q",
+  "E5-36": "6C",
+  "C0-20": "1P",
+  "C2-8": "5P",
+  "C3-6": "7D",
+  "C6-26": "2I",
+  "C7-29": "6I",
+  "S0-0": "1C",
+  "S2-6": "5I",
+  "J0-0": "1C"
+};
+var BAND_DIGIT = {
+  0: "1",
+  // L1 / E1 / B1C (1575.42)
+  1: "2",
+  // L2 / G2 (1227.6)
+  2: "5",
+  // L5 / E5a (1176.45)
+  3: "7",
+  // E5b / B2b (1207.14)
+  4: "8",
+  // E5 AltBOC (1191.795)
+  5: "6",
+  // E6 (1278.75)
+  6: "2",
+  // B1I (1561.098)
+  7: "6",
+  // B3I (1268.52)
+  8: "1",
+  // BeiDou E1 (1589.742)
+  9: "3",
+  // GLONASS G3 CDMA (1202.025)
+  11: "9"
+  // NavIC S-band (2492.028)
+};
+function rinexCode(sys, band, track) {
+  const key = `${sys}${band}-${track}`;
+  return CODES[key] ?? `${BAND_DIGIT[band] ?? "?"}X`;
+}
+function rinexPrn(sys, svid) {
+  let n = svid;
+  if (sys === "S")
+    n = svid - 100;
+  else if (sys === "J" && svid >= 193) n = svid - 192;
+  return `${sys}${String(n).padStart(2, "0")}`;
+}
+var readS16 = (dv, o) => dv.getInt16(o, false);
+function readS24(d, o) {
+  const v = d[o] << 16 | d[o + 1] << 8 | d[o + 2];
+  return v & 8388608 ? v - 16777216 : v;
+}
+function readS48(d, dv, o) {
+  const hi = d[o] << 8 | d[o + 1];
+  const lo = dv.getUint32(o + 2, false);
+  const v = hi * 4294967296 + lo;
+  return hi & 32768 ? v - 281474976710656 : v;
+}
+function parseTrimble(data) {
+  const stats = { badChecksum: 0 };
+  const epochs = [];
+  const obsCodes = {};
+  const recordCounts = {};
+  const retsvCounts = {};
+  let pages = [];
+  let msgReply = -1;
+  let lastPage = 0;
+  let msgPages = 0;
+  let epoch = null;
+  const finalize = () => {
+    if (epoch) decodeEpoch(epoch, epochs, obsCodes);
+    epoch = null;
+  };
+  const onMessage = (content, reply) => {
+    if (epoch && reply === epoch.reply) {
+      epoch.chunks.push(content);
+      return;
+    }
+    finalize();
+    if (content.length < 12) return;
+    const dv = new DataView(
+      content.buffer,
+      content.byteOffset,
+      content.byteLength
+    );
+    const headerLen = content[0];
+    const week = dv.getUint16(1, false);
+    const towMs = dv.getUint32(3, false);
+    const nsv = content[10];
+    const eflags = content[11];
+    if (week < 1024 || week > 4096 || headerLen < 12 || nsv === 0) return;
+    let svStart = headerLen;
+    if (eflags & 32 && svStart < content.length) svStart += content[svStart];
+    epoch = {
+      reply,
+      week,
+      towMs,
+      nsv,
+      chunks: [content.subarray(svStart)]
+    };
+  };
+  const completeMessage = () => {
+    const total = pages.reduce((a, c) => a + c.length, 0);
+    const content = new Uint8Array(total);
+    let o = 0;
+    for (const c of pages) {
+      content.set(c, o);
+      o += c.length;
+    }
+    const reply = msgReply;
+    pages = [];
+    onMessage(content, reply);
+  };
+  for (const f of trimbleFrames(data, stats)) {
+    if (f.type !== RAWDATA) {
+      if (f.type === 85) {
+        const sub = data[f.payload];
+        retsvCounts[sub] = (retsvCounts[sub] ?? 0) + 1;
+      }
+      continue;
+    }
+    const recType = data[f.payload];
+    recordCounts[recType] = (recordCounts[recType] ?? 0) + 1;
+    if (recType !== REC_TYPE_27) continue;
+    const pageByte = data[f.payload + 1];
+    const page = pageByte >> 4;
+    const pagesTotal = pageByte & 15;
+    const reply = data[f.payload + 2];
+    const content = data.subarray(f.payload + 4, f.payload + f.len);
+    if (page === 1) {
+      pages = [content];
+      msgReply = reply;
+      msgPages = pagesTotal;
+      lastPage = 1;
+    } else if (reply === msgReply && page === lastPage + 1) {
+      pages.push(content);
+      lastPage = page;
+    } else {
+      pages = [];
+      continue;
+    }
+    if (lastPage === msgPages) completeMessage();
+  }
+  finalize();
+  return {
+    epochs,
+    obsCodes,
+    recordCounts,
+    retsvCounts,
+    badChecksum: stats.badChecksum
+  };
+}
+function decodeEpoch(epoch, epochs, obsCodes) {
+  const total = epoch.chunks.reduce((a, c) => a + c.length, 0);
+  const M = new Uint8Array(total);
+  let o = 0;
+  for (const c of epoch.chunks) {
+    M.set(c, o);
+    o += c.length;
+  }
+  const dv = new DataView(M.buffer, M.byteOffset, M.byteLength);
+  const timeMs = GPS_EPOCH_MS2 + epoch.week * MS_PER_WEEK2 + epoch.towMs;
+  const meas = [];
+  let q = 0;
+  for (let sv = 0; sv < epoch.nsv && q + 8 <= M.length; sv++) {
+    const svLen = M[q];
+    if (svLen < 8 || q + svLen > M.length) break;
+    const svid = M[q + 1];
+    const sys = SYS[M[q + 2] & 63];
+    const chan = M[q + 3] << 24 >> 24;
+    const nblk = M[q + 4];
+    const prn = sys ? rinexPrn(sys, svid) : null;
+    const gloChannel = sys === "R" ? chan : null;
+    let r = q + svLen;
+    let pr0 = 0;
+    for (let b = 0; b < nblk; b++) {
+      if (r + 5 > M.length) break;
+      const mLen = M[r];
+      if (mLen < 5 || r + mLen > M.length) break;
+      const band = M[r + 1];
+      const track = M[r + 2];
+      const cn0 = dv.getUint16(r + 3, false) / 10;
+      let p = r + 5;
+      let pr;
+      if (b === 0) {
+        const scale = sys === "S" || sys === "J" ? 64 : 128;
+        pr = dv.getUint32(p, false) / scale;
+        p += 4;
+        pr0 = pr;
+      } else {
+        pr = pr0 + readS16(dv, p) / 256;
+        p += 2;
+      }
+      const cp = -readS48(M, dv, p) / 32768;
+      p += 6;
+      const cycleSlipCount = M[p];
+      p += 1;
+      const mf = M[p];
+      p += 1;
+      let mf2 = 0;
+      if (mf & 128) {
+        mf2 = M[p];
+        p += 1;
+        let f = mf2;
+        while (f & 128) {
+          f = M[p];
+          p += 1;
+        }
+      }
+      if (mf2 & 1) p += 1;
+      let doppler = null;
+      if (mf & 4) doppler = readS24(M, p) / 256;
+      if (prn && sys) {
+        const code = rinexCode(sys, band, track);
+        meas.push({
+          prn,
+          code,
+          pr: mf & 2 ? pr : null,
+          cp: mf & 1 ? cp : null,
+          doppler: mf & 1 ? doppler : null,
+          cn0,
+          cycleSlipCount,
+          cycleSlip: (mf & 8) !== 0,
+          gloChannel
+        });
+        const list = obsCodes[sys] ??= [];
+        if (!list.includes(code)) list.push(code);
+      }
+      r += mLen;
+    }
+    q = r;
+  }
+  epochs.push({ timeMs, meas });
+}
+// Annotate the CommonJS export names for ESM import in node:
+0 && (module.exports = {
+  parseTrimble,
+  parseTrimbleNav,
+  trimbleFrames
+});