gnss-js 2.13.0 → 2.14.0

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

@@ -2935,6 +2935,103 @@ function solidEarthTide(rcv, sun, moon) {
   return disp;
 }
 
+// src/positioning/ocean-loading.ts
+var BLQ_CONSTITUENTS = [
+  "M2",
+  "S2",
+  "N2",
+  "K2",
+  "K1",
+  "O1",
+  "P1",
+  "Q1",
+  "Mf",
+  "Mm",
+  "Ssa"
+];
+var ARGS = [
+  [140519e-9, 2, -2, 0, 0],
+  // M2
+  [145444e-9, 0, 0, 0, 0],
+  // S2
+  [13788e-8, 2, -3, 1, 0],
+  // N2
+  [145842e-9, 2, 0, 0, 0],
+  // K2
+  [72921e-9, 1, 0, 0, 0.25],
+  // K1
+  [67598e-9, 1, -2, 0, -0.25],
+  // O1
+  [72523e-9, -1, 0, 0, -0.25],
+  // P1
+  [64959e-9, 1, -3, 1, -0.25],
+  // Q1
+  [53234e-10, 0, 2, 0, 0],
+  // Mf
+  [26392e-10, 0, 1, -1, 0],
+  // Mm
+  [3982e-10, 2, 0, 0, 0]
+  // Ssa
+];
+var D2R2 = Math.PI / 180;
+function oceanLoadingEnu(c, utcMs) {
+  const day = Math.floor(utcMs / 864e5);
+  const fday = (utcMs - day * 864e5) / 1e3;
+  const days = day - Date.UTC(1975, 0, 1) / 864e5 + 1;
+  const t = (27392.500528 + 1.000000035 * days) / 36525;
+  const t2 = t * t;
+  const t3 = t2 * t;
+  const a = [
+    fday,
+    (279.69668 + 36000.768930485 * t + 303e-6 * t2) * D2R2,
+    // H0
+    (270.434358 + 481267.88314137 * t - 1133e-6 * t2 + 19e-7 * t3) * D2R2,
+    // S0
+    (334.329653 + 4069.0340329577 * t - 0.010325 * t2 - 12e-6 * t3) * D2R2,
+    // P0
+    2 * Math.PI
+  ];
+  const dp = [0, 0, 0];
+  for (let i = 0; i < ARGS.length; i++) {
+    let ang = 0;
+    for (let j = 0; j < 5; j++) ang += a[j] * ARGS[i][j];
+    for (let k = 0; k < 3; k++)
+      dp[k] += c.amplitude[i][k] * Math.cos(ang - c.phase[i][k] * D2R2);
+  }
+  return [-dp[1], -dp[2], dp[0]];
+}
+function parseBlq(text) {
+  const out = /* @__PURE__ */ new Map();
+  const lines = text.split(/\r?\n/);
+  for (let i = 0; i < lines.length; i++) {
+    const l = lines[i];
+    if (!l.trim() || l.startsWith("$$")) continue;
+    const name = l.trim().split(/\s+/)[0].toUpperCase();
+    const rows = [];
+    let j = i + 1;
+    while (j < lines.length && rows.length < 6) {
+      const r = lines[j];
+      j++;
+      if (r.startsWith("$$") || !r.trim()) continue;
+      const v = r.trim().split(/\s+/).map(Number);
+      if (v.length !== 11 || v.some((x) => !Number.isFinite(x))) break;
+      rows.push(v);
+    }
+    if (rows.length === 6) {
+      out.set(name, {
+        amplitude: BLQ_CONSTITUENTS.map(
+          (_, k) => [rows[0][k], rows[1][k], rows[2][k]]
+        ),
+        phase: BLQ_CONSTITUENTS.map(
+          (_, k) => [rows[3][k], rows[4][k], rows[5][k]]
+        )
+      });
+      i = j - 1;
+    }
+  }
+  return out;
+}
+
 // src/positioning/ppp-corrections.ts
 function applyCorrections(corrections, geom, ctx) {
   if (corrections.evaluate) return corrections.evaluate(geom, ctx);
@@ -2952,6 +3049,14 @@ function createPppCorrections(cfg) {
       const sun = sunEcef(timeMs);
       const moon = moonEcef(timeMs);
       const tideDisp = useTides ? solidEarthTide(rcvEcef, sun, moon) : [0, 0, 0];
+      if (cfg.oceanLoading) {
+        const [e, n, u] = oceanLoadingEnu(cfg.oceanLoading, timeMs);
+        const [lat, lon] = ecefToGeodetic(rcvEcef[0], rcvEcef[1], rcvEcef[2]);
+        const d = enuToEcef([e, n, u], lat, lon);
+        tideDisp[0] += d[0];
+        tideDisp[1] += d[1];
+        tideDisp[2] += d[2];
+      }
       return { timeMs, sun, moon, tideDisp };
     },
     evaluate(geom, ctx) {
@@ -3006,7 +3111,7 @@ function createPppCorrections(cfg) {
           rangeM += interpPcv(off, zenithDeg);
         }
       }
-      if (useTides && ctx) {
+      if (ctx) {
         rangeM -= dot(e, ctx.tideDisp);
       }
       if (useWindup && frame) {
@@ -3144,7 +3249,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;
@@ -3367,7 +3497,7 @@ var PppEngine = class {
     for (const ci of clkIdx.values()) P[ci][ci] = CLK_VAR;
     P[IDX_ZWD][IDX_ZWD] += ztdQ * (ei > 0 ? 1 : 0);
     const rcv = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
-    const [latRad, , heightM] = ecefToGeodetic(rcv[0], rcv[1], rcv[2]);
+    const [latRad, lonRad, heightM] = ecefToGeodetic(rcv[0], rcv[1], rcv[2]);
     const epochCtx = corrections?.epochContext?.(epoch.timeMs, rcv) ?? null;
     const vis = [];
     let epochZhd = 0;
@@ -3377,7 +3507,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,
@@ -3400,7 +3530,17 @@ var PppEngine = class {
       const pIf = g * o.c1 - (g - 1) * o.c2;
       const lIf = g * (o.l1 * lam1) - (g - 1) * (o.l2 * lam2);
       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 trop = opts.troposphere ?? true ? opts.troposphereModel ? opts.troposphereModel.delays(
+        epoch.timeMs,
+        latRad,
+        lonRad,
+        heightM,
+        elRad,
+        azRad
+      ) : {
+        ...niellMapping(elRad, latRad, heightM, epoch.timeMs),
+        gradient: 0
+      } : { zhd: 0, mh: 0, mw: 0, gradient: 0 };
       const { zhd, mh, mw } = trop;
       epochZhd = zhd;
       const geom = {
@@ -3425,7 +3565,7 @@ var PppEngine = class {
         lIf,
         satClkM: sat.state.clkM,
         range,
-        tropoFixed: zhd * mh + corr.rangeM,
+        tropoFixed: zhd * mh + trop.gradient + corr.rangeM,
         mw,
         windupM: corr.phaseWindupM,
         slip,
@@ -4380,6 +4520,9 @@ export {
   sunEcef,
   moonEcef,
   solidEarthTide,
+  BLQ_CONSTITUENTS,
+  oceanLoadingEnu,
+  parseBlq,
   createPppCorrections,
   Sp3EphemerisSource,
   PppEngine,