gnss-js 1.45.1 → 1.47.0

package/dist/{chunk-ST22ULND.js → chunk-75OYCNTY.js}

@@ -1150,31 +1150,50 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
       unmatched++;
       continue;
     }
-    const sol = engine.process(
-      toRtkEpoch(rEp.meas),
-      toRtkEpoch(bEp.meas),
-      rEp.timeMs
-    );
-    if (!sol) continue;
+    let point;
+    if (opts.codeOnly) {
+      const dg = solveDgnss(
+        toRtkEpoch(rEp.meas),
+        toRtkEpoch(bEp.meas),
+        baseEcef,
+        ephemerides,
+        rEp.timeMs,
+        {
+          elevationMaskDeg: opts.elevationMaskDeg,
+          troposphere: opts.troposphere,
+          codeSigmaM: opts.codeSigmaM
+        }
+      );
+      point = dg ? {
+        timeMs: rEp.timeMs,
+        position: dg.position,
+        status: "dgnss",
+        nSats: dg.nSats
+      } : null;
+    } else {
+      const sol = engine.process(
+        toRtkEpoch(rEp.meas),
+        toRtkEpoch(bEp.meas),
+        rEp.timeMs
+      );
+      point = sol ? {
+        timeMs: rEp.timeMs,
+        position: sol.position,
+        status: sol.status,
+        ratio: sol.ratio,
+        nSats: sol.nSats,
+        nFixed: sol.nFixed
+      } : null;
+    }
+    if (!point) continue;
     solved++;
-    if (sol.status === "fixed") fixed++;
+    if (point.status === "fixed") fixed++;
     const [lat, lon, height] = ecefToGeodetic(
-      sol.position[0],
-      sol.position[1],
-      sol.position[2]
+      point.position[0],
+      point.position[1],
+      point.position[2]
     );
-    track.push({
-      timeMs: rEp.timeMs,
-      position: sol.position,
-      lat,
-      lon,
-      height,
-      enu: enuOf(sol.position),
-      status: sol.status,
-      ratio: sol.ratio,
-      nSats: sol.nSats,
-      nFixed: sol.nFixed
-    });
+    track.push({ ...point, lat, lon, height, enu: enuOf(point.position) });
   }
   return { track, solved, unmatched, fixRate: solved ? fixed / solved : 0 };
 }
@@ -1810,6 +1829,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);
@@ -1838,6 +1858,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],
@@ -2035,6 +2056,7 @@ function solvePpp(epochs, sp3, opts) {
       position: pos,
       nSats,
       ztdWet: x[IDX_ZWD],
+      ztdHydrostatic: epochZhd,
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
@@ -2073,6 +2095,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

@@ -13083,31 +13083,50 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
       unmatched++;
       continue;
     }
-    const sol = engine.process(
-      toRtkEpoch(rEp.meas),
-      toRtkEpoch(bEp.meas),
-      rEp.timeMs
-    );
-    if (!sol) continue;
+    let point;
+    if (opts.codeOnly) {
+      const dg = solveDgnss(
+        toRtkEpoch(rEp.meas),
+        toRtkEpoch(bEp.meas),
+        baseEcef,
+        ephemerides,
+        rEp.timeMs,
+        {
+          elevationMaskDeg: opts.elevationMaskDeg,
+          troposphere: opts.troposphere,
+          codeSigmaM: opts.codeSigmaM
+        }
+      );
+      point = dg ? {
+        timeMs: rEp.timeMs,
+        position: dg.position,
+        status: "dgnss",
+        nSats: dg.nSats
+      } : null;
+    } else {
+      const sol = engine.process(
+        toRtkEpoch(rEp.meas),
+        toRtkEpoch(bEp.meas),
+        rEp.timeMs
+      );
+      point = sol ? {
+        timeMs: rEp.timeMs,
+        position: sol.position,
+        status: sol.status,
+        ratio: sol.ratio,
+        nSats: sol.nSats,
+        nFixed: sol.nFixed
+      } : null;
+    }
+    if (!point) continue;
     solved++;
-    if (sol.status === "fixed") fixed++;
+    if (point.status === "fixed") fixed++;
     const [lat, lon, height] = ecefToGeodetic(
-      sol.position[0],
-      sol.position[1],
-      sol.position[2]
+      point.position[0],
+      point.position[1],
+      point.position[2]
     );
-    track.push({
-      timeMs: rEp.timeMs,
-      position: sol.position,
-      lat,
-      lon,
-      height,
-      enu: enuOf(sol.position),
-      status: sol.status,
-      ratio: sol.ratio,
-      nSats: sol.nSats,
-      nFixed: sol.nFixed
-    });
+    track.push({ ...point, lat, lon, height, enu: enuOf(point.position) });
   }
   return { track, solved, unmatched, fixRate: solved ? fixed / solved : 0 };
 }
@@ -13743,6 +13762,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);
@@ -13771,6 +13791,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],
@@ -13968,6 +13989,7 @@ function solvePpp(epochs, sp3, opts) {
       position: pos,
       nSats,
       ztdWet: x[IDX_ZWD],
+      ztdHydrostatic: epochZhd,
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
@@ -14006,6 +14028,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-ST22ULND.js";
+} from "./chunk-75OYCNTY.js";
 import {
   almanacEpochMs,
   almanacSatPosition,

package/dist/positioning.cjs

@@ -1606,31 +1606,50 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
       unmatched++;
       continue;
     }
-    const sol = engine.process(
-      toRtkEpoch(rEp.meas),
-      toRtkEpoch(bEp.meas),
-      rEp.timeMs
-    );
-    if (!sol) continue;
+    let point;
+    if (opts.codeOnly) {
+      const dg = solveDgnss(
+        toRtkEpoch(rEp.meas),
+        toRtkEpoch(bEp.meas),
+        baseEcef,
+        ephemerides,
+        rEp.timeMs,
+        {
+          elevationMaskDeg: opts.elevationMaskDeg,
+          troposphere: opts.troposphere,
+          codeSigmaM: opts.codeSigmaM
+        }
+      );
+      point = dg ? {
+        timeMs: rEp.timeMs,
+        position: dg.position,
+        status: "dgnss",
+        nSats: dg.nSats
+      } : null;
+    } else {
+      const sol = engine.process(
+        toRtkEpoch(rEp.meas),
+        toRtkEpoch(bEp.meas),
+        rEp.timeMs
+      );
+      point = sol ? {
+        timeMs: rEp.timeMs,
+        position: sol.position,
+        status: sol.status,
+        ratio: sol.ratio,
+        nSats: sol.nSats,
+        nFixed: sol.nFixed
+      } : null;
+    }
+    if (!point) continue;
     solved++;
-    if (sol.status === "fixed") fixed++;
+    if (point.status === "fixed") fixed++;
     const [lat, lon, height] = ecefToGeodetic(
-      sol.position[0],
-      sol.position[1],
-      sol.position[2]
+      point.position[0],
+      point.position[1],
+      point.position[2]
     );
-    track.push({
-      timeMs: rEp.timeMs,
-      position: sol.position,
-      lat,
-      lon,
-      height,
-      enu: enuOf(sol.position),
-      status: sol.status,
-      ratio: sol.ratio,
-      nSats: sol.nSats,
-      nFixed: sol.nFixed
-    });
+    track.push({ ...point, lat, lon, height, enu: enuOf(point.position) });
   }
   return { track, solved, unmatched, fixRate: solved ? fixed / solved : 0 };
 }
@@ -2338,6 +2357,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);
@@ -2366,6 +2386,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],
@@ -2563,6 +2584,7 @@ function solvePpp(epochs, sp3, opts) {
       position: pos,
       nSats,
       ztdWet: x[IDX_ZWD],
+      ztdHydrostatic: epochZhd,
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
@@ -2601,6 +2623,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

@@ -437,6 +437,14 @@ interface RtkPostProcessOptions extends RtkFloatOptions {
      * mainly rejects a base gap rather than sub-second offsets.
      */
     pairToleranceMs?: number;
+    /**
+     * Code-only differential (DGNSS): solve each epoch from double-differenced
+     * pseudoranges via {@link solveDgnss} instead of the carrier-phase float
+     * engine. Every epoch is independent (no filter), so there is no ambiguity
+     * fixing — a robust metre→sub-metre fallback that needs no carrier phase.
+     * Track points come back with `status: 'dgnss'`. Default false.
+     */
+    codeOnly?: boolean;
 }
 interface RtkPostProcessResult {
     /** Solved rover epochs in time order. */
@@ -711,6 +719,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 +733,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

@@ -437,6 +437,14 @@ interface RtkPostProcessOptions extends RtkFloatOptions {
      * mainly rejects a base gap rather than sub-second offsets.
      */
     pairToleranceMs?: number;
+    /**
+     * Code-only differential (DGNSS): solve each epoch from double-differenced
+     * pseudoranges via {@link solveDgnss} instead of the carrier-phase float
+     * engine. Every epoch is independent (no filter), so there is no ambiguity
+     * fixing — a robust metre→sub-metre fallback that needs no carrier phase.
+     * Track points come back with `status: 'dgnss'`. Default false.
+     */
+    codeOnly?: boolean;
 }
 interface RtkPostProcessResult {
     /** Solved rover epochs in time order. */
@@ -711,6 +719,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 +733,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-ST22ULND.js";
+} from "./chunk-75OYCNTY.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.1",
+  "version": "1.47.0",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,