gnss-js 2.9.0 → 2.10.0

package/dist/{chunk-6SWJTN6H.js → chunk-56HM5B6T.js}

@@ -1566,6 +1566,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();
@@ -1612,11 +1627,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(
@@ -1888,6 +1901,7 @@ export {
   createGzipLineSink,
   stationHeaderLines,
   writeRinexObsBlob,
+  modalIntervalS,
   writeModernObsBlob,
   writeRinex2ObsBlob,
   writeRinex4ObsBlob

package/dist/{chunk-NOGA5KD7.js → chunk-K3GLUWZI.js}

@@ -1621,6 +1621,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/positioning/sbas.ts
 var P2_11 = 2 ** -11;
 var P2_31 = 2 ** -31;
@@ -4011,6 +4071,7 @@ export {
   RtkFloatEngine,
   postProcessRtk,
   ZeroBaselineEngine,
+  resolveMsmRanges,
   sbasLongTermDeg,
   sbasIonoDeg,
   SbasProcessor,

package/dist/{chunk-7GNRNKBR.js → chunk-N6GQUHM4.js}

@@ -1138,12 +1138,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 {
@@ -1188,7 +1190,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 {
@@ -1214,7 +1218,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);
@@ -1290,8 +1305,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++;
       }
@@ -1303,7 +1320,8 @@ function decodeMsmFull(frame) {
       messageType: type,
       epochMs,
       system: sys,
-      observations
+      observations,
+      ...compact ? { rangeModulo1ms: true } : {}
     };
   } catch (err) {
     reportDecodeError("msm", err);
@@ -1457,8 +1475,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/index.cjs

@@ -264,6 +264,7 @@ __export(src_exports, {
   looksLikeRtcm2: () => looksLikeRtcm2,
   maskRadForAzimuth: () => maskRadForAzimuth,
   meridianArc: () => meridianArc,
+  modalIntervalS: () => modalIntervalS,
   moonEcef: () => moonEcef,
   navTimesFromEph: () => navTimesFromEph,
   nearestValidEphemeris: () => nearestValidEphemeris,
@@ -323,6 +324,7 @@ __export(src_exports, {
   readString: () => readString,
   reportDecodeError: () => reportDecodeError,
   resetGloFreqCache: () => resetGloFreqCache,
+  resolveMsmRanges: () => resolveMsmRanges,
   resolvePppAmbiguities: () => resolvePppAmbiguities,
   rhumbLine: () => rhumbLine,
   rtcm2Dgps: () => rtcm2Dgps,
@@ -4259,6 +4261,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();
@@ -4305,11 +4322,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(
@@ -5944,12 +5959,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 {
@@ -5994,7 +6011,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 {
@@ -6020,7 +6039,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);
@@ -6096,8 +6126,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++;
       }
@@ -6109,7 +6141,8 @@ function decodeMsmFull(frame) {
       messageType: type,
       epochMs,
       system: sys,
-      observations
+      observations,
+      ...compact ? { rangeModulo1ms: true } : {}
     };
   } catch (err) {
     reportDecodeError("msm", err);
@@ -6263,8 +6296,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
@@ -14898,6 +14931,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/positioning/sbas.ts
 var P2_11 = 2 ** -11;
 var P2_31 = 2 ** -31;
@@ -19778,6 +19871,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   looksLikeRtcm2,
   maskRadForAzimuth,
   meridianArc,
+  modalIntervalS,
   moonEcef,
   navTimesFromEph,
   nearestValidEphemeris,
@@ -19837,6 +19931,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolveMsmRanges,
   resolvePppAmbiguities,
   rhumbLine,
   rtcm2Dgps,

package/dist/index.d.cts

@@ -5,13 +5,14 @@ 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, meridianArc, rad2deg, rhumbLine, vincenty, vincentyDirect } from './coordinates.cjs';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.cjs';
-export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.cjs';
+export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseBiasSinex, parseCrxDataLine, 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 { C as ClkFile, I as IonexGrid, S as Sp3File, a as Sp3Sample, c as clkBias, p as parseClk, b as parseIonex, d as parseSp3, s as sp3Position } from './clk-DGVBGRzA.cjs';
 export { RTCM2_MESSAGE_NAMES, Rtcm2Correction, Rtcm2Decoder, Rtcm2Frame, Rtcm2ObsRecord, Rtcm2Station, looksLikeRtcm2, rtcm2Dgps, rtcm2Observation, rtcm2Station, rtcm2Text, rtcm2Time } from './rtcm2.cjs';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.cjs';
-export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, ObsEpoch, ObsSatObs, ObsSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
+export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, rtcm3Constellation, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
+export { O as ObsDecodeOptions, a as ObsEpoch, b as ObsSatObs, c as ObsSignal, d as decodeMsmFull, o as obsEpochToDate, r as resetGloFreqCache, s as setGloFreqNumber } from './msm-DWtoTvx2.cjs';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxPvtResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxPvt, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.cjs';
 export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavicResult, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, extractSbfNavicSubframe, feedGeoBlock, feedGeoL5Block, feedNavicBlock, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfNavic, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.cjs';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.cjs';
@@ -19,7 +20,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';
 export { AllPositionsData, AlmanacPosition, DopValues, EPHEMERIS_VALIDITY_MS, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, isEphemerisValid, keplerPosition, maskRadForAzimuth, navTimesFromEph, nearestValidEphemeris, rangeRate, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
-export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, 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 } from './positioning.cjs';
+export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, MsmResolveOptions, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, 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 } from './positioning.cjs';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';

package/dist/index.d.ts

@@ -5,13 +5,14 @@ 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, meridianArc, rad2deg, rhumbLine, vincenty, vincentyDirect } from './coordinates.js';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.js';
-export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.js';
+export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, modalIntervalS, padL, padR, parseBiasSinex, parseCrxDataLine, 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 { C as ClkFile, I as IonexGrid, S as Sp3File, a as Sp3Sample, c as clkBias, p as parseClk, b as parseIonex, d as parseSp3, s as sp3Position } from './clk-DGVBGRzA.js';
 export { RTCM2_MESSAGE_NAMES, Rtcm2Correction, Rtcm2Decoder, Rtcm2Frame, Rtcm2ObsRecord, Rtcm2Station, looksLikeRtcm2, rtcm2Dgps, rtcm2Observation, rtcm2Station, rtcm2Text, rtcm2Time } from './rtcm2.js';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.js';
-export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, ObsEpoch, ObsSatObs, ObsSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
+export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, rtcm3Constellation, updateStationMeta, updateStreamStats } from './rtcm3.js';
+export { O as ObsDecodeOptions, a as ObsEpoch, b as ObsSatObs, c as ObsSignal, d as decodeMsmFull, o as obsEpochToDate, r as resetGloFreqCache, s as setGloFreqNumber } from './msm-CcQJV3PC.js';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxPvtResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxPvt, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.js';
 export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavicResult, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, extractSbfNavicSubframe, feedGeoBlock, feedGeoL5Block, feedNavicBlock, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfNavic, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.js';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.js';
@@ -19,7 +20,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';
 export { AllPositionsData, AlmanacPosition, DopValues, EPHEMERIS_VALIDITY_MS, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, isEphemerisValid, keplerPosition, maskRadForAzimuth, navTimesFromEph, nearestValidEphemeris, rangeRate, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
-export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, 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 } from './positioning.js';
+export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, MsmResolveOptions, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, 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 } from './positioning.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';

package/dist/index.js

@@ -77,7 +77,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-7GNRNKBR.js";
+} from "./chunk-N6GQUHM4.js";
 import {
   parseUbxCnav,
   parseUbxIonoUtc,
@@ -201,6 +201,7 @@ import {
   nlWavelength,
   postProcessRtk,
   rcvAntennaRangeM,
+  resolveMsmRanges,
   resolvePppAmbiguities,
   satClockCorrection,
   sbasIonoDeg,
@@ -213,7 +214,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-NOGA5KD7.js";
+} from "./chunk-K3GLUWZI.js";
 import "./chunk-WQX354KD.js";
 import {
   EPHEMERIS_VALIDITY_MS,
@@ -337,6 +338,7 @@ import {
   fmtD,
   fmtF,
   hdrLine,
+  modalIntervalS,
   padL,
   padR,
   parseBiasSinex,
@@ -350,7 +352,7 @@ import {
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-6SWJTN6H.js";
+} from "./chunk-56HM5B6T.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -646,6 +648,7 @@ export {
   looksLikeRtcm2,
   maskRadForAzimuth,
   meridianArc,
+  modalIntervalS,
   moonEcef,
   navTimesFromEph,
   nearestValidEphemeris,
@@ -705,6 +708,7 @@ export {
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolveMsmRanges,
   resolvePppAmbiguities,
   rhumbLine,
   rtcm2Dgps,

package/dist/msm-CcQJV3PC.d.ts

@@ -0,0 +1,108 @@
+import { b as Rtcm3Frame } from './ephemeris-BAnBLLB9.js';
+
+/**
+ * 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 };