gnss-js 2.9.0 → 2.10.0

package/dist/msm-DWtoTvx2.d.cts

@@ -0,0 +1,108 @@
+import { b as Rtcm3Frame } from './ephemeris-BAnBLLB9.cjs';
+
+/**
+ * Full RTCM3 MSM decoder — extracts pseudorange, carrier phase,
+ * Doppler, and C/N0 for every satellite × signal cell.
+ *
+ * Reference: BNC 2.13.4 RTCM3Decoder.cpp (BKG), RTKLIB rtcm3.c and
+ * RTCM 10403.3 §3.5.
+ *
+ * MSM4–7 carry full observations. MSM1–3 (compact: MSM1 code, MSM2 phase,
+ * MSM3 both; no C/N0, no Doppler) carry each satellite's rough range only
+ * modulo 1 ms (DF398, no DF397 integer milliseconds): their epochs come back
+ * with `rangeModulo1ms` set, ranges complete except for a whole number of
+ * light-milliseconds per satellite — see `resolveMsmRanges` (positioning).
+ */
+
+/** Single signal observation for one satellite */
+interface ObsSignal {
+    /** RINEX 2-char obs code, e.g. "1C", "5X" */
+    rinexCode: string;
+    /** Pseudorange in meters (undefined if invalid) */
+    pseudorange?: number;
+    /** Carrier phase in cycles (undefined if invalid) */
+    phase?: number;
+    /** Doppler in Hz (undefined if invalid) */
+    doppler?: number;
+    /** Signal-to-noise ratio in dB-Hz */
+    cn0?: number;
+    /** Lock time in seconds */
+    lockTime?: number;
+    /** Half-cycle ambiguity indicator */
+    halfCycle?: boolean;
+    /** Wavelength used (meters) — needed for RINEX phase */
+    wavelength: number;
+}
+/** All observations for one satellite in one epoch */
+interface ObsSatObs {
+    /** PRN string, e.g. "G01", "R15", "E04" */
+    prn: string;
+    /** System letter: G, R, E, C, J, I, S */
+    system: string;
+    /** Signal observations */
+    signals: ObsSignal[];
+}
+/**
+ * One decoded observation epoch. Originally the MSM4–7 output shape, now the
+ * neutral shape every RTCM3 obs source fills: MSM (`decodeMsmFull`), legacy
+ * 1001–1012 (`decodeLegacyObs`), via the `decodeObs` dispatcher.
+ */
+interface ObsEpoch {
+    /** Message type (1074, 1077, etc.) */
+    messageType: number;
+    /** Epoch time as GPS milliseconds of week (or GLONASS day-of-week ms, BDS SOW) */
+    epochMs: number;
+    /** System letter */
+    system: string;
+    /** Observations per satellite */
+    observations: ObsSatObs[];
+    /**
+     * Set for MSM1–3: every pseudorange and phase is short by an unknown whole
+     * number of light-milliseconds (~299.79 km) per satellite. Resolve with
+     * `resolveMsmRanges` before using the ranges.
+     */
+    rangeModulo1ms?: true;
+}
+/** Options for the RTCM3 observation decoders. */
+interface ObsDecodeOptions {
+    /**
+     * Also decode compact MSM1–3. Their epochs carry `rangeModulo1ms`: every
+     * range is short by whole light-milliseconds until `resolveMsmRanges`
+     * completes it, so a consumer must opt in (default false → null, as
+     * before).
+     */
+    compactMsm?: boolean;
+}
+/**
+ * Decode a single RTCM3 MSM frame into observations: MSM4–7, and MSM1–3
+ * with `opts.compactMsm` (flagged `rangeModulo1ms`). Null otherwise.
+ */
+declare function decodeMsmFull(frame: Rtcm3Frame, opts?: ObsDecodeOptions): ObsEpoch | null;
+/**
+ * Convert an observation epoch timestamp to a Date.
+ * For GPS/Galileo/BDS/QZSS: epochMs is ms of GPS week.
+ * For GLONASS: epochMs is day-of-week × 86400000 + ms of day (UTC+3).
+ *
+ * Requires an approximate reference time to resolve the current week.
+ */
+/**
+ * NOTE: returns a GPS-scale Date — the clock-face convention used by
+ * the RINEX parser and every consumer of epoch milliseconds in this
+ * library (satellite propagation, SPP, RINEX epoch tags). It is NOT
+ * UTC: converting to UTC here put every downstream satellite position
+ * 18 s (the current leap-second count) off — a ±14 km per-satellite
+ * range error that scattered live SPP solutions by kilometres.
+ */
+declare function obsEpochToDate(sys: string, epochMs: number, refTime?: Date): Date;
+/**
+ * Set a GLONASS FDMA frequency number (k = -7..+13) for a slot (1-24),
+ * e.g. from a decoded 1020 ephemeris. MSM4/6 carry no DF419, so without
+ * this the wavelength — and hence carrier phase — stays unknown.
+ */
+declare function setGloFreqNumber(slot: number, k: number): void;
+/**
+ * Reset the GLONASS frequency number cache (e.g., when switching mountpoints).
+ */
+declare function resetGloFreqCache(): void;
+
+export { type ObsDecodeOptions as O, type ObsEpoch as a, type ObsSatObs as b, type ObsSignal as c, decodeMsmFull as d, obsEpochToDate as o, resetGloFreqCache as r, setGloFreqNumber as s };

package/dist/positioning.cjs

@@ -49,6 +49,7 @@ __export(positioning_exports, {
   nlWavelength: () => nlWavelength,
   postProcessRtk: () => postProcessRtk,
   rcvAntennaRangeM: () => rcvAntennaRangeM,
+  resolveMsmRanges: () => resolveMsmRanges,
   resolvePppAmbiguities: () => resolvePppAmbiguities,
   satClockCorrection: () => satClockCorrection,
   sbasIonoDeg: () => sbasIonoDeg,
@@ -2134,6 +2135,66 @@ var ZeroBaselineEngine = class {
   }
 };
 
+// src/positioning/msm-resolve.ts
+var LIGHT_MS = C_LIGHT / 1e3;
+var MAX_FRACTION_MISFIT = 0.25;
+function predictedRange(eph, timeMs, rx) {
+  let tau = 0.075;
+  for (let i = 0; i < 2; i++) {
+    const tEmit = timeMs - tau * 1e3;
+    const p = computeSatPosition(eph, tEmit);
+    if (!Number.isFinite(p.x)) return null;
+    const rho = Math.hypot(p.x - rx[0], p.y - rx[1], p.z - rx[2]);
+    tau = rho / C_LIGHT;
+    if (i === 1) return rho - C_LIGHT * satClockCorrection(eph, tEmit);
+  }
+  return null;
+}
+function moduloRange(sat) {
+  for (const s of sat.signals)
+    if (s.pseudorange !== void 0) return s.pseudorange;
+  for (const s of sat.signals)
+    if (s.phase !== void 0 && s.wavelength > 0)
+      return s.phase * s.wavelength;
+  return null;
+}
+function resolveMsmRanges(epoch, opts) {
+  if (!epoch.rangeModulo1ms) return epoch;
+  const cand = [];
+  for (const sat of epoch.observations) {
+    const eph = opts.ephFor(sat.prn);
+    const mod = moduloRange(sat);
+    if (!eph || mod === null) continue;
+    const pred = predictedRange(eph, opts.timeMs, opts.rxEcef);
+    if (pred === null) continue;
+    cand.push({ sat, k: (pred - mod) / LIGHT_MS });
+  }
+  const { rangeModulo1ms: _drop, ...rest } = epoch;
+  if (cand.length === 0) return { ...rest, observations: [] };
+  let sx = 0;
+  let sy = 0;
+  for (const c of cand) {
+    sx += Math.cos(2 * Math.PI * c.k);
+    sy += Math.sin(2 * Math.PI * c.k);
+  }
+  const b = Math.atan2(sy, sx) / (2 * Math.PI);
+  const observations = [];
+  for (const { sat, k } of cand) {
+    const n = Math.round(k - b);
+    if (Math.abs(k - b - n) > MAX_FRACTION_MISFIT) continue;
+    const add = n * LIGHT_MS;
+    observations.push({
+      ...sat,
+      signals: sat.signals.map((s) => ({
+        ...s,
+        ...s.pseudorange !== void 0 ? { pseudorange: s.pseudorange + add } : {},
+        ...s.phase !== void 0 && s.wavelength > 0 ? { phase: s.phase + add / s.wavelength } : {}
+      }))
+    });
+  }
+  return { ...rest, observations };
+}
+
 // src/navbits/index.ts
 var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK = 7 * 86400;
@@ -4633,6 +4694,7 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   nlWavelength,
   postProcessRtk,
   rcvAntennaRangeM,
+  resolveMsmRanges,
   resolvePppAmbiguities,
   satClockCorrection,
   sbasIonoDeg,

package/dist/positioning.d.cts

@@ -2,6 +2,7 @@ import { E as Ephemeris } from './nav-DImXUdbp.cjs';
 import { DopValues } from './orbit.cjs';
 import { I as IonexGrid, C as ClkFile, S as Sp3File } from './clk-DGVBGRzA.cjs';
 import { AntexFile } from './antex.cjs';
+import { a as ObsEpoch } from './msm-DWtoTvx2.cjs';
 import { a as HasOrbitCorrection, b as HasClockCorrection } from './has-BaXnQarR.cjs';
 import './ephemeris-BAnBLLB9.cjs';
 import './ecef-CF0uAysr.cjs';
@@ -613,6 +614,23 @@ declare class ZeroBaselineEngine {
     private processReceiver;
 }
 
+interface MsmResolveOptions {
+    /** Epoch time, GPS-scale ms (as `obsEpochToDate` returns). */
+    timeMs: number;
+    /** Approximate receiver ECEF (m) — the station's 1005/1006, a sourcetable
+     *  position, or a previous fix. Tens of km of error are tolerated. */
+    rxEcef: readonly [number, number, number];
+    /** Broadcast ephemeris for a PRN at this time, or null. */
+    ephFor: (prn: string) => Ephemeris | null;
+}
+/**
+ * Add the missing whole light-milliseconds to an MSM1–3 epoch's ranges.
+ * Returns the epoch with `rangeModulo1ms` cleared, keeping only satellites
+ * whose range could be resolved (an ephemeris, and a prediction consistent
+ * with the rest of the epoch). Epochs that aren't modulo-1-ms pass through.
+ */
+declare function resolveMsmRanges(epoch: ObsEpoch, opts: MsmResolveOptions): ObsEpoch;
+
 /**
  * Integer ambiguity resolution: LAMBDA decorrelation + MLAMBDA
  * integer least-squares search.
@@ -1714,4 +1732,4 @@ declare function sbasProtectionLevels(sats: readonly {
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type MsmResolveOptions, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -2,6 +2,7 @@ import { E as Ephemeris } from './nav-DImXUdbp.js';
 import { DopValues } from './orbit.js';
 import { I as IonexGrid, C as ClkFile, S as Sp3File } from './clk-DGVBGRzA.js';
 import { AntexFile } from './antex.js';
+import { a as ObsEpoch } from './msm-CcQJV3PC.js';
 import { a as HasOrbitCorrection, b as HasClockCorrection } from './has-BaXnQarR.js';
 import './ephemeris-BAnBLLB9.js';
 import './ecef-CF0uAysr.js';
@@ -613,6 +614,23 @@ declare class ZeroBaselineEngine {
     private processReceiver;
 }
 
+interface MsmResolveOptions {
+    /** Epoch time, GPS-scale ms (as `obsEpochToDate` returns). */
+    timeMs: number;
+    /** Approximate receiver ECEF (m) — the station's 1005/1006, a sourcetable
+     *  position, or a previous fix. Tens of km of error are tolerated. */
+    rxEcef: readonly [number, number, number];
+    /** Broadcast ephemeris for a PRN at this time, or null. */
+    ephFor: (prn: string) => Ephemeris | null;
+}
+/**
+ * Add the missing whole light-milliseconds to an MSM1–3 epoch's ranges.
+ * Returns the epoch with `rangeModulo1ms` cleared, keeping only satellites
+ * whose range could be resolved (an ephemeris, and a prediction consistent
+ * with the rest of the epoch). Epochs that aren't modulo-1-ms pass through.
+ */
+declare function resolveMsmRanges(epoch: ObsEpoch, opts: MsmResolveOptions): ObsEpoch;
+
 /**
  * Integer ambiguity resolution: LAMBDA decorrelation + MLAMBDA
  * integer least-squares search.
@@ -1714,4 +1732,4 @@ declare function sbasProtectionLevels(sats: readonly {
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type MsmResolveOptions, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -28,6 +28,7 @@ import {
   nlWavelength,
   postProcessRtk,
   rcvAntennaRangeM,
+  resolveMsmRanges,
   resolvePppAmbiguities,
   satClockCorrection,
   sbasIonoDeg,
@@ -40,7 +41,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-NOGA5KD7.js";
+} from "./chunk-K3GLUWZI.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-7T5DLVZL.js";
 import "./chunk-HVXYFUCB.js";
@@ -82,6 +83,7 @@ export {
   nlWavelength,
   postProcessRtk,
   rcvAntennaRangeM,
+  resolveMsmRanges,
   resolvePppAmbiguities,
   satClockCorrection,
   sbasIonoDeg,

package/dist/rinex.cjs

@@ -32,6 +32,7 @@ __export(rinex_exports, {
   fmtD: () => fmtD,
   fmtF: () => fmtF,
   hdrLine: () => hdrLine,
+  modalIntervalS: () => modalIntervalS,
   padL: () => padL,
   padR: () => padR,
   parseBiasSinex: () => parseBiasSinex,
@@ -2600,6 +2601,21 @@ function writeRinexObsBlob(header, epochs, obsTypes) {
     comment: "Merged from multiple RINEX files"
   });
 }
+function modalIntervalS(epochs) {
+  const counts = /* @__PURE__ */ new Map();
+  for (let i = 1; i < epochs.length; i++) {
+    const d = Math.round(epochs[i].time - epochs[i - 1].time);
+    if (d > 0) counts.set(d, (counts.get(d) ?? 0) + 1);
+  }
+  let best = null;
+  let bestN = 0;
+  for (const [d, n] of counts)
+    if (n > bestN || n === bestN && best !== null && d < best) {
+      best = d;
+      bestN = n;
+    }
+  return best === null ? null : best / 1e3;
+}
 async function writeModernObsBlob(header, epochs, obsTypes, opts) {
   if (epochs.length === 0) return new Blob([], { type: "application/gzip" });
   const sink = createGzipLineSink();
@@ -2646,11 +2662,9 @@ async function writeModernObsBlob(header, epochs, obsTypes, opts) {
       "TIME OF LAST OBS"
     )
   );
-  if (epochs.length >= 2) {
-    const interval = (epochs[1].time - epochs[0].time) / 1e3;
-    if (interval > 0 && interval < 3600) {
-      sink.push(hdrLine(fmtF(interval, 10, 3), "INTERVAL"));
-    }
+  const interval = modalIntervalS(epochs);
+  if (interval !== null && interval < 3600) {
+    sink.push(hdrLine(fmtF(interval, 10, 3), "INTERVAL"));
   }
   if (header.glonassSlots && Object.keys(header.glonassSlots).length > 0) {
     const entries = Object.entries(header.glonassSlots).sort(
@@ -2916,6 +2930,7 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
   fmtD,
   fmtF,
   hdrLine,
+  modalIntervalS,
   padL,
   padR,
   parseBiasSinex,

package/dist/rinex.d.cts

@@ -258,6 +258,15 @@ interface CompactEpoch {
  * the Float64Array index order for each system.
  */
 declare function writeRinexObsBlob(header: RinexHeader, epochs: CompactEpoch[], obsTypes: Map<string, string[]>): Promise<Blob>;
+/**
+ * The nominal observation interval (s): the most common spacing between
+ * consecutive epochs. The first difference alone mis-states it whenever the
+ * file opens with an irregular gap (a late second epoch, a burst). Null with
+ * fewer than two epochs.
+ */
+declare function modalIntervalS(epochs: readonly {
+    time: number;
+}[]): number | null;
 /** Shared RINEX 3/4 implementation (also used by obs-writer-v4.ts). */
 declare function writeModernObsBlob(header: RinexHeader, epochs: CompactEpoch[], obsTypes: Map<string, string[]>, opts: {
     version: string;
@@ -311,4 +320,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.d.ts

@@ -258,6 +258,15 @@ interface CompactEpoch {
  * the Float64Array index order for each system.
  */
 declare function writeRinexObsBlob(header: RinexHeader, epochs: CompactEpoch[], obsTypes: Map<string, string[]>): Promise<Blob>;
+/**
+ * The nominal observation interval (s): the most common spacing between
+ * consecutive epochs. The first difference alone mis-states it whenever the
+ * file opens with an irregular gap (a late second epoch, a burst). Null with
+ * fewer than two epochs.
+ */
+declare function modalIntervalS(epochs: readonly {
+    time: number;
+}[]): number | null;
 /** Shared RINEX 3/4 implementation (also used by obs-writer-v4.ts). */
 declare function writeModernObsBlob(header: RinexHeader, epochs: CompactEpoch[], obsTypes: Map<string, string[]>, opts: {
     version: string;
@@ -311,4 +320,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.js

@@ -7,6 +7,7 @@ import {
   fmtD,
   fmtF,
   hdrLine,
+  modalIntervalS,
   padL,
   padR,
   parseBiasSinex,
@@ -20,7 +21,7 @@ import {
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-6SWJTN6H.js";
+} from "./chunk-56HM5B6T.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -52,6 +53,7 @@ export {
   fmtD,
   fmtF,
   hdrLine,
+  modalIntervalS,
   padL,
   padR,
   parseBiasSinex,

package/dist/rtcm3.cjs

@@ -1287,12 +1287,14 @@ function satPrn(sys, satIdx1) {
   const num = sys === "S" ? satIdx1 + 119 : satIdx1;
   return `${prefix}${String(num).padStart(2, "0")}`;
 }
-function decodeMsmFull(frame) {
+function decodeMsmFull(frame, opts = {}) {
   const type = frame.messageType;
   const sys = msmSystem(type);
   if (!sys) return null;
   const variant = type % 10;
-  if (variant < 4 || variant > 7) return null;
+  if (variant < 1 || variant > 7) return null;
+  const compact = variant <= 3;
+  if (compact && !opts.compactMsm) return null;
   const payload = frame.payload;
   if (!payload || payload.length < 10) return null;
   try {
@@ -1337,7 +1339,9 @@ function decodeMsmFull(frame) {
     const rrmod = new Float64Array(numSat);
     const extsat = new Int8Array(numSat);
     const rdop = new Float64Array(numSat);
-    if (variant === 4 || variant === 6) {
+    if (compact) {
+      for (let j = 0; j < numSat; j++) rrmod[j] = bits.readU(10) / 1024;
+    } else if (variant === 4 || variant === 6) {
       for (let j = 0; j < numSat; j++) rrint[j] = bits.readU(8);
       for (let j = 0; j < numSat; j++) rrmod[j] = bits.readU(10) / 1024;
     } else {
@@ -1363,7 +1367,18 @@ function decodeMsmFull(frame) {
     psr.fill(-1e30);
     cp.fill(-1e30);
     dop.fill(-1e30);
-    if (variant === 4) {
+    if (variant === 1) {
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
+    } else if (variant === 2) {
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+    } else if (variant === 3) {
+      for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
+      for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
+      for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
+      for (let i = 0; i < activeCells; i++) hc[i] = bits.readU(1);
+    } else if (variant === 4) {
       for (let i = 0; i < activeCells; i++) psr[i] = bits.readS(15) / 2 ** 24;
       for (let i = 0; i < activeCells; i++) cp[i] = bits.readS(22) / 2 ** 29;
       for (let i = 0; i < activeCells; i++) ll[i] = bits.readU(4);
@@ -1439,8 +1454,10 @@ function decodeMsmFull(frame) {
         if (cnrVal > 0) {
           signal.cn0 = cnrVal;
         }
-        signal.lockTime = lockTimeSec(type, ll[cellIdx]);
-        signal.halfCycle = hc[cellIdx] === 1;
+        if (variant !== 1) {
+          signal.lockTime = lockTimeSec(type, ll[cellIdx]);
+          signal.halfCycle = hc[cellIdx] === 1;
+        }
         signals.push(signal);
         cellIdx++;
       }
@@ -1452,7 +1469,8 @@ function decodeMsmFull(frame) {
       messageType: type,
       epochMs,
       system: sys,
-      observations
+      observations,
+      ...compact ? { rangeModulo1ms: true } : {}
     };
   } catch (err) {
     reportDecodeError("msm", err);
@@ -1606,8 +1624,8 @@ function decodeLegacyObs(frame) {
   }
   return { messageType: type, epochMs, system, observations };
 }
-function decodeObs(frame) {
-  return decodeMsmFull(frame) ?? decodeLegacyObs(frame);
+function decodeObs(frame, opts = {}) {
+  return decodeMsmFull(frame, opts) ?? decodeLegacyObs(frame);
 }
 
 // src/rtcm3/ssr.ts

package/dist/rtcm3.d.cts

@@ -1,5 +1,7 @@
 import { b as Rtcm3Frame, E as EphemerisInfo } from './ephemeris-BAnBLLB9.cjs';
 export { B as BitReader, R as Rtcm3DebugHandler, a as Rtcm3Decoder, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.cjs';
+import { a as ObsEpoch, O as ObsDecodeOptions } from './msm-DWtoTvx2.cjs';
+export { b as ObsSatObs, c as ObsSignal, d as decodeMsmFull, o as obsEpochToDate, r as resetGloFreqCache, s as setGloFreqNumber } from './msm-DWtoTvx2.cjs';
 
 /** One code bias for a specific signal (DF380/DF381 + DF383). */
 interface SsrCodeBias {
@@ -116,91 +118,6 @@ interface IgsSsrIonoLayer {
  */
 declare function decodeIgsSsr(frame: Rtcm3Frame): IgsSsrMessage | null;
 
-/**
- * Full RTCM3 MSM4-7 decoder — extracts pseudorange, carrier phase,
- * Doppler, and C/N0 for every satellite × signal cell.
- *
- * Reference: BNC 2.13.4 RTCM3Decoder.cpp (BKG) and RTCM 10403.3 §3.5.
- *
- * Only MSM types 4-7 are decoded (full data). Types 1-3 carry partial
- * data and are skipped.
- */
-
-/** Single signal observation for one satellite */
-interface ObsSignal {
-    /** RINEX 2-char obs code, e.g. "1C", "5X" */
-    rinexCode: string;
-    /** Pseudorange in meters (undefined if invalid) */
-    pseudorange?: number;
-    /** Carrier phase in cycles (undefined if invalid) */
-    phase?: number;
-    /** Doppler in Hz (undefined if invalid) */
-    doppler?: number;
-    /** Signal-to-noise ratio in dB-Hz */
-    cn0?: number;
-    /** Lock time in seconds */
-    lockTime?: number;
-    /** Half-cycle ambiguity indicator */
-    halfCycle?: boolean;
-    /** Wavelength used (meters) — needed for RINEX phase */
-    wavelength: number;
-}
-/** All observations for one satellite in one epoch */
-interface ObsSatObs {
-    /** PRN string, e.g. "G01", "R15", "E04" */
-    prn: string;
-    /** System letter: G, R, E, C, J, I, S */
-    system: string;
-    /** Signal observations */
-    signals: ObsSignal[];
-}
-/**
- * One decoded observation epoch. Originally the MSM4–7 output shape, now the
- * neutral shape every RTCM3 obs source fills: MSM (`decodeMsmFull`), legacy
- * 1001–1012 (`decodeLegacyObs`), via the `decodeObs` dispatcher.
- */
-interface ObsEpoch {
-    /** Message type (1074, 1077, etc.) */
-    messageType: number;
-    /** Epoch time as GPS milliseconds of week (or GLONASS day-of-week ms, BDS SOW) */
-    epochMs: number;
-    /** System letter */
-    system: string;
-    /** Observations per satellite */
-    observations: ObsSatObs[];
-}
-/**
- * Decode a single RTCM3 MSM4-7 frame into full observations.
- * Returns null for non-MSM or MSM1-3 frames.
- */
-declare function decodeMsmFull(frame: Rtcm3Frame): ObsEpoch | null;
-/**
- * Convert an observation epoch timestamp to a Date.
- * For GPS/Galileo/BDS/QZSS: epochMs is ms of GPS week.
- * For GLONASS: epochMs is day-of-week × 86400000 + ms of day (UTC+3).
- *
- * Requires an approximate reference time to resolve the current week.
- */
-/**
- * NOTE: returns a GPS-scale Date — the clock-face convention used by
- * the RINEX parser and every consumer of epoch milliseconds in this
- * library (satellite propagation, SPP, RINEX epoch tags). It is NOT
- * UTC: converting to UTC here put every downstream satellite position
- * 18 s (the current leap-second count) off — a ±14 km per-satellite
- * range error that scattered live SPP solutions by kilometres.
- */
-declare function obsEpochToDate(sys: string, epochMs: number, refTime?: Date): Date;
-/**
- * Set a GLONASS FDMA frequency number (k = -7..+13) for a slot (1-24),
- * e.g. from a decoded 1020 ephemeris. MSM4/6 carry no DF419, so without
- * this the wavelength — and hence carrier phase — stays unknown.
- */
-declare function setGloFreqNumber(slot: number, k: number): void;
-/**
- * Reset the GLONASS frequency number cache (e.g., when switching mountpoints).
- */
-declare function resetGloFreqCache(): void;
-
 /**
  * Decode a legacy RTCM3 observation frame (1001–1004 GPS, 1009–1012 GLONASS)
  * into an {@link ObsEpoch}, or null for any other message type.
@@ -208,13 +125,14 @@ declare function resetGloFreqCache(): void;
 declare function decodeLegacyObs(frame: Rtcm3Frame): ObsEpoch | null;
 /**
  * Decode any RTCM3 observation frame into an {@link ObsEpoch}: MSM4–7 first
- * ({@link decodeMsmFull}), then the legacy fixed-layout obs (1001–1004 GPS,
+ * ({@link decodeMsmFull}; MSM1–3 too with `opts.compactMsm` — resolve those
+ * with `resolveMsmRanges` before use), then the legacy fixed-layout obs (1001–1004 GPS,
  * 1009–1012 GLONASS) via {@link decodeLegacyObs}. Returns null for any
  * non-observation message. This is the single obs entry point every consumer
  * (stream stats, live SPP/QC/RTK, station monitor) should call, so a new obs
  * source is added here once rather than at every call site.
  */
-declare function decodeObs(frame: Rtcm3Frame): ObsEpoch | null;
+declare function decodeObs(frame: Rtcm3Frame, opts?: ObsDecodeOptions): ObsEpoch | null;
 
 /** One announced message in the station's schedule (DF055/DF056/DF057). */
 interface SystemParamsMessage {
@@ -319,4 +237,4 @@ declare function rtcm3Constellation(msgType: number): string | null;
 declare function createStreamStats(): StreamStats;
 declare function updateStreamStats(stats: StreamStats, frames: Rtcm3Frame[], rawBytes: number): void;
 
-export { EphemerisInfo, type IgsSsrIonoLayer, type IgsSsrMessage, type MessageTypeStats, type ObsEpoch, type ObsSatObs, type ObsSignal, RTCM3_MESSAGE_NAMES, Rtcm3Frame, type SatCn0, type SignalCn0, type SsrCodeBias, type SsrKind, type SsrMessage, type SsrPhaseBias, type SsrSatCorrection, type StationMeta, type StreamStats, type SystemParams, type SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats };
+export { EphemerisInfo, type IgsSsrIonoLayer, type IgsSsrMessage, type MessageTypeStats, ObsDecodeOptions, ObsEpoch, RTCM3_MESSAGE_NAMES, Rtcm3Frame, type SatCn0, type SignalCn0, type SsrCodeBias, type SsrKind, type SsrMessage, type SsrPhaseBias, type SsrSatCorrection, type StationMeta, type StreamStats, type SystemParams, type SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, rtcm3Constellation, updateStationMeta, updateStreamStats };