gnss-js 2.13.0 → 2.13.1

package/dist/{chunk-RWFJEVLJ.js → chunk-PQLPPFIK.js}

@@ -3144,7 +3144,32 @@ function sagnacRotate(x, y, z, travelS) {
   const s = Math.sin(theta);
   return [x * c + y * s, -x * s + y * c, z];
 }
-function satStateAtEmission(source, prn, recvTimeMs, rcv) {
+function satStateAtEmission(source, prn, recvTimeMs, rcv, pseudorangeM) {
+  if (pseudorangeM !== void 0 && pseudorangeM > 1e6) {
+    let tEmit = recvTimeMs - pseudorangeM / C_LIGHT * 1e3;
+    const s0 = source.satState(prn, tEmit);
+    if (!s0) return null;
+    tEmit -= s0.clkS * 1e3;
+    const s = source.satState(prn, tEmit);
+    if (!s) return null;
+    let travelS2 = 0.075;
+    let pos = [s.x, s.y, s.z];
+    for (let iter = 0; iter < 3; iter++) {
+      pos = sagnacRotate(s.x, s.y, s.z, travelS2);
+      travelS2 = Math.hypot(pos[0] - rcv[0], pos[1] - rcv[1], pos[2] - rcv[2]) / C_LIGHT;
+    }
+    const velDotPos2 = s.x * s.vx + s.y * s.vy + s.z * s.vz;
+    const relS2 = -2 * velDotPos2 / (C_LIGHT * C_LIGHT);
+    return {
+      state: {
+        x: pos[0],
+        y: pos[1],
+        z: pos[2],
+        clkM: C_LIGHT * (s.clkS + relS2)
+      },
+      travelS: travelS2
+    };
+  }
   let travelS = 0.075;
   let sx = 0;
   let sy = 0;
@@ -3377,7 +3402,7 @@ var PppEngine = class {
       const o = applySatBiases(source, raw, epoch.timeMs);
       if (!o) continue;
       const slip = detectSlip(o, ei);
-      const sat = satStateAtEmission(source, o.prn, epoch.timeMs, rcv);
+      const sat = satStateAtEmission(source, o.prn, epoch.timeMs, rcv, o.c1);
       if (!sat) continue;
       const [elRad, azRad] = getAer(
         sat.state.x,

package/dist/index.cjs

@@ -17013,7 +17013,32 @@ function sagnacRotate(x, y, z, travelS) {
   const s = Math.sin(theta);
   return [x * c + y * s, -x * s + y * c, z];
 }
-function satStateAtEmission(source, prn, recvTimeMs, rcv) {
+function satStateAtEmission(source, prn, recvTimeMs, rcv, pseudorangeM) {
+  if (pseudorangeM !== void 0 && pseudorangeM > 1e6) {
+    let tEmit = recvTimeMs - pseudorangeM / C_LIGHT * 1e3;
+    const s0 = source.satState(prn, tEmit);
+    if (!s0) return null;
+    tEmit -= s0.clkS * 1e3;
+    const s = source.satState(prn, tEmit);
+    if (!s) return null;
+    let travelS2 = 0.075;
+    let pos = [s.x, s.y, s.z];
+    for (let iter = 0; iter < 3; iter++) {
+      pos = sagnacRotate(s.x, s.y, s.z, travelS2);
+      travelS2 = Math.hypot(pos[0] - rcv[0], pos[1] - rcv[1], pos[2] - rcv[2]) / C_LIGHT;
+    }
+    const velDotPos2 = s.x * s.vx + s.y * s.vy + s.z * s.vz;
+    const relS2 = -2 * velDotPos2 / (C_LIGHT * C_LIGHT);
+    return {
+      state: {
+        x: pos[0],
+        y: pos[1],
+        z: pos[2],
+        clkM: C_LIGHT * (s.clkS + relS2)
+      },
+      travelS: travelS2
+    };
+  }
   let travelS = 0.075;
   let sx = 0;
   let sy = 0;
@@ -17246,7 +17271,7 @@ var PppEngine = class {
       const o = applySatBiases(source, raw, epoch.timeMs);
       if (!o) continue;
       const slip = detectSlip(o, ei);
-      const sat = satStateAtEmission(source, o.prn, epoch.timeMs, rcv);
+      const sat = satStateAtEmission(source, o.prn, epoch.timeMs, rcv, o.c1);
       if (!sat) continue;
       const [elRad, azRad] = getAer(
         sat.state.x,

package/dist/index.js

@@ -218,7 +218,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-RWFJEVLJ.js";
+} from "./chunk-PQLPPFIK.js";
 import "./chunk-H3AKFD25.js";
 import {
   EPHEMERIS_VALIDITY_MS,

package/dist/positioning.cjs

@@ -3753,7 +3753,32 @@ function sagnacRotate(x, y, z, travelS) {
   const s = Math.sin(theta);
   return [x * c + y * s, -x * s + y * c, z];
 }
-function satStateAtEmission(source, prn, recvTimeMs, rcv) {
+function satStateAtEmission(source, prn, recvTimeMs, rcv, pseudorangeM) {
+  if (pseudorangeM !== void 0 && pseudorangeM > 1e6) {
+    let tEmit = recvTimeMs - pseudorangeM / C_LIGHT * 1e3;
+    const s0 = source.satState(prn, tEmit);
+    if (!s0) return null;
+    tEmit -= s0.clkS * 1e3;
+    const s = source.satState(prn, tEmit);
+    if (!s) return null;
+    let travelS2 = 0.075;
+    let pos = [s.x, s.y, s.z];
+    for (let iter = 0; iter < 3; iter++) {
+      pos = sagnacRotate(s.x, s.y, s.z, travelS2);
+      travelS2 = Math.hypot(pos[0] - rcv[0], pos[1] - rcv[1], pos[2] - rcv[2]) / C_LIGHT;
+    }
+    const velDotPos2 = s.x * s.vx + s.y * s.vy + s.z * s.vz;
+    const relS2 = -2 * velDotPos2 / (C_LIGHT * C_LIGHT);
+    return {
+      state: {
+        x: pos[0],
+        y: pos[1],
+        z: pos[2],
+        clkM: C_LIGHT * (s.clkS + relS2)
+      },
+      travelS: travelS2
+    };
+  }
   let travelS = 0.075;
   let sx = 0;
   let sy = 0;
@@ -3986,7 +4011,7 @@ var PppEngine = class {
       const o = applySatBiases(source, raw, epoch.timeMs);
       if (!o) continue;
       const slip = detectSlip(o, ei);
-      const sat = satStateAtEmission(source, o.prn, epoch.timeMs, rcv);
+      const sat = satStateAtEmission(source, o.prn, epoch.timeMs, rcv, o.c1);
       if (!sat) continue;
       const [elRad, azRad] = getAer(
         sat.state.x,

package/dist/positioning.js

@@ -41,7 +41,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-RWFJEVLJ.js";
+} from "./chunk-PQLPPFIK.js";
 import "./chunk-H3AKFD25.js";
 import "./chunk-7T5DLVZL.js";
 import "./chunk-HVXYFUCB.js";

package/package.json

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