gnss-js 1.31.0 → 1.32.1

package/dist/{chunk-EVFD4UNZ.js → chunk-KD2KJGLQ.js}

@@ -1788,8 +1788,6 @@ function solvePpp(epochs, sp3, opts) {
       P[ci][ci] = CLK_VAR;
     }
     let nSats = 0;
-    let sumPhaseResSq = 0;
-    let nPhaseRes = 0;
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
@@ -1840,15 +1838,37 @@ function solvePpp(epochs, sp3, opts) {
           h[ci] = 1;
           h[IDX_ZWD] = v.mw;
           h[ai] = 1;
-          if (ekfUpdateScalar(x, P, h, innov, rPhase, 4)) {
-            sumPhaseResSq += innov * innov;
-            nPhaseRes++;
-          }
+          ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
+    let sumPhaseResSq = 0;
+    let nPhaseRes = 0;
+    {
+      const bySysRes = /* @__PURE__ */ new Map();
+      for (const v of vis) {
+        const ai = ambIdx.get(v.prn);
+        if (ai === void 0) continue;
+        const sys = v.prn[0];
+        const ci = clockOf(sys);
+        const commonPh = v.range + x[ci] - v.satClkM + v.tropoFixed + x[IDX_ZWD] * v.mw;
+        const innov = v.lIf - (commonPh + v.windupM + x[ai]);
+        if (Math.abs(innov) >= 0.15) continue;
+        const arr = bySysRes.get(sys);
+        if (arr) arr.push(innov);
+        else bySysRes.set(sys, [innov]);
+      }
+      for (const res of bySysRes.values()) {
+        const mean = res.reduce((a, b) => a + b, 0) / res.length;
+        for (const r of res) {
+          const d = r - mean;
+          sumPhaseResSq += d * d;
+          nPhaseRes++;
+        }
+      }
+    }
     for (const [prn, seen] of [...lastSeen]) {
       if (ei - seen > 20) {
         const idx = ambIdx.get(prn);
@@ -2001,6 +2021,122 @@ function resolvePppAmbiguities(sats, Qaif, opts = {}) {
   };
 }
 
+// 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;
@@ -2250,6 +2386,8 @@ export {
   wlWavelength,
   nlWavelength,
   resolvePppAmbiguities,
+  extractWidelaneArcs,
+  estimateWidelaneFcb,
   satClockCorrection,
   ionoFree,
   solveSpp

package/dist/index.cjs

@@ -149,7 +149,9 @@ __export(src_exports, {
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
   ephInfoToEphemeris: () => ephInfoToEphemeris,
+  estimateWidelaneFcb: () => estimateWidelaneFcb,
   euclidean3D: () => euclidean3D,
+  extractWidelaneArcs: () => extractWidelaneArcs,
   fetchSourcetable: () => fetchSourcetable,
   findSatBias: () => findSatBias,
   fmtD: () => fmtD,
@@ -13182,8 +13184,6 @@ function solvePpp(epochs, sp3, opts) {
       P[ci][ci] = CLK_VAR;
     }
     let nSats = 0;
-    let sumPhaseResSq = 0;
-    let nPhaseRes = 0;
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
@@ -13234,15 +13234,37 @@ function solvePpp(epochs, sp3, opts) {
           h[ci] = 1;
           h[IDX_ZWD] = v.mw;
           h[ai] = 1;
-          if (ekfUpdateScalar(x, P, h, innov, rPhase, 4)) {
-            sumPhaseResSq += innov * innov;
-            nPhaseRes++;
-          }
+          ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
+    let sumPhaseResSq = 0;
+    let nPhaseRes = 0;
+    {
+      const bySysRes = /* @__PURE__ */ new Map();
+      for (const v of vis) {
+        const ai = ambIdx.get(v.prn);
+        if (ai === void 0) continue;
+        const sys = v.prn[0];
+        const ci = clockOf(sys);
+        const commonPh = v.range + x[ci] - v.satClkM + v.tropoFixed + x[IDX_ZWD] * v.mw;
+        const innov = v.lIf - (commonPh + v.windupM + x[ai]);
+        if (Math.abs(innov) >= 0.15) continue;
+        const arr = bySysRes.get(sys);
+        if (arr) arr.push(innov);
+        else bySysRes.set(sys, [innov]);
+      }
+      for (const res of bySysRes.values()) {
+        const mean = res.reduce((a, b) => a + b, 0) / res.length;
+        for (const r of res) {
+          const d = r - mean;
+          sumPhaseResSq += d * d;
+          nPhaseRes++;
+        }
+      }
+    }
     for (const [prn, seen] of [...lastSeen]) {
       if (ei - seen > 20) {
         const idx = ambIdx.get(prn);
@@ -13395,6 +13417,122 @@ function resolvePppAmbiguities(sats, Qaif, opts = {}) {
   };
 }
 
+// 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;
@@ -16012,7 +16150,9 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateWidelaneFcb,
   euclidean3D,
+  extractWidelaneArcs,
   fetchSourcetable,
   findSatBias,
   fmtD,

package/dist/index.d.cts

@@ -17,7 +17,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, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.cjs';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateWidelaneFcb, extractWidelaneArcs, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, 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

@@ -17,7 +17,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, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.js';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateWidelaneFcb, extractWidelaneArcs, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, 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

@@ -146,6 +146,8 @@ import {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateWidelaneFcb,
+  extractWidelaneArcs,
   gimSlantIonoDelayL1,
   gimVerticalTec,
   gmst,
@@ -165,7 +167,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-EVFD4UNZ.js";
+} from "./chunk-KD2KJGLQ.js";
 import {
   almanacEpochMs,
   almanacSatPosition,
@@ -472,7 +474,9 @@ export {
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateWidelaneFcb,
   euclidean3D,
+  extractWidelaneArcs,
   fetchSourcetable,
   findSatBias,
   fmtD,

package/dist/positioning.cjs

@@ -24,6 +24,8 @@ __export(positioning_exports, {
   RtkFloatEngine: () => RtkFloatEngine,
   buildPppAntenna: () => buildPppAntenna,
   createPppCorrections: () => createPppCorrections,
+  estimateWidelaneFcb: () => estimateWidelaneFcb,
+  extractWidelaneArcs: () => extractWidelaneArcs,
   gimSlantIonoDelayL1: () => gimSlantIonoDelayL1,
   gimVerticalTec: () => gimVerticalTec,
   gmst: () => gmst,
@@ -2310,8 +2312,6 @@ function solvePpp(epochs, sp3, opts) {
       P[ci][ci] = CLK_VAR;
     }
     let nSats = 0;
-    let sumPhaseResSq = 0;
-    let nPhaseRes = 0;
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
@@ -2362,15 +2362,37 @@ function solvePpp(epochs, sp3, opts) {
           h[ci] = 1;
           h[IDX_ZWD] = v.mw;
           h[ai] = 1;
-          if (ekfUpdateScalar(x, P, h, innov, rPhase, 4)) {
-            sumPhaseResSq += innov * innov;
-            nPhaseRes++;
-          }
+          ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
+    let sumPhaseResSq = 0;
+    let nPhaseRes = 0;
+    {
+      const bySysRes = /* @__PURE__ */ new Map();
+      for (const v of vis) {
+        const ai = ambIdx.get(v.prn);
+        if (ai === void 0) continue;
+        const sys = v.prn[0];
+        const ci = clockOf(sys);
+        const commonPh = v.range + x[ci] - v.satClkM + v.tropoFixed + x[IDX_ZWD] * v.mw;
+        const innov = v.lIf - (commonPh + v.windupM + x[ai]);
+        if (Math.abs(innov) >= 0.15) continue;
+        const arr = bySysRes.get(sys);
+        if (arr) arr.push(innov);
+        else bySysRes.set(sys, [innov]);
+      }
+      for (const res of bySysRes.values()) {
+        const mean = res.reduce((a, b) => a + b, 0) / res.length;
+        for (const r of res) {
+          const d = r - mean;
+          sumPhaseResSq += d * d;
+          nPhaseRes++;
+        }
+      }
+    }
     for (const [prn, seen] of [...lastSeen]) {
       if (ei - seen > 20) {
         const idx = ambIdx.get(prn);
@@ -2523,6 +2545,122 @@ function resolvePppAmbiguities(sats, Qaif, opts = {}) {
   };
 }
 
+// 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_GPS2 = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -2756,6 +2894,8 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateWidelaneFcb,
+  extractWidelaneArcs,
   gimSlantIonoDelayL1,
   gimVerticalTec,
   gmst,

package/dist/positioning.d.cts

@@ -793,6 +793,93 @@ interface ArResult {
  */
 declare function resolvePppAmbiguities(sats: ArSat[], Qaif: Float64Array, opts?: ArOptions): ArResult;
 
+/**
+ * Wide-lane fractional-cycle bias (FCB) estimation for PPP-AR — the
+ * self-contained route to integer ambiguities that needs no external
+ * phase-bias product, only a small network of stations.
+ *
+ * A single receiver cannot resolve ambiguities because its Melbourne–Wübbena
+ * wide-lane float is `Ñ_wl = N (integer) + b_rcv + b_sat`, and the integer
+ * and the satellite fractional bias `b_sat` are entangled. Across a *network*
+ * they separate: every station sees the same `b_sat` for a given satellite
+ * but its own integers and receiver bias, so the fractional part common to
+ * all stations *is* the satellite FCB. This is the classic Ge/Gendt method.
+ *
+ * The estimator alternates:
+ *   1. round each arc to its integer  N = round(Ñ_wl − b_rcv − b_sat);
+ *   2. update each satellite FCB      = circular-mean over stations of
+ *      (Ñ_wl − N − b_rcv);
+ *   3. update each receiver bias       = circular-mean over satellites of
+ *      (Ñ_wl − N − b_sat);
+ *   4. pin the datum (b_rcv of a reference station ≡ 0) — receiver and
+ *      satellite biases otherwise trade off by an arbitrary constant.
+ *
+ * The output `satFcb` feeds `resolvePppAmbiguities`' `wlBiasCyc`.
+ */
+/** Minimal per-satellite observation for wide-lane extraction. */
+interface WlObs {
+    prn: string;
+    f1: number;
+    f2: number;
+    /** Code on f1/f2 (m) and carrier phase on f1/f2 (cycles). */
+    c1: number;
+    c2: number;
+    l1: number;
+    l2: number;
+}
+/**
+ * Extract arc-averaged Melbourne–Wübbena wide-lanes for one station from its
+ * epoch stream. The MW combination
+ *   MW = (L1 − L2) − [(f1·C1 + f2·C2)/(f1+f2)] / (c/(f1−f2))   [cycles]
+ * is geometry-, clock- and (first-order) ionosphere-free, so it holds the
+ * wide-lane ambiguity plus biases and is near-constant within an arc. Arcs
+ * are split on a running-mean jump (cycle slip) or a time gap.
+ */
+declare function extractWidelaneArcs(station: string, epochs: {
+    timeMs: number;
+    obs: WlObs[];
+}[], opts?: {
+    slipCyc?: number;
+    gapMs?: number;
+    minObs?: number;
+}): WlArc[];
+/** One continuous wide-lane arc at a station. */
+interface WlArc {
+    station: string;
+    prn: string;
+    /** Arc-averaged Melbourne–Wübbena wide-lane float (cycles). */
+    wlFloat: number;
+    /** Number of epochs averaged (used as the weight). */
+    nObs: number;
+}
+interface WlFcbOptions {
+    /** Max iterations. Default 50. */
+    maxIter?: number;
+    /** Convergence threshold on the largest bias change (cycles). Default 1e-4. */
+    tol?: number;
+    /** Fix-rate acceptance window |resid| (cycles). Default 0.25. */
+    fixWindow?: number;
+}
+interface WlFcbResult {
+    /** Per-satellite wide-lane FCB (cycles, wrapped to [−0.5, 0.5)). */
+    satFcb: Map<string, number>;
+    /** Per-station wide-lane bias (cycles); the reference station is 0. */
+    rcvBias: Map<string, number>;
+    /** RMS of the post-fit arc residuals Ñ_wl − N − b_rcv − b_sat (cycles). */
+    residRms: number;
+    /** Fraction of arcs whose residual lies within `fixWindow` of an integer. */
+    fixRate: number;
+    /** Reference station pinning the datum. */
+    refStation: string;
+    iterations: number;
+}
+/**
+ * Estimate per-satellite wide-lane FCBs from a network of stations' arc
+ * wide-lane floats. Needs ≥2 stations sharing satellites to break the
+ * integer/bias entanglement.
+ */
+declare function estimateWidelaneFcb(arcs: WlArc[], opts?: WlFcbOptions): WlFcbResult;
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -882,4 +969,4 @@ declare function ionoFree(p1: number, p2: number, f1: number, f2: number): numbe
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateWidelaneFcb, extractWidelaneArcs, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -793,6 +793,93 @@ interface ArResult {
  */
 declare function resolvePppAmbiguities(sats: ArSat[], Qaif: Float64Array, opts?: ArOptions): ArResult;
 
+/**
+ * Wide-lane fractional-cycle bias (FCB) estimation for PPP-AR — the
+ * self-contained route to integer ambiguities that needs no external
+ * phase-bias product, only a small network of stations.
+ *
+ * A single receiver cannot resolve ambiguities because its Melbourne–Wübbena
+ * wide-lane float is `Ñ_wl = N (integer) + b_rcv + b_sat`, and the integer
+ * and the satellite fractional bias `b_sat` are entangled. Across a *network*
+ * they separate: every station sees the same `b_sat` for a given satellite
+ * but its own integers and receiver bias, so the fractional part common to
+ * all stations *is* the satellite FCB. This is the classic Ge/Gendt method.
+ *
+ * The estimator alternates:
+ *   1. round each arc to its integer  N = round(Ñ_wl − b_rcv − b_sat);
+ *   2. update each satellite FCB      = circular-mean over stations of
+ *      (Ñ_wl − N − b_rcv);
+ *   3. update each receiver bias       = circular-mean over satellites of
+ *      (Ñ_wl − N − b_sat);
+ *   4. pin the datum (b_rcv of a reference station ≡ 0) — receiver and
+ *      satellite biases otherwise trade off by an arbitrary constant.
+ *
+ * The output `satFcb` feeds `resolvePppAmbiguities`' `wlBiasCyc`.
+ */
+/** Minimal per-satellite observation for wide-lane extraction. */
+interface WlObs {
+    prn: string;
+    f1: number;
+    f2: number;
+    /** Code on f1/f2 (m) and carrier phase on f1/f2 (cycles). */
+    c1: number;
+    c2: number;
+    l1: number;
+    l2: number;
+}
+/**
+ * Extract arc-averaged Melbourne–Wübbena wide-lanes for one station from its
+ * epoch stream. The MW combination
+ *   MW = (L1 − L2) − [(f1·C1 + f2·C2)/(f1+f2)] / (c/(f1−f2))   [cycles]
+ * is geometry-, clock- and (first-order) ionosphere-free, so it holds the
+ * wide-lane ambiguity plus biases and is near-constant within an arc. Arcs
+ * are split on a running-mean jump (cycle slip) or a time gap.
+ */
+declare function extractWidelaneArcs(station: string, epochs: {
+    timeMs: number;
+    obs: WlObs[];
+}[], opts?: {
+    slipCyc?: number;
+    gapMs?: number;
+    minObs?: number;
+}): WlArc[];
+/** One continuous wide-lane arc at a station. */
+interface WlArc {
+    station: string;
+    prn: string;
+    /** Arc-averaged Melbourne–Wübbena wide-lane float (cycles). */
+    wlFloat: number;
+    /** Number of epochs averaged (used as the weight). */
+    nObs: number;
+}
+interface WlFcbOptions {
+    /** Max iterations. Default 50. */
+    maxIter?: number;
+    /** Convergence threshold on the largest bias change (cycles). Default 1e-4. */
+    tol?: number;
+    /** Fix-rate acceptance window |resid| (cycles). Default 0.25. */
+    fixWindow?: number;
+}
+interface WlFcbResult {
+    /** Per-satellite wide-lane FCB (cycles, wrapped to [−0.5, 0.5)). */
+    satFcb: Map<string, number>;
+    /** Per-station wide-lane bias (cycles); the reference station is 0. */
+    rcvBias: Map<string, number>;
+    /** RMS of the post-fit arc residuals Ñ_wl − N − b_rcv − b_sat (cycles). */
+    residRms: number;
+    /** Fraction of arcs whose residual lies within `fixWindow` of an integer. */
+    fixRate: number;
+    /** Reference station pinning the datum. */
+    refStation: string;
+    iterations: number;
+}
+/**
+ * Estimate per-satellite wide-lane FCBs from a network of stations' arc
+ * wide-lane floats. Needs ≥2 stations sharing satellites to break the
+ * integer/bias entanglement.
+ */
+declare function estimateWidelaneFcb(arcs: WlArc[], opts?: WlFcbOptions): WlFcbResult;
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -882,4 +969,4 @@ declare function ionoFree(p1: number, p2: number, f1: number, f2: number): numbe
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateWidelaneFcb, extractWidelaneArcs, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -3,6 +3,8 @@ import {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateWidelaneFcb,
+  extractWidelaneArcs,
   gimSlantIonoDelayL1,
   gimVerticalTec,
   gmst,
@@ -22,7 +24,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-EVFD4UNZ.js";
+} from "./chunk-KD2KJGLQ.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
@@ -35,6 +37,8 @@ export {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateWidelaneFcb,
+  extractWidelaneArcs,
   gimSlantIonoDelayL1,
   gimVerticalTec,
   gmst,

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.31.0",
+  "version": "1.32.1",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,