gnss-js 1.33.0 → 1.34.0

package/dist/{chunk-PVOF4IGX.js → chunk-H73PZB5C.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,55 @@ 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 };
+}
 
 // src/positioning/index.ts
 var GM_GPS = 3986005e8;
@@ -2388,6 +2473,7 @@ export {
   resolvePppAmbiguities,
   extractWidelaneArcs,
   estimateWidelaneFcb,
+  estimateNarrowlaneFcb,
   satClockCorrection,
   ionoFree,
   solveSpp

package/dist/index.cjs

@@ -149,6 +149,7 @@ __export(src_exports, {
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
   ephInfoToEphemeris: () => ephInfoToEphemeris,
+  estimateNarrowlaneFcb: () => estimateNarrowlaneFcb,
   estimateWidelaneFcb: () => estimateWidelaneFcb,
   euclidean3D: () => euclidean3D,
   extractWidelaneArcs: () => extractWidelaneArcs,
@@ -1067,7 +1068,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 +1084,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 +1190,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 +11858,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 +13046,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 +13181,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 +13208,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 +13242,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 +13260,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 +13316,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 +13356,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 +13371,8 @@ function solvePpp(epochs, sp3, opts) {
     ztdWet: x[IDX_ZWD],
     epochsUsed: epochs.length,
     convergenceSec,
-    finalError3d
+    finalError3d,
+    arcs
   };
 }
 
@@ -13464,34 +13513,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 +13550,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 +13569,55 @@ 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 };
+}
 
 // src/positioning/index.ts
 var GM_GPS3 = 3986005e8;
@@ -13943,9 +14029,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 +14243,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 +14513,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 +14523,7 @@ var IonoAccumulator = class {
           C: cVal,
           L: lVal,
           f: freq,
-          code: codes[C] ?? `C${band}`
+          code: codes[C2] ?? `C${band}`
         });
       }
     }
@@ -16150,6 +16236,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateNarrowlaneFcb,
   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, 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, 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, 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, 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,7 @@ import {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateNarrowlaneFcb,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,
@@ -167,7 +168,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PVOF4IGX.js";
+} from "./chunk-H73PZB5C.js";
 import {
   almanacEpochMs,
   almanacSatPosition,
@@ -474,6 +475,7 @@ export {
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateNarrowlaneFcb,
   estimateWidelaneFcb,
   euclidean3D,
   extractWidelaneArcs,

package/dist/positioning.cjs

@@ -24,6 +24,7 @@ __export(positioning_exports, {
   RtkFloatEngine: () => RtkFloatEngine,
   buildPppAntenna: () => buildPppAntenna,
   createPppCorrections: () => createPppCorrections,
+  estimateNarrowlaneFcb: () => estimateNarrowlaneFcb,
   estimateWidelaneFcb: () => estimateWidelaneFcb,
   extractWidelaneArcs: () => extractWidelaneArcs,
   gimSlantIonoDelayL1: () => gimSlantIonoDelayL1,
@@ -913,15 +914,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 +2174,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 +2309,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 +2336,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 +2370,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 +2388,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 +2444,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 +2484,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 +2499,8 @@ function solvePpp(epochs, sp3, opts) {
     ztdWet: x[IDX_ZWD],
     epochsUsed: epochs.length,
     convergenceSec,
-    finalError3d
+    finalError3d,
+    arcs
   };
 }
 
@@ -2592,34 +2641,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 +2678,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 +2697,55 @@ 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 };
+}
 
 // src/positioning/index.ts
 var GM_GPS2 = 3986005e8;
@@ -2894,6 +2980,7 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateNarrowlaneFcb,
   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,39 @@ 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;
 
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
@@ -969,4 +1022,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 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, 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,39 @@ 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;
 
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
@@ -969,4 +1022,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 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, 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,7 @@ import {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateNarrowlaneFcb,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,
@@ -24,7 +25,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PVOF4IGX.js";
+} from "./chunk-H73PZB5C.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
@@ -37,6 +38,7 @@ export {
   RtkFloatEngine,
   buildPppAntenna,
   createPppCorrections,
+  estimateNarrowlaneFcb,
   estimateWidelaneFcb,
   extractWidelaneArcs,
   gimSlantIonoDelayL1,

package/package.json

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