gnss-js 1.33.0 → 1.35.0

package/dist/{chunk-PVOF4IGX.js → chunk-BLGBMRX7.js}

@@ -461,15 +461,15 @@ function matMul(A, B) {
   const n = A.length;
   const k = B.length;
   const m = B[0]?.length ?? 0;
-  const C = zeros(n, m);
+  const C2 = zeros(n, m);
   for (let i = 0; i < n; i++) {
     for (let l = 0; l < k; l++) {
       const a = A[i][l];
       if (a === 0) continue;
-      for (let j = 0; j < m; j++) C[i][j] += a * B[l][j];
+      for (let j = 0; j < m; j++) C2[i][j] += a * B[l][j];
     }
   }
-  return C;
+  return C2;
 }
 function transpose(A) {
   const n = A.length;
@@ -1649,6 +1649,24 @@ function solvePpp(epochs, sp3, opts) {
   const clockSeeded = /* @__PURE__ */ new Set();
   const ambIdx = /* @__PURE__ */ new Map();
   const lastSeen = /* @__PURE__ */ new Map();
+  const arcs = [];
+  const arcSnap = /* @__PURE__ */ new Map();
+  const flushArc = (prn, minEpochs = 10) => {
+    const s = arcSnap.get(prn);
+    arcSnap.delete(prn);
+    if (!s || s.nEpochs < minEpochs) return;
+    arcs.push({
+      prn,
+      aIF: s.aIF,
+      mwCyc: s.mwCyc,
+      f1: s.f1,
+      f2: s.f2,
+      meanElevDeg: s.elevSum / s.nEpochs,
+      nEpochs: s.nEpochs,
+      startMs: s.startMs,
+      endMs: s.endMs
+    });
+  };
   const mwState = /* @__PURE__ */ new Map();
   const detectSlip = (o, ei) => {
     const lamW = C_LIGHT / (o.f1 - o.f2);
@@ -1766,7 +1784,9 @@ function solvePpp(epochs, sp3, opts) {
         tropoFixed: zhd * mh + corr.rangeM,
         mw,
         windupM: corr.phaseWindupM,
-        slip
+        slip,
+        f1: o.f1,
+        f2: o.f2
       });
     }
     const zwd0 = x[IDX_ZWD];
@@ -1791,6 +1811,7 @@ function solvePpp(epochs, sp3, opts) {
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
+      if (v.slip) flushArc(v.prn);
       let ai = ambIdx.get(v.prn);
       if (ai === void 0 || v.slip) {
         const bInit = v.lIf - (common + v.windupM);
@@ -1824,6 +1845,7 @@ function solvePpp(epochs, sp3, opts) {
         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) {
+          flushArc(v.prn);
           x[ai] = v.lIf - (commonPh + v.windupM);
           for (let j = 0; j < P.length; j++) {
             P[ai][j] = 0;
@@ -1841,6 +1863,29 @@ function solvePpp(epochs, sp3, opts) {
           ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
+      {
+        const elevDeg = v.elRad * 180 / Math.PI;
+        const mwMean = mwState.get(v.prn)?.mean ?? 0;
+        const s = arcSnap.get(v.prn);
+        if (s) {
+          s.aIF = x[ai];
+          s.mwCyc = mwMean;
+          s.elevSum += elevDeg;
+          s.nEpochs++;
+          s.endMs = epoch.timeMs;
+        } else {
+          arcSnap.set(v.prn, {
+            aIF: x[ai],
+            mwCyc: mwMean,
+            f1: v.f1,
+            f2: v.f2,
+            elevSum: elevDeg,
+            nEpochs: 1,
+            startMs: epoch.timeMs,
+            endMs: epoch.timeMs
+          });
+        }
+      }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
@@ -1874,6 +1919,7 @@ function solvePpp(epochs, sp3, opts) {
         const idx = ambIdx.get(prn);
         if (idx !== void 0) dropState(idx);
         lastSeen.delete(prn);
+        flushArc(prn);
       }
     }
     const pos = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
@@ -1913,6 +1959,7 @@ function solvePpp(epochs, sp3, opts) {
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
+  for (const prn of [...arcSnap.keys()]) flushArc(prn);
   const finalPos = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
   const llh = ecefToGeodetic(finalPos[0], finalPos[1], finalPos[2]);
   const finalError3d = truth ? Math.hypot(
@@ -1927,7 +1974,8 @@ function solvePpp(epochs, sp3, opts) {
     ztdWet: x[IDX_ZWD],
     epochsUsed: epochs.length,
     convergenceSec,
-    finalError3d
+    finalError3d,
+    arcs
   };
 }
 
@@ -2068,34 +2116,34 @@ function circMean(vals, weights) {
   if (s === 0 && c === 0) return 0;
   return wrapHalf(Math.atan2(s, c) / (2 * Math.PI));
 }
-function estimateWidelaneFcb(arcs, opts = {}) {
+function decomposeFcb(entries, 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 sats = [...new Set(entries.map((e) => e.sat))].sort();
+  const stations = [...new Set(entries.map((e) => e.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);
+  const nInt = new Array(entries.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];
+    for (let i = 0; i < entries.length; i++) {
+      const e = entries[i];
       nInt[i] = Math.round(
-        a.wlFloat - rcvBias.get(a.station) - satFcb.get(a.prn)
+        e.value - rcvBias.get(e.station) - satFcb.get(e.sat)
       );
     }
     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);
+      for (let i = 0; i < entries.length; i++) {
+        const e = entries[i];
+        if (e.sat !== s) continue;
+        vals.push(e.value - nInt[i] - rcvBias.get(e.station));
+        w.push(e.weight);
       }
       if (!vals.length) continue;
       const nv = circMean(vals, w);
@@ -2105,11 +2153,11 @@ function estimateWidelaneFcb(arcs, opts = {}) {
     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);
+      for (let i = 0; i < entries.length; i++) {
+        const e = entries[i];
+        if (e.station !== r) continue;
+        vals.push(e.value - nInt[i] - satFcb.get(e.sat));
+        w.push(e.weight);
       }
       if (!vals.length) continue;
       rcvBias.set(r, circMean(vals, w));
@@ -2124,18 +2172,88 @@ function estimateWidelaneFcb(arcs, opts = {}) {
   }
   let sumSq = 0;
   let fixed = 0;
-  for (let i = 0; i < arcs.length; i++) {
-    const a = arcs[i];
+  for (let i = 0; i < entries.length; i++) {
+    const e = entries[i];
     const resid = wrapHalf(
-      a.wlFloat - nInt[i] - rcvBias.get(a.station) - satFcb.get(a.prn)
+      e.value - nInt[i] - rcvBias.get(e.station) - satFcb.get(e.sat)
     );
     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;
+  const residRms = entries.length ? Math.sqrt(sumSq / entries.length) : 0;
+  const fixRate = entries.length ? fixed / entries.length : 0;
   return { satFcb, rcvBias, residRms, fixRate, refStation, iterations };
 }
+function estimateWidelaneFcb(arcs, opts = {}) {
+  return decomposeFcb(
+    arcs.map((a) => ({
+      station: a.station,
+      sat: a.prn,
+      value: a.wlFloat,
+      weight: a.nObs
+    })),
+    opts
+  );
+}
+var C = 299792458;
+function estimateNarrowlaneFcb(arcs, wl, opts = {}) {
+  const entries = [];
+  let wlRejected = 0;
+  for (const a of arcs) {
+    const satWl = wl.satFcb.get(a.prn);
+    const rcvWl = wl.rcvBias.get(a.station);
+    if (satWl === void 0 || rcvWl === void 0) continue;
+    const wlFloat = a.mwCyc - satWl - rcvWl;
+    const nWl = Math.round(wlFloat);
+    if (Math.abs(wlFloat - nWl) > (opts.fixWindow ?? 0.25)) {
+      wlRejected++;
+      continue;
+    }
+    const lamNl = C / (a.f1 + a.f2);
+    const n1Float = a.aIF / lamNl - a.f2 / (a.f1 - a.f2) * nWl;
+    entries.push({
+      station: a.station,
+      sat: a.prn,
+      value: n1Float,
+      weight: a.nEpochs
+    });
+  }
+  const res = decomposeFcb(entries, { fixWindow: 0.15, ...opts });
+  return { ...res, usedArcs: entries.length, wlRejected };
+}
+function estimateNetworkFcbs(arcs, opts = {}) {
+  const minArcEpochs = opts.minArcEpochs ?? 120;
+  const minWlObs = opts.minWlObs ?? 40;
+  const fixWindow = opts.fixWindow ?? 0.15;
+  const satWlFcb = /* @__PURE__ */ new Map();
+  const satNlFcb = /* @__PURE__ */ new Map();
+  const perSystem = {};
+  const systems = [...new Set(arcs.map((a) => a.prn[0]))].sort();
+  for (const sys of systems) {
+    const sysArcs = arcs.filter((a) => a.prn[0] === sys);
+    const wlArcs = sysArcs.filter((a) => a.nEpochs >= minWlObs).map((a) => ({
+      station: a.station,
+      prn: a.prn,
+      wlFloat: a.mwCyc,
+      nObs: a.nEpochs
+    }));
+    if (wlArcs.length < 3) continue;
+    const wl = estimateWidelaneFcb(wlArcs, { fixWindow });
+    const nlInput = sysArcs.filter((a) => a.nEpochs >= minArcEpochs);
+    const nl = estimateNarrowlaneFcb(nlInput, wl, { fixWindow });
+    for (const [prn, v] of wl.satFcb) satWlFcb.set(prn, v);
+    for (const [prn, v] of nl.satFcb) satNlFcb.set(prn, v);
+    perSystem[sys] = {
+      wlFixRate: wl.fixRate,
+      wlResidRms: wl.residRms,
+      nlUsedArcs: nl.usedArcs,
+      nlFixRate: nl.fixRate,
+      nlResidRms: nl.residRms,
+      wlRejected: nl.wlRejected
+    };
+  }
+  return { satWlFcb, satNlFcb, perSystem };
+}
 
 // src/positioning/index.ts
 var GM_GPS = 3986005e8;
@@ -2388,6 +2506,8 @@ export {
   resolvePppAmbiguities,
   extractWidelaneArcs,
   estimateWidelaneFcb,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   satClockCorrection,
   ionoFree,
   solveSpp

package/dist/index.cjs

@@ -149,6 +149,8 @@ __export(src_exports, {
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
   ephInfoToEphemeris: () => ephInfoToEphemeris,
+  estimateNarrowlaneFcb: () => estimateNarrowlaneFcb,
+  estimateNetworkFcbs: () => estimateNetworkFcbs,
   estimateWidelaneFcb: () => estimateWidelaneFcb,
   euclidean3D: () => euclidean3D,
   extractWidelaneArcs: () => extractWidelaneArcs,
@@ -1067,7 +1069,7 @@ function vincenty(lat1, lon1, lat2, lon2) {
   let lambdaPrev;
   let sinSigma, cosSigma, sigma;
   let sinAlpha, cos2Alpha, cos2SigmaM;
-  let C;
+  let C2;
   let iter = 0;
   do {
     const sinLambda2 = Math.sin(lambda);
@@ -1083,9 +1085,9 @@ function vincenty(lat1, lon1, lat2, lon2) {
     sinAlpha = cosU1 * cosU2 * sinLambda2 / sinSigma;
     cos2Alpha = 1 - sinAlpha ** 2;
     cos2SigmaM = cos2Alpha !== 0 ? cosSigma - 2 * sinU1 * sinU2 / cos2Alpha : 0;
-    C = WGS84_FLATTENING / 16 * cos2Alpha * (4 + WGS84_FLATTENING * (4 - 3 * cos2Alpha));
+    C2 = WGS84_FLATTENING / 16 * cos2Alpha * (4 + WGS84_FLATTENING * (4 - 3 * cos2Alpha));
     lambdaPrev = lambda;
-    lambda = lon2 - lon1 + (1 - C) * WGS84_FLATTENING * sinAlpha * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM ** 2)));
+    lambda = lon2 - lon1 + (1 - C2) * WGS84_FLATTENING * sinAlpha * (sigma + C2 * sinSigma * (cos2SigmaM + C2 * cosSigma * (-1 + 2 * cos2SigmaM ** 2)));
   } while (Math.abs(lambda - lambdaPrev) > 1e-12 && ++iter < 200);
   if (iter >= 200) {
     const cosLat1 = Math.cos(lat1), sinLat1 = Math.sin(lat1);
@@ -1189,11 +1191,11 @@ function geodeticToUtm(lat, lon) {
   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 C = ep2 * Math.cos(lat) ** 2;
+  const C2 = ep2 * Math.cos(lat) ** 2;
   const A = Math.cos(lat) * (lon - lon0);
   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 + C) * A ** 3 / 6 + (5 - 18 * T + T ** 2 + 72 * C - 58 * ep2) * A ** 5 / 120) + 5e5;
-  let northing = k0 * (M + N * Math.tan(lat) * (A ** 2 / 2 + (5 - T + 9 * C + 4 * C ** 2) * A ** 4 / 24 + (61 - 58 * T + T ** 2 + 600 * C - 330 * ep2) * A ** 6 / 720));
+  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));
   const hemisphere = latDeg >= 0 ? "N" : "S";
   if (hemisphere === "S") {
     northing += 1e7;
@@ -11857,15 +11859,15 @@ function matMul(A, B) {
   const n = A.length;
   const k = B.length;
   const m = B[0]?.length ?? 0;
-  const C = zeros(n, m);
+  const C2 = zeros(n, m);
   for (let i = 0; i < n; i++) {
     for (let l = 0; l < k; l++) {
       const a = A[i][l];
       if (a === 0) continue;
-      for (let j = 0; j < m; j++) C[i][j] += a * B[l][j];
+      for (let j = 0; j < m; j++) C2[i][j] += a * B[l][j];
     }
   }
-  return C;
+  return C2;
 }
 function transpose(A) {
   const n = A.length;
@@ -13045,6 +13047,24 @@ function solvePpp(epochs, sp3, opts) {
   const clockSeeded = /* @__PURE__ */ new Set();
   const ambIdx = /* @__PURE__ */ new Map();
   const lastSeen = /* @__PURE__ */ new Map();
+  const arcs = [];
+  const arcSnap = /* @__PURE__ */ new Map();
+  const flushArc = (prn, minEpochs = 10) => {
+    const s = arcSnap.get(prn);
+    arcSnap.delete(prn);
+    if (!s || s.nEpochs < minEpochs) return;
+    arcs.push({
+      prn,
+      aIF: s.aIF,
+      mwCyc: s.mwCyc,
+      f1: s.f1,
+      f2: s.f2,
+      meanElevDeg: s.elevSum / s.nEpochs,
+      nEpochs: s.nEpochs,
+      startMs: s.startMs,
+      endMs: s.endMs
+    });
+  };
   const mwState = /* @__PURE__ */ new Map();
   const detectSlip = (o, ei) => {
     const lamW = C_LIGHT / (o.f1 - o.f2);
@@ -13162,7 +13182,9 @@ function solvePpp(epochs, sp3, opts) {
         tropoFixed: zhd * mh + corr.rangeM,
         mw,
         windupM: corr.phaseWindupM,
-        slip
+        slip,
+        f1: o.f1,
+        f2: o.f2
       });
     }
     const zwd0 = x[IDX_ZWD];
@@ -13187,6 +13209,7 @@ function solvePpp(epochs, sp3, opts) {
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
+      if (v.slip) flushArc(v.prn);
       let ai = ambIdx.get(v.prn);
       if (ai === void 0 || v.slip) {
         const bInit = v.lIf - (common + v.windupM);
@@ -13220,6 +13243,7 @@ function solvePpp(epochs, sp3, opts) {
         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) {
+          flushArc(v.prn);
           x[ai] = v.lIf - (commonPh + v.windupM);
           for (let j = 0; j < P.length; j++) {
             P[ai][j] = 0;
@@ -13237,6 +13261,29 @@ function solvePpp(epochs, sp3, opts) {
           ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
+      {
+        const elevDeg = v.elRad * 180 / Math.PI;
+        const mwMean = mwState.get(v.prn)?.mean ?? 0;
+        const s = arcSnap.get(v.prn);
+        if (s) {
+          s.aIF = x[ai];
+          s.mwCyc = mwMean;
+          s.elevSum += elevDeg;
+          s.nEpochs++;
+          s.endMs = epoch.timeMs;
+        } else {
+          arcSnap.set(v.prn, {
+            aIF: x[ai],
+            mwCyc: mwMean,
+            f1: v.f1,
+            f2: v.f2,
+            elevSum: elevDeg,
+            nEpochs: 1,
+            startMs: epoch.timeMs,
+            endMs: epoch.timeMs
+          });
+        }
+      }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
@@ -13270,6 +13317,7 @@ function solvePpp(epochs, sp3, opts) {
         const idx = ambIdx.get(prn);
         if (idx !== void 0) dropState(idx);
         lastSeen.delete(prn);
+        flushArc(prn);
       }
     }
     const pos = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
@@ -13309,6 +13357,7 @@ function solvePpp(epochs, sp3, opts) {
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
+  for (const prn of [...arcSnap.keys()]) flushArc(prn);
   const finalPos = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
   const llh = ecefToGeodetic(finalPos[0], finalPos[1], finalPos[2]);
   const finalError3d = truth ? Math.hypot(
@@ -13323,7 +13372,8 @@ function solvePpp(epochs, sp3, opts) {
     ztdWet: x[IDX_ZWD],
     epochsUsed: epochs.length,
     convergenceSec,
-    finalError3d
+    finalError3d,
+    arcs
   };
 }
 
@@ -13464,34 +13514,34 @@ function circMean(vals, weights) {
   if (s === 0 && c === 0) return 0;
   return wrapHalf(Math.atan2(s, c) / (2 * Math.PI));
 }
-function estimateWidelaneFcb(arcs, opts = {}) {
+function decomposeFcb(entries, 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 sats = [...new Set(entries.map((e) => e.sat))].sort();
+  const stations = [...new Set(entries.map((e) => e.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);
+  const nInt = new Array(entries.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];
+    for (let i = 0; i < entries.length; i++) {
+      const e = entries[i];
       nInt[i] = Math.round(
-        a.wlFloat - rcvBias.get(a.station) - satFcb.get(a.prn)
+        e.value - rcvBias.get(e.station) - satFcb.get(e.sat)
       );
     }
     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);
+      for (let i = 0; i < entries.length; i++) {
+        const e = entries[i];
+        if (e.sat !== s) continue;
+        vals.push(e.value - nInt[i] - rcvBias.get(e.station));
+        w.push(e.weight);
       }
       if (!vals.length) continue;
       const nv = circMean(vals, w);
@@ -13501,11 +13551,11 @@ function estimateWidelaneFcb(arcs, opts = {}) {
     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);
+      for (let i = 0; i < entries.length; i++) {
+        const e = entries[i];
+        if (e.station !== r) continue;
+        vals.push(e.value - nInt[i] - satFcb.get(e.sat));
+        w.push(e.weight);
       }
       if (!vals.length) continue;
       rcvBias.set(r, circMean(vals, w));
@@ -13520,18 +13570,88 @@ function estimateWidelaneFcb(arcs, opts = {}) {
   }
   let sumSq = 0;
   let fixed = 0;
-  for (let i = 0; i < arcs.length; i++) {
-    const a = arcs[i];
+  for (let i = 0; i < entries.length; i++) {
+    const e = entries[i];
     const resid = wrapHalf(
-      a.wlFloat - nInt[i] - rcvBias.get(a.station) - satFcb.get(a.prn)
+      e.value - nInt[i] - rcvBias.get(e.station) - satFcb.get(e.sat)
     );
     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;
+  const residRms = entries.length ? Math.sqrt(sumSq / entries.length) : 0;
+  const fixRate = entries.length ? fixed / entries.length : 0;
   return { satFcb, rcvBias, residRms, fixRate, refStation, iterations };
 }
+function estimateWidelaneFcb(arcs, opts = {}) {
+  return decomposeFcb(
+    arcs.map((a) => ({
+      station: a.station,
+      sat: a.prn,
+      value: a.wlFloat,
+      weight: a.nObs
+    })),
+    opts
+  );
+}
+var C = 299792458;
+function estimateNarrowlaneFcb(arcs, wl, opts = {}) {
+  const entries = [];
+  let wlRejected = 0;
+  for (const a of arcs) {
+    const satWl = wl.satFcb.get(a.prn);
+    const rcvWl = wl.rcvBias.get(a.station);
+    if (satWl === void 0 || rcvWl === void 0) continue;
+    const wlFloat = a.mwCyc - satWl - rcvWl;
+    const nWl = Math.round(wlFloat);
+    if (Math.abs(wlFloat - nWl) > (opts.fixWindow ?? 0.25)) {
+      wlRejected++;
+      continue;
+    }
+    const lamNl = C / (a.f1 + a.f2);
+    const n1Float = a.aIF / lamNl - a.f2 / (a.f1 - a.f2) * nWl;
+    entries.push({
+      station: a.station,
+      sat: a.prn,
+      value: n1Float,
+      weight: a.nEpochs
+    });
+  }
+  const res = decomposeFcb(entries, { fixWindow: 0.15, ...opts });
+  return { ...res, usedArcs: entries.length, wlRejected };
+}
+function estimateNetworkFcbs(arcs, opts = {}) {
+  const minArcEpochs = opts.minArcEpochs ?? 120;
+  const minWlObs = opts.minWlObs ?? 40;
+  const fixWindow = opts.fixWindow ?? 0.15;
+  const satWlFcb = /* @__PURE__ */ new Map();
+  const satNlFcb = /* @__PURE__ */ new Map();
+  const perSystem = {};
+  const systems = [...new Set(arcs.map((a) => a.prn[0]))].sort();
+  for (const sys of systems) {
+    const sysArcs = arcs.filter((a) => a.prn[0] === sys);
+    const wlArcs = sysArcs.filter((a) => a.nEpochs >= minWlObs).map((a) => ({
+      station: a.station,
+      prn: a.prn,
+      wlFloat: a.mwCyc,
+      nObs: a.nEpochs
+    }));
+    if (wlArcs.length < 3) continue;
+    const wl = estimateWidelaneFcb(wlArcs, { fixWindow });
+    const nlInput = sysArcs.filter((a) => a.nEpochs >= minArcEpochs);
+    const nl = estimateNarrowlaneFcb(nlInput, wl, { fixWindow });
+    for (const [prn, v] of wl.satFcb) satWlFcb.set(prn, v);
+    for (const [prn, v] of nl.satFcb) satNlFcb.set(prn, v);
+    perSystem[sys] = {
+      wlFixRate: wl.fixRate,
+      wlResidRms: wl.residRms,
+      nlUsedArcs: nl.usedArcs,
+      nlFixRate: nl.fixRate,
+      nlResidRms: nl.residRms,
+      wlRejected: nl.wlRejected
+    };
+  }
+  return { satWlFcb, satNlFcb, perSystem };
+}
 
 // src/positioning/index.ts
 var GM_GPS3 = 3986005e8;
@@ -13943,9 +14063,9 @@ var MultipathAccumulator = class {
     const bandMap = this.obsIndices.get(sys) ?? this.obsIndices.get("_v2");
     if (!bandMap) return;
     const bandData = /* @__PURE__ */ new Map();
-    for (const [band, { C, L }] of bandMap) {
-      if (C === null) continue;
-      const cVal = values[C];
+    for (const [band, { C: C2, L }] of bandMap) {
+      if (C2 === null) continue;
+      const cVal = values[C2];
       const lVal = values[L];
       const freq = getFreq(this.gloChannels, prn, band);
       if (cVal != null && cVal !== 0 && lVal != null && lVal !== 0 && freq) {
@@ -14157,13 +14277,13 @@ var CycleSlipAccumulator = class {
     if (sys !== "R" && !FREQ[sys]) return;
     const phaseM = /* @__PURE__ */ new Map();
     const codeM = /* @__PURE__ */ new Map();
-    for (const [band, { L, C }] of idxMap) {
+    for (const [band, { L, C: C2 }] of idxMap) {
       const freq = this.getFreqForPrn(prn, band);
       if (!freq) continue;
       const lVal = values[L];
       if (lVal != null && lVal !== 0) phaseM.set(band, lVal * (C_LIGHT / freq));
-      if (C !== null) {
-        const cVal = values[C];
+      if (C2 !== null) {
+        const cVal = values[C2];
         if (cVal != null && cVal !== 0) codeM.set(band, cVal);
       }
     }
@@ -14427,9 +14547,9 @@ var IonoAccumulator = class {
     const bandMap = this.obsIndices.get(sys) ?? this.obsIndices.get("_v2");
     if (!bandMap) return;
     const bandData = /* @__PURE__ */ new Map();
-    for (const [band, { C, L }] of bandMap) {
-      if (C === null) continue;
-      const cVal = values[C];
+    for (const [band, { C: C2, L }] of bandMap) {
+      if (C2 === null) continue;
+      const cVal = values[C2];
       const lVal = values[L];
       const freq = getFreq(this.gloChannels, prn, band);
       if (cVal != null && cVal !== 0 && lVal != null && lVal !== 0 && freq) {
@@ -14437,7 +14557,7 @@ var IonoAccumulator = class {
           C: cVal,
           L: lVal,
           f: freq,
-          code: codes[C] ?? `C${band}`
+          code: codes[C2] ?? `C${band}`
         });
       }
     }
@@ -16150,6 +16270,8 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   estimateWidelaneFcb,
   euclidean3D,
   extractWidelaneArcs,

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, 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 { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, 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, 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 { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, 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,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,
@@ -167,7 +169,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PVOF4IGX.js";
+} from "./chunk-BLGBMRX7.js";
 import {
   almanacEpochMs,
   almanacSatPosition,
@@ -474,6 +476,8 @@ export {
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   estimateWidelaneFcb,
   euclidean3D,
   extractWidelaneArcs,

package/dist/positioning.cjs

@@ -24,6 +24,8 @@ __export(positioning_exports, {
   RtkFloatEngine: () => RtkFloatEngine,
   buildPppAntenna: () => buildPppAntenna,
   createPppCorrections: () => createPppCorrections,
+  estimateNarrowlaneFcb: () => estimateNarrowlaneFcb,
+  estimateNetworkFcbs: () => estimateNetworkFcbs,
   estimateWidelaneFcb: () => estimateWidelaneFcb,
   extractWidelaneArcs: () => extractWidelaneArcs,
   gimSlantIonoDelayL1: () => gimSlantIonoDelayL1,
@@ -913,15 +915,15 @@ function matMul(A, B) {
   const n = A.length;
   const k = B.length;
   const m = B[0]?.length ?? 0;
-  const C = zeros(n, m);
+  const C2 = zeros(n, m);
   for (let i = 0; i < n; i++) {
     for (let l = 0; l < k; l++) {
       const a = A[i][l];
       if (a === 0) continue;
-      for (let j = 0; j < m; j++) C[i][j] += a * B[l][j];
+      for (let j = 0; j < m; j++) C2[i][j] += a * B[l][j];
     }
   }
-  return C;
+  return C2;
 }
 function transpose(A) {
   const n = A.length;
@@ -2173,6 +2175,24 @@ function solvePpp(epochs, sp3, opts) {
   const clockSeeded = /* @__PURE__ */ new Set();
   const ambIdx = /* @__PURE__ */ new Map();
   const lastSeen = /* @__PURE__ */ new Map();
+  const arcs = [];
+  const arcSnap = /* @__PURE__ */ new Map();
+  const flushArc = (prn, minEpochs = 10) => {
+    const s = arcSnap.get(prn);
+    arcSnap.delete(prn);
+    if (!s || s.nEpochs < minEpochs) return;
+    arcs.push({
+      prn,
+      aIF: s.aIF,
+      mwCyc: s.mwCyc,
+      f1: s.f1,
+      f2: s.f2,
+      meanElevDeg: s.elevSum / s.nEpochs,
+      nEpochs: s.nEpochs,
+      startMs: s.startMs,
+      endMs: s.endMs
+    });
+  };
   const mwState = /* @__PURE__ */ new Map();
   const detectSlip = (o, ei) => {
     const lamW = C_LIGHT / (o.f1 - o.f2);
@@ -2290,7 +2310,9 @@ function solvePpp(epochs, sp3, opts) {
         tropoFixed: zhd * mh + corr.rangeM,
         mw,
         windupM: corr.phaseWindupM,
-        slip
+        slip,
+        f1: o.f1,
+        f2: o.f2
       });
     }
     const zwd0 = x[IDX_ZWD];
@@ -2315,6 +2337,7 @@ function solvePpp(epochs, sp3, opts) {
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
+      if (v.slip) flushArc(v.prn);
       let ai = ambIdx.get(v.prn);
       if (ai === void 0 || v.slip) {
         const bInit = v.lIf - (common + v.windupM);
@@ -2348,6 +2371,7 @@ function solvePpp(epochs, sp3, opts) {
         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) {
+          flushArc(v.prn);
           x[ai] = v.lIf - (commonPh + v.windupM);
           for (let j = 0; j < P.length; j++) {
             P[ai][j] = 0;
@@ -2365,6 +2389,29 @@ function solvePpp(epochs, sp3, opts) {
           ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
+      {
+        const elevDeg = v.elRad * 180 / Math.PI;
+        const mwMean = mwState.get(v.prn)?.mean ?? 0;
+        const s = arcSnap.get(v.prn);
+        if (s) {
+          s.aIF = x[ai];
+          s.mwCyc = mwMean;
+          s.elevSum += elevDeg;
+          s.nEpochs++;
+          s.endMs = epoch.timeMs;
+        } else {
+          arcSnap.set(v.prn, {
+            aIF: x[ai],
+            mwCyc: mwMean,
+            f1: v.f1,
+            f2: v.f2,
+            elevSum: elevDeg,
+            nEpochs: 1,
+            startMs: epoch.timeMs,
+            endMs: epoch.timeMs
+          });
+        }
+      }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
@@ -2398,6 +2445,7 @@ function solvePpp(epochs, sp3, opts) {
         const idx = ambIdx.get(prn);
         if (idx !== void 0) dropState(idx);
         lastSeen.delete(prn);
+        flushArc(prn);
       }
     }
     const pos = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
@@ -2437,6 +2485,7 @@ function solvePpp(epochs, sp3, opts) {
       phaseResRms: nPhaseRes > 0 ? Math.sqrt(sumPhaseResSq / nPhaseRes) : 0
     });
   }
+  for (const prn of [...arcSnap.keys()]) flushArc(prn);
   const finalPos = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
   const llh = ecefToGeodetic(finalPos[0], finalPos[1], finalPos[2]);
   const finalError3d = truth ? Math.hypot(
@@ -2451,7 +2500,8 @@ function solvePpp(epochs, sp3, opts) {
     ztdWet: x[IDX_ZWD],
     epochsUsed: epochs.length,
     convergenceSec,
-    finalError3d
+    finalError3d,
+    arcs
   };
 }
 
@@ -2592,34 +2642,34 @@ function circMean(vals, weights) {
   if (s === 0 && c === 0) return 0;
   return wrapHalf(Math.atan2(s, c) / (2 * Math.PI));
 }
-function estimateWidelaneFcb(arcs, opts = {}) {
+function decomposeFcb(entries, 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 sats = [...new Set(entries.map((e) => e.sat))].sort();
+  const stations = [...new Set(entries.map((e) => e.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);
+  const nInt = new Array(entries.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];
+    for (let i = 0; i < entries.length; i++) {
+      const e = entries[i];
       nInt[i] = Math.round(
-        a.wlFloat - rcvBias.get(a.station) - satFcb.get(a.prn)
+        e.value - rcvBias.get(e.station) - satFcb.get(e.sat)
       );
     }
     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);
+      for (let i = 0; i < entries.length; i++) {
+        const e = entries[i];
+        if (e.sat !== s) continue;
+        vals.push(e.value - nInt[i] - rcvBias.get(e.station));
+        w.push(e.weight);
       }
       if (!vals.length) continue;
       const nv = circMean(vals, w);
@@ -2629,11 +2679,11 @@ function estimateWidelaneFcb(arcs, opts = {}) {
     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);
+      for (let i = 0; i < entries.length; i++) {
+        const e = entries[i];
+        if (e.station !== r) continue;
+        vals.push(e.value - nInt[i] - satFcb.get(e.sat));
+        w.push(e.weight);
       }
       if (!vals.length) continue;
       rcvBias.set(r, circMean(vals, w));
@@ -2648,18 +2698,88 @@ function estimateWidelaneFcb(arcs, opts = {}) {
   }
   let sumSq = 0;
   let fixed = 0;
-  for (let i = 0; i < arcs.length; i++) {
-    const a = arcs[i];
+  for (let i = 0; i < entries.length; i++) {
+    const e = entries[i];
     const resid = wrapHalf(
-      a.wlFloat - nInt[i] - rcvBias.get(a.station) - satFcb.get(a.prn)
+      e.value - nInt[i] - rcvBias.get(e.station) - satFcb.get(e.sat)
     );
     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;
+  const residRms = entries.length ? Math.sqrt(sumSq / entries.length) : 0;
+  const fixRate = entries.length ? fixed / entries.length : 0;
   return { satFcb, rcvBias, residRms, fixRate, refStation, iterations };
 }
+function estimateWidelaneFcb(arcs, opts = {}) {
+  return decomposeFcb(
+    arcs.map((a) => ({
+      station: a.station,
+      sat: a.prn,
+      value: a.wlFloat,
+      weight: a.nObs
+    })),
+    opts
+  );
+}
+var C = 299792458;
+function estimateNarrowlaneFcb(arcs, wl, opts = {}) {
+  const entries = [];
+  let wlRejected = 0;
+  for (const a of arcs) {
+    const satWl = wl.satFcb.get(a.prn);
+    const rcvWl = wl.rcvBias.get(a.station);
+    if (satWl === void 0 || rcvWl === void 0) continue;
+    const wlFloat = a.mwCyc - satWl - rcvWl;
+    const nWl = Math.round(wlFloat);
+    if (Math.abs(wlFloat - nWl) > (opts.fixWindow ?? 0.25)) {
+      wlRejected++;
+      continue;
+    }
+    const lamNl = C / (a.f1 + a.f2);
+    const n1Float = a.aIF / lamNl - a.f2 / (a.f1 - a.f2) * nWl;
+    entries.push({
+      station: a.station,
+      sat: a.prn,
+      value: n1Float,
+      weight: a.nEpochs
+    });
+  }
+  const res = decomposeFcb(entries, { fixWindow: 0.15, ...opts });
+  return { ...res, usedArcs: entries.length, wlRejected };
+}
+function estimateNetworkFcbs(arcs, opts = {}) {
+  const minArcEpochs = opts.minArcEpochs ?? 120;
+  const minWlObs = opts.minWlObs ?? 40;
+  const fixWindow = opts.fixWindow ?? 0.15;
+  const satWlFcb = /* @__PURE__ */ new Map();
+  const satNlFcb = /* @__PURE__ */ new Map();
+  const perSystem = {};
+  const systems = [...new Set(arcs.map((a) => a.prn[0]))].sort();
+  for (const sys of systems) {
+    const sysArcs = arcs.filter((a) => a.prn[0] === sys);
+    const wlArcs = sysArcs.filter((a) => a.nEpochs >= minWlObs).map((a) => ({
+      station: a.station,
+      prn: a.prn,
+      wlFloat: a.mwCyc,
+      nObs: a.nEpochs
+    }));
+    if (wlArcs.length < 3) continue;
+    const wl = estimateWidelaneFcb(wlArcs, { fixWindow });
+    const nlInput = sysArcs.filter((a) => a.nEpochs >= minArcEpochs);
+    const nl = estimateNarrowlaneFcb(nlInput, wl, { fixWindow });
+    for (const [prn, v] of wl.satFcb) satWlFcb.set(prn, v);
+    for (const [prn, v] of nl.satFcb) satNlFcb.set(prn, v);
+    perSystem[sys] = {
+      wlFixRate: wl.fixRate,
+      wlResidRms: wl.residRms,
+      nlUsedArcs: nl.usedArcs,
+      nlFixRate: nl.fixRate,
+      nlResidRms: nl.residRms,
+      wlRejected: nl.wlRejected
+    };
+  }
+  return { satWlFcb, satNlFcb, perSystem };
+}
 
 // src/positioning/index.ts
 var GM_GPS2 = 3986005e8;
@@ -2894,6 +3014,8 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,

package/dist/positioning.d.cts

@@ -666,6 +666,26 @@ interface PppSolution {
     convergenceSec: number | null;
     /** Final 3D error vs ground truth (m), if given. */
     finalError3d: number | null;
+    /** Per-arc converged float ambiguities — one per continuous satellite
+     * tracking arc — the raw material for ambiguity resolution (PPP-AR). */
+    arcs: PppArc[];
+}
+/** A converged carrier-phase arc: the float ionosphere-free ambiguity plus
+ * the wide-lane, for integer resolution / FCB estimation. */
+interface PppArc {
+    prn: string;
+    /** Final float ionosphere-free ambiguity (m). */
+    aIF: number;
+    /** Arc-averaged Melbourne–Wübbena wide-lane (cycles). */
+    mwCyc: number;
+    f1: number;
+    f2: number;
+    /** Mean satellite elevation over the arc (deg). */
+    meanElevDeg: number;
+    /** Epochs in the arc. */
+    nEpochs: number;
+    startMs: number;
+    endMs: number;
 }
 declare function solvePpp(epochs: PppEpoch[], sp3: Sp3File, opts: PppOptions): PppSolution;
 
@@ -879,6 +899,73 @@ interface WlFcbResult {
  * integer/bias entanglement.
  */
 declare function estimateWidelaneFcb(arcs: WlArc[], opts?: WlFcbOptions): WlFcbResult;
+/** One converged arc (from `solvePpp`) tagged with its station, for NL FCB. */
+interface NlArc {
+    station: string;
+    prn: string;
+    /** Float ionosphere-free ambiguity (m). */
+    aIF: number;
+    /** Arc-averaged Melbourne–Wübbena wide-lane (cycles). */
+    mwCyc: number;
+    f1: number;
+    f2: number;
+    nEpochs: number;
+}
+interface NlFcbResult extends WlFcbResult {
+    /** Arcs whose wide-lane fixed and were used in the narrow-lane estimation. */
+    usedArcs: number;
+    /** Arcs dropped because the wide-lane did not round confidently. */
+    wlRejected: number;
+}
+/**
+ * Estimate per-satellite **narrow-lane** FCBs — the centimetre half of
+ * PPP-AR. Needs the wide-lane FCBs already solved (`wl`, from the same
+ * network) to fix each arc's N_WL, from which the float narrow-lane
+ * ambiguity is recovered:
+ *
+ *   A_IF = λ_NL·N1 + λ_NL·(f2/(f1−f2))·N_WL
+ *   ⇒ N1_float = A_IF/λ_NL − (f2/(f1−f2))·N_WL,   λ_NL = c/(f1+f2)
+ *
+ * then decomposed across the network exactly like the wide-lane. The
+ * narrow-lane wavelength is ~10.7 cm, so this is far more sensitive to
+ * residual orbit/clock/troposphere/position error than the ~86 cm wide-lane
+ * — the fix rate is the honest measure of how close the network is to cm.
+ */
+declare function estimateNarrowlaneFcb(arcs: NlArc[], wl: WlFcbResult, opts?: WlFcbOptions): NlFcbResult;
+interface NetworkFcbSummary {
+    wlFixRate: number;
+    wlResidRms: number;
+    nlUsedArcs: number;
+    nlFixRate: number;
+    nlResidRms: number;
+    wlRejected: number;
+}
+interface NetworkFcbResult {
+    /** Per-satellite wide-lane FCB (cycles), all constellations merged. */
+    satWlFcb: Map<string, number>;
+    /** Per-satellite narrow-lane FCB (cycles), all constellations merged. */
+    satNlFcb: Map<string, number>;
+    /** Per-constellation fix-rate / residual summary, keyed by system letter. */
+    perSystem: Record<string, NetworkFcbSummary>;
+}
+/**
+ * Full network FCB calibration for PPP-AR from a set of `solvePpp` arcs (each
+ * tagged with its station), **per constellation** — the receiver wide- and
+ * narrow-lane biases are system-specific, so lumping GPS/Galileo/BeiDou under
+ * one datum corrupts the fit. Runs the wide-lane then narrow-lane
+ * decomposition for each system and merges the per-satellite biases.
+ *
+ * The satellite FCBs (`satWlFcb`/`satNlFcb`) are the reusable product: a rover
+ * applies them and solves only its own receiver bias. Multi-GNSS input
+ * multiplies the fixable arcs (Galileo roughly doubles GPS-only) without
+ * changing per-arc behaviour — each satellite's arc length is independent of
+ * how many constellations are present.
+ */
+declare function estimateNetworkFcbs(arcs: NlArc[], opts?: {
+    minArcEpochs?: number;
+    minWlObs?: number;
+    fixWindow?: number;
+}): NetworkFcbResult;
 
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
@@ -969,4 +1056,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, 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 };
+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 NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, 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, estimateNarrowlaneFcb, estimateNetworkFcbs, 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

@@ -666,6 +666,26 @@ interface PppSolution {
     convergenceSec: number | null;
     /** Final 3D error vs ground truth (m), if given. */
     finalError3d: number | null;
+    /** Per-arc converged float ambiguities — one per continuous satellite
+     * tracking arc — the raw material for ambiguity resolution (PPP-AR). */
+    arcs: PppArc[];
+}
+/** A converged carrier-phase arc: the float ionosphere-free ambiguity plus
+ * the wide-lane, for integer resolution / FCB estimation. */
+interface PppArc {
+    prn: string;
+    /** Final float ionosphere-free ambiguity (m). */
+    aIF: number;
+    /** Arc-averaged Melbourne–Wübbena wide-lane (cycles). */
+    mwCyc: number;
+    f1: number;
+    f2: number;
+    /** Mean satellite elevation over the arc (deg). */
+    meanElevDeg: number;
+    /** Epochs in the arc. */
+    nEpochs: number;
+    startMs: number;
+    endMs: number;
 }
 declare function solvePpp(epochs: PppEpoch[], sp3: Sp3File, opts: PppOptions): PppSolution;
 
@@ -879,6 +899,73 @@ interface WlFcbResult {
  * integer/bias entanglement.
  */
 declare function estimateWidelaneFcb(arcs: WlArc[], opts?: WlFcbOptions): WlFcbResult;
+/** One converged arc (from `solvePpp`) tagged with its station, for NL FCB. */
+interface NlArc {
+    station: string;
+    prn: string;
+    /** Float ionosphere-free ambiguity (m). */
+    aIF: number;
+    /** Arc-averaged Melbourne–Wübbena wide-lane (cycles). */
+    mwCyc: number;
+    f1: number;
+    f2: number;
+    nEpochs: number;
+}
+interface NlFcbResult extends WlFcbResult {
+    /** Arcs whose wide-lane fixed and were used in the narrow-lane estimation. */
+    usedArcs: number;
+    /** Arcs dropped because the wide-lane did not round confidently. */
+    wlRejected: number;
+}
+/**
+ * Estimate per-satellite **narrow-lane** FCBs — the centimetre half of
+ * PPP-AR. Needs the wide-lane FCBs already solved (`wl`, from the same
+ * network) to fix each arc's N_WL, from which the float narrow-lane
+ * ambiguity is recovered:
+ *
+ *   A_IF = λ_NL·N1 + λ_NL·(f2/(f1−f2))·N_WL
+ *   ⇒ N1_float = A_IF/λ_NL − (f2/(f1−f2))·N_WL,   λ_NL = c/(f1+f2)
+ *
+ * then decomposed across the network exactly like the wide-lane. The
+ * narrow-lane wavelength is ~10.7 cm, so this is far more sensitive to
+ * residual orbit/clock/troposphere/position error than the ~86 cm wide-lane
+ * — the fix rate is the honest measure of how close the network is to cm.
+ */
+declare function estimateNarrowlaneFcb(arcs: NlArc[], wl: WlFcbResult, opts?: WlFcbOptions): NlFcbResult;
+interface NetworkFcbSummary {
+    wlFixRate: number;
+    wlResidRms: number;
+    nlUsedArcs: number;
+    nlFixRate: number;
+    nlResidRms: number;
+    wlRejected: number;
+}
+interface NetworkFcbResult {
+    /** Per-satellite wide-lane FCB (cycles), all constellations merged. */
+    satWlFcb: Map<string, number>;
+    /** Per-satellite narrow-lane FCB (cycles), all constellations merged. */
+    satNlFcb: Map<string, number>;
+    /** Per-constellation fix-rate / residual summary, keyed by system letter. */
+    perSystem: Record<string, NetworkFcbSummary>;
+}
+/**
+ * Full network FCB calibration for PPP-AR from a set of `solvePpp` arcs (each
+ * tagged with its station), **per constellation** — the receiver wide- and
+ * narrow-lane biases are system-specific, so lumping GPS/Galileo/BeiDou under
+ * one datum corrupts the fit. Runs the wide-lane then narrow-lane
+ * decomposition for each system and merges the per-satellite biases.
+ *
+ * The satellite FCBs (`satWlFcb`/`satNlFcb`) are the reusable product: a rover
+ * applies them and solves only its own receiver bias. Multi-GNSS input
+ * multiplies the fixable arcs (Galileo roughly doubles GPS-only) without
+ * changing per-arc behaviour — each satellite's arc length is independent of
+ * how many constellations are present.
+ */
+declare function estimateNetworkFcbs(arcs: NlArc[], opts?: {
+    minArcEpochs?: number;
+    minWlObs?: number;
+    fixWindow?: number;
+}): NetworkFcbResult;
 
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
@@ -969,4 +1056,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, 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 };
+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 NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, 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, estimateNarrowlaneFcb, estimateNetworkFcbs, 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,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,
@@ -24,7 +26,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PVOF4IGX.js";
+} from "./chunk-BLGBMRX7.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
@@ -37,6 +39,8 @@ export {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateNarrowlaneFcb,
+  estimateNetworkFcbs,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,

package/package.json

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