gnss-js 1.30.0 → 1.32.0

package/dist/{chunk-KDM2NVBO.js → chunk-5RNDQQ6J.js}

@@ -5,6 +5,9 @@ import {
 import {
   GPS_PI
 } from "./chunk-IQ5OIMQD.js";
+import {
+  C_LIGHT
+} from "./chunk-FIEWO4J4.js";
 
 // src/rinex/crx-writer.ts
 var ARC_ORDER = 3;
@@ -1417,6 +1420,81 @@ function parseIonex(text) {
   return { epochs, lats, lons, maps };
 }
 
+// src/rinex/bias-sinex.ts
+function parseEpoch(s, openTo) {
+  const m = /^(\d{4}):(\d{3}):(\d{5})$/.exec(s.trim());
+  if (!m) return openTo;
+  const y = +m[1];
+  if (y === 0) return openTo;
+  const doy = +m[2];
+  const sod = +m[3];
+  return Date.UTC(y, 0, 1) + (doy - 1) * 864e5 + sod * 1e3;
+}
+function parseBiasSinex(text) {
+  const records = [];
+  let inBlock = false;
+  for (const line of text.split("\n")) {
+    if (line.startsWith("+BIAS/SOLUTION")) {
+      inBlock = true;
+      continue;
+    }
+    if (line.startsWith("-BIAS/SOLUTION")) {
+      inBlock = false;
+      continue;
+    }
+    if (!inBlock || line.length < 65 || line.startsWith("*")) continue;
+    const type = line.slice(1, 4).trim();
+    if (type !== "DSB" && type !== "OSB" && type !== "ISB") continue;
+    const prn = line.slice(11, 14).trim();
+    const station = line.slice(15, 24).trim();
+    const obs1 = line.slice(25, 29).trim();
+    const obs2 = line.slice(30, 34).trim();
+    const startMs = parseEpoch(line.slice(35, 49), -Infinity);
+    const endMs = parseEpoch(line.slice(50, 64), Infinity);
+    const tail = line.slice(65).trim().split(/\s+/);
+    if (tail.length < 2) continue;
+    const unit = tail[0];
+    const value = Number(tail[1]);
+    if (!Number.isFinite(value)) continue;
+    records.push({
+      type,
+      prn,
+      station,
+      obs1,
+      obs2,
+      startMs,
+      endMs,
+      value,
+      unit
+    });
+  }
+  const bySat = /* @__PURE__ */ new Map();
+  for (const r of records) {
+    if (!r.prn || r.station) continue;
+    let a = bySat.get(r.prn);
+    if (!a) {
+      a = [];
+      bySat.set(r.prn, a);
+    }
+    a.push(r);
+  }
+  return { records, bySat };
+}
+function biasMetres(r) {
+  return r.unit === "ns" ? r.value * 1e-9 * C_LIGHT : r.value;
+}
+function findSatBias(bias, prn, obs1, obs2, tMs) {
+  const recs = bias.bySat.get(prn);
+  if (!recs) return null;
+  for (const r of recs) {
+    if (r.obs1 !== obs1) continue;
+    if (obs2 === void 0 ? r.obs2 !== "" : r.obs2 !== obs2) continue;
+    if (tMs != null && (tMs < r.startMs || tMs > r.endMs)) continue;
+    return r;
+  }
+  return null;
+}
+
 // src/rinex/obs-writer-common.ts
 function createGzipLineSink() {
   const encoder = new TextEncoder();
@@ -1804,6 +1882,9 @@ export {
   writeRinexNav,
   writeRinexNav4,
   parseIonex,
+  parseBiasSinex,
+  biasMetres,
+  findSatBias,
   createGzipLineSink,
   stationHeaderLines,
   writeRinexObsBlob,

package/dist/{chunk-CRFFUU5Y.js → chunk-HOH27N42.js}

@@ -1911,6 +1911,212 @@ function solvePpp(epochs, sp3, opts) {
   };
 }
 
+// src/positioning/ppp-ar.ts
+function wlWavelength(f1, f2) {
+  return C_LIGHT / (f1 - f2);
+}
+function nlWavelength(f1, f2) {
+  return C_LIGHT / (f1 + f2);
+}
+function resolvePppAmbiguities(sats, Qaif, opts = {}) {
+  const wlThreshold = opts.wlThreshold ?? 0.25;
+  const ratioThreshold = opts.ratioThreshold ?? 3;
+  const kept = [];
+  const nWl = [];
+  const rejectedWl = [];
+  for (let i = 0; i < sats.length; i++) {
+    const s = sats[i];
+    const wlFloat = s.mwCyc - (s.wlBiasCyc ?? 0);
+    const r = Math.round(wlFloat);
+    if (Math.abs(wlFloat - r) <= wlThreshold) {
+      kept.push(i);
+      nWl.push(r);
+    } else {
+      rejectedWl.push(s.prn);
+    }
+  }
+  if (kept.length < 2) {
+    return {
+      fixed: false,
+      ratio: 0,
+      refPrn: null,
+      sats: [],
+      rejectedWl
+    };
+  }
+  const n1Float = [];
+  const lamNl = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    const lNl = nlWavelength(s.f1, s.f2);
+    lamNl.push(lNl);
+    const n1 = s.aIF / lNl - s.f2 / (s.f1 - s.f2) * nWl[k] - (s.nlBiasCyc ?? 0);
+    n1Float.push(n1);
+  }
+  let ref = 0;
+  for (let k = 1; k < kept.length; k++) {
+    const e = sats[kept[k]].elevDeg ?? -1;
+    if (e > (sats[kept[ref]].elevDeg ?? -1)) ref = k;
+  }
+  const m = kept.length - 1;
+  const sd = new Float64Array(m);
+  const Qsd = new Float64Array(m * m);
+  const others = [];
+  for (let k = 0; k < kept.length; k++) if (k !== ref) others.push(k);
+  const n = sats.length;
+  const qN1 = (ka, kb) => Qaif[kept[ka] * n + kept[kb]] / (lamNl[ka] * lamNl[kb]);
+  for (let a = 0; a < m; a++) {
+    const ka = others[a];
+    sd[a] = n1Float[ka] - n1Float[ref];
+    for (let b = 0; b < m; b++) {
+      const kb = others[b];
+      Qsd[a * m + b] = qN1(ka, kb) - qN1(ka, ref) - qN1(ref, kb) + qN1(ref, ref);
+    }
+  }
+  const res = lambdaSearch(sd, Qsd, m, 2);
+  const ratio = res ? Math.min(res.ratio, 999.9) : 0;
+  const fixed = !!res && ratio >= ratioThreshold;
+  const refN1 = Math.round(n1Float[ref]);
+  const sdFix = res?.candidates[0];
+  const outSats = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    let n1;
+    if (k === ref) {
+      n1 = refN1;
+    } else {
+      const a = others.indexOf(k);
+      n1 = fixed && sdFix ? refN1 + Math.round(sdFix[a]) : Math.round(n1Float[k]);
+    }
+    const lNl = lamNl[k];
+    const aIF = lNl * n1 + lNl * (s.f2 / (s.f1 - s.f2)) * nWl[k];
+    outSats.push({ prn: s.prn, nWl: nWl[k], n1, aIF });
+  }
+  return {
+    fixed,
+    ratio,
+    refPrn: sats[kept[ref]].prn,
+    sats: outSats,
+    rejectedWl
+  };
+}
+
+// src/positioning/ppp-fcb.ts
+function extractWidelaneArcs(station, epochs, opts = {}) {
+  const slipCyc = opts.slipCyc ?? 4;
+  const gapMs = opts.gapMs ?? 12e4;
+  const minObs = opts.minObs ?? 20;
+  const state = /* @__PURE__ */ new Map();
+  const out = [];
+  const flush = (prn) => {
+    const s = state.get(prn);
+    if (s && s.n >= minObs)
+      out.push({ station, prn, wlFloat: s.sum / s.n, nObs: s.n });
+  };
+  for (const ep of epochs) {
+    for (const o of ep.obs) {
+      if (!o.c1 || !o.c2 || !o.l1 || !o.l2) continue;
+      const lamW = C_LIGHT / (o.f1 - o.f2);
+      const codeNl = (o.f1 * o.c1 + o.f2 * o.c2) / (o.f1 + o.f2);
+      const mw = o.l1 - o.l2 - codeNl / lamW;
+      const s = state.get(o.prn);
+      if (!s || ep.timeMs - s.lastT > gapMs || Math.abs(mw - s.mean) > slipCyc) {
+        if (s) flush(o.prn);
+        state.set(o.prn, { sum: mw, n: 1, mean: mw, lastT: ep.timeMs });
+      } else {
+        s.sum += mw;
+        s.n += 1;
+        s.mean = s.sum / s.n;
+        s.lastT = ep.timeMs;
+      }
+    }
+  }
+  for (const prn of state.keys()) flush(prn);
+  return out;
+}
+function wrapHalf(x) {
+  return x - Math.round(x);
+}
+function circMean(vals, weights) {
+  let s = 0;
+  let c = 0;
+  for (let i = 0; i < vals.length; i++) {
+    const a = 2 * Math.PI * vals[i];
+    s += weights[i] * Math.sin(a);
+    c += weights[i] * Math.cos(a);
+  }
+  if (s === 0 && c === 0) return 0;
+  return wrapHalf(Math.atan2(s, c) / (2 * Math.PI));
+}
+function estimateWidelaneFcb(arcs, opts = {}) {
+  const maxIter = opts.maxIter ?? 50;
+  const tol = opts.tol ?? 1e-4;
+  const fixWindow = opts.fixWindow ?? 0.25;
+  const sats = [...new Set(arcs.map((a) => a.prn))].sort();
+  const stations = [...new Set(arcs.map((a) => a.station))].sort();
+  const refStation = stations[0] ?? "";
+  const satFcb = new Map(sats.map((s) => [s, 0]));
+  const rcvBias = new Map(stations.map((s) => [s, 0]));
+  const nInt = new Array(arcs.length).fill(0);
+  let iterations = 0;
+  for (let it = 0; it < maxIter; it++) {
+    iterations = it + 1;
+    for (let i = 0; i < arcs.length; i++) {
+      const a = arcs[i];
+      nInt[i] = Math.round(
+        a.wlFloat - rcvBias.get(a.station) - satFcb.get(a.prn)
+      );
+    }
+    let maxChange = 0;
+    for (const s of sats) {
+      const vals = [];
+      const w = [];
+      for (let i = 0; i < arcs.length; i++) {
+        const a = arcs[i];
+        if (a.prn !== s) continue;
+        vals.push(a.wlFloat - nInt[i] - rcvBias.get(a.station));
+        w.push(a.nObs);
+      }
+      if (!vals.length) continue;
+      const nv = circMean(vals, w);
+      maxChange = Math.max(maxChange, Math.abs(wrapHalf(nv - satFcb.get(s))));
+      satFcb.set(s, nv);
+    }
+    for (const r of stations) {
+      const vals = [];
+      const w = [];
+      for (let i = 0; i < arcs.length; i++) {
+        const a = arcs[i];
+        if (a.station !== r) continue;
+        vals.push(a.wlFloat - nInt[i] - satFcb.get(a.prn));
+        w.push(a.nObs);
+      }
+      if (!vals.length) continue;
+      rcvBias.set(r, circMean(vals, w));
+    }
+    const shift = rcvBias.get(refStation);
+    if (shift !== 0) {
+      for (const r of stations)
+        rcvBias.set(r, wrapHalf(rcvBias.get(r) - shift));
+      for (const s of sats) satFcb.set(s, wrapHalf(satFcb.get(s) + shift));
+    }
+    if (maxChange < tol) break;
+  }
+  let sumSq = 0;
+  let fixed = 0;
+  for (let i = 0; i < arcs.length; i++) {
+    const a = arcs[i];
+    const resid = wrapHalf(
+      a.wlFloat - nInt[i] - rcvBias.get(a.station) - satFcb.get(a.prn)
+    );
+    sumSq += resid * resid;
+    if (Math.abs(resid) <= fixWindow) fixed++;
+  }
+  const residRms = arcs.length ? Math.sqrt(sumSq / arcs.length) : 0;
+  const fixRate = arcs.length ? fixed / arcs.length : 0;
+  return { satFcb, rcvBias, residRms, fixRate, refStation, iterations };
+}
+
 // src/positioning/index.ts
 var GM_GPS = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -2157,6 +2363,11 @@ export {
   buildPppAntenna,
   createPppCorrections,
   solvePpp,
+  wlWavelength,
+  nlWavelength,
+  resolvePppAmbiguities,
+  extractWidelaneArcs,
+  estimateWidelaneFcb,
   satClockCorrection,
   ionoFree,
   solveSpp

package/dist/index.cjs

@@ -111,6 +111,7 @@ __export(src_exports, {
   analyzeQuality: () => analyzeQuality,
   applyHelmert: () => applyHelmert,
   applyIonoDcb: () => applyIonoDcb,
+  biasMetres: () => biasMetres,
   binexCrc16: () => binexCrc16,
   binexCrc32: () => binexCrc32,
   binexCsum8: () => binexCsum8,
@@ -148,8 +149,11 @@ __export(src_exports, {
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
   ephInfoToEphemeris: () => ephInfoToEphemeris,
+  estimateWidelaneFcb: () => estimateWidelaneFcb,
   euclidean3D: () => euclidean3D,
+  extractWidelaneArcs: () => extractWidelaneArcs,
   fetchSourcetable: () => fetchSourcetable,
+  findSatBias: () => findSatBias,
   fmtD: () => fmtD,
   fmtF: () => fmtF,
   formatUTCTime: () => formatUTCTime,
@@ -232,11 +236,13 @@ __export(src_exports, {
   msmEpochToDate: () => msmEpochToDate,
   navTimesFromEph: () => navTimesFromEph,
   niellMapping: () => niellMapping,
+  nlWavelength: () => nlWavelength,
   nmeaCoordToDecimal: () => nmeaCoordToDecimal,
   oem4Frames: () => oem4Frames,
   padL: () => padL,
   padR: () => padR,
   parseAntex: () => parseAntex,
+  parseBiasSinex: () => parseBiasSinex,
   parseBinex: () => parseBinex,
   parseClk: () => parseClk,
   parseCrxDataLine: () => parseCrxDataLine,
@@ -274,6 +280,7 @@ __export(src_exports, {
   readString: () => readString,
   reportDecodeError: () => reportDecodeError,
   resetGloFreqCache: () => resetGloFreqCache,
+  resolvePppAmbiguities: () => resolvePppAmbiguities,
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   satClockCorrection: () => satClockCorrection,
@@ -299,6 +306,7 @@ __export(src_exports, {
   updateStreamStats: () => updateStreamStats,
   verifyChecksum: () => verifyChecksum,
   vincenty: () => vincenty,
+  wlWavelength: () => wlWavelength,
   writeCrx: () => writeCrx,
   writeModernObsBlob: () => writeModernObsBlob,
   writeRinex2ObsBlob: () => writeRinex2ObsBlob,
@@ -3958,6 +3966,81 @@ function clkBias(clk, prn, tMs) {
   return c[lo] + f * (c[lo + 1] - c[lo]);
 }
 
+// src/rinex/bias-sinex.ts
+function parseEpoch(s, openTo) {
+  const m = /^(\d{4}):(\d{3}):(\d{5})$/.exec(s.trim());
+  if (!m) return openTo;
+  const y = +m[1];
+  if (y === 0) return openTo;
+  const doy = +m[2];
+  const sod = +m[3];
+  return Date.UTC(y, 0, 1) + (doy - 1) * 864e5 + sod * 1e3;
+}
+function parseBiasSinex(text) {
+  const records = [];
+  let inBlock = false;
+  for (const line of text.split("\n")) {
+    if (line.startsWith("+BIAS/SOLUTION")) {
+      inBlock = true;
+      continue;
+    }
+    if (line.startsWith("-BIAS/SOLUTION")) {
+      inBlock = false;
+      continue;
+    }
+    if (!inBlock || line.length < 65 || line.startsWith("*")) continue;
+    const type = line.slice(1, 4).trim();
+    if (type !== "DSB" && type !== "OSB" && type !== "ISB") continue;
+    const prn = line.slice(11, 14).trim();
+    const station = line.slice(15, 24).trim();
+    const obs1 = line.slice(25, 29).trim();
+    const obs2 = line.slice(30, 34).trim();
+    const startMs = parseEpoch(line.slice(35, 49), -Infinity);
+    const endMs = parseEpoch(line.slice(50, 64), Infinity);
+    const tail = line.slice(65).trim().split(/\s+/);
+    if (tail.length < 2) continue;
+    const unit = tail[0];
+    const value = Number(tail[1]);
+    if (!Number.isFinite(value)) continue;
+    records.push({
+      type,
+      prn,
+      station,
+      obs1,
+      obs2,
+      startMs,
+      endMs,
+      value,
+      unit
+    });
+  }
+  const bySat = /* @__PURE__ */ new Map();
+  for (const r of records) {
+    if (!r.prn || r.station) continue;
+    let a = bySat.get(r.prn);
+    if (!a) {
+      a = [];
+      bySat.set(r.prn, a);
+    }
+    a.push(r);
+  }
+  return { records, bySat };
+}
+function biasMetres(r) {
+  return r.unit === "ns" ? r.value * 1e-9 * C_LIGHT : r.value;
+}
+function findSatBias(bias, prn, obs1, obs2, tMs) {
+  const recs = bias.bySat.get(prn);
+  if (!recs) return null;
+  for (const r of recs) {
+    if (r.obs1 !== obs1) continue;
+    if (obs2 === void 0 ? r.obs2 !== "" : r.obs2 !== obs2) continue;
+    if (tMs != null && (tMs < r.startMs || tMs > r.endMs)) continue;
+    return r;
+  }
+  return null;
+}
+
 // src/rinex/obs-writer-common.ts
 function createGzipLineSink() {
   const encoder = new TextEncoder();
@@ -13224,6 +13307,212 @@ function solvePpp(epochs, sp3, opts) {
   };
 }
 
+// src/positioning/ppp-ar.ts
+function wlWavelength(f1, f2) {
+  return C_LIGHT / (f1 - f2);
+}
+function nlWavelength(f1, f2) {
+  return C_LIGHT / (f1 + f2);
+}
+function resolvePppAmbiguities(sats, Qaif, opts = {}) {
+  const wlThreshold = opts.wlThreshold ?? 0.25;
+  const ratioThreshold = opts.ratioThreshold ?? 3;
+  const kept = [];
+  const nWl = [];
+  const rejectedWl = [];
+  for (let i = 0; i < sats.length; i++) {
+    const s = sats[i];
+    const wlFloat = s.mwCyc - (s.wlBiasCyc ?? 0);
+    const r = Math.round(wlFloat);
+    if (Math.abs(wlFloat - r) <= wlThreshold) {
+      kept.push(i);
+      nWl.push(r);
+    } else {
+      rejectedWl.push(s.prn);
+    }
+  }
+  if (kept.length < 2) {
+    return {
+      fixed: false,
+      ratio: 0,
+      refPrn: null,
+      sats: [],
+      rejectedWl
+    };
+  }
+  const n1Float = [];
+  const lamNl = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    const lNl = nlWavelength(s.f1, s.f2);
+    lamNl.push(lNl);
+    const n1 = s.aIF / lNl - s.f2 / (s.f1 - s.f2) * nWl[k] - (s.nlBiasCyc ?? 0);
+    n1Float.push(n1);
+  }
+  let ref = 0;
+  for (let k = 1; k < kept.length; k++) {
+    const e = sats[kept[k]].elevDeg ?? -1;
+    if (e > (sats[kept[ref]].elevDeg ?? -1)) ref = k;
+  }
+  const m = kept.length - 1;
+  const sd2 = new Float64Array(m);
+  const Qsd = new Float64Array(m * m);
+  const others = [];
+  for (let k = 0; k < kept.length; k++) if (k !== ref) others.push(k);
+  const n = sats.length;
+  const qN1 = (ka, kb) => Qaif[kept[ka] * n + kept[kb]] / (lamNl[ka] * lamNl[kb]);
+  for (let a = 0; a < m; a++) {
+    const ka = others[a];
+    sd2[a] = n1Float[ka] - n1Float[ref];
+    for (let b = 0; b < m; b++) {
+      const kb = others[b];
+      Qsd[a * m + b] = qN1(ka, kb) - qN1(ka, ref) - qN1(ref, kb) + qN1(ref, ref);
+    }
+  }
+  const res = lambdaSearch(sd2, Qsd, m, 2);
+  const ratio = res ? Math.min(res.ratio, 999.9) : 0;
+  const fixed = !!res && ratio >= ratioThreshold;
+  const refN1 = Math.round(n1Float[ref]);
+  const sdFix = res?.candidates[0];
+  const outSats = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    let n1;
+    if (k === ref) {
+      n1 = refN1;
+    } else {
+      const a = others.indexOf(k);
+      n1 = fixed && sdFix ? refN1 + Math.round(sdFix[a]) : Math.round(n1Float[k]);
+    }
+    const lNl = lamNl[k];
+    const aIF = lNl * n1 + lNl * (s.f2 / (s.f1 - s.f2)) * nWl[k];
+    outSats.push({ prn: s.prn, nWl: nWl[k], n1, aIF });
+  }
+  return {
+    fixed,
+    ratio,
+    refPrn: sats[kept[ref]].prn,
+    sats: outSats,
+    rejectedWl
+  };
+}
+
+// src/positioning/ppp-fcb.ts
+function extractWidelaneArcs(station, epochs, opts = {}) {
+  const slipCyc = opts.slipCyc ?? 4;
+  const gapMs = opts.gapMs ?? 12e4;
+  const minObs = opts.minObs ?? 20;
+  const state = /* @__PURE__ */ new Map();
+  const out = [];
+  const flush = (prn) => {
+    const s = state.get(prn);
+    if (s && s.n >= minObs)
+      out.push({ station, prn, wlFloat: s.sum / s.n, nObs: s.n });
+  };
+  for (const ep of epochs) {
+    for (const o of ep.obs) {
+      if (!o.c1 || !o.c2 || !o.l1 || !o.l2) continue;
+      const lamW = C_LIGHT / (o.f1 - o.f2);
+      const codeNl = (o.f1 * o.c1 + o.f2 * o.c2) / (o.f1 + o.f2);
+      const mw = o.l1 - o.l2 - codeNl / lamW;
+      const s = state.get(o.prn);
+      if (!s || ep.timeMs - s.lastT > gapMs || Math.abs(mw - s.mean) > slipCyc) {
+        if (s) flush(o.prn);
+        state.set(o.prn, { sum: mw, n: 1, mean: mw, lastT: ep.timeMs });
+      } else {
+        s.sum += mw;
+        s.n += 1;
+        s.mean = s.sum / s.n;
+        s.lastT = ep.timeMs;
+      }
+    }
+  }
+  for (const prn of state.keys()) flush(prn);
+  return out;
+}
+function wrapHalf(x) {
+  return x - Math.round(x);
+}
+function circMean(vals, weights) {
+  let s = 0;
+  let c = 0;
+  for (let i = 0; i < vals.length; i++) {
+    const a = 2 * Math.PI * vals[i];
+    s += weights[i] * Math.sin(a);
+    c += weights[i] * Math.cos(a);
+  }
+  if (s === 0 && c === 0) return 0;
+  return wrapHalf(Math.atan2(s, c) / (2 * Math.PI));
+}
+function estimateWidelaneFcb(arcs, opts = {}) {
+  const maxIter = opts.maxIter ?? 50;
+  const tol = opts.tol ?? 1e-4;
+  const fixWindow = opts.fixWindow ?? 0.25;
+  const sats = [...new Set(arcs.map((a) => a.prn))].sort();
+  const stations = [...new Set(arcs.map((a) => a.station))].sort();
+  const refStation = stations[0] ?? "";
+  const satFcb = new Map(sats.map((s) => [s, 0]));
+  const rcvBias = new Map(stations.map((s) => [s, 0]));
+  const nInt = new Array(arcs.length).fill(0);
+  let iterations = 0;
+  for (let it = 0; it < maxIter; it++) {
+    iterations = it + 1;
+    for (let i = 0; i < arcs.length; i++) {
+      const a = arcs[i];
+      nInt[i] = Math.round(
+        a.wlFloat - rcvBias.get(a.station) - satFcb.get(a.prn)
+      );
+    }
+    let maxChange = 0;
+    for (const s of sats) {
+      const vals = [];
+      const w = [];
+      for (let i = 0; i < arcs.length; i++) {
+        const a = arcs[i];
+        if (a.prn !== s) continue;
+        vals.push(a.wlFloat - nInt[i] - rcvBias.get(a.station));
+        w.push(a.nObs);
+      }
+      if (!vals.length) continue;
+      const nv = circMean(vals, w);
+      maxChange = Math.max(maxChange, Math.abs(wrapHalf(nv - satFcb.get(s))));
+      satFcb.set(s, nv);
+    }
+    for (const r of stations) {
+      const vals = [];
+      const w = [];
+      for (let i = 0; i < arcs.length; i++) {
+        const a = arcs[i];
+        if (a.station !== r) continue;
+        vals.push(a.wlFloat - nInt[i] - satFcb.get(a.prn));
+        w.push(a.nObs);
+      }
+      if (!vals.length) continue;
+      rcvBias.set(r, circMean(vals, w));
+    }
+    const shift = rcvBias.get(refStation);
+    if (shift !== 0) {
+      for (const r of stations)
+        rcvBias.set(r, wrapHalf(rcvBias.get(r) - shift));
+      for (const s of sats) satFcb.set(s, wrapHalf(satFcb.get(s) + shift));
+    }
+    if (maxChange < tol) break;
+  }
+  let sumSq = 0;
+  let fixed = 0;
+  for (let i = 0; i < arcs.length; i++) {
+    const a = arcs[i];
+    const resid = wrapHalf(
+      a.wlFloat - nInt[i] - rcvBias.get(a.station) - satFcb.get(a.prn)
+    );
+    sumSq += resid * resid;
+    if (Math.abs(resid) <= fixWindow) fixed++;
+  }
+  const residRms = arcs.length ? Math.sqrt(sumSq / arcs.length) : 0;
+  const fixRate = arcs.length ? fixed / arcs.length : 0;
+  return { satFcb, rcvBias, residRms, fixRate, refStation, iterations };
+}
+
 // src/positioning/index.ts
 var GM_GPS3 = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -15803,6 +16092,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   analyzeQuality,
   applyHelmert,
   applyIonoDcb,
+  biasMetres,
   binexCrc16,
   binexCrc32,
   binexCsum8,
@@ -15840,8 +16130,11 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateWidelaneFcb,
   euclidean3D,
+  extractWidelaneArcs,
   fetchSourcetable,
+  findSatBias,
   fmtD,
   fmtF,
   formatUTCTime,
@@ -15924,11 +16217,13 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   msmEpochToDate,
   navTimesFromEph,
   niellMapping,
+  nlWavelength,
   nmeaCoordToDecimal,
   oem4Frames,
   padL,
   padR,
   parseAntex,
+  parseBiasSinex,
   parseBinex,
   parseClk,
   parseCrxDataLine,
@@ -15966,6 +16261,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolvePppAmbiguities,
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
@@ -15991,6 +16287,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  wlWavelength,
   writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,