gnss-js 1.45.0 → 1.46.0

package/dist/{chunk-XYCECWWS.js → chunk-VNDMPBKK.js}

@@ -1102,12 +1102,23 @@ var RtkFloatEngine = class {
 
 // src/positioning/rtk-postprocess.ts
 var epochsOf = (o) => Array.isArray(o) ? o : o.epochs;
+function medianIntervalMs(epochs) {
+  const times = [...new Set(epochs.map((e) => e.timeMs))].sort((a, b) => a - b);
+  if (times.length < 2) return null;
+  const gaps = [];
+  for (let i = 1; i < times.length; i++) gaps.push(times[i] - times[i - 1]);
+  gaps.sort((a, b) => a - b);
+  return gaps[Math.floor(gaps.length / 2)];
+}
 function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   const tol = opts.pairToleranceMs ?? 500;
   const round = (ms) => Math.round(ms / 1e3) * 1e3;
+  const roverEpochs = epochsOf(rover);
+  const interval = medianIntervalMs(roverEpochs);
+  const engineOpts = opts.maxGapMs != null || interval == null ? opts : { ...opts, maxGapMs: Math.max(1e4, Math.round(interval * 2.5)) };
   const baseByT = /* @__PURE__ */ new Map();
   for (const e of epochsOf(base)) baseByT.set(round(e.timeMs), e);
-  const engine = new RtkFloatEngine(baseEcef, ephemerides, opts);
+  const engine = new RtkFloatEngine(baseEcef, ephemerides, engineOpts);
   const [bLat, bLon] = ecefToGeodetic(baseEcef[0], baseEcef[1], baseEcef[2]);
   const [sLat, cLat, sLon, cLon] = [
     Math.sin(bLat),
@@ -1130,7 +1141,7 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   let solved = 0;
   let unmatched = 0;
   let fixed = 0;
-  const sorted = [...epochsOf(rover)].sort((a, b) => a.timeMs - b.timeMs);
+  const sorted = [...roverEpochs].sort((a, b) => a.timeMs - b.timeMs);
   for (const rEp of sorted) {
     if (seen.has(rEp.timeMs)) continue;
     seen.add(rEp.timeMs);
@@ -1799,6 +1810,7 @@ function solvePpp(epochs, sp3, opts) {
     const latRad = latDeg * Math.PI / 180;
     const epochCtx = corrections?.epochContext?.(epoch.timeMs, rcv) ?? null;
     const vis = [];
+    let epochZhd = 0;
     for (const o of epoch.obs) {
       if (o.c1 === 0 || o.c2 === 0 || o.l1 === 0 || o.l2 === 0) continue;
       const slip = detectSlip(o, ei);
@@ -1827,6 +1839,7 @@ function solvePpp(epochs, sp3, opts) {
       if (!Number.isFinite(pIf) || !Number.isFinite(lIf)) continue;
       const trop = opts.troposphere ?? true ? niellMapping(elRad, latRad, heightM, epoch.timeMs) : { zhd: 0, mh: 0, mw: 0 };
       const { zhd, mh, mw } = trop;
+      epochZhd = zhd;
       const geom = {
         prn: o.prn,
         satEcef: [sat.state.x, sat.state.y, sat.state.z],
@@ -2024,6 +2037,7 @@ function solvePpp(epochs, sp3, opts) {
       position: pos,
       nSats,
       ztdWet: x[IDX_ZWD],
+      ztdHydrostatic: epochZhd,
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
@@ -2062,6 +2076,7 @@ function solvePpp(epochs, sp3, opts) {
     llh,
     series,
     ztdWet: x[IDX_ZWD],
+    ztdHydrostatic: series.length ? series[series.length - 1].ztdHydrostatic : 0,
     epochsUsed: epochs.length,
     convergenceSec,
     finalError3d,

package/dist/index.cjs

@@ -13035,12 +13035,23 @@ var RtkFloatEngine = class {
 
 // src/positioning/rtk-postprocess.ts
 var epochsOf = (o) => Array.isArray(o) ? o : o.epochs;
+function medianIntervalMs(epochs) {
+  const times = [...new Set(epochs.map((e) => e.timeMs))].sort((a, b) => a - b);
+  if (times.length < 2) return null;
+  const gaps = [];
+  for (let i = 1; i < times.length; i++) gaps.push(times[i] - times[i - 1]);
+  gaps.sort((a, b) => a - b);
+  return gaps[Math.floor(gaps.length / 2)];
+}
 function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   const tol = opts.pairToleranceMs ?? 500;
   const round = (ms) => Math.round(ms / 1e3) * 1e3;
+  const roverEpochs = epochsOf(rover);
+  const interval = medianIntervalMs(roverEpochs);
+  const engineOpts = opts.maxGapMs != null || interval == null ? opts : { ...opts, maxGapMs: Math.max(1e4, Math.round(interval * 2.5)) };
   const baseByT = /* @__PURE__ */ new Map();
   for (const e of epochsOf(base)) baseByT.set(round(e.timeMs), e);
-  const engine = new RtkFloatEngine(baseEcef, ephemerides, opts);
+  const engine = new RtkFloatEngine(baseEcef, ephemerides, engineOpts);
   const [bLat, bLon] = ecefToGeodetic(baseEcef[0], baseEcef[1], baseEcef[2]);
   const [sLat, cLat, sLon, cLon] = [
     Math.sin(bLat),
@@ -13063,7 +13074,7 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   let solved = 0;
   let unmatched = 0;
   let fixed = 0;
-  const sorted = [...epochsOf(rover)].sort((a, b) => a.timeMs - b.timeMs);
+  const sorted = [...roverEpochs].sort((a, b) => a.timeMs - b.timeMs);
   for (const rEp of sorted) {
     if (seen.has(rEp.timeMs)) continue;
     seen.add(rEp.timeMs);
@@ -13732,6 +13743,7 @@ function solvePpp(epochs, sp3, opts) {
     const latRad = latDeg * Math.PI / 180;
     const epochCtx = corrections?.epochContext?.(epoch.timeMs, rcv) ?? null;
     const vis = [];
+    let epochZhd = 0;
     for (const o of epoch.obs) {
       if (o.c1 === 0 || o.c2 === 0 || o.l1 === 0 || o.l2 === 0) continue;
       const slip = detectSlip(o, ei);
@@ -13760,6 +13772,7 @@ function solvePpp(epochs, sp3, opts) {
       if (!Number.isFinite(pIf) || !Number.isFinite(lIf)) continue;
       const trop = opts.troposphere ?? true ? niellMapping(elRad, latRad, heightM, epoch.timeMs) : { zhd: 0, mh: 0, mw: 0 };
       const { zhd, mh, mw } = trop;
+      epochZhd = zhd;
       const geom = {
         prn: o.prn,
         satEcef: [sat.state.x, sat.state.y, sat.state.z],
@@ -13957,6 +13970,7 @@ function solvePpp(epochs, sp3, opts) {
       position: pos,
       nSats,
       ztdWet: x[IDX_ZWD],
+      ztdHydrostatic: epochZhd,
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
@@ -13995,6 +14009,7 @@ function solvePpp(epochs, sp3, opts) {
     llh,
     series,
     ztdWet: x[IDX_ZWD],
+    ztdHydrostatic: series.length ? series[series.length - 1].ztdHydrostatic : 0,
     epochsUsed: epochs.length,
     convergenceSec,
     finalError3d,

package/dist/index.js

@@ -164,7 +164,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-XYCECWWS.js";
+} from "./chunk-VNDMPBKK.js";
 import {
   almanacEpochMs,
   almanacSatPosition,

package/dist/positioning.cjs

@@ -1558,12 +1558,23 @@ var RtkFloatEngine = class {
 
 // src/positioning/rtk-postprocess.ts
 var epochsOf = (o) => Array.isArray(o) ? o : o.epochs;
+function medianIntervalMs(epochs) {
+  const times = [...new Set(epochs.map((e) => e.timeMs))].sort((a, b) => a - b);
+  if (times.length < 2) return null;
+  const gaps = [];
+  for (let i = 1; i < times.length; i++) gaps.push(times[i] - times[i - 1]);
+  gaps.sort((a, b) => a - b);
+  return gaps[Math.floor(gaps.length / 2)];
+}
 function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   const tol = opts.pairToleranceMs ?? 500;
   const round = (ms) => Math.round(ms / 1e3) * 1e3;
+  const roverEpochs = epochsOf(rover);
+  const interval = medianIntervalMs(roverEpochs);
+  const engineOpts = opts.maxGapMs != null || interval == null ? opts : { ...opts, maxGapMs: Math.max(1e4, Math.round(interval * 2.5)) };
   const baseByT = /* @__PURE__ */ new Map();
   for (const e of epochsOf(base)) baseByT.set(round(e.timeMs), e);
-  const engine = new RtkFloatEngine(baseEcef, ephemerides, opts);
+  const engine = new RtkFloatEngine(baseEcef, ephemerides, engineOpts);
   const [bLat, bLon] = ecefToGeodetic(baseEcef[0], baseEcef[1], baseEcef[2]);
   const [sLat, cLat, sLon, cLon] = [
     Math.sin(bLat),
@@ -1586,7 +1597,7 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   let solved = 0;
   let unmatched = 0;
   let fixed = 0;
-  const sorted = [...epochsOf(rover)].sort((a, b) => a.timeMs - b.timeMs);
+  const sorted = [...roverEpochs].sort((a, b) => a.timeMs - b.timeMs);
   for (const rEp of sorted) {
     if (seen.has(rEp.timeMs)) continue;
     seen.add(rEp.timeMs);
@@ -2327,6 +2338,7 @@ function solvePpp(epochs, sp3, opts) {
     const latRad = latDeg * Math.PI / 180;
     const epochCtx = corrections?.epochContext?.(epoch.timeMs, rcv) ?? null;
     const vis = [];
+    let epochZhd = 0;
     for (const o of epoch.obs) {
       if (o.c1 === 0 || o.c2 === 0 || o.l1 === 0 || o.l2 === 0) continue;
       const slip = detectSlip(o, ei);
@@ -2355,6 +2367,7 @@ function solvePpp(epochs, sp3, opts) {
       if (!Number.isFinite(pIf) || !Number.isFinite(lIf)) continue;
       const trop = opts.troposphere ?? true ? niellMapping(elRad, latRad, heightM, epoch.timeMs) : { zhd: 0, mh: 0, mw: 0 };
       const { zhd, mh, mw } = trop;
+      epochZhd = zhd;
       const geom = {
         prn: o.prn,
         satEcef: [sat.state.x, sat.state.y, sat.state.z],
@@ -2552,6 +2565,7 @@ function solvePpp(epochs, sp3, opts) {
       position: pos,
       nSats,
       ztdWet: x[IDX_ZWD],
+      ztdHydrostatic: epochZhd,
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
@@ -2590,6 +2604,7 @@ function solvePpp(epochs, sp3, opts) {
     llh,
     series,
     ztdWet: x[IDX_ZWD],
+    ztdHydrostatic: series.length ? series[series.length - 1].ztdHydrostatic : 0,
     epochsUsed: epochs.length,
     convergenceSec,
     finalError3d,

package/dist/positioning.d.cts

@@ -711,6 +711,8 @@ interface PppEpochResult {
     nSats: number;
     /** Estimated zenith wet delay (m). */
     ztdWet: number;
+    /** A priori zenith hydrostatic delay (m, Saastamoinen). ZTD = this + ztdWet. */
+    ztdHydrostatic: number;
     /** RMS of post-fit phase residuals this epoch (m). */
     phaseResRms: number;
 }
@@ -723,6 +725,8 @@ interface PppSolution {
     series: PppEpochResult[];
     /** Estimated zenith wet delay at the end (m). */
     ztdWet: number;
+    /** A priori zenith hydrostatic delay at the end (m). ZTD = this + ztdWet. */
+    ztdHydrostatic: number;
     /** Epochs processed. */
     epochsUsed: number;
     /** Seconds from first epoch until 3D error first stays < 0.1 m

package/dist/positioning.d.ts

@@ -711,6 +711,8 @@ interface PppEpochResult {
     nSats: number;
     /** Estimated zenith wet delay (m). */
     ztdWet: number;
+    /** A priori zenith hydrostatic delay (m, Saastamoinen). ZTD = this + ztdWet. */
+    ztdHydrostatic: number;
     /** RMS of post-fit phase residuals this epoch (m). */
     phaseResRms: number;
 }
@@ -723,6 +725,8 @@ interface PppSolution {
     series: PppEpochResult[];
     /** Estimated zenith wet delay at the end (m). */
     ztdWet: number;
+    /** A priori zenith hydrostatic delay at the end (m). ZTD = this + ztdWet. */
+    ztdHydrostatic: number;
     /** Epochs processed. */
     epochsUsed: number;
     /** Seconds from first epoch until 3D error first stays < 0.1 m

package/dist/positioning.js

@@ -28,7 +28,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-XYCECWWS.js";
+} from "./chunk-VNDMPBKK.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-MIIM4LDY.js";
 import "./chunk-37QNKGTC.js";

package/package.json

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