gnss-js 1.25.0 → 1.27.0

package/dist/index.d.cts

@@ -4,7 +4,7 @@ export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSat
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.cjs';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.cjs';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.cjs';
-export { CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, Sp3File, Sp3Sample, WarningAccumulator, WarningSeverity, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.cjs';
+export { CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, Sp3File, Sp3Sample, WarningAccumulator, WarningSeverity, WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.cjs';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, 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';
 export { C as CnavEphemeris, E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, R as RinexCnavEphemeris, p as parseNavFile } from './nav-DImXUdbp.cjs';
 export { I as IonexGrid, p as parseIonex } from './ionex-DaW5x9nq.cjs';
@@ -13,6 +13,8 @@ export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES,
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.cjs';
 export { HasMessage, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGloNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfParseResult, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGloNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, sbfCrc16, scanSbfFrames } from './sbf.cjs';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelRange } from './novatel.cjs';
+export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.cjs';
+export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
 export { DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SppOptions, SppSolution, gimSlantIonoDelayL1, gimVerticalTec, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, satClockCorrection, solveDgnss, solveSpp, toRtkEpoch } from './positioning.cjs';

package/dist/index.d.ts

@@ -4,7 +4,7 @@ export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSat
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.js';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.js';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.js';
-export { CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, Sp3File, Sp3Sample, WarningAccumulator, WarningSeverity, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.js';
+export { CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, Sp3File, Sp3Sample, WarningAccumulator, WarningSeverity, WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.js';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, 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';
 export { C as CnavEphemeris, E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, R as RinexCnavEphemeris, p as parseNavFile } from './nav-DImXUdbp.js';
 export { I as IonexGrid, p as parseIonex } from './ionex-DaW5x9nq.js';
@@ -13,6 +13,8 @@ export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES,
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.js';
 export { HasMessage, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGloNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfParseResult, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGloNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, sbfCrc16, scanSbfFrames } from './sbf.js';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelRange } from './novatel.js';
+export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.js';
+export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
 export { DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SppOptions, SppSolution, gimSlantIonoDelayL1, gimVerticalTec, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, satClockCorrection, solveDgnss, solveSpp, toRtkEpoch } from './positioning.js';

package/dist/index.js

@@ -1,3 +1,13 @@
+import {
+  frequencyLabel,
+  parseAntex
+} from "./chunk-5CMSYIYM.js";
+import {
+  computeStats,
+  nmeaCoordToDecimal,
+  parseNmeaFile,
+  verifyChecksum
+} from "./chunk-YDMYO3YK.js";
 import {
   BAND_PRESETS,
   CONSTELLATIONS,
@@ -40,7 +50,7 @@ import {
   parseUbxRawNav,
   parseUbxRawx,
   ubxFrames
-} from "./chunk-MCAVU3XL.js";
+} from "./chunk-FEOMYZEU.js";
 import {
   crc16,
   parseSbfAlmanac,
@@ -53,15 +63,29 @@ import {
   parseSbfMeas,
   parseSbfNav,
   scanSbfFrames
-} from "./chunk-PZNVFCKJ.js";
-import "./chunk-VTX3VCL4.js";
+} from "./chunk-AFEDNWXI.js";
+import "./chunk-5AUK42AL.js";
 import {
   OEM4_HLEN,
   crc32,
   oem4Frames,
   parseNovatelNav,
   parseNovatelRange
-} from "./chunk-7RKLWE6R.js";
+} from "./chunk-F4CM6Q4F.js";
+import {
+  parseTrimble,
+  parseTrimbleNav,
+  trimbleFrames
+} from "./chunk-CKRNXZHT.js";
+import {
+  binexCrc16,
+  binexCrc32,
+  binexCsum8,
+  binexRecords,
+  decodeBinexEph,
+  getBnxi,
+  parseBinex
+} from "./chunk-OK4GGNOO.js";
 import {
   connectToMountpoint,
   fetchSourcetable,
@@ -78,16 +102,6 @@ import {
   detrendIonoArcs,
   parseSinexBiasDcb
 } from "./chunk-Z5UZUNCB.js";
-import {
-  frequencyLabel,
-  parseAntex
-} from "./chunk-5CMSYIYM.js";
-import {
-  computeStats,
-  nmeaCoordToDecimal,
-  parseNmeaFile,
-  verifyChecksum
-} from "./chunk-YDMYO3YK.js";
 import {
   Scale,
   getDateFromDayOfWeek,
@@ -190,14 +204,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,
@@ -400,6 +416,10 @@ export {
   analyzeQuality,
   applyHelmert,
   applyIonoDcb,
+  binexCrc16,
+  binexCrc32,
+  binexCsum8,
+  binexRecords,
   buildGloChannelMap,
   buildObsIndices,
   clampUnit,
@@ -420,6 +440,7 @@ export {
   crxDecompress,
   crxRepair,
   dateToEpoch,
+  decodeBinexEph,
   decodeEphemeris,
   decodeMsmFull,
   deg2dms,
@@ -441,6 +462,7 @@ export {
   geodeticToUtm,
   getAer,
   getBdsTime,
+  getBnxi,
   getDateFromBdsTime,
   getDateFromDayOfWeek,
   getDateFromDayOfYear,
@@ -514,6 +536,7 @@ export {
   padL,
   padR,
   parseAntex,
+  parseBinex,
   parseCrxDataLine,
   parseIonex,
   parseNavFile,
@@ -533,6 +556,8 @@ export {
   parseSinexBiasDcb,
   parseSourcetable,
   parseSp3,
+  parseTrimble,
+  parseTrimbleNav,
   parseUbxCnav,
   parseUbxIonoUtc,
   parseUbxNav,
@@ -563,11 +588,13 @@ export {
   systemName,
   toRtkEpoch,
   transformFrame,
+  trimbleFrames,
   ubxFrames,
   updateStationMeta,
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,

package/dist/novatel.js

@@ -4,8 +4,9 @@ import {
   oem4Frames,
   parseNovatelNav,
   parseNovatelRange
-} from "./chunk-7RKLWE6R.js";
-import "./chunk-WZDFZ76K.js";
+} from "./chunk-F4CM6Q4F.js";
+import "./chunk-726LQBGM.js";
+import "./chunk-IQ5OIMQD.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";
 export {

package/dist/rinex.cjs

@@ -40,6 +40,7 @@ __export(rinex_exports, {
   stationHeaderLines: () => stationHeaderLines,
   systemCmp: () => systemCmp,
   systemName: () => systemName,
+  writeCrx: () => writeCrx,
   writeModernObsBlob: () => writeModernObsBlob,
   writeRinex2ObsBlob: () => writeRinex2ObsBlob,
   writeRinex4ObsBlob: () => writeRinex4ObsBlob,
@@ -134,6 +135,442 @@ function crxDecompress(prev, field) {
   };
 }
 
+// src/rinex/crx-writer.ts
+var ARC_ORDER = 3;
+function emptyData() {
+  return { u: [0, 0, 0, 0], l: [0, 0, 0, 0], order: -1 };
+}
+function cAtol(s) {
+  let i = 0;
+  while (i < s.length && (s[i] === " " || s[i] === "	")) i++;
+  let sign = 1;
+  if (s[i] === "+" || s[i] === "-") {
+    if (s[i] === "-") sign = -1;
+    i++;
+  }
+  let n = 0;
+  let any = false;
+  while (i < s.length && s[i] >= "0" && s[i] <= "9") {
+    n = n * 10 + (s.charCodeAt(i) - 48);
+    i++;
+    any = true;
+  }
+  return any ? sign * n : 0;
+}
+function strdiff(s1, s2) {
+  let ds = "";
+  const n = Math.min(s1.length, s2.length);
+  let i = 0;
+  for (; i < n; i++) {
+    if (s2[i] === s1[i]) ds += " ";
+    else if (s2[i] === " ") ds += "&";
+    else ds += s2[i];
+  }
+  if (s1.length > s2.length) {
+    for (let k = i; k < s1.length; k++) ds += s1[k] === " " ? " " : "&";
+  } else if (s2.length > s1.length) {
+    ds += s2.slice(i);
+  }
+  let end = ds.length;
+  while (end > 0 && ds[end - 1] === " ") end--;
+  return ds.slice(0, end);
+}
+function readValue(field14) {
+  const c = field14.padEnd(14, " ").split("");
+  const p7 = c[7];
+  c[10] = c[9];
+  c[9] = c[8];
+  let l = cAtol(c.slice(9, 14).join(""));
+  let u;
+  if (p7 === " ") {
+    u = 0;
+  } else if (p7 === "-") {
+    u = 0;
+    l = -l;
+  } else {
+    c[8] = ".";
+    u = cAtol(c.slice(0, 14).join(""));
+    if (u < 0) l = -l;
+  }
+  return { u, l };
+}
+function putdiff(dddu, dddl) {
+  dddu += Math.trunc(dddl / 1e5);
+  dddl %= 1e5;
+  if (dddu < 0 && dddl > 0) {
+    dddu++;
+    dddl -= 1e5;
+  } else if (dddu > 0 && dddl < 0) {
+    dddu--;
+    dddl += 1e5;
+  }
+  if (dddu === 0) return String(dddl);
+  return `${dddu}${String(Math.abs(dddl)).padStart(5, "0")}`;
+}
+function putClock(du, dl, cOrder) {
+  du += Math.trunc(dl / 1e8);
+  dl %= 1e8;
+  if (du < 0 && dl > 0) {
+    du++;
+    dl -= 1e8;
+  } else if (du > 0 && dl < 0) {
+    du--;
+    dl += 1e8;
+  }
+  let out = "";
+  if (cOrder === 0) out += `${ARC_ORDER}&`;
+  if (du === 0) out += String(dl);
+  else out += `${du}${String(Math.abs(dl)).padStart(8, "0")}`;
+  return out;
+}
+function chop(line) {
+  let end = line.length;
+  if (end > 0 && line[end - 1] === "\r") end--;
+  while (end > 0 && line[end - 1] === " ") end--;
+  return line.slice(0, end);
+}
+function readSat(lines, cursor, rinexVersion, ntype, ntypeGnss) {
+  let line = chop(lines[cursor.li++] ?? "");
+  let maxField;
+  let ntypeRec;
+  let firstRec;
+  let satId = "";
+  if (rinexVersion === 2) {
+    maxField = 5;
+    ntypeRec = ntype;
+    firstRec = 0;
+  } else {
+    satId = line.slice(0, 3);
+    const sys = line[0];
+    maxField = ntypeRec = ntypeGnss[sys] ?? -1;
+    if (maxField < 0) return null;
+    firstRec = 3;
+  }
+  const data = [];
+  let flag = "";
+  for (let i = 0; i < ntypeRec; i += maxField) {
+    const nfield = Math.min(ntypeRec - i, maxField);
+    const need = firstRec + 16 * nfield;
+    let body = line;
+    if (body.length < need) body = body.padEnd(need, " ");
+    for (let j = 0; j < nfield; j++) {
+      const p = firstRec + j * 16;
+      const rec = body.slice(p, p + 16);
+      const d = emptyData();
+      if (rec[10] === ".") {
+        flag += rec[14] ?? " ";
+        flag += rec[15] ?? " ";
+        const v = readValue(rec.slice(0, 14));
+        d.u[0] = v.u;
+        d.l[0] = v.l;
+        d.order = 0;
+      } else if (rec.slice(0, 14).trim() === "") {
+        flag += rec[14] ?? " ";
+        flag += rec[15] ?? " ";
+        d.order = -1;
+      } else {
+        return null;
+      }
+      data.push(d);
+    }
+    if (i + maxField < ntypeRec) {
+      line = chop(lines[cursor.li++] ?? "");
+    }
+  }
+  return { sat: { data, flag, ntypeRec }, satId };
+}
+var MONTHS = [
+  "Jan",
+  "Feb",
+  "Mar",
+  "Apr",
+  "May",
+  "Jun",
+  "Jul",
+  "Aug",
+  "Sep",
+  "Oct",
+  "Nov",
+  "Dec"
+];
+function crinexDate() {
+  const d = /* @__PURE__ */ new Date();
+  const dd = String(d.getUTCDate()).padStart(2, "0");
+  const mon = MONTHS[d.getUTCMonth()];
+  const yy = String(d.getUTCFullYear() % 100).padStart(2, "0");
+  const hh = String(d.getUTCHours()).padStart(2, "0");
+  const mm = String(d.getUTCMinutes()).padStart(2, "0");
+  return `${dd}-${mon}-${yy} ${hh}:${mm}`;
+}
+function writeCrx(rinexText, opts = {}) {
+  const rawLines = rinexText.split("\n");
+  const cursor = { li: 0 };
+  const out = [];
+  const firstLine = chop(rawLines[cursor.li++] ?? "");
+  if (firstLine.slice(60, 80) !== "RINEX VERSION / TYPE" || firstLine[20] !== "O") {
+    throw new Error("writeCrx: not a RINEX observation file");
+  }
+  const rinexVersion = parseInt(firstLine, 10);
+  const rinexVersionFull = readVersionFull(firstLine);
+  if (rinexVersionFull < 0) throw new Error("writeCrx: invalid RINEX version");
+  const crxVersion = rinexVersionFull >= 402 ? "3.1" : rinexVersion === 2 ? "1.0" : "3.0";
+  out.push(
+    `${crxVersion.padEnd(20)}${"COMPACT RINEX FORMAT".padEnd(40)}${"CRINEX VERS   / TYPE".padEnd(20)}`
+  );
+  const prog = (opts.prog ?? "gnss-js writeCrx").slice(0, 40);
+  const date = opts.date ?? crinexDate();
+  out.push(`${prog.padEnd(40)}${date.padEnd(20)}CRINEX PROG / DATE`);
+  out.push(firstLine);
+  let ntype = 0;
+  const ntypeGnss = {};
+  let line;
+  do {
+    line = chop(rawLines[cursor.li++] ?? "");
+    out.push(line);
+    const label = line.slice(60).trimStart();
+    if (label.startsWith("# / TYPES OF OBSERV") && line[5] !== " ") {
+      ntype = parseInt(line, 10);
+    } else if (label.startsWith("SYS / # / OBS TYPES") && line[0] !== " ") {
+      ntypeGnss[line[0]] = parseInt(line.slice(3), 10);
+    }
+    if (cursor.li >= rawLines.length) break;
+  } while (!line.slice(60).trimStart().startsWith("END OF HEADER"));
+  const v2 = rinexVersion === 2;
+  const off = v2 ? { event: 28, nsat: 29, satlst: 32, clock: 68, satOld: 32, shiftClk: 1 } : {
+    event: 31,
+    nsat: 32,
+    satlst: 41,
+    clock: 41,
+    satOld: 41,
+    shiftClk: 4,
+    psec: 57
+  };
+  let oldline = "&";
+  let nsatOld = 0;
+  let clk0 = { u: [0, 0, 0, 0], l: [0, 0, 0, 0] };
+  let dy0 = [];
+  let flag0 = [];
+  let oldpsec = "";
+  let clkOrderPrev = -1;
+  const resetArcs = () => {
+    oldline = "&";
+    nsatOld = 0;
+    clkOrderPrev = -1;
+  };
+  for (; ; ) {
+    if (cursor.li >= rawLines.length) break;
+    let epoch = rawLines[cursor.li];
+    if (epoch === void 0) break;
+    if (epoch === "" && cursor.li === rawLines.length - 1) break;
+    cursor.li++;
+    epoch = chop(epoch);
+    if (epoch === "") continue;
+    let newline;
+    if (v2) {
+      newline = epoch;
+    } else {
+      newline = epoch.length < 41 ? epoch.padEnd(41, " ") : epoch;
+    }
+    const eventFlag = parseInt(newline[off.event] ?? "0", 10) || 0;
+    if (eventFlag > 1) {
+      putEventData(rawLines, cursor, v2, newline, out, ntypeGnss, (n) => {
+        ntype = n;
+      });
+      resetArcs();
+      clk0 = { u: [0, 0, 0, 0], l: [0, 0, 0, 0] };
+      dy0 = [];
+      flag0 = [];
+      continue;
+    }
+    let newpsec = "";
+    if (rinexVersionFull >= 402 && newline.length === 62) {
+      const psec = newline.slice(off.psec, off.psec + 5);
+      if (/^\d{5}$/.test(psec)) newpsec = psec;
+      newline = newline.slice(0, off.psec).replace(/\s+$/, "");
+    }
+    let clkOrder = -1;
+    const clk1 = { u: [0, 0, 0, 0], l: [0, 0, 0, 0] };
+    if (newline.length > off.clock) {
+      clkOrder = readClock(
+        newline.slice(off.clock),
+        off.shiftClk,
+        clk1,
+        clkOrderPrev
+      );
+      if (clkOrder > -1) {
+        newline = newline.slice(0, off.clock).replace(/\s+$/, "");
+        if (!v2 && newline.length < 41) newline = newline.padEnd(41, " ");
+      }
+    }
+    clkOrderPrev = clkOrder;
+    const nsat = parseInt((newline.slice(off.nsat) || "0").trim(), 10) || 0;
+    if (v2 && nsat > 12) {
+      newline = readMoreSat(rawLines, cursor, nsat, newline);
+    }
+    const sats = [];
+    const satIds = [];
+    let bad = false;
+    let satlist = v2 ? newline.slice(off.satlst, off.satlst + nsat * 3) : "";
+    for (let i = 0; i < nsat; i++) {
+      const r = readSat(rawLines, cursor, rinexVersion, ntype, ntypeGnss);
+      if (!r) {
+        bad = true;
+        break;
+      }
+      sats.push(r.sat);
+      satIds.push(v2 ? satlist.slice(i * 3, i * 3 + 3) : r.satId);
+    }
+    if (bad) continue;
+    if (!v2) satlist = satIds.join("");
+    const fullNew = v2 ? newline : newline.slice(0, 41) + satlist;
+    const oldSatList = v2 ? oldline.slice(off.satOld, off.satOld + nsatOld * 3) : oldline.slice(off.satOld);
+    const sattbl = [];
+    const dup = /* @__PURE__ */ new Set();
+    for (let i = 0; i < nsat; i++) {
+      const id = satIds[i];
+      if (dup.has(id)) {
+        bad = true;
+        break;
+      }
+      dup.add(id);
+      let found = -1;
+      for (let j = 0; j < nsatOld; j++) {
+        if (oldSatList.slice(j * 3, j * 3 + 3) === id) {
+          found = j;
+          break;
+        }
+      }
+      sattbl.push(found);
+    }
+    if (bad) continue;
+    out.push(strdiff(oldline, fullNew));
+    let clockLine = "";
+    if (clkOrder > -1) {
+      if (clkOrder > 0) processClock(clk1, clk0, clkOrder);
+      clockLine += putClock(clk1.u[clkOrder], clk1.l[clkOrder], clkOrder);
+    }
+    if (newpsec.length === 5) {
+      clockLine += " " + strdiff(oldpsec, newpsec);
+    }
+    out.push(clockLine);
+    const dataLines = buildData(sats, sattbl, dy0, flag0, rinexVersion);
+    for (const dl of dataLines) out.push(dl);
+    oldline = fullNew;
+    nsatOld = nsat;
+    clk0 = clk1;
+    oldpsec = newpsec;
+    dy0 = sats.map((s) => s.data);
+    flag0 = sats.map((s) => s.flag);
+  }
+  return out.join("\n") + "\n";
+}
+function readClock(field, shiftClk, clk1, prevOrder) {
+  const c = field.split("");
+  if (c[2] !== ".") return -1;
+  for (let k = 0; k < shiftClk; k++) c[2 + k] = c[3 + k];
+  c[2 + shiftClk] = ".";
+  const s = c.join("");
+  const m = /(-?\d+)\.(\d+)/.exec(s);
+  if (!m) return -1;
+  clk1.u[0] = parseInt(m[1], 10);
+  clk1.l[0] = parseInt(m[2], 10);
+  if (c[0] === "-" || c[1] === "-") clk1.l[0] = -clk1.l[0];
+  return Math.min(prevOrder + 1, ARC_ORDER);
+}
+function processClock(clk1, clk0, clkOrder) {
+  for (let i = 0; i < clkOrder; i++) {
+    clk1.u[i + 1] = clk1.u[i] - clk0.u[i];
+    clk1.l[i + 1] = clk1.l[i] - clk0.l[i];
+  }
+}
+function takeDiff(py1, py0) {
+  py1.order = py0.order;
+  if (py1.order < ARC_ORDER) py1.order++;
+  for (let k = 0; k < py1.order; k++) {
+    py1.u[k + 1] = py1.u[k] - py0.u[k];
+    py1.l[k + 1] = py1.l[k] - py0.l[k];
+  }
+}
+function buildData(sats, sattbl, dy0, flag0, rinexVersion) {
+  const lines = [];
+  for (let i = 0; i < sats.length; i++) {
+    const sat = sats[i];
+    const i0 = sattbl[i];
+    const prevFlag = (i0 >= 0 ? flag0[i0] ?? "" : "").split("");
+    let line = "";
+    for (let j = 0; j < sat.ntypeRec; j++) {
+      const py1 = sat.data[j];
+      if (py1.order >= 0) {
+        if (i0 < 0 || dy0[i0]?.[j]?.order === -1) {
+          py1.order = 0;
+          line += `${ARC_ORDER}&`;
+        } else {
+          takeDiff(py1, dy0[i0][j]);
+          if (Math.abs(py1.u[py1.order]) > 1e5) {
+            py1.order = 0;
+            line += `${ARC_ORDER}&`;
+          }
+        }
+        line += putdiff(py1.u[py1.order], py1.l[py1.order]);
+      } else if (i0 >= 0 && rinexVersion === 2) {
+        prevFlag[j * 2] = " ";
+        prevFlag[j * 2 + 1] = " ";
+      }
+      line += " ";
+    }
+    if (i0 < 0 && rinexVersion !== 2) {
+      for (const ch of sat.flag) line += ch === " " ? "&" : ch;
+    } else {
+      line += i0 < 0 ? strdiff("", sat.flag) : strdiff(prevFlag.join(""), sat.flag);
+      line = line.replace(/ +$/, "");
+    }
+    lines.push(line);
+  }
+  return lines;
+}
+function putEventData(lines, cursor, v2, epochLine, out, ntypeGnss, setNtype) {
+  if (v2) {
+    out.push("&" + epochLine.slice(1));
+    if (epochLine.length > 29) {
+      const n = parseInt(epochLine.slice(29), 10) || 0;
+      for (let i = 0; i < n; i++) {
+        const l = chop(lines[cursor.li++] ?? "");
+        out.push(l);
+        if (l.slice(60).trimStart().startsWith("# / TYPES OF OBSERV") && l[5] !== " ") {
+          setNtype(parseInt(l, 10));
+        }
+      }
+    }
+  } else {
+    out.push(epochLine.replace(/\s+$/, ""));
+    const n = parseInt(epochLine.slice(32), 10) || 0;
+    for (let i = 0; i < n; i++) {
+      const l = chop(lines[cursor.li++] ?? "");
+      out.push(l);
+      if (l.slice(60).trimStart().startsWith("SYS / # / OBS TYPES") && l[0] !== " ") {
+        ntypeGnss[l[0]] = parseInt(l.slice(3), 10);
+      }
+    }
+  }
+}
+function readMoreSat(lines, cursor, n, epochLine) {
+  let result = epochLine;
+  let remaining = n;
+  while (remaining > 12) {
+    const l = lines[cursor.li++] ?? "";
+    const chopped = chop(l);
+    result += chopped[2] === " " ? chopped.slice(32) : chopped;
+    remaining -= 12;
+  }
+  return result;
+}
+function readVersionFull(line) {
+  if (line[5] === "2" && line[6] === " ") return 200;
+  const d = (c) => c !== void 0 && c >= "0" && c <= "9";
+  if (!d(line[5]) || line[6] !== "." || !d(line[7]) || !d(line[8])) return -1;
+  return (line.charCodeAt(5) - 48) * 100 + (line.charCodeAt(7) - 48) * 10 + (line.charCodeAt(8) - 48);
+}
+
 // src/rinex/format.ts
 function padL(s, w) {
   return s.length >= w ? s.slice(0, w) : s + " ".repeat(w - s.length);
@@ -2346,6 +2783,7 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
   stationHeaderLines,
   systemCmp,
   systemName,
+  writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,
   writeRinex4ObsBlob,