gnss-js 2.8.1 → 2.8.2

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

@@ -409,12 +409,18 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  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 clock = dopClock(prn);

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

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

package/dist/index.cjs

@@ -12684,12 +12684,18 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  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 clock = dopClock(prn);

package/dist/index.js

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

package/dist/orbit.cjs

@@ -570,12 +570,18 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  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 clock = dopClock(prn);

package/dist/orbit.d.cts

@@ -164,8 +164,11 @@ declare function ecefToAzEl(rxX: number, rxY: number, rxZ: number, satX: number,
  * 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. Untagged
- * satellites share a single clock (the classic single-system model).
+ * 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

@@ -164,8 +164,11 @@ declare function ecefToAzEl(rxX: number, rxY: number, rxZ: number, satX: number,
  * 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. Untagged
- * satellites share a single clock (the classic single-system model).
+ * 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-T5M3LVBJ.js";
+} from "./chunk-5OWRCF56.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,

package/dist/positioning.cjs

@@ -442,12 +442,18 @@ function dopClock(prn) {
   return sys === "J" || sys === "S" ? "G" : sys;
 }
 function computeDop(satAzEls) {
-  const n = satAzEls.length;
-  const clocks = [...new Set(satAzEls.map((s) => dopClock(s.prn)))];
+  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 clock = dopClock(prn);

package/dist/positioning.js

@@ -40,9 +40,9 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-F6KOFOXN.js";
+} from "./chunk-CLNLOYBX.js";
 import "./chunk-WQX354KD.js";
-import "./chunk-T5M3LVBJ.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.1",
+  "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,