gnss-js 2.8.0 → 2.8.2

package/dist/{chunk-GLMQRZDL.js → chunk-5OWRCF56.js}

@@ -404,17 +404,27 @@ function ecefToAzEl(rxX, rxY, rxZ, satX, satY, satZ) {
   const el = Math.atan2(up, Math.sqrt(east * east + north * north));
   return { az: (az + TWO_PI) % TWO_PI, el };
 }
+function dopClock(prn) {
+  const sys = prn?.[0] ?? "";
+  return sys === "J" || sys === "S" ? "G" : sys;
+}
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => s.prn?.[0] ?? ""))];
+  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) {
+  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 sys = prn?.[0] ?? "";
-    for (const c of clocks) row.push(c === sys ? 1 : 0);
+    const clock = dopClock(prn);
+    for (const c of clocks) row.push(c === clock ? 1 : 0);
     H.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));

package/dist/{chunk-PHS7OFHN.js → chunk-CLNLOYBX.js}

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

package/dist/index.cjs

@@ -12679,17 +12679,27 @@ function ecefToAzEl(rxX, rxY, rxZ, satX, satY, satZ) {
   const el = Math.atan2(up, Math.sqrt(east * east + north * north));
   return { az: (az + TWO_PI) % TWO_PI, el };
 }
+function dopClock(prn) {
+  const sys = prn?.[0] ?? "";
+  return sys === "J" || sys === "S" ? "G" : sys;
+}
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => s.prn?.[0] ?? ""))];
+  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 sys = prn?.[0] ?? "";
-    for (const c of clocks) row.push(c === sys ? 1 : 0);
+    const clock = dopClock(prn);
+    for (const c of clocks) row.push(c === clock ? 1 : 0);
     H3.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));

package/dist/index.js

@@ -213,7 +213,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PHS7OFHN.js";
+} from "./chunk-CLNLOYBX.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,
@@ -235,7 +235,7 @@ import {
   navTimesFromEph,
   rangeRate,
   selectEphemeris
-} from "./chunk-GLMQRZDL.js";
+} from "./chunk-5OWRCF56.js";
 import {
   getBdsTime,
   getDateFromBdsTime,

package/dist/orbit.cjs

@@ -565,17 +565,27 @@ function ecefToAzEl(rxX, rxY, rxZ, satX, satY, satZ) {
   const el = Math.atan2(up, Math.sqrt(east * east + north * north));
   return { az: (az + TWO_PI) % TWO_PI, el };
 }
+function dopClock(prn) {
+  const sys = prn?.[0] ?? "";
+  return sys === "J" || sys === "S" ? "G" : sys;
+}
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => s.prn?.[0] ?? ""))];
+  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) {
+  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 sys = prn?.[0] ?? "";
-    for (const c of clocks) row.push(c === sys ? 1 : 0);
+    const clock = dopClock(prn);
+    for (const c of clocks) row.push(c === clock ? 1 : 0);
     H.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));

package/dist/orbit.d.cts

@@ -161,10 +161,14 @@ declare function ecefToAzEl(rxX: number, rxY: number, rxZ: number, satX: number,
  * Compute DOP values from a set of satellite az/el.
  *
  * Satellites tagged with a `prn` get one receiver-clock unknown per
- * constellation (the PRN's system letter), matching how a multi-GNSS
- * receiver solves and how `solveSpp` does: each system runs on its own
- * time scale, so every extra constellation costs one observation. Untagged
- * satellites share a single clock (the classic single-system model).
+ * 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. 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.

package/dist/orbit.d.ts

@@ -161,10 +161,14 @@ declare function ecefToAzEl(rxX: number, rxY: number, rxZ: number, satX: number,
  * Compute DOP values from a set of satellite az/el.
  *
  * Satellites tagged with a `prn` get one receiver-clock unknown per
- * constellation (the PRN's system letter), matching how a multi-GNSS
- * receiver solves and how `solveSpp` does: each system runs on its own
- * time scale, so every extra constellation costs one observation. Untagged
- * satellites share a single clock (the classic single-system model).
+ * 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. 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.

package/dist/orbit.js

@@ -18,7 +18,7 @@ import {
   navTimesFromEph,
   rangeRate,
   selectEphemeris
-} from "./chunk-GLMQRZDL.js";
+} from "./chunk-5OWRCF56.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,

package/dist/positioning.cjs

@@ -437,17 +437,27 @@ function ecefToAzEl(rxX, rxY, rxZ, satX, satY, satZ) {
   const el = Math.atan2(up, Math.sqrt(east * east + north * north));
   return { az: (az + TWO_PI) % TWO_PI, el };
 }
+function dopClock(prn) {
+  const sys = prn?.[0] ?? "";
+  return sys === "J" || sys === "S" ? "G" : sys;
+}
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => s.prn?.[0] ?? ""))];
+  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) {
+  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 sys = prn?.[0] ?? "";
-    for (const c of clocks) row.push(c === sys ? 1 : 0);
+    const clock = dopClock(prn);
+    for (const c of clocks) row.push(c === clock ? 1 : 0);
     H.push(row);
   }
   const HtH = Array.from({ length: m }, () => new Array(m).fill(0));

package/dist/positioning.js

@@ -40,9 +40,9 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PHS7OFHN.js";
+} from "./chunk-CLNLOYBX.js";
 import "./chunk-WQX354KD.js";
-import "./chunk-GLMQRZDL.js";
+import "./chunk-5OWRCF56.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.0",
+  "version": "2.8.2",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,