gnss-js 2.5.0 → 2.7.0

package/dist/coordinates.d.ts

@@ -14,13 +14,29 @@ declare function deg2dms(deg: number): [number, number, number];
  * Vincenty's inverse formula — orthodromic (great-circle) distance,
  * initial bearing, and final bearing on the WGS84 ellipsoid.
  * All angles in radians.
- * @returns { distance (m), initialBearing (rad), finalBearing (rad) }
+ * @returns { distance (m), initialBearing (rad), finalBearing (rad),
+ *   converged } — `converged` is false for near-antipodal pairs, where the
+ *   iteration fails and a spherical approximation is returned instead
+ *   (distance good to ~0.5 %, bearings unreliable).
  */
 declare function vincenty(lat1: number, lon1: number, lat2: number, lon2: number): {
     distance: number;
     initialBearing: number;
     finalBearing: number;
+    converged: boolean;
 };
+/**
+ * Vincenty's direct formula — the point reached from (lat1, lon1) after
+ * `distance` metres along the geodesic leaving at `bearing` (radians,
+ * clockwise from north) on the WGS84 ellipsoid.
+ * @returns [lat, lon] in radians, lon wrapped to (-π, π]
+ */
+declare function vincentyDirect(lat1: number, lon1: number, bearing: number, distance: number): [number, number];
+/**
+ * Meridian arc length from the equator to geodetic latitude `lat` (radians)
+ * on WGS84, via the rectifying latitude (Helmert series in n, sub-mm).
+ */
+declare function meridianArc(lat: number): number;
 /**
  * Loxodromic (rhumb line) distance and constant bearing on the WGS84 ellipsoid.
  * Uses the ellipsoidal rhumb-line formula with isometric latitudes.
@@ -41,7 +57,8 @@ declare function rhumbLine(lat1: number, lon1: number, lat2: number, lon2: numbe
  */
 declare function euclidean3D(x1: number, y1: number, z1: number, x2: number, y2: number, z2: number): number;
 /**
- * Geodetic midpoint along the great-circle arc.
+ * Midpoint of the geodesic (shortest path on the WGS84 ellipsoid) between
+ * two points: half the Vincenty distance along the initial bearing.
  * @returns [lat, lon] in radians
  */
 declare function greatCircleMidpoint(lat1: number, lon1: number, lat2: number, lon2: number): [number, number];
@@ -57,12 +74,18 @@ declare function horizonDistance(heightMeters: number): number;
  * Convert geodetic coordinates to UTM.
  * @param lat Latitude in radians
  * @param lon Longitude in radians
+ * @returns easting/northing (m), zone (1–60, with the Norway/Svalbard
+ *   exceptions), hemisphere, MGRS latitude `band` letter (C–X), and
+ *   `inRange` — false outside UTM's 80°S–84°N, where the values are
+ *   computed at the clamped latitude and are not a valid UTM position.
  */
 declare function geodeticToUtm(lat: number, lon: number): {
     easting: number;
     northing: number;
     zone: number;
     hemisphere: 'N' | 'S';
+    band: string;
+    inRange: boolean;
 };
 /**
  * Maidenhead grid locator (6 characters) from geodetic coordinates.
@@ -78,4 +101,4 @@ declare function geodeticToMaidenhead(lat: number, lon: number): string;
  */
 declare function geodeticToGeohash(lat: number, lon: number, precision?: number): string;
 
-export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty };
+export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, meridianArc, rad2deg, rhumbLine, vincenty, vincentyDirect };

package/dist/coordinates.js

@@ -7,10 +7,12 @@ import {
   geodeticToUtm,
   greatCircleMidpoint,
   horizonDistance,
+  meridianArc,
   rad2deg,
   rhumbLine,
-  vincenty
-} from "./chunk-MIIM4LDY.js";
+  vincenty,
+  vincentyDirect
+} from "./chunk-CVCTRIWU.js";
 import {
   clampUnit,
   ecefToGeodetic,
@@ -33,7 +35,9 @@ export {
   getEnuDifference,
   greatCircleMidpoint,
   horizonDistance,
+  meridianArc,
   rad2deg,
   rhumbLine,
-  vincenty
+  vincenty,
+  vincentyDirect
 };

package/dist/index.cjs

@@ -113,6 +113,7 @@ __export(src_exports, {
   WGS84_SEMI_MAJOR_AXIS: () => WGS84_SEMI_MAJOR_AXIS,
   WGS84_SEMI_MINOR_AXIS: () => WGS84_SEMI_MINOR_AXIS,
   WarningAccumulator: () => WarningAccumulator,
+  ZeroBaselineEngine: () => ZeroBaselineEngine,
   almanacEpochMs: () => almanacEpochMs,
   almanacSatPosition: () => almanacSatPosition,
   analyzeQuality: () => analyzeQuality,
@@ -260,6 +261,7 @@ __export(src_exports, {
   lambdaSearch: () => lambdaSearch,
   looksLikeRtcm2: () => looksLikeRtcm2,
   maskRadForAzimuth: () => maskRadForAzimuth,
+  meridianArc: () => meridianArc,
   moonEcef: () => moonEcef,
   navTimesFromEph: () => navTimesFromEph,
   niellMapping: () => niellMapping,
@@ -354,6 +356,7 @@ __export(src_exports, {
   updateStreamStats: () => updateStreamStats,
   verifyChecksum: () => verifyChecksum,
   vincenty: () => vincenty,
+  vincentyDirect: () => vincentyDirect,
   wlWavelength: () => wlWavelength,
   writeCrx: () => writeCrx,
   writeModernObsBlob: () => writeModernObsBlob,
@@ -1124,7 +1127,12 @@ function vincenty(lat1, lon1, lat2, lon2) {
       (cosU2 * sinLambda2) ** 2 + (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda2) ** 2
     );
     if (sinSigma === 0) {
-      return { distance: 0, initialBearing: 0, finalBearing: 0 };
+      return {
+        distance: 0,
+        initialBearing: 0,
+        finalBearing: 0,
+        converged: true
+      };
     }
     cosSigma = sinU1 * sinU2 + cosU1 * cosU2 * cosLambda2;
     sigma = Math.atan2(sinSigma, cosSigma);
@@ -1154,7 +1162,7 @@ function vincenty(lat1, lon1, lat2, lon2) {
       cosLat1 * Math.sin(-dLon),
       -sinLat1 * cosLat2 + cosLat1 * sinLat2 * Math.cos(-dLon)
     );
-    return { distance: distance2, initialBearing: ib, finalBearing: fb };
+    return { distance: distance2, initialBearing: ib, finalBearing: fb, converged: false };
   }
   const uSq = cos2Alpha * (WGS84_SEMI_MAJOR_AXIS ** 2 - WGS84_SEMI_MINOR_AXIS ** 2) / WGS84_SEMI_MINOR_AXIS ** 2;
   const A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));
@@ -1171,10 +1179,53 @@ function vincenty(lat1, lon1, lat2, lon2) {
     cosU1 * sinLambda,
     -sinU1 * cosU2 + cosU1 * sinU2 * cosLambda
   );
-  return { distance, initialBearing, finalBearing };
+  return { distance, initialBearing, finalBearing, converged: true };
+}
+function vincentyDirect(lat1, lon1, bearing, distance) {
+  const a = WGS84_SEMI_MAJOR_AXIS, b = WGS84_SEMI_MINOR_AXIS, f = WGS84_FLATTENING;
+  const sinA1 = Math.sin(bearing), cosA1 = Math.cos(bearing);
+  const tanU1 = (1 - f) * Math.tan(lat1);
+  const cosU1 = 1 / Math.sqrt(1 + tanU1 ** 2), sinU1 = tanU1 * cosU1;
+  const sigma1 = Math.atan2(tanU1, cosA1);
+  const sinAlpha = cosU1 * sinA1;
+  const cos2Alpha = 1 - sinAlpha ** 2;
+  const uSq = cos2Alpha * (a ** 2 - b ** 2) / b ** 2;
+  const A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));
+  const B = uSq / 1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq)));
+  let sigma = distance / (b * A);
+  let sigmaPrev;
+  let cos2SigmaM, sinSigma, cosSigma;
+  let iter = 0;
+  do {
+    cos2SigmaM = Math.cos(2 * sigma1 + sigma);
+    sinSigma = Math.sin(sigma);
+    cosSigma = Math.cos(sigma);
+    const deltaSigma = B * sinSigma * (cos2SigmaM + B / 4 * (cosSigma * (-1 + 2 * cos2SigmaM ** 2) - B / 6 * cos2SigmaM * (-3 + 4 * sinSigma ** 2) * (-3 + 4 * cos2SigmaM ** 2)));
+    sigmaPrev = sigma;
+    sigma = distance / (b * A) + deltaSigma;
+  } while (Math.abs(sigma - sigmaPrev) > 1e-12 && ++iter < 200);
+  const tmp = sinU1 * sinSigma - cosU1 * cosSigma * cosA1;
+  const lat2 = Math.atan2(
+    sinU1 * cosSigma + cosU1 * sinSigma * cosA1,
+    (1 - f) * Math.sqrt(sinAlpha ** 2 + tmp ** 2)
+  );
+  const lambda = Math.atan2(
+    sinSigma * sinA1,
+    cosU1 * cosSigma - sinU1 * sinSigma * cosA1
+  );
+  const C2 = f / 16 * cos2Alpha * (4 + f * (4 - 3 * cos2Alpha));
+  const L = lambda - (1 - C2) * f * sinAlpha * (sigma + C2 * sinSigma * (cos2SigmaM + C2 * cosSigma * (-1 + 2 * cos2SigmaM ** 2)));
+  let lon2 = lon1 + L;
+  lon2 = Math.atan2(Math.sin(lon2), Math.cos(lon2));
+  return [lat2, lon2];
+}
+function meridianArc(lat) {
+  const n = WGS84_FLATTENING / (2 - WGS84_FLATTENING);
+  const n2 = n * n, n3 = n2 * n, n4 = n3 * n;
+  const mu = lat + (-3 * n / 2 + 9 * n3 / 16) * Math.sin(2 * lat) + (15 * n2 / 16 - 15 * n4 / 32) * Math.sin(4 * lat) + -35 * n3 / 48 * Math.sin(6 * lat) + 315 * n4 / 512 * Math.sin(8 * lat);
+  return WGS84_SEMI_MAJOR_AXIS / (1 + n) * (1 + n2 / 4 + n4 / 64) * mu;
 }
 function rhumbLine(lat1, lon1, lat2, lon2) {
-  const dLat = lat2 - lat1;
   let dLon = lon2 - lon1;
   const MAX_LAT = Math.PI / 2 - 1e-10;
   const isometricLat = (lat) => {
@@ -1188,33 +1239,28 @@ function rhumbLine(lat1, lon1, lat2, lon2) {
   const psi1 = isometricLat(lat1);
   const psi2 = isometricLat(lat2);
   const dPsi = psi2 - psi1;
-  const q = Math.abs(dPsi) > 1e-12 ? dLat / dPsi : Math.cos(lat1);
   if (Math.abs(dLon) > Math.PI) {
     dLon = dLon > 0 ? -(2 * Math.PI - dLon) : 2 * Math.PI + dLon;
   }
   const bearing = Math.atan2(dLon, dPsi);
-  const e22 = WGS84_ECCENTRICITY_SQUARED;
-  const e4 = e22 * e22;
-  const e6 = e4 * e22;
-  const meridionalRadius = WGS84_SEMI_MAJOR_AXIS * (1 - e22) * (1 + 3 / 4 * e22 + 45 / 64 * e4 + 175 / 256 * e6);
-  const distance = Math.sqrt(dLat ** 2 + q ** 2 * dLon ** 2) * meridionalRadius;
+  const dM = meridianArc(lat2) - meridianArc(lat1);
+  let q;
+  if (Math.abs(dPsi) > 1e-6) {
+    q = dM / dPsi;
+  } else {
+    const latM = (lat1 + lat2) / 2;
+    q = WGS84_SEMI_MAJOR_AXIS * Math.cos(latM) / Math.sqrt(1 - WGS84_ECCENTRICITY_SQUARED * Math.sin(latM) ** 2);
+  }
+  const distance = Math.sqrt(dM ** 2 + (q * dLon) ** 2);
   return { distance, bearing };
 }
 function euclidean3D(x12, y1, z1, x22, y2, z2) {
   return Math.sqrt((x22 - x12) ** 2 + (y2 - y1) ** 2 + (z2 - z1) ** 2);
 }
 function greatCircleMidpoint(lat1, lon1, lat2, lon2) {
-  const dLon = lon2 - lon1;
-  const Bx = Math.cos(lat2) * Math.cos(dLon);
-  const By = Math.cos(lat2) * Math.sin(dLon);
-  const sumX = Math.cos(lat1) + Bx;
-  const denom = Math.sqrt(sumX ** 2 + By ** 2);
-  if (denom < 1e-15) {
-    return [(lat1 + lat2) / 2, lon1];
-  }
-  const latMid = Math.atan2(Math.sin(lat1) + Math.sin(lat2), denom);
-  const lonMid = lon1 + Math.atan2(By, sumX);
-  return [latMid, lonMid];
+  const { distance, initialBearing } = vincenty(lat1, lon1, lat2, lon2);
+  if (distance === 0) return [lat1, lon1];
+  return vincentyDirect(lat1, lon1, initialBearing, distance / 2);
 }
 function horizonDistance(heightMeters) {
   if (heightMeters <= 0) return 0;
@@ -1227,18 +1273,36 @@ function horizonDistance(heightMeters) {
 // src/coordinates/projections.ts
 var e2 = WGS84_ECCENTRICITY_SQUARED;
 var k0 = 0.9996;
+var UTM_BANDS = "CDEFGHJKLMNPQRSTUVWX";
+function utmZone(latDeg, lonDeg) {
+  let zone = Math.floor((lonDeg + 180) / 6) + 1;
+  if (zone > 60) zone = 1;
+  if (latDeg >= 56 && latDeg < 64 && lonDeg >= 3 && lonDeg < 12) return 32;
+  if (latDeg >= 72 && latDeg <= 84) {
+    if (lonDeg >= 0 && lonDeg < 9) return 31;
+    if (lonDeg >= 9 && lonDeg < 21) return 33;
+    if (lonDeg >= 21 && lonDeg < 33) return 35;
+    if (lonDeg >= 33 && lonDeg < 42) return 37;
+  }
+  return zone;
+}
 function geodeticToUtm(lat, lon) {
+  const rawLatDeg = lat * 180 / Math.PI;
+  const inRange = rawLatDeg >= -80 && rawLatDeg <= 84;
   const MAX_UTM_LAT = 84 * (Math.PI / 180);
   lat = Math.max(-MAX_UTM_LAT, Math.min(MAX_UTM_LAT, lat));
   const latDeg = lat * 180 / Math.PI;
+  lon = Math.atan2(Math.sin(lon), Math.cos(lon));
   const lonDeg = lon * 180 / Math.PI;
-  const zone = Math.floor((lonDeg + 180) / 6) + 1;
+  const zone = utmZone(latDeg, lonDeg);
+  const band = UTM_BANDS[Math.max(0, Math.min(19, Math.floor((latDeg + 80) / 8)))];
   const lon0 = ((zone - 1) * 6 - 180 + 3) * (Math.PI / 180);
   const ep2 = e2 / (1 - e2);
   const N = WGS84_SEMI_MAJOR_AXIS / Math.sqrt(1 - e2 * Math.sin(lat) ** 2);
   const T = Math.tan(lat) ** 2;
   const C2 = ep2 * Math.cos(lat) ** 2;
-  const A = Math.cos(lat) * (lon - lon0);
+  const dLon0 = Math.atan2(Math.sin(lon - lon0), Math.cos(lon - lon0));
+  const A = Math.cos(lat) * dLon0;
   const M = WGS84_SEMI_MAJOR_AXIS * ((1 - e2 / 4 - 3 * e2 ** 2 / 64 - 5 * e2 ** 3 / 256) * lat - (3 * e2 / 8 + 3 * e2 ** 2 / 32 + 45 * e2 ** 3 / 1024) * Math.sin(2 * lat) + (15 * e2 ** 2 / 256 + 45 * e2 ** 3 / 1024) * Math.sin(4 * lat) - 35 * e2 ** 3 / 3072 * Math.sin(6 * lat));
   const easting = k0 * N * (A + (1 - T + C2) * A ** 3 / 6 + (5 - 18 * T + T ** 2 + 72 * C2 - 58 * ep2) * A ** 5 / 120) + 5e5;
   let northing = k0 * (M + N * Math.tan(lat) * (A ** 2 / 2 + (5 - T + 9 * C2 + 4 * C2 ** 2) * A ** 4 / 24 + (61 - 58 * T + T ** 2 + 600 * C2 - 330 * ep2) * A ** 6 / 720));
@@ -1246,7 +1310,7 @@ function geodeticToUtm(lat, lon) {
   if (hemisphere === "S") {
     northing += 1e7;
   }
-  return { easting, northing, zone, hemisphere };
+  return { easting, northing, zone, hemisphere, band, inRange };
 }
 function geodeticToMaidenhead(lat, lon) {
   const latDeg = Math.max(0, Math.min(179.999999, lat * 180 / Math.PI + 90));
@@ -14636,6 +14700,162 @@ function postProcessRtk(base, rover, ephemerides, baseEcef, opts = {}) {
   return { track, solved, unmatched, fixRate: solved ? fixed / solved : 0 };
 }
 
+// src/positioning/zero-baseline.ts
+function carrierFreqHz2(sys, band, gloChannel) {
+  if (sys === "R") {
+    if (band === "1") return 1602e6 + (gloChannel ?? 0) * 562500;
+    if (band === "2") return 1246e6 + (gloChannel ?? 0) * 437500;
+    return FREQ.R?.[band] ?? 0;
+  }
+  return FREQ[sys]?.[band] ?? 0;
+}
+function median(a) {
+  if (a.length === 0) return 0;
+  const s = [...a].sort((x, y) => x - y);
+  const m = s.length >> 1;
+  return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
+}
+var ZeroBaselineEngine = class {
+  ref;
+  refSys;
+  slip;
+  /** Per (receiver|prn|code) continuous phase-SD arc. */
+  arcs = /* @__PURE__ */ new Map();
+  /** Per receiver: integrated phase-smoothed clock (m). */
+  phaseClock = /* @__PURE__ */ new Map();
+  /** Per (receiver|system): accumulated ISB (m) + count. */
+  isb = /* @__PURE__ */ new Map();
+  constructor(opts) {
+    this.ref = opts.reference;
+    this.refSys = opts.referenceSystem ?? "G";
+    this.slip = opts.slipThreshM ?? 0.05;
+  }
+  reset() {
+    this.arcs.clear();
+    this.phaseClock.clear();
+    this.isb.clear();
+  }
+  /** ISB/IFB estimates accumulated so far, per receiver × system. */
+  biases() {
+    const out = [];
+    for (const [key, { sum, n }] of this.isb) {
+      const [receiver, system] = key.split("|");
+      out.push({ receiver, system, biasM: sum / n, n });
+    }
+    return out;
+  }
+  /**
+   * Ingest one time-aligned instant: observations per receiver id. Returns one
+   * {@link ReceiverClockSample} per non-reference receiver present this epoch.
+   */
+  process(epoch) {
+    const refObs = epoch.obs[this.ref];
+    if (!refObs) return [];
+    const refByKey = /* @__PURE__ */ new Map();
+    for (const o of refObs)
+      if (o.pr != null && Number.isFinite(o.pr))
+        refByKey.set(`${o.prn}|${o.code}`, o);
+    const out = [];
+    for (const rx of Object.keys(epoch.obs)) {
+      if (rx === this.ref) continue;
+      out.push(
+        this.processReceiver(rx, epoch.obs[rx], refByKey, epoch.timeMs)
+      );
+    }
+    return out;
+  }
+  processReceiver(rx, obs, refByKey, t) {
+    const codeBySys = /* @__PURE__ */ new Map();
+    const phase = [];
+    for (const o of obs) {
+      if (o.pr == null || !Number.isFinite(o.pr)) continue;
+      const r = refByKey.get(`${o.prn}|${o.code}`);
+      if (!r || r.pr == null) continue;
+      const sys = o.prn[0];
+      (codeBySys.get(sys) ?? codeBySys.set(sys, []).get(sys)).push(
+        o.pr - r.pr
+      );
+      if (o.cp != null && r.cp != null) {
+        const f = carrierFreqHz2(
+          sys,
+          o.code[0],
+          o.gloChannel ?? r.gloChannel ?? null
+        );
+        if (f > 0)
+          phase.push({
+            key: `${rx}|${o.prn}|${o.code}`,
+            sdM: (o.cp - r.cp) * (C_LIGHT / f),
+            lock: o.lockTimeS
+          });
+      }
+    }
+    if (codeBySys.size === 0)
+      return {
+        t,
+        receiver: rx,
+        clockOffsetM: null,
+        clockOffsetPhaseM: this.phaseClock.get(rx) ?? null,
+        used: 0,
+        rejected: 0,
+        residualRmsM: 0
+      };
+    const sysMedian = /* @__PURE__ */ new Map();
+    for (const [sys, v] of codeBySys) sysMedian.set(sys, median(v));
+    const allSd = [];
+    for (const v of codeBySys.values()) allSd.push(...v);
+    const clk = sysMedian.get(this.refSys) ?? median(allSd);
+    for (const [sys, med] of sysMedian) {
+      if (sys === this.refSys) continue;
+      const cur = this.isb.get(`${rx}|${sys}`) ?? { sum: 0, n: 0 };
+      cur.sum += med - clk;
+      cur.n += 1;
+      this.isb.set(`${rx}|${sys}`, cur);
+    }
+    let ss = 0, used = 0, rejected = 0;
+    for (const [sys, v] of codeBySys) {
+      const isb = sys === this.refSys ? 0 : sysMedian.get(sys) - clk;
+      for (const sd2 of v) {
+        const res = sd2 - clk - isb;
+        if (Math.abs(res) > 5) {
+          rejected++;
+          continue;
+        }
+        ss += res * res;
+        used++;
+      }
+    }
+    const residualRmsM = used ? Math.sqrt(ss / used) : 0;
+    const deltas = [];
+    for (const p of phase) {
+      const arc = this.arcs.get(p.key);
+      this.arcs.set(p.key, { lastPhaseSdM: p.sdM, lastLock: p.lock });
+      if (!arc) continue;
+      const lockReset = p.lock != null && arc.lastLock != null && p.lock < arc.lastLock;
+      if (lockReset) continue;
+      deltas.push(p.sdM - arc.lastPhaseSdM);
+    }
+    let phaseClockM = this.phaseClock.get(rx) ?? null;
+    if (phaseClockM == null) {
+      phaseClockM = clk;
+    } else if (deltas.length) {
+      const md = median(deltas);
+      const good = deltas.filter((d) => Math.abs(d - md) < this.slip);
+      const step = good.length ? good.reduce((a, b) => a + b, 0) / good.length : md;
+      phaseClockM += step;
+    }
+    this.phaseClock.set(rx, phaseClockM);
+    return {
+      t,
+      receiver: rx,
+      clockOffsetM: clk,
+      clockOffsetPhaseM: phaseClockM,
+      used,
+      rejected,
+      residualRmsM
+    };
+  }
+};
+
 // src/positioning/sbas.ts
 var P2_11 = 2 ** -11;
 var P2_31 = 2 ** -31;
@@ -17183,7 +17403,7 @@ async function connectToMountpoint(proxyUrl, info, signal) {
 
 // src/analysis/stats-util.ts
 var MIN_ARC_LENGTH = 10;
-function median(values) {
+function median2(values) {
   const s = [...values].sort((a, b) => a - b);
   const m = s.length >> 1;
   return s.length % 2 ? s[m] : (s[m - 1] + s[m]) / 2;
@@ -17262,7 +17482,7 @@ var MultipathAccumulator = class {
       this.closeArc(prn, band, refBand, bandState);
     } else if (bandState.arc.rawMp.length > 0) {
       const recent = bandState.arc.rawMp.slice(-MP_JUMP_WINDOW);
-      if (Math.abs(mp - median(recent)) > MP_JUMP_M) {
+      if (Math.abs(mp - median2(recent)) > MP_JUMP_M) {
         this.closeArc(prn, band, refBand, bandState);
       }
     }
@@ -17783,7 +18003,7 @@ var IonoAccumulator = class {
   closeArc(prn, pairKey, ps) {
     const arc = ps.arc;
     if (arc.times.length >= MIN_ARC_LENGTH) {
-      const level = median(arc.l4.map((v, k) => v - arc.p4[k]));
+      const level = median2(arc.l4.map((v, k) => v - arc.p4[k]));
       const toTecu = arc.fi * arc.fi / TEC_FACTOR / (arc.gamma - 1);
       let satArcs = this.closed.get(prn);
       if (!satArcs) {
@@ -19387,6 +19607,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   WGS84_SEMI_MAJOR_AXIS,
   WGS84_SEMI_MINOR_AXIS,
   WarningAccumulator,
+  ZeroBaselineEngine,
   almanacEpochMs,
   almanacSatPosition,
   analyzeQuality,
@@ -19534,6 +19755,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   lambdaSearch,
   looksLikeRtcm2,
   maskRadForAzimuth,
+  meridianArc,
   moonEcef,
   navTimesFromEph,
   niellMapping,
@@ -19628,6 +19850,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  vincentyDirect,
   wlWavelength,
   writeCrx,
   writeModernObsBlob,

package/dist/index.d.cts

@@ -3,7 +3,7 @@ export { D as DAYS_MJD2000_MJD, a as DAYS_MJD_JULIAN, H as HourCode, M as MILLIS
 export { WGS84_ECCENTRICITY_SQUARED, WGS84_FLATTENING, WGS84_SEMI_MAJOR_AXIS, WGS84_SEMI_MINOR_AXIS } from './constants.cjs';
 export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OMEGA_E, i as OnSlipDetected, S as SYSTEM_NAMES, j as SYS_SHORT, k as buildGloChannelMap, l as buildObsIndices, m as formatUTCTime, n as getFreq, o as gloFreq } from './gnss-BtXrG3zH.cjs';
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.cjs';
-export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.cjs';
+export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, meridianArc, rad2deg, rhumbLine, vincenty, vincentyDirect } from './coordinates.cjs';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.cjs';
 export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.cjs';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.cjs';
@@ -19,7 +19,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
-export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.cjs';
+export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.cjs';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';

package/dist/index.d.ts

@@ -3,7 +3,7 @@ export { D as DAYS_MJD2000_MJD, a as DAYS_MJD_JULIAN, H as HourCode, M as MILLIS
 export { WGS84_ECCENTRICITY_SQUARED, WGS84_FLATTENING, WGS84_SEMI_MAJOR_AXIS, WGS84_SEMI_MINOR_AXIS } from './constants.js';
 export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OMEGA_E, i as OnSlipDetected, S as SYSTEM_NAMES, j as SYS_SHORT, k as buildGloChannelMap, l as buildObsIndices, m as formatUTCTime, n as getFreq, o as gloFreq } from './gnss-DhnEr0Dy.js';
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.js';
-export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.js';
+export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, meridianArc, rad2deg, rhumbLine, vincenty, vincentyDirect } from './coordinates.js';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.js';
 export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.js';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.js';
@@ -19,7 +19,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
-export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.js';
+export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';

package/dist/index.js

@@ -177,6 +177,7 @@ import {
   RtkFloatEngine,
   SbasProcessor,
   Sp3EphemerisSource,
+  ZeroBaselineEngine,
   antexFreqOfGroup,
   applyOrbitClock,
   buildPppAntenna,
@@ -212,7 +213,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-B2XWO4RP.js";
+} from "./chunk-BWLLHLLS.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,
@@ -279,10 +280,12 @@ import {
   geodeticToUtm,
   greatCircleMidpoint,
   horizonDistance,
+  meridianArc,
   rad2deg,
   rhumbLine,
-  vincenty
-} from "./chunk-MIIM4LDY.js";
+  vincenty,
+  vincentyDirect
+} from "./chunk-CVCTRIWU.js";
 import {
   clampUnit,
   ecefToGeodetic,
@@ -489,6 +492,7 @@ export {
   WGS84_SEMI_MAJOR_AXIS,
   WGS84_SEMI_MINOR_AXIS,
   WarningAccumulator,
+  ZeroBaselineEngine,
   almanacEpochMs,
   almanacSatPosition,
   analyzeQuality,
@@ -636,6 +640,7 @@ export {
   lambdaSearch,
   looksLikeRtcm2,
   maskRadForAzimuth,
+  meridianArc,
   moonEcef,
   navTimesFromEph,
   niellMapping,
@@ -730,6 +735,7 @@ export {
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  vincentyDirect,
   wlWavelength,
   writeCrx,
   writeModernObsBlob,