gnss-js 1.3.1 → 1.4.0

package/dist/index.cjs

@@ -67,6 +67,7 @@ __export(src_exports, {
   GLO_F2_BASE: () => GLO_F2_BASE,
   GLO_F2_STEP: () => GLO_F2_STEP,
   GLO_F3: () => GLO_F3,
+  IonoAccumulator: () => IonoAccumulator,
   MILLISECONDS_GPS_TAI: () => MILLISECONDS_GPS_TAI,
   MILLISECONDS_IN_DAY: () => MILLISECONDS_IN_DAY,
   MILLISECONDS_IN_HOUR: () => MILLISECONDS_IN_HOUR,
@@ -75,6 +76,7 @@ __export(src_exports, {
   MILLISECONDS_IN_WEEK: () => MILLISECONDS_IN_WEEK,
   MILLISECONDS_TT_TAI: () => MILLISECONDS_TT_TAI,
   MultipathAccumulator: () => MultipathAccumulator,
+  OMEGA_E: () => OMEGA_E,
   REFERENCE_FRAMES: () => REFERENCE_FRAMES,
   RINEX_CODES: () => RINEX_CODES,
   RTCM3_MESSAGE_NAMES: () => RTCM3_MESSAGE_NAMES,
@@ -308,6 +310,7 @@ var SYS_SHORT = {
   S: "SBS"
 };
 var C_LIGHT = 299792458;
+var OMEGA_E = 72921151467e-15;
 var FREQ = {
   G: { "1": 157542e4, "2": 12276e5, "5": 117645e4 },
   // R bands 4/6 are the CDMA L1OC/L2OC center frequencies
@@ -3814,6 +3817,7 @@ var RTCM3_MESSAGE_NAMES = {
   1031: "GLONASS Network RTK Residual",
   1032: "Physical Reference Station",
   1033: "Receiver+Antenna Descriptor",
+  1041: "NavIC Ephemeris",
   1042: "BeiDou Ephemeris",
   1043: "SBAS Ephemeris",
   1044: "QZSS Ephemeris",
@@ -3878,6 +3882,8 @@ var RTCM3_MESSAGE_NAMES = {
   1064: "GLONASS SSR Clock",
   1065: "GLONASS SSR Code Bias",
   1066: "GLONASS SSR Orbit+Clock",
+  // Biases
+  1230: "GLONASS Code-Phase Biases",
   // IGS SSR
   4076: "IGS SSR"
 };
@@ -3886,6 +3892,7 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1009 && msgType <= 1012) return "GLONASS";
   if (msgType === 1019) return "GPS";
   if (msgType === 1020) return "GLONASS";
+  if (msgType === 1041) return "NavIC";
   if (msgType === 1042) return "BeiDou";
   if (msgType === 1043) return "SBAS";
   if (msgType === 1044) return "QZSS";
@@ -3897,6 +3904,7 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1111 && msgType <= 1117) return "QZSS";
   if (msgType >= 1121 && msgType <= 1127) return "BeiDou";
   if (msgType >= 1131 && msgType <= 1137) return "NavIC";
+  if (msgType === 1230) return "GLONASS";
   return null;
 }
 function createStreamStats() {
@@ -3991,7 +3999,6 @@ var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
-var OMEGA_E = 72921151467e-15;
 var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
 var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
@@ -4726,8 +4733,7 @@ function tropoDelay(elevationRad) {
   return 2.47 / (sinEl + 0.0121);
 }
 function sagnac(pos, travelTimeS) {
-  const OMEGA_E2 = 72921151467e-15;
-  const theta = OMEGA_E2 * travelTimeS;
+  const theta = OMEGA_E * travelTimeS;
   const c = Math.cos(theta);
   const s = Math.sin(theta);
   return [pos.x * c + pos.y * s, -pos.x * s + pos.y * c, pos.z];
@@ -5038,6 +5044,14 @@ async function connectToMountpoint(proxyUrl, info, signal) {
 
 // src/analysis/multipath.ts
 var MIN_ARC_LENGTH = 10;
+var MP_JUMP_M = 1.25;
+var MP_JUMP_WINDOW = 5;
+var MP_EDIT_SIGMA = 3;
+function median(values) {
+  const s = [...values].sort((a, b) => a - b);
+  const m = s.length >> 1;
+  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
+}
 var MultipathAccumulator = class {
   state = /* @__PURE__ */ new Map();
   results = [];
@@ -5094,6 +5108,11 @@ var MultipathAccumulator = class {
     const gap = bandState.lastTime > 0 ? (time - bandState.lastTime) / 1e3 : 0;
     if (bandState.lastTime > 0 && gap > this.interval * ARC_GAP_FACTOR) {
       this.closeArc(prn, band, refBand, bandState);
+    } else if (bandState.arc.rawMp.length > 0) {
+      const recent = bandState.arc.rawMp.slice(-MP_JUMP_WINDOW);
+      if (Math.abs(mp - median(recent)) > MP_JUMP_M) {
+        this.closeArc(prn, band, refBand, bandState);
+      }
     }
     bandState.arc.times.push(time);
     bandState.arc.rawMp.push(mp);
@@ -5113,13 +5132,25 @@ var MultipathAccumulator = class {
   closeArc(prn, band, refBand, state) {
     const arc = state.arc;
     if (arc.times.length >= MIN_ARC_LENGTH) {
-      const mean = arc.rawMp.reduce((a, b) => a + b, 0) / arc.rawMp.length;
+      let keep = arc.times.map((_, i) => i);
+      let mean = 0;
+      for (let pass = 0; pass < 3; pass++) {
+        mean = keep.reduce((s, i) => s + arc.rawMp[i], 0) / keep.length;
+        const sigma = Math.sqrt(
+          keep.reduce((s, i) => s + (arc.rawMp[i] - mean) ** 2, 0) / keep.length
+        );
+        const kept = keep.filter(
+          (i) => Math.abs(arc.rawMp[i] - mean) <= MP_EDIT_SIGMA * sigma
+        );
+        if (kept.length === keep.length || kept.length < MIN_ARC_LENGTH) break;
+        keep = kept;
+      }
       const sys = prn[0];
       const bLabel = BAND_LABELS[sys]?.[band] ?? band;
       const rLabel = BAND_LABELS[sys]?.[refBand] ?? refBand;
       const label2 = `${prn} MP ${bLabel}-${rLabel}`;
-      const points = arc.times.map((t, i) => ({
-        time: t,
+      const points = keep.map((i) => ({
+        time: arc.times[i],
         mp: arc.rawMp[i] - mean
       }));
       const rms = Math.sqrt(
@@ -5173,7 +5204,7 @@ var MultipathAccumulator = class {
         refBand,
         rms,
         count,
-        satellites: seriesList.length
+        satellites: new Set(seriesList.map((s) => s.prn)).size
       });
     }
     const sysOrder = "GRECIJS";
@@ -5491,11 +5522,178 @@ var CompletenessAccumulator = class {
   }
 };
 
+// src/analysis/ionosphere.ts
+var TEC_FACTOR = 403e15;
+var GF_JUMP_M = 0.15;
+var MIN_ARC_LENGTH2 = 10;
+function median2(values) {
+  const s = [...values].sort((a, b) => a - b);
+  const m = s.length >> 1;
+  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
+}
+var IonoAccumulator = class {
+  state = /* @__PURE__ */ new Map();
+  closed = /* @__PURE__ */ new Map();
+  interval;
+  obsIndices;
+  gloChannels;
+  pairLabel = /* @__PURE__ */ new Map();
+  constructor(header) {
+    this.interval = header.interval ?? 30;
+    this.obsIndices = buildObsIndices(header);
+    this.gloChannels = buildGloChannelMap(header.glonassSlots);
+  }
+  /** Observation callback — wire this into parseRinexStream. */
+  onObservation = (time, prn, _codes, values) => {
+    const sys = prn[0];
+    const bandMap = this.obsIndices.get(sys) ?? this.obsIndices.get("_v2");
+    if (!bandMap) return;
+    const bandData = /* @__PURE__ */ new Map();
+    for (const [band, { C, L }] of bandMap) {
+      if (C === null) continue;
+      const cVal = values[C];
+      const lVal = values[L];
+      const freq = getFreq(this.gloChannels, prn, band);
+      if (cVal != null && cVal !== 0 && lVal != null && lVal !== 0 && freq) {
+        bandData.set(band, { C: cVal, L: lVal, f: freq });
+      }
+    }
+    if (bandData.size < 2) return;
+    const pairs = DUAL_FREQ_PAIRS[sys] ?? [];
+    for (const [bi, bj] of pairs) {
+      const di = bandData.get(bi);
+      const dj = bandData.get(bj);
+      if (!di || !dj) continue;
+      const \u03BBi = C_LIGHT / di.f;
+      const \u03BBj = C_LIGHT / dj.f;
+      const gamma = di.f * di.f / (dj.f * dj.f);
+      const l4 = di.L * \u03BBi - dj.L * \u03BBj;
+      const p4 = dj.C - di.C;
+      const pairKey = `${bi}-${bj}`;
+      if (!this.pairLabel.has(`${sys}:${pairKey}`)) {
+        const bLabel = BAND_LABELS[sys]?.[bi] ?? bi;
+        const rLabel = BAND_LABELS[sys]?.[bj] ?? bj;
+        this.pairLabel.set(`${sys}:${pairKey}`, `${bLabel}-${rLabel}`);
+      }
+      this.push(prn, pairKey, time, l4, p4, di.f, gamma);
+      break;
+    }
+  };
+  push(prn, pairKey, time, l4, p4, fi, gamma) {
+    if (!isFinite(l4) || !isFinite(p4)) return;
+    let satStates = this.state.get(prn);
+    if (!satStates) {
+      satStates = /* @__PURE__ */ new Map();
+      this.state.set(prn, satStates);
+    }
+    let ps = satStates.get(pairKey);
+    if (!ps) {
+      ps = {
+        arc: { times: [], l4: [], p4: [], fi, gamma },
+        lastTime: 0,
+        lastL4: 0
+      };
+      satStates.set(pairKey, ps);
+    }
+    const gap = ps.lastTime > 0 ? (time - ps.lastTime) / 1e3 : 0;
+    if (ps.lastTime > 0 && gap > this.interval * ARC_GAP_FACTOR) {
+      this.closeArc(prn, pairKey, ps);
+    } else if (ps.arc.times.length > 0 && Math.abs(l4 - ps.lastL4) > GF_JUMP_M) {
+      this.closeArc(prn, pairKey, ps);
+    }
+    if (ps.arc.times.length === 0) {
+      ps.arc.fi = fi;
+      ps.arc.gamma = gamma;
+    }
+    ps.arc.times.push(time);
+    ps.arc.l4.push(l4);
+    ps.arc.p4.push(p4);
+    ps.lastTime = time;
+    ps.lastL4 = l4;
+  }
+  /** External cycle-slip notification — close affected arcs. */
+  notifySlip(_time, prn, bands) {
+    const satStates = this.state.get(prn);
+    if (!satStates) return;
+    for (const [pairKey, ps] of satStates) {
+      const [bi, bj] = pairKey.split("-");
+      if (bands.has(bi) || bands.has(bj)) {
+        this.closeArc(prn, pairKey, ps);
+      }
+    }
+  }
+  closeArc(prn, pairKey, ps) {
+    const arc = ps.arc;
+    if (arc.times.length >= MIN_ARC_LENGTH2) {
+      const level = median2(arc.l4.map((v, k) => v - arc.p4[k]));
+      const toTecu = arc.fi * arc.fi / TEC_FACTOR / (arc.gamma - 1);
+      let satArcs = this.closed.get(prn);
+      if (!satArcs) {
+        satArcs = /* @__PURE__ */ new Map();
+        this.closed.set(prn, satArcs);
+      }
+      let points = satArcs.get(pairKey);
+      if (!points) {
+        points = [];
+        satArcs.set(pairKey, points);
+      }
+      for (let k = 0; k < arc.times.length; k++) {
+        points.push({
+          time: arc.times[k],
+          stec: (arc.l4[k] - level) * toTecu
+        });
+      }
+    }
+    ps.arc = { times: [], l4: [], p4: [], fi: arc.fi, gamma: arc.gamma };
+  }
+  /** Finalize: close remaining arcs, keep one pair per satellite. */
+  finalize() {
+    for (const [prn, satStates] of this.state) {
+      for (const [pairKey, ps] of satStates) {
+        this.closeArc(prn, pairKey, ps);
+      }
+    }
+    const series = [];
+    let sum = 0;
+    let count = 0;
+    let maxStec = 0;
+    for (const [prn, satArcs] of this.closed) {
+      let bestKey = null;
+      let bestLen = 0;
+      for (const [pairKey, points2] of satArcs) {
+        if (points2.length > bestLen) {
+          bestLen = points2.length;
+          bestKey = pairKey;
+        }
+      }
+      if (!bestKey) continue;
+      const points = satArcs.get(bestKey);
+      points.sort((a, b) => a.time - b.time);
+      const sys = prn[0];
+      for (const p of points) {
+        sum += p.stec;
+        count++;
+        if (p.stec > maxStec) maxStec = p.stec;
+      }
+      series.push({
+        prn,
+        system: sys,
+        label: this.pairLabel.get(`${sys}:${bestKey}`) ?? bestKey,
+        points
+      });
+    }
+    series.sort((a, b) => a.prn.localeCompare(b.prn));
+    return { series, maxStec, meanStec: count > 0 ? sum / count : 0 };
+  }
+};
+
 // src/analysis/quality-analysis.ts
 async function analyzeQuality(file, header, onProgress, signal) {
   const mpAccum = new MultipathAccumulator(header);
+  const ionoAccum = new IonoAccumulator(header);
   const csAccum = new CycleSlipAccumulator(header, (time, prn, bands) => {
     mpAccum.notifySlip(time, prn, bands);
+    ionoAccum.notifySlip(time, prn, bands);
   });
   const compAccum = new CompletenessAccumulator(header);
   await parseRinexStream(
@@ -5505,6 +5703,7 @@ async function analyzeQuality(file, header, onProgress, signal) {
     (time, prn, codes, values) => {
       csAccum.onObservation(time, prn, codes, values);
       mpAccum.onObservation(time, prn, codes, values);
+      ionoAccum.onObservation(time, prn, codes, values);
       compAccum.onObservation(time, prn, codes, values);
     },
     true
@@ -5513,7 +5712,8 @@ async function analyzeQuality(file, header, onProgress, signal) {
   return {
     cycleSlips: csAccum.finalize(),
     completeness: compAccum.finalize(),
-    multipath: mpAccum.finalize()
+    multipath: mpAccum.finalize(),
+    iono: ionoAccum.finalize()
   };
 }
 
@@ -6755,6 +6955,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   GLO_F2_BASE,
   GLO_F2_STEP,
   GLO_F3,
+  IonoAccumulator,
   MILLISECONDS_GPS_TAI,
   MILLISECONDS_IN_DAY,
   MILLISECONDS_IN_HOUR,
@@ -6763,6 +6964,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   MILLISECONDS_IN_WEEK,
   MILLISECONDS_TT_TAI,
   MultipathAccumulator,
+  OMEGA_E,
   REFERENCE_FRAMES,
   RINEX_CODES,
   RTCM3_MESSAGE_NAMES,

package/dist/index.d.cts

@@ -1,6 +1,6 @@
 export { D as DAYS_MJD2000_MJD, a as DAYS_MJD_JULIAN, H as HourCode, M as MILLISECONDS_GPS_TAI, b as MILLISECONDS_IN_DAY, c as MILLISECONDS_IN_HOUR, d as MILLISECONDS_IN_MINUTE, e as MILLISECONDS_IN_SECOND, f as MILLISECONDS_IN_WEEK, g as MILLISECONDS_TT_TAI, R as RINEX_CODES, S as SECONDS_IN_DAY, h as SECONDS_IN_HOUR, i as SECONDS_IN_MINUTE, j as SECONDS_IN_WEEK, k as START_BDS_TIME, l as START_GAL_TIME, m as START_GLO_LEAP, n as START_GPS_TIME, o as START_JULIAN_TAI, p as START_MJD_UNIX_SECONDS, q as START_NTP_TIME, r as START_TAI_TIME, s as START_UNIX_TIME, T as TimeDifference } from './time-DnI1VpE8.cjs';
 export { WGS84_ECCENTRICITY_SQUARED, WGS84_FLATTENING, WGS84_SEMI_MAJOR_AXIS, WGS84_SEMI_MINOR_AXIS } from './constants.cjs';
-export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OnSlipDetected, S as SYSTEM_NAMES, i as SYS_SHORT, j as buildGloChannelMap, k as buildObsIndices, l as formatUTCTime, m as getFreq, n as gloFreq } from './gnss-BT6ulR17.cjs';
+export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OMEGA_E, i as OnSlipDetected, S as SYSTEM_NAMES, j as SYS_SHORT, k as buildGloChannelMap, l as buildObsIndices, m as formatUTCTime, n as getFreq, o as gloFreq } from './gnss-BtXrG3zH.cjs';
 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';
@@ -13,7 +13,7 @@ export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPositi
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
 export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } 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, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.cjs';
+export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';
 export { NmeaFix, NmeaStats, NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum } from './nmea.cjs';
 export { BAND_PRESETS, BandPreset, CONSTELLATIONS, ConstellationRow, DELTA_GLO_L1, DELTA_GLO_L2, FREQ_BDS_B1A, FREQ_BDS_B1C, FREQ_BDS_B1I, FREQ_BDS_B2I, FREQ_BDS_B3A, FREQ_BDS_B3I, FREQ_GAL_E1, FREQ_GAL_E5, FREQ_GAL_E5a, FREQ_GAL_E5b, FREQ_GAL_E6, FREQ_GLO_L1, FREQ_GLO_L1OC, FREQ_GLO_L2, FREQ_GLO_L2OC, FREQ_GLO_L3OC, FREQ_GPS_L1, FREQ_GPS_L2, FREQ_GPS_L5, FREQ_NAVIC_L5, FREQ_QZS_L1, FREQ_QZS_L2, FREQ_QZS_L5, FREQ_QZS_L6, Modulation, SignalDef, computePsdDb, phiAltBOC, phiBOCc, phiBOCs, phiBPSK } from './signals.cjs';

package/dist/index.d.ts

@@ -1,6 +1,6 @@
 export { D as DAYS_MJD2000_MJD, a as DAYS_MJD_JULIAN, H as HourCode, M as MILLISECONDS_GPS_TAI, b as MILLISECONDS_IN_DAY, c as MILLISECONDS_IN_HOUR, d as MILLISECONDS_IN_MINUTE, e as MILLISECONDS_IN_SECOND, f as MILLISECONDS_IN_WEEK, g as MILLISECONDS_TT_TAI, R as RINEX_CODES, S as SECONDS_IN_DAY, h as SECONDS_IN_HOUR, i as SECONDS_IN_MINUTE, j as SECONDS_IN_WEEK, k as START_BDS_TIME, l as START_GAL_TIME, m as START_GLO_LEAP, n as START_GPS_TIME, o as START_JULIAN_TAI, p as START_MJD_UNIX_SECONDS, q as START_NTP_TIME, r as START_TAI_TIME, s as START_UNIX_TIME, T as TimeDifference } from './time-DnI1VpE8.js';
 export { WGS84_ECCENTRICITY_SQUARED, WGS84_FLATTENING, WGS84_SEMI_MAJOR_AXIS, WGS84_SEMI_MINOR_AXIS } from './constants.js';
-export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OnSlipDetected, S as SYSTEM_NAMES, i as SYS_SHORT, j as buildGloChannelMap, k as buildObsIndices, l as formatUTCTime, m as getFreq, n as gloFreq } from './gnss-C-tgoYNa.js';
+export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OMEGA_E, i as OnSlipDetected, S as SYSTEM_NAMES, j as SYS_SHORT, k as buildGloChannelMap, l as buildObsIndices, m as formatUTCTime, n as getFreq, o as gloFreq } from './gnss-DhnEr0Dy.js';
 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';
@@ -13,7 +13,7 @@ export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPositi
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
 export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } 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, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.js';
+export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoPoint, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';
 export { NmeaFix, NmeaStats, NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum } from './nmea.js';
 export { BAND_PRESETS, BandPreset, CONSTELLATIONS, ConstellationRow, DELTA_GLO_L1, DELTA_GLO_L2, FREQ_BDS_B1A, FREQ_BDS_B1C, FREQ_BDS_B1I, FREQ_BDS_B2I, FREQ_BDS_B3A, FREQ_BDS_B3I, FREQ_GAL_E1, FREQ_GAL_E5, FREQ_GAL_E5a, FREQ_GAL_E5b, FREQ_GAL_E6, FREQ_GLO_L1, FREQ_GLO_L1OC, FREQ_GLO_L2, FREQ_GLO_L2OC, FREQ_GLO_L3OC, FREQ_GPS_L1, FREQ_GPS_L2, FREQ_GPS_L5, FREQ_NAVIC_L5, FREQ_QZS_L1, FREQ_QZS_L2, FREQ_QZS_L5, FREQ_QZS_L6, Modulation, SignalDef, computePsdDb, phiAltBOC, phiBOCc, phiBOCs, phiBPSK } from './signals.js';

package/dist/index.js

@@ -6,9 +6,10 @@ import {
 import {
   CompletenessAccumulator,
   CycleSlipAccumulator,
+  IonoAccumulator,
   MultipathAccumulator,
   analyzeQuality
-} from "./chunk-G3N4S3DM.js";
+} from "./chunk-3ZAFWOT6.js";
 import {
   frequencyLabel,
   parseAntex
@@ -53,7 +54,7 @@ import {
   phiBOCc,
   phiBOCs,
   phiBPSK
-} from "./chunk-5SPJH4MG.js";
+} from "./chunk-2K3FCJBX.js";
 import {
   Scale,
   getDateFromDayOfWeek,
@@ -97,7 +98,7 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-USHP4IND.js";
+} from "./chunk-PX7TPPTQ.js";
 import {
   computeAllPositions,
   computeDop,
@@ -110,7 +111,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-CC2R4JGS.js";
+} from "./chunk-PQYBTE42.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -178,7 +179,7 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-OZCYOM5D.js";
+} from "./chunk-3U5AX7PY.js";
 import {
   BitReader,
   RTCM3_MESSAGE_NAMES,
@@ -196,7 +197,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-7EEWQ5DU.js";
+} from "./chunk-65CQINSB.js";
 import {
   DAYS_MJD2000_MJD,
   DAYS_MJD_JULIAN,
@@ -236,6 +237,7 @@ import {
   GLO_F2_BASE,
   GLO_F2_STEP,
   GLO_F3,
+  OMEGA_E,
   SYSTEM_NAMES,
   SYS_SHORT,
   buildGloChannelMap,
@@ -243,7 +245,7 @@ import {
   formatUTCTime,
   getFreq,
   gloFreq
-} from "./chunk-W5WKEV7U.js";
+} from "./chunk-FIEWO4J4.js";
 export {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -292,6 +294,7 @@ export {
   GLO_F2_BASE,
   GLO_F2_STEP,
   GLO_F3,
+  IonoAccumulator,
   MILLISECONDS_GPS_TAI,
   MILLISECONDS_IN_DAY,
   MILLISECONDS_IN_HOUR,
@@ -300,6 +303,7 @@ export {
   MILLISECONDS_IN_WEEK,
   MILLISECONDS_TT_TAI,
   MultipathAccumulator,
+  OMEGA_E,
   REFERENCE_FRAMES,
   RINEX_CODES,
   RTCM3_MESSAGE_NAMES,

package/dist/orbit.cjs

@@ -80,6 +80,9 @@ var START_GPS_TIME = /* @__PURE__ */ new Date("1980-01-06T00:00:00Z");
 var START_BDS_TIME = /* @__PURE__ */ new Date("2006-01-01T00:00:14Z");
 var START_NTP_TIME = /* @__PURE__ */ new Date("1900-01-01T00:00:00Z");
 
+// src/constants/gnss.ts
+var OMEGA_E = 72921151467e-15;
+
 // src/time/gnss.ts
 function getTaiDate(date) {
   return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
@@ -167,7 +170,6 @@ var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
-var OMEGA_E = 72921151467e-15;
 var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
 var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;

package/dist/orbit.js

@@ -10,7 +10,7 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-CC2R4JGS.js";
+} from "./chunk-PQYBTE42.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
@@ -18,6 +18,7 @@ import {
 } from "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-LEEU5OIO.js";
+import "./chunk-FIEWO4J4.js";
 export {
   computeAllPositions,
   computeDop,

package/dist/positioning.cjs

@@ -63,6 +63,10 @@ var START_GPS_TIME = /* @__PURE__ */ new Date("1980-01-06T00:00:00Z");
 var START_BDS_TIME = /* @__PURE__ */ new Date("2006-01-01T00:00:14Z");
 var START_NTP_TIME = /* @__PURE__ */ new Date("1900-01-01T00:00:00Z");
 
+// src/constants/gnss.ts
+var C_LIGHT = 299792458;
+var OMEGA_E = 72921151467e-15;
+
 // src/time/gnss.ts
 function getTaiDate(date) {
   return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
@@ -150,7 +154,6 @@ var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
-var OMEGA_E = 72921151467e-15;
 var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
 var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
@@ -366,9 +369,6 @@ function computeSatPosition(eph, timeMs) {
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
 
-// src/constants/gnss.ts
-var C_LIGHT = 299792458;
-
 // src/positioning/index.ts
 var GM_GPS2 = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -400,8 +400,7 @@ function tropoDelay(elevationRad) {
   return 2.47 / (sinEl + 0.0121);
 }
 function sagnac(pos, travelTimeS) {
-  const OMEGA_E2 = 72921151467e-15;
-  const theta = OMEGA_E2 * travelTimeS;
+  const theta = OMEGA_E * travelTimeS;
   const c = Math.cos(theta);
   const s = Math.sin(theta);
   return [pos.x * c + pos.y * s, -pos.x * s + pos.y * c, pos.z];

package/dist/positioning.js

@@ -2,13 +2,13 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-USHP4IND.js";
-import "./chunk-CC2R4JGS.js";
+} from "./chunk-PX7TPPTQ.js";
+import "./chunk-PQYBTE42.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-LEEU5OIO.js";
-import "./chunk-W5WKEV7U.js";
+import "./chunk-FIEWO4J4.js";
 export {
   ionoFree,
   satClockCorrection,

package/dist/rinex.js

@@ -16,8 +16,8 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-OZCYOM5D.js";
-import "./chunk-W5WKEV7U.js";
+} from "./chunk-3U5AX7PY.js";
+import "./chunk-FIEWO4J4.js";
 export {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,

package/dist/rtcm3.cjs

@@ -1319,6 +1319,7 @@ var RTCM3_MESSAGE_NAMES = {
   1031: "GLONASS Network RTK Residual",
   1032: "Physical Reference Station",
   1033: "Receiver+Antenna Descriptor",
+  1041: "NavIC Ephemeris",
   1042: "BeiDou Ephemeris",
   1043: "SBAS Ephemeris",
   1044: "QZSS Ephemeris",
@@ -1383,6 +1384,8 @@ var RTCM3_MESSAGE_NAMES = {
   1064: "GLONASS SSR Clock",
   1065: "GLONASS SSR Code Bias",
   1066: "GLONASS SSR Orbit+Clock",
+  // Biases
+  1230: "GLONASS Code-Phase Biases",
   // IGS SSR
   4076: "IGS SSR"
 };
@@ -1391,6 +1394,7 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1009 && msgType <= 1012) return "GLONASS";
   if (msgType === 1019) return "GPS";
   if (msgType === 1020) return "GLONASS";
+  if (msgType === 1041) return "NavIC";
   if (msgType === 1042) return "BeiDou";
   if (msgType === 1043) return "SBAS";
   if (msgType === 1044) return "QZSS";
@@ -1402,6 +1406,7 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1111 && msgType <= 1117) return "QZSS";
   if (msgType >= 1121 && msgType <= 1127) return "BeiDou";
   if (msgType >= 1131 && msgType <= 1137) return "NavIC";
+  if (msgType === 1230) return "GLONASS";
   return null;
 }
 function createStreamStats() {

package/dist/rtcm3.js

@@ -15,9 +15,9 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-7EEWQ5DU.js";
+} from "./chunk-65CQINSB.js";
 import "./chunk-LEEU5OIO.js";
-import "./chunk-W5WKEV7U.js";
+import "./chunk-FIEWO4J4.js";
 export {
   BitReader,
   RTCM3_MESSAGE_NAMES,

package/dist/signals.js

@@ -32,8 +32,8 @@ import {
   phiBOCc,
   phiBOCs,
   phiBPSK
-} from "./chunk-5SPJH4MG.js";
-import "./chunk-W5WKEV7U.js";
+} from "./chunk-2K3FCJBX.js";
+import "./chunk-FIEWO4J4.js";
 export {
   BAND_PRESETS,
   CONSTELLATIONS,

package/package.json

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