gnss-js 2.8.1 → 2.9.0

package/dist/{chunk-T5M3LVBJ.js → chunk-7T5DLVZL.js}

@@ -37,6 +37,42 @@ function invertMatrix(A) {
   return M.map((row) => row.slice(n));
 }
 
+// src/orbit/validity.ts
+var H = 36e5;
+var MIN = 6e4;
+var EPHEMERIS_VALIDITY_MS = {
+  G: 2 * H,
+  J: 2 * H,
+  E: 2.5 * H,
+  C: 90 * MIN,
+  R: 15 * MIN,
+  I: 2 * H,
+  S: 6 * MIN
+};
+function tocDistanceMs(eph, tGpsMs) {
+  const t = eph.system === "R" || eph.system === "S" ? tGpsMs - getGpsLeap(new Date(tGpsMs)) * 1e3 : tGpsMs;
+  return Math.abs(eph.tocDate.getTime() - t);
+}
+function isEphemerisValid(eph, tGpsMs) {
+  const w = EPHEMERIS_VALIDITY_MS[eph.system];
+  return w !== void 0 && tocDistanceMs(eph, tGpsMs) <= w;
+}
+function nearestValidEphemeris(ephs, tGpsMs) {
+  if (!ephs) return null;
+  let best = null;
+  let bestDt = Infinity;
+  for (const eph of ephs) {
+    const dt = tocDistanceMs(eph, tGpsMs);
+    const w = EPHEMERIS_VALIDITY_MS[eph.system];
+    if (w === void 0 || dt > w) continue;
+    if (dt < bestDt) {
+      bestDt = dt;
+      best = eph;
+    }
+  }
+  return best;
+}
+
 // src/orbit/almanac.ts
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
@@ -409,24 +445,30 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  const perClock = /* @__PURE__ */ new Map();
+  for (const s of satAzEls) {
+    const c = dopClock(s.prn);
+    perClock.set(c, (perClock.get(c) ?? 0) + 1);
+  }
+  const sats = perClock.size > 1 ? satAzEls.filter((s) => perClock.get(dopClock(s.prn)) > 1) : satAzEls;
+  const n = sats.length;
+  const clocks = [...new Set(sats.map((s) => dopClock(s.prn)))];
   const m = 3 + clocks.length;
   if (n < Math.max(4, m)) return null;
-  const H = [];
-  for (const { az, el, prn } of satAzEls) {
+  const H2 = [];
+  for (const { az, el, prn } of sats) {
     const cosEl = Math.cos(el);
     const row = [cosEl * Math.sin(az), cosEl * Math.cos(az), Math.sin(el)];
     const clock = dopClock(prn);
     for (const c of clocks) row.push(c === clock ? 1 : 0);
-    H.push(row);
+    H2.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));
   for (let i = 0; i < m; i++) {
     for (let j = 0; j < m; j++) {
       let sum = 0;
       for (let k = 0; k < n; k++) {
-        sum += H[k][i] * H[k][j];
+        sum += H2[k][i] * H2[k][j];
       }
       HtH[i][j] = sum;
     }
@@ -527,18 +569,7 @@ function computeAllPositions(ephemerides, times, rxPos) {
   return { prns, times, positions };
 }
 function selectBest(ephs, timeMs) {
-  let best = null;
-  let bestDt = Infinity;
-  for (const eph of ephs) {
-    const dt = Math.abs(timeMs - eph.tocDate.getTime());
-    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
-    if (dt > maxAge) continue;
-    if (dt < bestDt) {
-      bestDt = dt;
-      best = eph;
-    }
-  }
-  return best;
+  return nearestValidEphemeris(ephs, timeMs);
 }
 var GPS_EPOCH_MS2 = START_GPS_TIME.getTime();
 var MS_PER_WEEK = 7 * 864e5;
@@ -780,6 +811,9 @@ function computeVisibilityFromPositions(all, elevationMaskDeg = 10) {
 
 export {
   invertMatrix,
+  EPHEMERIS_VALIDITY_MS,
+  isEphemerisValid,
+  nearestValidEphemeris,
   almanacEpochMs,
   almanacSatPosition,
   glonassAlmanacEpochMs,

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

@@ -7,7 +7,7 @@ import {
   ecefToAzEl,
   invertMatrix,
   selectEphemeris
-} from "./chunk-T5M3LVBJ.js";
+} from "./chunk-7T5DLVZL.js";
 import {
   ecefToGeodetic,
   getAer,

package/dist/index.cjs

@@ -36,6 +36,7 @@ __export(src_exports, {
   DELTA_GLO_L2: () => DELTA_GLO_L2,
   DUAL_FREQ_PAIRS: () => DUAL_FREQ_PAIRS,
   EMPTY_WARNINGS: () => EMPTY_WARNINGS,
+  EPHEMERIS_VALIDITY_MS: () => EPHEMERIS_VALIDITY_MS,
   FREQ: () => FREQ,
   FREQ_BDS_B1A: () => FREQ_BDS_B1A,
   FREQ_BDS_B1C: () => FREQ_BDS_B1C,
@@ -252,6 +253,7 @@ __export(src_exports, {
   hdrLine: () => hdrLine,
   horizonDistance: () => horizonDistance,
   ionoFree: () => ionoFree,
+  isEphemerisValid: () => isEphemerisValid,
   isSbasL1Preamble: () => isSbasL1Preamble,
   isSbasL5Correction: () => isSbasL5Correction,
   isSsrMessage: () => isSsrMessage,
@@ -264,6 +266,7 @@ __export(src_exports, {
   meridianArc: () => meridianArc,
   moonEcef: () => moonEcef,
   navTimesFromEph: () => navTimesFromEph,
+  nearestValidEphemeris: () => nearestValidEphemeris,
   niellMapping: () => niellMapping,
   nlWavelength: () => nlWavelength,
   nmeaCoordToDecimal: () => nmeaCoordToDecimal,
@@ -12312,6 +12315,42 @@ function invertMatrix(A) {
   return M.map((row) => row.slice(n));
 }
 
+// src/orbit/validity.ts
+var H = 36e5;
+var MIN = 6e4;
+var EPHEMERIS_VALIDITY_MS = {
+  G: 2 * H,
+  J: 2 * H,
+  E: 2.5 * H,
+  C: 90 * MIN,
+  R: 15 * MIN,
+  I: 2 * H,
+  S: 6 * MIN
+};
+function tocDistanceMs(eph, tGpsMs) {
+  const t = eph.system === "R" || eph.system === "S" ? tGpsMs - getGpsLeap(new Date(tGpsMs)) * 1e3 : tGpsMs;
+  return Math.abs(eph.tocDate.getTime() - t);
+}
+function isEphemerisValid(eph, tGpsMs) {
+  const w = EPHEMERIS_VALIDITY_MS[eph.system];
+  return w !== void 0 && tocDistanceMs(eph, tGpsMs) <= w;
+}
+function nearestValidEphemeris(ephs, tGpsMs) {
+  if (!ephs) return null;
+  let best = null;
+  let bestDt = Infinity;
+  for (const eph of ephs) {
+    const dt = tocDistanceMs(eph, tGpsMs);
+    const w = EPHEMERIS_VALIDITY_MS[eph.system];
+    if (w === void 0 || dt > w) continue;
+    if (dt < bestDt) {
+      bestDt = dt;
+      best = eph;
+    }
+  }
+  return best;
+}
+
 // src/orbit/almanac.ts
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
@@ -12684,12 +12723,18 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  const perClock = /* @__PURE__ */ new Map();
+  for (const s of satAzEls) {
+    const c = dopClock(s.prn);
+    perClock.set(c, (perClock.get(c) ?? 0) + 1);
+  }
+  const sats = perClock.size > 1 ? satAzEls.filter((s) => perClock.get(dopClock(s.prn)) > 1) : satAzEls;
+  const n = sats.length;
+  const clocks = [...new Set(sats.map((s) => dopClock(s.prn)))];
   const m = 3 + clocks.length;
   if (n < Math.max(4, m)) return null;
   const H3 = [];
-  for (const { az, el, prn } of satAzEls) {
+  for (const { az, el, prn } of sats) {
     const cosEl = Math.cos(el);
     const row = [cosEl * Math.sin(az), cosEl * Math.cos(az), Math.sin(el)];
     const clock = dopClock(prn);
@@ -12802,18 +12847,7 @@ function computeAllPositions(ephemerides, times, rxPos) {
   return { prns, times, positions };
 }
 function selectBest(ephs, timeMs) {
-  let best = null;
-  let bestDt = Infinity;
-  for (const eph of ephs) {
-    const dt = Math.abs(timeMs - eph.tocDate.getTime());
-    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
-    if (dt > maxAge) continue;
-    if (dt < bestDt) {
-      bestDt = dt;
-      best = eph;
-    }
-  }
-  return best;
+  return nearestValidEphemeris(ephs, timeMs);
 }
 var GPS_EPOCH_MS24 = START_GPS_TIME.getTime();
 var MS_PER_WEEK13 = 7 * 864e5;
@@ -13054,14 +13088,14 @@ function computeVisibilityFromPositions(all, elevationMaskDeg = 10) {
 }
 
 // src/frames/index.ts
-var H = (t, d, r, tDot, dDot, rDot, epoch) => ({ t, d, r, tDot, dDot, rDot, epoch });
+var H2 = (t, d, r, tDot, dDot, rDot, epoch) => ({ t, d, r, tDot, dDot, rDot, epoch });
 var Z = [0, 0, 0];
 var FROM_ITRF2020 = {
-  ITRF2014: H([-1.4, -0.9, 1.4], -0.42, Z, [0, -0.1, 0.2], 0, Z, 2015),
-  ITRF2008: H([0.2, 1, 3.3], -0.29, Z, [0, -0.1, 0.1], 0.03, Z, 2015),
-  ITRF2005: H([2.7, 0.1, -1.4], 0.65, Z, [0.3, -0.1, 0.1], 0.03, Z, 2015),
-  ITRF2000: H([-0.2, 0.8, -34.2], 2.25, Z, [0.1, 0, -1.7], 0.11, Z, 2015),
-  ITRF97: H(
+  ITRF2014: H2([-1.4, -0.9, 1.4], -0.42, Z, [0, -0.1, 0.2], 0, Z, 2015),
+  ITRF2008: H2([0.2, 1, 3.3], -0.29, Z, [0, -0.1, 0.1], 0.03, Z, 2015),
+  ITRF2005: H2([2.7, 0.1, -1.4], 0.65, Z, [0.3, -0.1, 0.1], 0.03, Z, 2015),
+  ITRF2000: H2([-0.2, 0.8, -34.2], 2.25, Z, [0.1, 0, -1.7], 0.11, Z, 2015),
+  ITRF97: H2(
     [6.5, -3.9, -77.9],
     3.98,
     [0, 0, 0.36],
@@ -13070,7 +13104,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF96: H(
+  ITRF96: H2(
     [6.5, -3.9, -77.9],
     3.98,
     [0, 0, 0.36],
@@ -13079,7 +13113,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF94: H(
+  ITRF94: H2(
     [6.5, -3.9, -77.9],
     3.98,
     [0, 0, 0.36],
@@ -13088,7 +13122,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF93: H(
+  ITRF93: H2(
     [-65.8, 1.9, -71.3],
     4.47,
     [-3.36, -4.33, 0.75],
@@ -13097,7 +13131,7 @@ var FROM_ITRF2020 = {
     [-0.11, -0.19, 0.07],
     2015
   ),
-  ITRF92: H(
+  ITRF92: H2(
     [14.5, -1.9, -85.9],
     3.27,
     [0, 0, 0.36],
@@ -13106,7 +13140,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF91: H(
+  ITRF91: H2(
     [26.5, 12.1, -91.9],
     4.67,
     [0, 0, 0.36],
@@ -13115,7 +13149,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF90: H(
+  ITRF90: H2(
     [24.5, 8.1, -107.9],
     4.97,
     [0, 0, 0.36],
@@ -13124,7 +13158,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF89: H(
+  ITRF89: H2(
     [29.5, 32.1, -145.9],
     8.37,
     [0, 0, 0.36],
@@ -13133,7 +13167,7 @@ var FROM_ITRF2020 = {
     [0, 0, 0.02],
     2015
   ),
-  ITRF88: H(
+  ITRF88: H2(
     [24.5, -3.9, -169.9],
     11.47,
     [0.1, 0, 0.36],
@@ -13143,10 +13177,10 @@ var FROM_ITRF2020 = {
     2015
   ),
   /** ITRF2020 → ETRF2020, EPSG:10572 (epoch-1989 variant) */
-  ETRF2020: H(Z, 0, Z, Z, 0, [0.086, 0.519, -0.753], 1989)
+  ETRF2020: H2(Z, 0, Z, Z, 0, [0.086, 0.519, -0.753], 1989)
 };
-var ITRF2014_TO_ETRF2014 = H(Z, 0, Z, Z, 0, [0.085, 0.531, -0.77], 1989);
-var ITRF2014_TO_ETRF2000 = H(
+var ITRF2014_TO_ETRF2014 = H2(Z, 0, Z, Z, 0, [0.085, 0.531, -0.77], 1989);
+var ITRF2014_TO_ETRF2000 = H2(
   [54.7, 52.2, -74.1],
   2.12,
   [1.701, 10.29, -16.632],
@@ -13155,7 +13189,7 @@ var ITRF2014_TO_ETRF2000 = H(
   [0.081, 0.49, -0.792],
   2010
 );
-var ITRF2014_TO_NAD83_2011 = H(
+var ITRF2014_TO_NAD83_2011 = H2(
   [1005.3, -1909.21, -541.57],
   0.36891,
   [-26.78138, 0.42027, -10.93206],
@@ -15831,7 +15865,7 @@ function moonEcef(timeMs) {
   ];
   return eciToEcef(eci, gmst(timeMs));
 }
-var H2 = 0.6078;
+var H22 = 0.6078;
 var L2 = 0.0847;
 var RE = 6378137;
 var GM_SUN_RATIO = 332946.0482;
@@ -15856,7 +15890,7 @@ function solidEarthTide(rcv, sun, moon) {
     ];
     const dot2 = rHat[0] * jHat[0] + rHat[1] * jHat[1] + rHat[2] * jHat[2];
     const factor = ratio * RE * RE * RE * RE / (rj * rj * rj);
-    const radialScale = H2 * (1.5 * dot2 * dot2 - 0.5);
+    const radialScale = H22 * (1.5 * dot2 * dot2 - 0.5);
     const tanScale = 3 * L2 * dot2;
     for (let k = 0; k < 3; k++) {
       disp[k] += factor * (radialScale * rHat[k] + tanScale * (jHat[k] - dot2 * rHat[k]));
@@ -19516,6 +19550,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   DELTA_GLO_L2,
   DUAL_FREQ_PAIRS,
   EMPTY_WARNINGS,
+  EPHEMERIS_VALIDITY_MS,
   FREQ,
   FREQ_BDS_B1A,
   FREQ_BDS_B1C,
@@ -19732,6 +19767,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   hdrLine,
   horizonDistance,
   ionoFree,
+  isEphemerisValid,
   isSbasL1Preamble,
   isSbasL5Correction,
   isSsrMessage,
@@ -19744,6 +19780,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   meridianArc,
   moonEcef,
   navTimesFromEph,
+  nearestValidEphemeris,
   niellMapping,
   nlWavelength,
   nmeaCoordToDecimal,

package/dist/index.d.cts

@@ -17,7 +17,7 @@ export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAM
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.cjs';
 export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.cjs';
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';
-export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.cjs';
+export { 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 { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';

package/dist/index.d.ts

@@ -17,7 +17,7 @@ export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAM
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.js';
 export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.js';
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';
-export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.js';
+export { 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 { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';

package/dist/index.js

@@ -213,9 +213,10 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-F6KOFOXN.js";
+} from "./chunk-NOGA5KD7.js";
 import "./chunk-WQX354KD.js";
 import {
+  EPHEMERIS_VALIDITY_MS,
   almanacEpochMs,
   almanacSatPosition,
   computeAllPositions,
@@ -230,12 +231,14 @@ import {
   glonassAlmanacEpochMs,
   glonassAlmanacPosition,
   glonassPosition,
+  isEphemerisValid,
   keplerPosition,
   maskRadForAzimuth,
   navTimesFromEph,
+  nearestValidEphemeris,
   rangeRate,
   selectEphemeris
-} from "./chunk-T5M3LVBJ.js";
+} from "./chunk-7T5DLVZL.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -415,6 +418,7 @@ export {
   DELTA_GLO_L2,
   DUAL_FREQ_PAIRS,
   EMPTY_WARNINGS,
+  EPHEMERIS_VALIDITY_MS,
   FREQ,
   FREQ_BDS_B1A,
   FREQ_BDS_B1C,
@@ -631,6 +635,7 @@ export {
   hdrLine,
   horizonDistance,
   ionoFree,
+  isEphemerisValid,
   isSbasL1Preamble,
   isSbasL5Correction,
   isSsrMessage,
@@ -643,6 +648,7 @@ export {
   meridianArc,
   moonEcef,
   navTimesFromEph,
+  nearestValidEphemeris,
   niellMapping,
   nlWavelength,
   nmeaCoordToDecimal,

package/dist/orbit.cjs

@@ -20,6 +20,7 @@ var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: tru
 // src/orbit/index.ts
 var orbit_exports = {};
 __export(orbit_exports, {
+  EPHEMERIS_VALIDITY_MS: () => EPHEMERIS_VALIDITY_MS,
   almanacEpochMs: () => almanacEpochMs,
   almanacSatPosition: () => almanacSatPosition,
   computeAllPositions: () => computeAllPositions,
@@ -36,9 +37,11 @@ __export(orbit_exports, {
   glonassAlmanacEpochMs: () => glonassAlmanacEpochMs,
   glonassAlmanacPosition: () => glonassAlmanacPosition,
   glonassPosition: () => glonassPosition,
+  isEphemerisValid: () => isEphemerisValid,
   keplerPosition: () => keplerPosition,
   maskRadForAzimuth: () => maskRadForAzimuth,
   navTimesFromEph: () => navTimesFromEph,
+  nearestValidEphemeris: () => nearestValidEphemeris,
   rangeRate: () => rangeRate,
   selectEphemeris: () => selectEphemeris
 });
@@ -198,6 +201,42 @@ function invertMatrix(A) {
   return M.map((row) => row.slice(n));
 }
 
+// src/orbit/validity.ts
+var H = 36e5;
+var MIN = 6e4;
+var EPHEMERIS_VALIDITY_MS = {
+  G: 2 * H,
+  J: 2 * H,
+  E: 2.5 * H,
+  C: 90 * MIN,
+  R: 15 * MIN,
+  I: 2 * H,
+  S: 6 * MIN
+};
+function tocDistanceMs(eph, tGpsMs) {
+  const t = eph.system === "R" || eph.system === "S" ? tGpsMs - getGpsLeap(new Date(tGpsMs)) * 1e3 : tGpsMs;
+  return Math.abs(eph.tocDate.getTime() - t);
+}
+function isEphemerisValid(eph, tGpsMs) {
+  const w = EPHEMERIS_VALIDITY_MS[eph.system];
+  return w !== void 0 && tocDistanceMs(eph, tGpsMs) <= w;
+}
+function nearestValidEphemeris(ephs, tGpsMs) {
+  if (!ephs) return null;
+  let best = null;
+  let bestDt = Infinity;
+  for (const eph of ephs) {
+    const dt = tocDistanceMs(eph, tGpsMs);
+    const w = EPHEMERIS_VALIDITY_MS[eph.system];
+    if (w === void 0 || dt > w) continue;
+    if (dt < bestDt) {
+      bestDt = dt;
+      best = eph;
+    }
+  }
+  return best;
+}
+
 // src/orbit/almanac.ts
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
@@ -570,24 +609,30 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  const perClock = /* @__PURE__ */ new Map();
+  for (const s of satAzEls) {
+    const c = dopClock(s.prn);
+    perClock.set(c, (perClock.get(c) ?? 0) + 1);
+  }
+  const sats = perClock.size > 1 ? satAzEls.filter((s) => perClock.get(dopClock(s.prn)) > 1) : satAzEls;
+  const n = sats.length;
+  const clocks = [...new Set(sats.map((s) => dopClock(s.prn)))];
   const m = 3 + clocks.length;
   if (n < Math.max(4, m)) return null;
-  const H = [];
-  for (const { az, el, prn } of satAzEls) {
+  const H2 = [];
+  for (const { az, el, prn } of sats) {
     const cosEl = Math.cos(el);
     const row = [cosEl * Math.sin(az), cosEl * Math.cos(az), Math.sin(el)];
     const clock = dopClock(prn);
     for (const c of clocks) row.push(c === clock ? 1 : 0);
-    H.push(row);
+    H2.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));
   for (let i = 0; i < m; i++) {
     for (let j = 0; j < m; j++) {
       let sum = 0;
       for (let k = 0; k < n; k++) {
-        sum += H[k][i] * H[k][j];
+        sum += H2[k][i] * H2[k][j];
       }
       HtH[i][j] = sum;
     }
@@ -688,18 +733,7 @@ function computeAllPositions(ephemerides, times, rxPos) {
   return { prns, times, positions };
 }
 function selectBest(ephs, timeMs) {
-  let best = null;
-  let bestDt = Infinity;
-  for (const eph of ephs) {
-    const dt = Math.abs(timeMs - eph.tocDate.getTime());
-    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
-    if (dt > maxAge) continue;
-    if (dt < bestDt) {
-      bestDt = dt;
-      best = eph;
-    }
-  }
-  return best;
+  return nearestValidEphemeris(ephs, timeMs);
 }
 var GPS_EPOCH_MS2 = START_GPS_TIME.getTime();
 var MS_PER_WEEK = 7 * 864e5;
@@ -940,6 +974,7 @@ function computeVisibilityFromPositions(all, elevationMaskDeg = 10) {
 }
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
+  EPHEMERIS_VALIDITY_MS,
   almanacEpochMs,
   almanacSatPosition,
   computeAllPositions,
@@ -956,9 +991,11 @@ function computeVisibilityFromPositions(all, elevationMaskDeg = 10) {
   glonassAlmanacEpochMs,
   glonassAlmanacPosition,
   glonassPosition,
+  isEphemerisValid,
   keplerPosition,
   maskRadForAzimuth,
   navTimesFromEph,
+  nearestValidEphemeris,
   rangeRate,
   selectEphemeris
 });

package/dist/orbit.d.cts

@@ -91,6 +91,40 @@ declare function glonassAlmanacEpochMs(alm: SbfGlonassAlmanac): number;
  */
 declare function glonassAlmanacPosition(alm: SbfGlonassAlmanac, timeMs: number): AlmanacPosition;
 
+/**
+ * Broadcast-ephemeris validity — whether an ephemeris may position its
+ * satellite at a given epoch, and the nearest one that may.
+ *
+ * Each system's half-window |t − toc| follows its broadcast design, never
+ * wider than the documented validity:
+ *  - GPS / QZSS: the 4-h curve fit centred on toe → ±2 h (IS-GPS-200,
+ *    IS-QZSS; RTKLIB MAXDTOE).
+ *  - Galileo: 2.5 h, inside the 4-h nominal validity (RTKLIB MAXDTOE_GAL).
+ *  - BeiDou: 90 min — ephemerides are refreshed hourly.
+ *  - GLONASS: the immediate data cover a 30-min interval centred on tb →
+ *    ±15 min.
+ *  - NavIC: ±2 h (2-h update cycle).
+ *  - SBAS GEO (MT9 state vector): 6 min (RTKLIB MAXDTOE_SBS) — a GEO
+ *    extrapolated from an hours-old state vector is not an orbit.
+ * Any other system is never valid.
+ *
+ * Time: `tGpsMs` is GPS-scale epoch ms, like everywhere in this module.
+ * GLONASS and SBAS tocs are on the UTC axis (see {@link computeSatPosition}),
+ * so for those t is moved to UTC first — on a ±15-min window the 18-s
+ * leap-second offset matters.
+ */
+
+/** Validity half-window |t − toc| per system letter, in ms. */
+declare const EPHEMERIS_VALIDITY_MS: Readonly<Record<string, number>>;
+/** May `eph` position its satellite at GPS-scale epoch `tGpsMs`? */
+declare function isEphemerisValid(eph: Ephemeris, tGpsMs: number): boolean;
+/**
+ * The ephemeris in `ephs` whose toc is nearest to `tGpsMs` among those valid
+ * at `tGpsMs`, or null when none is. `ephs` may mix satellites; filter by PRN
+ * first if it does.
+ */
+declare function nearestValidEphemeris(ephs: readonly Ephemeris[] | undefined, tGpsMs: number): Ephemeris | null;
+
 /**
  * Satellite orbit computation from broadcast ephemerides.
  * Supports GPS, Galileo, BeiDou (Keplerian) and GLONASS (Runge-Kutta).
@@ -164,8 +198,11 @@ declare function ecefToAzEl(rxX: number, rxY: number, rxZ: number, satX: number,
  * constellation time scale, as a multi-GNSS receiver solves: every extra
  * time scale costs one observation. QZSS and SBAS share the GPS clock —
  * QZSST and SBAS network time are steered to GPST (RTKLIB pntpos does the
- * same) — so a lone SBAS GEO doesn't inflate TDOP/GDOP. Untagged
- * satellites share a single clock (the classic single-system model).
+ * same) — so a lone SBAS GEO doesn't inflate TDOP/GDOP. A time scale
+ * with a single satellite in view is left out: that satellite only fixes
+ * its own clock, adding nothing to PDOP/HDOP/VDOP and an unrecoverable ~1
+ * to TDOP². Untagged satellites share a single clock (the classic
+ * single-system model).
  *
  * @param satAzEls Satellite azimuth/elevation pairs (radians), optionally
  *   tagged with their PRN.
@@ -178,11 +215,12 @@ declare function computeDop(satAzEls: {
     prn?: string;
 }[]): DopValues | null;
 /**
- * Find the best ephemeris for a given PRN at a given time.
- * Selects the one with the closest epoch (toe/toc) to the target time.
+ * Find the best ephemeris for a given PRN at a given time: the one with the
+ * closest epoch (toe/toc) among those within their broadcast validity at
+ * that time ({@link isEphemerisValid}); null when none is valid.
  * @param ephemerides Array of candidate ephemerides
  * @param prn Satellite PRN identifier (e.g. "G01")
- * @param timeMs Target time in Unix milliseconds
+ * @param timeMs Target time, GPS-scale epoch milliseconds
  */
 declare function selectEphemeris(ephemerides: Ephemeris[], prn: string, timeMs: number): Ephemeris | null;
 /**
@@ -300,4 +338,4 @@ declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, num
  */
 declare function computeVisibilityFromPositions(all: AllPositionsData, elevationMaskDeg?: number | number[]): VisibilityResult;
 
-export { type AllPositionsData, type AlmanacPosition, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris };
+export { type AllPositionsData, type AlmanacPosition, type DopValues, EPHEMERIS_VALIDITY_MS, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, isEphemerisValid, keplerPosition, maskRadForAzimuth, navTimesFromEph, nearestValidEphemeris, rangeRate, selectEphemeris };

package/dist/orbit.d.ts

@@ -91,6 +91,40 @@ declare function glonassAlmanacEpochMs(alm: SbfGlonassAlmanac): number;
  */
 declare function glonassAlmanacPosition(alm: SbfGlonassAlmanac, timeMs: number): AlmanacPosition;
 
+/**
+ * Broadcast-ephemeris validity — whether an ephemeris may position its
+ * satellite at a given epoch, and the nearest one that may.
+ *
+ * Each system's half-window |t − toc| follows its broadcast design, never
+ * wider than the documented validity:
+ *  - GPS / QZSS: the 4-h curve fit centred on toe → ±2 h (IS-GPS-200,
+ *    IS-QZSS; RTKLIB MAXDTOE).
+ *  - Galileo: 2.5 h, inside the 4-h nominal validity (RTKLIB MAXDTOE_GAL).
+ *  - BeiDou: 90 min — ephemerides are refreshed hourly.
+ *  - GLONASS: the immediate data cover a 30-min interval centred on tb →
+ *    ±15 min.
+ *  - NavIC: ±2 h (2-h update cycle).
+ *  - SBAS GEO (MT9 state vector): 6 min (RTKLIB MAXDTOE_SBS) — a GEO
+ *    extrapolated from an hours-old state vector is not an orbit.
+ * Any other system is never valid.
+ *
+ * Time: `tGpsMs` is GPS-scale epoch ms, like everywhere in this module.
+ * GLONASS and SBAS tocs are on the UTC axis (see {@link computeSatPosition}),
+ * so for those t is moved to UTC first — on a ±15-min window the 18-s
+ * leap-second offset matters.
+ */
+
+/** Validity half-window |t − toc| per system letter, in ms. */
+declare const EPHEMERIS_VALIDITY_MS: Readonly<Record<string, number>>;
+/** May `eph` position its satellite at GPS-scale epoch `tGpsMs`? */
+declare function isEphemerisValid(eph: Ephemeris, tGpsMs: number): boolean;
+/**
+ * The ephemeris in `ephs` whose toc is nearest to `tGpsMs` among those valid
+ * at `tGpsMs`, or null when none is. `ephs` may mix satellites; filter by PRN
+ * first if it does.
+ */
+declare function nearestValidEphemeris(ephs: readonly Ephemeris[] | undefined, tGpsMs: number): Ephemeris | null;
+
 /**
  * Satellite orbit computation from broadcast ephemerides.
  * Supports GPS, Galileo, BeiDou (Keplerian) and GLONASS (Runge-Kutta).
@@ -164,8 +198,11 @@ declare function ecefToAzEl(rxX: number, rxY: number, rxZ: number, satX: number,
  * constellation time scale, as a multi-GNSS receiver solves: every extra
  * time scale costs one observation. QZSS and SBAS share the GPS clock —
  * QZSST and SBAS network time are steered to GPST (RTKLIB pntpos does the
- * same) — so a lone SBAS GEO doesn't inflate TDOP/GDOP. Untagged
- * satellites share a single clock (the classic single-system model).
+ * same) — so a lone SBAS GEO doesn't inflate TDOP/GDOP. A time scale
+ * with a single satellite in view is left out: that satellite only fixes
+ * its own clock, adding nothing to PDOP/HDOP/VDOP and an unrecoverable ~1
+ * to TDOP². Untagged satellites share a single clock (the classic
+ * single-system model).
  *
  * @param satAzEls Satellite azimuth/elevation pairs (radians), optionally
  *   tagged with their PRN.
@@ -178,11 +215,12 @@ declare function computeDop(satAzEls: {
     prn?: string;
 }[]): DopValues | null;
 /**
- * Find the best ephemeris for a given PRN at a given time.
- * Selects the one with the closest epoch (toe/toc) to the target time.
+ * Find the best ephemeris for a given PRN at a given time: the one with the
+ * closest epoch (toe/toc) among those within their broadcast validity at
+ * that time ({@link isEphemerisValid}); null when none is valid.
  * @param ephemerides Array of candidate ephemerides
  * @param prn Satellite PRN identifier (e.g. "G01")
- * @param timeMs Target time in Unix milliseconds
+ * @param timeMs Target time, GPS-scale epoch milliseconds
  */
 declare function selectEphemeris(ephemerides: Ephemeris[], prn: string, timeMs: number): Ephemeris | null;
 /**
@@ -300,4 +338,4 @@ declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, num
  */
 declare function computeVisibilityFromPositions(all: AllPositionsData, elevationMaskDeg?: number | number[]): VisibilityResult;
 
-export { type AllPositionsData, type AlmanacPosition, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris };
+export { type AllPositionsData, type AlmanacPosition, type DopValues, EPHEMERIS_VALIDITY_MS, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, isEphemerisValid, keplerPosition, maskRadForAzimuth, navTimesFromEph, nearestValidEphemeris, rangeRate, selectEphemeris };

package/dist/orbit.js

@@ -1,4 +1,5 @@
 import {
+  EPHEMERIS_VALIDITY_MS,
   almanacEpochMs,
   almanacSatPosition,
   computeAllPositions,
@@ -13,12 +14,14 @@ import {
   glonassAlmanacEpochMs,
   glonassAlmanacPosition,
   glonassPosition,
+  isEphemerisValid,
   keplerPosition,
   maskRadForAzimuth,
   navTimesFromEph,
+  nearestValidEphemeris,
   rangeRate,
   selectEphemeris
-} from "./chunk-T5M3LVBJ.js";
+} from "./chunk-7T5DLVZL.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
@@ -28,6 +31,7 @@ import "./chunk-LEEU5OIO.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-FIEWO4J4.js";
 export {
+  EPHEMERIS_VALIDITY_MS,
   almanacEpochMs,
   almanacSatPosition,
   computeAllPositions,
@@ -44,9 +48,11 @@ export {
   glonassAlmanacEpochMs,
   glonassAlmanacPosition,
   glonassPosition,
+  isEphemerisValid,
   keplerPosition,
   maskRadForAzimuth,
   navTimesFromEph,
+  nearestValidEphemeris,
   rangeRate,
   selectEphemeris
 };

package/dist/positioning.cjs

@@ -254,6 +254,38 @@ function invertMatrix(A) {
   return M.map((row) => row.slice(n));
 }
 
+// src/orbit/validity.ts
+var H = 36e5;
+var MIN = 6e4;
+var EPHEMERIS_VALIDITY_MS = {
+  G: 2 * H,
+  J: 2 * H,
+  E: 2.5 * H,
+  C: 90 * MIN,
+  R: 15 * MIN,
+  I: 2 * H,
+  S: 6 * MIN
+};
+function tocDistanceMs(eph, tGpsMs) {
+  const t = eph.system === "R" || eph.system === "S" ? tGpsMs - getGpsLeap(new Date(tGpsMs)) * 1e3 : tGpsMs;
+  return Math.abs(eph.tocDate.getTime() - t);
+}
+function nearestValidEphemeris(ephs, tGpsMs) {
+  if (!ephs) return null;
+  let best = null;
+  let bestDt = Infinity;
+  for (const eph of ephs) {
+    const dt = tocDistanceMs(eph, tGpsMs);
+    const w = EPHEMERIS_VALIDITY_MS[eph.system];
+    if (w === void 0 || dt > w) continue;
+    if (dt < bestDt) {
+      bestDt = dt;
+      best = eph;
+    }
+  }
+  return best;
+}
+
 // src/orbit/index.ts
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
@@ -442,24 +474,30 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  const perClock = /* @__PURE__ */ new Map();
+  for (const s of satAzEls) {
+    const c = dopClock(s.prn);
+    perClock.set(c, (perClock.get(c) ?? 0) + 1);
+  }
+  const sats = perClock.size > 1 ? satAzEls.filter((s) => perClock.get(dopClock(s.prn)) > 1) : satAzEls;
+  const n = sats.length;
+  const clocks = [...new Set(sats.map((s) => dopClock(s.prn)))];
   const m = 3 + clocks.length;
   if (n < Math.max(4, m)) return null;
-  const H = [];
-  for (const { az, el, prn } of satAzEls) {
+  const H3 = [];
+  for (const { az, el, prn } of sats) {
     const cosEl = Math.cos(el);
     const row = [cosEl * Math.sin(az), cosEl * Math.cos(az), Math.sin(el)];
     const clock = dopClock(prn);
     for (const c of clocks) row.push(c === clock ? 1 : 0);
-    H.push(row);
+    H3.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));
   for (let i = 0; i < m; i++) {
     for (let j = 0; j < m; j++) {
       let sum = 0;
       for (let k = 0; k < n; k++) {
-        sum += H[k][i] * H[k][j];
+        sum += H3[k][i] * H3[k][j];
       }
       HtH[i][j] = sum;
     }
@@ -495,18 +533,7 @@ function computeSatPosition(eph, timeMs) {
   return keplerPosition(eph, tow);
 }
 function selectBest(ephs, timeMs) {
-  let best = null;
-  let bestDt = Infinity;
-  for (const eph of ephs) {
-    const dt = Math.abs(timeMs - eph.tocDate.getTime());
-    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
-    if (dt > maxAge) continue;
-    if (dt < bestDt) {
-      bestDt = dt;
-      best = eph;
-    }
-  }
-  return best;
+  return nearestValidEphemeris(ephs, timeMs);
 }
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
 var MS_PER_WEEK = 7 * 864e5;
@@ -1302,7 +1329,7 @@ function solveDgnss(rover, base, basePos, ephemerides, timeMs, opts = {}) {
     let cov = null;
     for (let iter = 0; iter < maxIterations; iter++) {
       iterations = iter + 1;
-      const H = [];
+      const H3 = [];
       const v = [];
       const terms = /* @__PURE__ */ new Map();
       const termOf = (g) => {
@@ -1318,11 +1345,11 @@ function solveDgnss(rover, base, basePos, ephemerides, timeMs, opts = {}) {
         const tr = termOf(row.ref);
         const pred = ts.rho - row.g.rhoB + ts.dTropo - (tr.rho - row.ref.rhoB + tr.dTropo);
         v.push(row.z - pred);
-        H.push([ts.u[0] - tr.u[0], ts.u[1] - tr.u[1], ts.u[2] - tr.u[2]]);
+        H3.push([ts.u[0] - tr.u[0], ts.u[1] - tr.u[1], ts.u[2] - tr.u[2]]);
       }
-      const Ht = transpose(H);
+      const Ht = transpose(H3);
       const HtW = matMul(Ht, W);
-      const N = matMul(HtW, H);
+      const N = matMul(HtW, H3);
       const b = matMul(
         HtW,
         v.map((s) => [s])
@@ -1338,7 +1365,7 @@ function solveDgnss(rover, base, basePos, ephemerides, timeMs, opts = {}) {
       z += dx[2][0];
       cov = Ninv;
       rows.forEach((row, i) => {
-        residuals[row.g.prn] = v[i] - (H[i][0] * dx[0][0] + H[i][1] * dx[1][0] + H[i][2] * dx[2][0]);
+        residuals[row.g.prn] = v[i] - (H3[i][0] * dx[0][0] + H3[i][1] * dx[1][0] + H3[i][2] * dx[2][0]);
       });
       if (Math.hypot(dx[0][0], dx[1][0], dx[2][0]) < convergenceM) {
         converged = true;
@@ -1732,19 +1759,19 @@ var RtkFloatEngine = class {
     const PPrior = this.P;
     let xi = [...this.x];
     let K = [];
-    let H = [];
+    let H3 = [];
     for (let it = 0; it < Math.max(o.updateIterations, 1); it++) {
-      H = [];
+      H3 = [];
       const v = [];
       for (const row of rows) {
         const { pred, h } = predict(row, xi);
         let corr = 0;
         for (let j = 0; j < n; j++) corr += h[j] * (xPrior[j] - xi[j]);
         v.push(row.z - pred - corr);
-        H.push(h);
+        H3.push(h);
       }
-      const Ht = transpose(H);
-      const S = matMul(matMul(H, PPrior), Ht);
+      const Ht = transpose(H3);
+      const S = matMul(matMul(H3, PPrior), Ht);
       for (let i = 0; i < rows.length; i++)
         for (let j = 0; j < rows.length; j++) S[i][j] += R[i][j];
       const Sinv = matInv(S);
@@ -1760,7 +1787,7 @@ var RtkFloatEngine = class {
       xi = xPrior.map((s, i) => s + dx[i][0]);
     }
     this.x = xi;
-    const KH = matMul(K, H);
+    const KH = matMul(K, H3);
     const IKH = zeros(n, n);
     for (let i = 0; i < n; i++)
       for (let j = 0; j < n; j++) IKH[i][j] = (i === j ? 1 : 0) - KH[i][j];

package/dist/positioning.js

@@ -40,9 +40,9 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-F6KOFOXN.js";
+} from "./chunk-NOGA5KD7.js";
 import "./chunk-WQX354KD.js";
-import "./chunk-T5M3LVBJ.js";
+import "./chunk-7T5DLVZL.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-CVCTRIWU.js";
 import "./chunk-37QNKGTC.js";

package/package.json

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