gnss-js 1.18.0 → 1.19.0

package/dist/rinex.cjs

@@ -26,6 +26,7 @@ __export(rinex_exports, {
   createGzipLineSink: () => createGzipLineSink,
   crxDecompress: () => crxDecompress,
   crxRepair: () => crxRepair,
+  fmtD: () => fmtD,
   fmtF: () => fmtF,
   hdrLine: () => hdrLine,
   padL: () => padL,
@@ -43,6 +44,7 @@ __export(rinex_exports, {
   writeRinex2ObsBlob: () => writeRinex2ObsBlob,
   writeRinex4ObsBlob: () => writeRinex4ObsBlob,
   writeRinexNav: () => writeRinexNav,
+  writeRinexNav4: () => writeRinexNav4,
   writeRinexObsBlob: () => writeRinexObsBlob
 });
 module.exports = __toCommonJS(rinex_exports);
@@ -1038,6 +1040,17 @@ var WarningAccumulator = class {
   }
 };
 
+// src/navbits/index.ts
+var GPS_PI = 3.1415926535898;
+var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
+var SEC_PER_WEEK = 7 * 86400;
+
+// src/navbits/cnav.ts
+var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
+var SEC_PER_WEEK2 = 7 * 86400;
+var CNAV_A_REF = 26559710;
+var CNAV_OMEGA_DOT_REF = -26e-10;
+
 // src/rinex/nav.ts
 var DATA_LINES = {
   G: 7,
@@ -1061,19 +1074,19 @@ var R4_MSG_LINES = {
   SBAS: 4,
   FDMA: 5,
   // RINEX 4 adds orbit-4 (status flags) vs 3 in RINEX 3
-  // EPH types we skip (different field layout)
+  // EPH type parsed into RinexCnavEphemeris for G/J (skipped otherwise)
   CNAV: 9,
+  // EPH types we skip (different field layout)
   CNV1: 10,
   CNV2: 10,
   CNV3: 9,
   L1NV: 8,
   L1OC: 9,
-  L3OC: 9,
-  // Non-EPH records
-  STO: 2,
-  EOP: 3,
-  ION: 3
+  L3OC: 9
 };
+var GPS_EPOCH_MS3 = Date.UTC(1980, 0, 6);
+var SEC_PER_WEEK3 = 7 * 86400;
+var QZSS_CNAV_A_REF = 42164200;
 var SUPPORTED_EPH_MSGS = /* @__PURE__ */ new Set([
   "LNAV",
   "INAV",
@@ -1137,6 +1150,20 @@ function parseDataLine(line, colOffset = 4) {
   }
   return values;
 }
+function parseDataLineSlots(line) {
+  const values = [];
+  for (let i = 0; i < 4; i++) {
+    const start = 4 + i * 19;
+    const s = start < line.length ? line.substring(start, start + 19).trim() : "";
+    if (s.length === 0) {
+      values.push(null);
+      continue;
+    }
+    const v = parseFloat19(s);
+    values.push(Number.isFinite(v) ? v : null);
+  }
+  return values;
+}
 function buildKeplerEphemeris(sys, prn, date, epochVals, data) {
   const d = [];
   for (const line of data) {
@@ -1199,6 +1226,68 @@ function buildStateVectorEphemeris(sys, prn, date, epochVals, data) {
     zAcc: d[10] ?? 0
   };
 }
+function buildCnavEphemeris(sys, prn, date, epochVals, d) {
+  const num = (r, c) => d[r]?.[c] ?? 0;
+  const opt = (r, c) => d[r]?.[c] ?? null;
+  const t = (date.getTime() - GPS_EPOCH_MS3) / 1e3;
+  const week = Math.floor(t / SEC_PER_WEEK3);
+  const toc = t - week * SEC_PER_WEEK3;
+  const sqrtA = num(1, 3);
+  const a = sqrtA * sqrtA;
+  const aRef = sys === "J" ? QZSS_CNAV_A_REF : CNAV_A_REF;
+  const omegaDot = num(3, 3);
+  const wnOp = d[7]?.[1];
+  return {
+    system: sys,
+    prn,
+    week,
+    health: num(5, 1),
+    uraEd: num(5, 0),
+    uraNed0: num(4, 2),
+    uraNed1: num(4, 3),
+    uraNed2: num(5, 3),
+    top: num(2, 0),
+    toe: toc,
+    toeDate: date,
+    a,
+    deltaA: a - aRef,
+    aDot: num(0, 0),
+    deltaN0: num(0, 2),
+    deltaN0Dot: num(4, 1),
+    m0: num(0, 3),
+    e: num(1, 1),
+    omega: num(3, 2),
+    omega0: num(2, 2),
+    i0: num(3, 0),
+    omegaDot,
+    deltaOmegaDot: omegaDot - CNAV_OMEGA_DOT_REF * GPS_PI,
+    i0Dot: num(4, 0),
+    cis: num(2, 3),
+    cic: num(2, 1),
+    crs: num(0, 1),
+    crc: num(3, 1),
+    cus: num(1, 2),
+    cuc: num(1, 0),
+    toc,
+    tocDate: date,
+    af0: epochVals[0] ?? 0,
+    af1: epochVals[1] ?? 0,
+    af2: epochVals[2] ?? 0,
+    clockMsgType: 0,
+    // not recorded in RINEX
+    tgd: opt(5, 2),
+    iscL1ca: opt(6, 0),
+    iscL2c: opt(6, 1),
+    iscL5i5: opt(6, 2),
+    iscL5q5: opt(6, 3),
+    ...wnOp != null && { wnOp },
+    integrityFlag: false,
+    // not recorded in RINEX
+    l2cPhasing: false,
+    // not recorded in RINEX
+    tow: num(7, 0)
+  };
+}
 function parseNavFile(text) {
   const lines = text.split("\n");
   const header = {
@@ -1208,6 +1297,7 @@ function parseNavFile(text) {
     ionoCorrections: {}
   };
   const ephemerides = [];
+  const cnav = [];
   let inHeader = true;
   let i = 0;
   while (i < lines.length) {
@@ -1268,8 +1358,15 @@ function parseNavFile(text) {
     if (isV4 && line.charAt(0) === ">") {
       const parts = line.substring(2).trim().split(/\s+/);
       const recType = parts[0] ?? "";
+      const svSys = (parts[1] ?? "").charAt(0);
       const msgType = parts[2] ?? recType;
-      if (recType !== "EPH" || !SUPPORTED_EPH_MSGS.has(msgType)) {
+      if (recType !== "EPH") {
+        const totalLines = recType === "STO" ? 2 : recType === "EOP" ? 3 : recType === "ION" ? svSys === "E" ? 2 : 3 : 2;
+        i += totalLines + 1;
+        continue;
+      }
+      const isCnav = msgType === "CNAV" && (svSys === "G" || svSys === "J");
+      if (!SUPPORTED_EPH_MSGS.has(msgType) && !isCnav) {
         const totalLines = R4_MSG_LINES[msgType] ?? 2;
         i += totalLines + 1;
         continue;
@@ -1281,6 +1378,21 @@ function parseNavFile(text) {
       const sys2 = epochLine.charAt(0);
       const parsed = parseNavEpochV3(epochLine);
       const prn2 = `${sys2}${parsed.prn.substring(1)}`;
+      if (isCnav) {
+        const slotLines = [];
+        for (let j = 0; j < numDataLines2; j++) {
+          i++;
+          if (i >= lines.length) break;
+          slotLines.push(parseDataLineSlots(lines[i]));
+        }
+        if (slotLines.length === numDataLines2 && (sys2 === "G" || sys2 === "J")) {
+          cnav.push(
+            buildCnavEphemeris(sys2, prn2, parsed.date, parsed.values, slotLines)
+          );
+        }
+        i++;
+        continue;
+      }
       const dataLines2 = [];
       for (let j = 0; j < numDataLines2; j++) {
         i++;
@@ -1354,7 +1466,7 @@ function parseNavFile(text) {
     }
     i++;
   }
-  return { header, ephemerides };
+  return cnav.length > 0 ? { header, ephemerides, cnav } : { header, ephemerides };
 }
 
 // src/rinex/nav-writer.ts
@@ -1362,9 +1474,13 @@ function fmtD(val) {
   if (val === 0) return " 0.000000000000E+00";
   const sign = val < 0 ? "-" : " ";
   const abs = Math.abs(val);
-  const exp = Math.floor(Math.log10(abs));
+  let exp = Math.floor(Math.log10(abs));
   const mantissa = abs / 10 ** exp;
-  const mStr = mantissa.toFixed(12);
+  let mStr = mantissa.toFixed(12);
+  if (mStr.length > 14) {
+    exp += 1;
+    mStr = (mantissa / 10).toFixed(12);
+  }
   const expSign = exp >= 0 ? "+" : "-";
   const expStr = String(Math.abs(exp)).padStart(2, "0");
   return `${sign}${mStr}E${expSign}${expStr}`;
@@ -1439,7 +1555,7 @@ function writeKeplerRecord(eph) {
 function writeGlonassRecord(eph) {
   const lines = [];
   lines.push(
-    `${padL(eph.prn, 3)} ${fmtEpoch(eph.tocDate)}${fmtD(-eph.tauN)}${fmtD(eph.gammaN)}${fmtD(eph.messageFrameTime)}`
+    `${padL(eph.prn, 3)} ${fmtEpoch(eph.tocDate)}${fmtD(eph.tauN)}${fmtD(eph.gammaN)}${fmtD(eph.messageFrameTime)}`
   );
   lines.push(orbitLine([eph.x, eph.xDot, eph.xAcc, eph.health]));
   lines.push(orbitLine([eph.y, eph.yDot, eph.yAcc, eph.freqNum]));
@@ -1459,6 +1575,181 @@ function writeRinexNav(nav) {
   return header + records.join("\n") + "\n";
 }
 
+// src/rinex/nav-writer-v4.ts
+var BLANK19 = " ".repeat(19);
+function fmtEpoch4(d) {
+  const p2 = (v) => String(v).padStart(2, "0");
+  return `${d.getUTCFullYear()} ${p2(d.getUTCMonth() + 1)} ${p2(d.getUTCDate())} ${p2(d.getUTCHours())} ${p2(d.getUTCMinutes())} ${p2(d.getUTCSeconds())}`;
+}
+function orbitLine4(values) {
+  const line = "    " + values.map((v) => v == null ? BLANK19 : fmtD(v)).join("");
+  return line.replace(/ +$/, "");
+}
+function epochLine4(prn, date, vals) {
+  return `${padL(prn, 3)} ${fmtEpoch4(date)}${vals.map(fmtD).join("")}`;
+}
+var BDS_GEO = /* @__PURE__ */ new Set([
+  "C01",
+  "C02",
+  "C03",
+  "C04",
+  "C05",
+  "C59",
+  "C60",
+  "C61",
+  "C62",
+  "C63"
+]);
+function msgLabel(eph) {
+  switch (eph.system) {
+    case "E":
+      return "INAV";
+    case "C":
+      return BDS_GEO.has(eph.prn) ? "D2" : "D1";
+    case "R":
+      return "FDMA";
+    case "S":
+      return "SBAS";
+    default:
+      return "LNAV";
+  }
+}
+function isKepler2(eph) {
+  return "af0" in eph;
+}
+function keplerRecord(eph) {
+  return [
+    epochLine4(eph.prn, eph.tocDate, [eph.af0, eph.af1, eph.af2]),
+    orbitLine4([eph.iode, eph.crs, eph.deltaN, eph.m0]),
+    orbitLine4([eph.cuc, eph.e, eph.cus, eph.sqrtA]),
+    orbitLine4([eph.toe, eph.cic, eph.omega0, eph.cis]),
+    orbitLine4([eph.i0, eph.crc, eph.omega, eph.omegaDot]),
+    orbitLine4([eph.idot, 0, eph.week, 0]),
+    orbitLine4([0, eph.svHealth, eph.tgd, 0]),
+    orbitLine4([eph.toe, 0, 0, 0])
+  ].join("\n");
+}
+function stateVectorRecord(eph) {
+  const lines = [
+    epochLine4(eph.prn, eph.tocDate, [
+      eph.tauN,
+      eph.gammaN,
+      eph.messageFrameTime
+    ]),
+    orbitLine4([eph.x, eph.xDot, eph.xAcc, eph.health]),
+    orbitLine4([eph.y, eph.yDot, eph.yAcc, eph.freqNum]),
+    orbitLine4([eph.z, eph.zDot, eph.zAcc, 0])
+  ];
+  if (eph.system === "R") {
+    lines.push(orbitLine4([null, 999999999999e-3, 15, null]));
+  }
+  return lines.join("\n");
+}
+function resolveWnOp(week, wnOp) {
+  if (wnOp == null) return week;
+  return week - ((week - wnOp) % 256 + 256) % 256;
+}
+function cnavRecord(eph) {
+  return [
+    epochLine4(eph.prn, eph.tocDate, [eph.af0, eph.af1, eph.af2]),
+    orbitLine4([eph.aDot, eph.crs, eph.deltaN0, eph.m0]),
+    orbitLine4([eph.cuc, eph.e, eph.cus, Math.sqrt(eph.a)]),
+    orbitLine4([eph.top, eph.cic, eph.omega0, eph.cis]),
+    orbitLine4([eph.i0, eph.crc, eph.omega, eph.omegaDot]),
+    orbitLine4([eph.i0Dot, eph.deltaN0Dot, eph.uraNed0, eph.uraNed1]),
+    orbitLine4([eph.uraEd, eph.health, eph.tgd, eph.uraNed2]),
+    orbitLine4([eph.iscL1ca, eph.iscL2c, eph.iscL5i5, eph.iscL5q5]),
+    orbitLine4([eph.tow, resolveWnOp(eph.week, eph.wnOp)])
+  ].join("\n");
+}
+var KLOBUCHAR_SOURCES = [
+  ["GPSA", "GPSB", "G", "LNAV"],
+  ["QZSA", "QZSB", "J", "LNAV"],
+  ["BDSA", "BDSB", "C", "D1D2"],
+  ["IRNA", "IRNB", "I", "LNAV"]
+];
+function ionRecords(header, epoch) {
+  const records = [];
+  const iono = header.ionoCorrections;
+  for (const [alphaKey, betaKey, sys, msg] of KLOBUCHAR_SOURCES) {
+    const alpha = iono[alphaKey];
+    if (!alpha) continue;
+    const beta = iono[betaKey] ?? [0, 0, 0, 0];
+    records.push(
+      [
+        `> ION ${sys}   ${msg}`,
+        `    ${fmtEpoch4(epoch)}${[alpha[0] ?? 0, alpha[1] ?? 0, alpha[2] ?? 0].map(fmtD).join("")}`,
+        orbitLine4([alpha[3] ?? 0, beta[0] ?? 0, beta[1] ?? 0, beta[2] ?? 0]),
+        orbitLine4([beta[3] ?? 0, null])
+      ].join("\n")
+    );
+  }
+  const gal = iono["GAL"];
+  if (gal) {
+    records.push(
+      [
+        "> ION E   IFNV",
+        `    ${fmtEpoch4(epoch)}${[gal[0] ?? 0, gal[1] ?? 0, gal[2] ?? 0].map(fmtD).join("")}`,
+        orbitLine4([gal[3] ?? 0])
+      ].join("\n")
+    );
+  }
+  return records;
+}
+function writeHeader4(header) {
+  const lines = [];
+  lines.push(
+    hdrLine(
+      padL("     4.01", 9) + "           " + padL("NAVIGATION DATA", 20) + padL("M", 20),
+      "RINEX VERSION / TYPE"
+    )
+  );
+  const now = /* @__PURE__ */ new Date();
+  const dateStr = now.getUTCFullYear().toString() + String(now.getUTCMonth() + 1).padStart(2, "0") + String(now.getUTCDate()).padStart(2, "0") + " " + String(now.getUTCHours()).padStart(2, "0") + String(now.getUTCMinutes()).padStart(2, "0") + String(now.getUTCSeconds()).padStart(2, "0") + " UTC";
+  lines.push(
+    hdrLine(
+      padL("GNSSCalc", 20) + padL("", 20) + padL(dateStr, 20),
+      "PGM / RUN BY / DATE"
+    )
+  );
+  if (header.leapSeconds != null) {
+    lines.push(hdrLine(padR(String(header.leapSeconds), 6), "LEAP SECONDS"));
+  }
+  lines.push(hdrLine("", "END OF HEADER"));
+  return lines.join("\n") + "\n";
+}
+function writeRinexNav4(nav) {
+  const records = [];
+  for (const eph of nav.ephemerides) {
+    const label = msgLabel(eph);
+    records.push({
+      prn: eph.prn,
+      time: eph.tocDate.getTime(),
+      label,
+      body: `> EPH ${padL(eph.prn, 3)} ${label}
+` + (isKepler2(eph) ? keplerRecord(eph) : stateVectorRecord(eph))
+    });
+  }
+  for (const eph of nav.cnav ?? []) {
+    records.push({
+      prn: eph.prn,
+      time: eph.tocDate.getTime(),
+      label: "CNAV",
+      body: `> EPH ${padL(eph.prn, 3)} CNAV
+` + cnavRecord(eph)
+    });
+  }
+  records.sort(
+    (a, b) => a.prn.localeCompare(b.prn) || a.time - b.time || a.label.localeCompare(b.label)
+  );
+  const bodies = records.map((r) => r.body);
+  if (records.length > 0) {
+    const earliest = new Date(Math.min(...records.map((r) => r.time)));
+    bodies.push(...ionRecords(nav.header, earliest));
+  }
+  return writeHeader4(nav.header) + (bodies.length > 0 ? bodies.join("\n") + "\n" : "");
+}
+
 // src/rinex/ionex.ts
 function gridRange(l1, l2, dl) {
   const out = [];
@@ -2041,6 +2332,7 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
   createGzipLineSink,
   crxDecompress,
   crxRepair,
+  fmtD,
   fmtF,
   hdrLine,
   padL,
@@ -2058,5 +2350,6 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,
   writeRinexNav,
+  writeRinexNav4,
   writeRinexObsBlob
 });

package/dist/rinex.d.cts

@@ -1,7 +1,7 @@
 import { R as RinexHeader } from './parser-JPjjFgeP.cjs';
 export { E as EpochSummary, P as ParseOptions, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.cjs';
-import { a as NavResult } from './nav-BAI1a9vK.cjs';
-export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, p as parseNavFile } from './nav-BAI1a9vK.cjs';
+import { a as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris } from './nav-DImXUdbp.cjs';
+export { G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-DImXUdbp.cjs';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -106,6 +106,16 @@ declare class WarningAccumulator {
  * Produces a standards-compliant RINEX 3.04 nav file from parsed NavResult data.
  */
 
+/**
+ * Format a number in RINEX "D19.12" scientific notation.
+ *
+ * 19 characters wide, 12 decimal digits, explicit exponent sign.
+ * Examples:
+ *   0.0       → " 0.000000000000E+00"
+ *  -1.5e-8    → "-1.500000000000E-08"
+ *   5153.71875→ " 5.153718750000E+03"
+ */
+declare function fmtD(val: number): string;
 /**
  * Write a RINEX 3.04 mixed navigation file from parsed navigation data.
  *
@@ -113,6 +123,38 @@ declare class WarningAccumulator {
  */
 declare function writeRinexNav(nav: NavResult): string;
 
+/**
+ * RINEX 4.01 mixed navigation file writer.
+ *
+ * Produces `> EPH` data records for the classic ephemeris types
+ * (GPS/QZSS/NavIC LNAV, Galileo INAV, BeiDou D1/D2, GLONASS FDMA,
+ * SBAS), `> EPH ... CNAV` records for GPS/QZSS CNAV ephemerides, and
+ * `> ION` records for the header's Klobuchar/NeQuick-G coefficients,
+ * per the RINEX 4.01 record tables (A9-A18, A21, A26, A32, A33).
+ */
+
+/** Input for {@link writeRinexNav4}: a NavResult, optionally with CNAV. */
+interface Nav4Input {
+    header: NavHeader;
+    ephemerides: Ephemeris[];
+    cnav?: RinexCnavEphemeris[];
+}
+/**
+ * Write a RINEX 4.01 mixed navigation file.
+ *
+ * Emits `> EPH` records for every classic ephemeris (message labels
+ * per constellation: LNAV, INAV, D1/D2, FDMA, SBAS) and every CNAV
+ * ephemeris (`CNAV`), sorted by PRN, epoch, then label; followed by
+ * `> ION` records when the header carries Klobuchar (GPSA/GPSB,
+ * QZSA/QZSB, BDSA/BDSB, IRNA/IRNB) or NeQuick-G (GAL) coefficients.
+ * ION records need a transmit-time epoch, which RINEX 3-style headers
+ * do not store — the earliest record epoch is used, so a file with no
+ * ephemerides gets no ION records. There is no source for system
+ * time offsets in NavHeader, so no `> STO` records are produced; leap
+ * seconds go in the LEAP SECONDS header line.
+ */
+declare function writeRinexNav4(nav: Nav4Input): string;
+
 /**
  * IONEX 1.0 parser — global ionosphere maps (GIMs) as published by
  * the IGS analysis centres (e.g. ESA0OPSRAP_*_GIM.INX).
@@ -263,4 +305,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, type IonexGrid, NavResult, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, createGzipLineSink, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseIonex, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexObsBlob };
+export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type IonexGrid, 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, parseIonex, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.d.ts

@@ -1,7 +1,7 @@
 import { R as RinexHeader } from './parser-JPjjFgeP.js';
 export { E as EpochSummary, P as ParseOptions, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.js';
-import { a as NavResult } from './nav-BAI1a9vK.js';
-export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, p as parseNavFile } from './nav-BAI1a9vK.js';
+import { a as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris } from './nav-DImXUdbp.js';
+export { G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-DImXUdbp.js';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -106,6 +106,16 @@ declare class WarningAccumulator {
  * Produces a standards-compliant RINEX 3.04 nav file from parsed NavResult data.
  */
 
+/**
+ * Format a number in RINEX "D19.12" scientific notation.
+ *
+ * 19 characters wide, 12 decimal digits, explicit exponent sign.
+ * Examples:
+ *   0.0       → " 0.000000000000E+00"
+ *  -1.5e-8    → "-1.500000000000E-08"
+ *   5153.71875→ " 5.153718750000E+03"
+ */
+declare function fmtD(val: number): string;
 /**
  * Write a RINEX 3.04 mixed navigation file from parsed navigation data.
  *
@@ -113,6 +123,38 @@ declare class WarningAccumulator {
  */
 declare function writeRinexNav(nav: NavResult): string;
 
+/**
+ * RINEX 4.01 mixed navigation file writer.
+ *
+ * Produces `> EPH` data records for the classic ephemeris types
+ * (GPS/QZSS/NavIC LNAV, Galileo INAV, BeiDou D1/D2, GLONASS FDMA,
+ * SBAS), `> EPH ... CNAV` records for GPS/QZSS CNAV ephemerides, and
+ * `> ION` records for the header's Klobuchar/NeQuick-G coefficients,
+ * per the RINEX 4.01 record tables (A9-A18, A21, A26, A32, A33).
+ */
+
+/** Input for {@link writeRinexNav4}: a NavResult, optionally with CNAV. */
+interface Nav4Input {
+    header: NavHeader;
+    ephemerides: Ephemeris[];
+    cnav?: RinexCnavEphemeris[];
+}
+/**
+ * Write a RINEX 4.01 mixed navigation file.
+ *
+ * Emits `> EPH` records for every classic ephemeris (message labels
+ * per constellation: LNAV, INAV, D1/D2, FDMA, SBAS) and every CNAV
+ * ephemeris (`CNAV`), sorted by PRN, epoch, then label; followed by
+ * `> ION` records when the header carries Klobuchar (GPSA/GPSB,
+ * QZSA/QZSB, BDSA/BDSB, IRNA/IRNB) or NeQuick-G (GAL) coefficients.
+ * ION records need a transmit-time epoch, which RINEX 3-style headers
+ * do not store — the earliest record epoch is used, so a file with no
+ * ephemerides gets no ION records. There is no source for system
+ * time offsets in NavHeader, so no `> STO` records are produced; leap
+ * seconds go in the LEAP SECONDS header line.
+ */
+declare function writeRinexNav4(nav: Nav4Input): string;
+
 /**
  * IONEX 1.0 parser — global ionosphere maps (GIMs) as published by
  * the IGS analysis centres (e.g. ESA0OPSRAP_*_GIM.INX).
@@ -263,4 +305,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, type IonexGrid, NavResult, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, createGzipLineSink, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseIonex, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexObsBlob };
+export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type IonexGrid, 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, parseIonex, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.js

@@ -2,6 +2,7 @@ import {
   EMPTY_WARNINGS,
   WarningAccumulator,
   createGzipLineSink,
+  fmtD,
   fmtF,
   hdrLine,
   padL,
@@ -15,8 +16,11 @@ import {
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,
   writeRinexNav,
+  writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-HXF7XVW7.js";
+} from "./chunk-K73SERD5.js";
+import "./chunk-6NCMACDO.js";
+import "./chunk-IQ5OIMQD.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -34,6 +38,7 @@ export {
   createGzipLineSink,
   crxDecompress,
   crxRepair,
+  fmtD,
   fmtF,
   hdrLine,
   padL,
@@ -51,5 +56,6 @@ export {
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,
   writeRinexNav,
+  writeRinexNav4,
   writeRinexObsBlob
 };

package/dist/sbf.d.cts

@@ -1,5 +1,4 @@
-import { E as Ephemeris } from './nav-BAI1a9vK.cjs';
-import { C as CnavEphemeris } from './cnav-CRns0ECn.cjs';
+import { E as Ephemeris, C as CnavEphemeris } from './nav-DImXUdbp.cjs';
 
 /**
  * Septentrio SBF framing shared by the measurement and navigation

package/dist/sbf.d.ts

@@ -1,5 +1,4 @@
-import { E as Ephemeris } from './nav-BAI1a9vK.js';
-import { C as CnavEphemeris } from './cnav-CRns0ECn.js';
+import { E as Ephemeris, C as CnavEphemeris } from './nav-DImXUdbp.js';
 
 /**
  * Septentrio SBF framing shared by the measurement and navigation

package/dist/sbf.js

@@ -6,8 +6,8 @@ import {
   parseSbfMeas,
   parseSbfNav,
   scanSbfFrames
-} from "./chunk-CG7SF45F.js";
-import "./chunk-Y5IWUPJQ.js";
+} from "./chunk-5S5HMAKW.js";
+import "./chunk-6NCMACDO.js";
 import "./chunk-IQ5OIMQD.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";

package/dist/ubx.d.cts

@@ -1,5 +1,4 @@
-import { E as Ephemeris } from './nav-BAI1a9vK.cjs';
-import { C as CnavEphemeris } from './cnav-CRns0ECn.cjs';
+import { E as Ephemeris, C as CnavEphemeris } from './nav-DImXUdbp.cjs';
 
 /**
  * u-blox UBX framing shared by the measurement and navigation decoders:

package/dist/ubx.d.ts

@@ -1,5 +1,4 @@
-import { E as Ephemeris } from './nav-BAI1a9vK.js';
-import { C as CnavEphemeris } from './cnav-CRns0ECn.js';
+import { E as Ephemeris, C as CnavEphemeris } from './nav-DImXUdbp.js';
 
 /**
  * u-blox UBX framing shared by the measurement and navigation decoders:

package/dist/ubx.js

@@ -4,8 +4,8 @@ import {
   parseUbxNav,
   parseUbxRawx,
   ubxFrames
-} from "./chunk-FEGLXNWA.js";
-import "./chunk-Y5IWUPJQ.js";
+} from "./chunk-JW2ICDR7.js";
+import "./chunk-6NCMACDO.js";
 import "./chunk-IQ5OIMQD.js";
 export {
   parseUbxCnav,

package/package.json

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