gnss-js 1.30.0 → 1.31.0

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

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

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

@@ -1911,6 +1911,96 @@ function solvePpp(epochs, sp3, opts) {
   };
 }
 
+// src/positioning/ppp-ar.ts
+function wlWavelength(f1, f2) {
+  return C_LIGHT / (f1 - f2);
+}
+function nlWavelength(f1, f2) {
+  return C_LIGHT / (f1 + f2);
+}
+function resolvePppAmbiguities(sats, Qaif, opts = {}) {
+  const wlThreshold = opts.wlThreshold ?? 0.25;
+  const ratioThreshold = opts.ratioThreshold ?? 3;
+  const kept = [];
+  const nWl = [];
+  const rejectedWl = [];
+  for (let i = 0; i < sats.length; i++) {
+    const s = sats[i];
+    const wlFloat = s.mwCyc - (s.wlBiasCyc ?? 0);
+    const r = Math.round(wlFloat);
+    if (Math.abs(wlFloat - r) <= wlThreshold) {
+      kept.push(i);
+      nWl.push(r);
+    } else {
+      rejectedWl.push(s.prn);
+    }
+  }
+  if (kept.length < 2) {
+    return {
+      fixed: false,
+      ratio: 0,
+      refPrn: null,
+      sats: [],
+      rejectedWl
+    };
+  }
+  const n1Float = [];
+  const lamNl = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    const lNl = nlWavelength(s.f1, s.f2);
+    lamNl.push(lNl);
+    const n1 = s.aIF / lNl - s.f2 / (s.f1 - s.f2) * nWl[k] - (s.nlBiasCyc ?? 0);
+    n1Float.push(n1);
+  }
+  let ref = 0;
+  for (let k = 1; k < kept.length; k++) {
+    const e = sats[kept[k]].elevDeg ?? -1;
+    if (e > (sats[kept[ref]].elevDeg ?? -1)) ref = k;
+  }
+  const m = kept.length - 1;
+  const sd = new Float64Array(m);
+  const Qsd = new Float64Array(m * m);
+  const others = [];
+  for (let k = 0; k < kept.length; k++) if (k !== ref) others.push(k);
+  const n = sats.length;
+  const qN1 = (ka, kb) => Qaif[kept[ka] * n + kept[kb]] / (lamNl[ka] * lamNl[kb]);
+  for (let a = 0; a < m; a++) {
+    const ka = others[a];
+    sd[a] = n1Float[ka] - n1Float[ref];
+    for (let b = 0; b < m; b++) {
+      const kb = others[b];
+      Qsd[a * m + b] = qN1(ka, kb) - qN1(ka, ref) - qN1(ref, kb) + qN1(ref, ref);
+    }
+  }
+  const res = lambdaSearch(sd, Qsd, m, 2);
+  const ratio = res ? Math.min(res.ratio, 999.9) : 0;
+  const fixed = !!res && ratio >= ratioThreshold;
+  const refN1 = Math.round(n1Float[ref]);
+  const sdFix = res?.candidates[0];
+  const outSats = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    let n1;
+    if (k === ref) {
+      n1 = refN1;
+    } else {
+      const a = others.indexOf(k);
+      n1 = fixed && sdFix ? refN1 + Math.round(sdFix[a]) : Math.round(n1Float[k]);
+    }
+    const lNl = lamNl[k];
+    const aIF = lNl * n1 + lNl * (s.f2 / (s.f1 - s.f2)) * nWl[k];
+    outSats.push({ prn: s.prn, nWl: nWl[k], n1, aIF });
+  }
+  return {
+    fixed,
+    ratio,
+    refPrn: sats[kept[ref]].prn,
+    sats: outSats,
+    rejectedWl
+  };
+}
+
 // src/positioning/index.ts
 var GM_GPS = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -2157,6 +2247,9 @@ export {
   buildPppAntenna,
   createPppCorrections,
   solvePpp,
+  wlWavelength,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   ionoFree,
   solveSpp

package/dist/index.cjs

@@ -111,6 +111,7 @@ __export(src_exports, {
   analyzeQuality: () => analyzeQuality,
   applyHelmert: () => applyHelmert,
   applyIonoDcb: () => applyIonoDcb,
+  biasMetres: () => biasMetres,
   binexCrc16: () => binexCrc16,
   binexCrc32: () => binexCrc32,
   binexCsum8: () => binexCsum8,
@@ -150,6 +151,7 @@ __export(src_exports, {
   ephInfoToEphemeris: () => ephInfoToEphemeris,
   euclidean3D: () => euclidean3D,
   fetchSourcetable: () => fetchSourcetable,
+  findSatBias: () => findSatBias,
   fmtD: () => fmtD,
   fmtF: () => fmtF,
   formatUTCTime: () => formatUTCTime,
@@ -232,11 +234,13 @@ __export(src_exports, {
   msmEpochToDate: () => msmEpochToDate,
   navTimesFromEph: () => navTimesFromEph,
   niellMapping: () => niellMapping,
+  nlWavelength: () => nlWavelength,
   nmeaCoordToDecimal: () => nmeaCoordToDecimal,
   oem4Frames: () => oem4Frames,
   padL: () => padL,
   padR: () => padR,
   parseAntex: () => parseAntex,
+  parseBiasSinex: () => parseBiasSinex,
   parseBinex: () => parseBinex,
   parseClk: () => parseClk,
   parseCrxDataLine: () => parseCrxDataLine,
@@ -274,6 +278,7 @@ __export(src_exports, {
   readString: () => readString,
   reportDecodeError: () => reportDecodeError,
   resetGloFreqCache: () => resetGloFreqCache,
+  resolvePppAmbiguities: () => resolvePppAmbiguities,
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   satClockCorrection: () => satClockCorrection,
@@ -299,6 +304,7 @@ __export(src_exports, {
   updateStreamStats: () => updateStreamStats,
   verifyChecksum: () => verifyChecksum,
   vincenty: () => vincenty,
+  wlWavelength: () => wlWavelength,
   writeCrx: () => writeCrx,
   writeModernObsBlob: () => writeModernObsBlob,
   writeRinex2ObsBlob: () => writeRinex2ObsBlob,
@@ -3958,6 +3964,81 @@ function clkBias(clk, prn, tMs) {
   return c[lo] + f * (c[lo + 1] - c[lo]);
 }
 
+// src/rinex/bias-sinex.ts
+function parseEpoch(s, openTo) {
+  const m = /^(\d{4}):(\d{3}):(\d{5})$/.exec(s.trim());
+  if (!m) return openTo;
+  const y = +m[1];
+  if (y === 0) return openTo;
+  const doy = +m[2];
+  const sod = +m[3];
+  return Date.UTC(y, 0, 1) + (doy - 1) * 864e5 + sod * 1e3;
+}
+function parseBiasSinex(text) {
+  const records = [];
+  let inBlock = false;
+  for (const line of text.split("\n")) {
+    if (line.startsWith("+BIAS/SOLUTION")) {
+      inBlock = true;
+      continue;
+    }
+    if (line.startsWith("-BIAS/SOLUTION")) {
+      inBlock = false;
+      continue;
+    }
+    if (!inBlock || line.length < 65 || line.startsWith("*")) continue;
+    const type = line.slice(1, 4).trim();
+    if (type !== "DSB" && type !== "OSB" && type !== "ISB") continue;
+    const prn = line.slice(11, 14).trim();
+    const station = line.slice(15, 24).trim();
+    const obs1 = line.slice(25, 29).trim();
+    const obs2 = line.slice(30, 34).trim();
+    const startMs = parseEpoch(line.slice(35, 49), -Infinity);
+    const endMs = parseEpoch(line.slice(50, 64), Infinity);
+    const tail = line.slice(65).trim().split(/\s+/);
+    if (tail.length < 2) continue;
+    const unit = tail[0];
+    const value = Number(tail[1]);
+    if (!Number.isFinite(value)) continue;
+    records.push({
+      type,
+      prn,
+      station,
+      obs1,
+      obs2,
+      startMs,
+      endMs,
+      value,
+      unit
+    });
+  }
+  const bySat = /* @__PURE__ */ new Map();
+  for (const r of records) {
+    if (!r.prn || r.station) continue;
+    let a = bySat.get(r.prn);
+    if (!a) {
+      a = [];
+      bySat.set(r.prn, a);
+    }
+    a.push(r);
+  }
+  return { records, bySat };
+}
+function biasMetres(r) {
+  return r.unit === "ns" ? r.value * 1e-9 * C_LIGHT : r.value;
+}
+function findSatBias(bias, prn, obs1, obs2, tMs) {
+  const recs = bias.bySat.get(prn);
+  if (!recs) return null;
+  for (const r of recs) {
+    if (r.obs1 !== obs1) continue;
+    if (obs2 === void 0 ? r.obs2 !== "" : r.obs2 !== obs2) continue;
+    if (tMs != null && (tMs < r.startMs || tMs > r.endMs)) continue;
+    return r;
+  }
+  return null;
+}
+
 // src/rinex/obs-writer-common.ts
 function createGzipLineSink() {
   const encoder = new TextEncoder();
@@ -13224,6 +13305,96 @@ function solvePpp(epochs, sp3, opts) {
   };
 }
 
+// src/positioning/ppp-ar.ts
+function wlWavelength(f1, f2) {
+  return C_LIGHT / (f1 - f2);
+}
+function nlWavelength(f1, f2) {
+  return C_LIGHT / (f1 + f2);
+}
+function resolvePppAmbiguities(sats, Qaif, opts = {}) {
+  const wlThreshold = opts.wlThreshold ?? 0.25;
+  const ratioThreshold = opts.ratioThreshold ?? 3;
+  const kept = [];
+  const nWl = [];
+  const rejectedWl = [];
+  for (let i = 0; i < sats.length; i++) {
+    const s = sats[i];
+    const wlFloat = s.mwCyc - (s.wlBiasCyc ?? 0);
+    const r = Math.round(wlFloat);
+    if (Math.abs(wlFloat - r) <= wlThreshold) {
+      kept.push(i);
+      nWl.push(r);
+    } else {
+      rejectedWl.push(s.prn);
+    }
+  }
+  if (kept.length < 2) {
+    return {
+      fixed: false,
+      ratio: 0,
+      refPrn: null,
+      sats: [],
+      rejectedWl
+    };
+  }
+  const n1Float = [];
+  const lamNl = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    const lNl = nlWavelength(s.f1, s.f2);
+    lamNl.push(lNl);
+    const n1 = s.aIF / lNl - s.f2 / (s.f1 - s.f2) * nWl[k] - (s.nlBiasCyc ?? 0);
+    n1Float.push(n1);
+  }
+  let ref = 0;
+  for (let k = 1; k < kept.length; k++) {
+    const e = sats[kept[k]].elevDeg ?? -1;
+    if (e > (sats[kept[ref]].elevDeg ?? -1)) ref = k;
+  }
+  const m = kept.length - 1;
+  const sd2 = new Float64Array(m);
+  const Qsd = new Float64Array(m * m);
+  const others = [];
+  for (let k = 0; k < kept.length; k++) if (k !== ref) others.push(k);
+  const n = sats.length;
+  const qN1 = (ka, kb) => Qaif[kept[ka] * n + kept[kb]] / (lamNl[ka] * lamNl[kb]);
+  for (let a = 0; a < m; a++) {
+    const ka = others[a];
+    sd2[a] = n1Float[ka] - n1Float[ref];
+    for (let b = 0; b < m; b++) {
+      const kb = others[b];
+      Qsd[a * m + b] = qN1(ka, kb) - qN1(ka, ref) - qN1(ref, kb) + qN1(ref, ref);
+    }
+  }
+  const res = lambdaSearch(sd2, Qsd, m, 2);
+  const ratio = res ? Math.min(res.ratio, 999.9) : 0;
+  const fixed = !!res && ratio >= ratioThreshold;
+  const refN1 = Math.round(n1Float[ref]);
+  const sdFix = res?.candidates[0];
+  const outSats = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    let n1;
+    if (k === ref) {
+      n1 = refN1;
+    } else {
+      const a = others.indexOf(k);
+      n1 = fixed && sdFix ? refN1 + Math.round(sdFix[a]) : Math.round(n1Float[k]);
+    }
+    const lNl = lamNl[k];
+    const aIF = lNl * n1 + lNl * (s.f2 / (s.f1 - s.f2)) * nWl[k];
+    outSats.push({ prn: s.prn, nWl: nWl[k], n1, aIF });
+  }
+  return {
+    fixed,
+    ratio,
+    refPrn: sats[kept[ref]].prn,
+    sats: outSats,
+    rejectedWl
+  };
+}
+
 // src/positioning/index.ts
 var GM_GPS3 = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -15803,6 +15974,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   analyzeQuality,
   applyHelmert,
   applyIonoDcb,
+  biasMetres,
   binexCrc16,
   binexCrc32,
   binexCsum8,
@@ -15842,6 +16014,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ephInfoToEphemeris,
   euclidean3D,
   fetchSourcetable,
+  findSatBias,
   fmtD,
   fmtF,
   formatUTCTime,
@@ -15924,11 +16097,13 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   msmEpochToDate,
   navTimesFromEph,
   niellMapping,
+  nlWavelength,
   nmeaCoordToDecimal,
   oem4Frames,
   padL,
   padR,
   parseAntex,
+  parseBiasSinex,
   parseBinex,
   parseClk,
   parseCrxDataLine,
@@ -15966,6 +16141,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolvePppAmbiguities,
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
@@ -15991,6 +16167,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  wlWavelength,
   writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,

package/dist/index.d.cts

@@ -4,7 +4,7 @@ export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSat
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.cjs';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.cjs';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.cjs';
-export { CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.cjs';
+export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.cjs';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.cjs';
 export { C as CnavEphemeris, E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, R as RinexCnavEphemeris, p as parseNavFile } from './nav-DImXUdbp.cjs';
 export { C as ClkFile, I as IonexGrid, S as Sp3File, a as Sp3Sample, c as clkBias, p as parseClk, b as parseIonex, d as parseSp3, s as sp3Position } from './clk-DGVBGRzA.cjs';
@@ -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, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch } from './positioning.cjs';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.cjs';
 export { 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

@@ -4,7 +4,7 @@ export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSat
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.js';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.js';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.js';
-export { CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.js';
+export { BiasRecord, BiasSinex, CompactEpoch, CrxField, DiffState, EMPTY_WARNINGS, Nav4Input, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob } from './rinex.js';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.js';
 export { C as CnavEphemeris, E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, R as RinexCnavEphemeris, p as parseNavFile } from './nav-DImXUdbp.js';
 export { C as ClkFile, I as IonexGrid, S as Sp3File, a as Sp3Sample, c as clkBias, p as parseClk, b as parseIonex, d as parseSp3, s as sp3Position } from './clk-DGVBGRzA.js';
@@ -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, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch } from './positioning.js';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, PppAntennaModel, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, SatGeom, SppOptions, SppSolution, TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.js';
 export { 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

@@ -155,14 +155,17 @@ import {
   lambdaSearch,
   moonEcef,
   niellMapping,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   solidEarthTide,
   solveDgnss,
   solvePpp,
   solveSpp,
   sunEcef,
-  toRtkEpoch
-} from "./chunk-CRFFUU5Y.js";
+  toRtkEpoch,
+  wlWavelength
+} from "./chunk-EVFD4UNZ.js";
 import {
   almanacEpochMs,
   almanacSatPosition,
@@ -201,12 +204,15 @@ import {
 import {
   EMPTY_WARNINGS,
   WarningAccumulator,
+  biasMetres,
   createGzipLineSink,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
   parseIonex,
   parseNavFile,
   stationHeaderLines,
@@ -217,7 +223,7 @@ import {
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-KDM2NVBO.js";
+} from "./chunk-5RNDQQ6J.js";
 import "./chunk-726LQBGM.js";
 import "./chunk-IQ5OIMQD.js";
 import {
@@ -428,6 +434,7 @@ export {
   analyzeQuality,
   applyHelmert,
   applyIonoDcb,
+  biasMetres,
   binexCrc16,
   binexCrc32,
   binexCsum8,
@@ -467,6 +474,7 @@ export {
   ephInfoToEphemeris,
   euclidean3D,
   fetchSourcetable,
+  findSatBias,
   fmtD,
   fmtF,
   formatUTCTime,
@@ -549,11 +557,13 @@ export {
   msmEpochToDate,
   navTimesFromEph,
   niellMapping,
+  nlWavelength,
   nmeaCoordToDecimal,
   oem4Frames,
   padL,
   padR,
   parseAntex,
+  parseBiasSinex,
   parseBinex,
   parseClk,
   parseCrxDataLine,
@@ -591,6 +601,7 @@ export {
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolvePppAmbiguities,
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
@@ -616,6 +627,7 @@ export {
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  wlWavelength,
   writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,

package/dist/positioning.cjs

@@ -33,13 +33,16 @@ __export(positioning_exports, {
   lambdaSearch: () => lambdaSearch,
   moonEcef: () => moonEcef,
   niellMapping: () => niellMapping,
+  nlWavelength: () => nlWavelength,
+  resolvePppAmbiguities: () => resolvePppAmbiguities,
   satClockCorrection: () => satClockCorrection,
   solidEarthTide: () => solidEarthTide,
   solveDgnss: () => solveDgnss,
   solvePpp: () => solvePpp,
   solveSpp: () => solveSpp,
   sunEcef: () => sunEcef,
-  toRtkEpoch: () => toRtkEpoch
+  toRtkEpoch: () => toRtkEpoch,
+  wlWavelength: () => wlWavelength
 });
 module.exports = __toCommonJS(positioning_exports);
 
@@ -2430,6 +2433,96 @@ function solvePpp(epochs, sp3, opts) {
   };
 }
 
+// src/positioning/ppp-ar.ts
+function wlWavelength(f1, f2) {
+  return C_LIGHT / (f1 - f2);
+}
+function nlWavelength(f1, f2) {
+  return C_LIGHT / (f1 + f2);
+}
+function resolvePppAmbiguities(sats, Qaif, opts = {}) {
+  const wlThreshold = opts.wlThreshold ?? 0.25;
+  const ratioThreshold = opts.ratioThreshold ?? 3;
+  const kept = [];
+  const nWl = [];
+  const rejectedWl = [];
+  for (let i = 0; i < sats.length; i++) {
+    const s = sats[i];
+    const wlFloat = s.mwCyc - (s.wlBiasCyc ?? 0);
+    const r = Math.round(wlFloat);
+    if (Math.abs(wlFloat - r) <= wlThreshold) {
+      kept.push(i);
+      nWl.push(r);
+    } else {
+      rejectedWl.push(s.prn);
+    }
+  }
+  if (kept.length < 2) {
+    return {
+      fixed: false,
+      ratio: 0,
+      refPrn: null,
+      sats: [],
+      rejectedWl
+    };
+  }
+  const n1Float = [];
+  const lamNl = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    const lNl = nlWavelength(s.f1, s.f2);
+    lamNl.push(lNl);
+    const n1 = s.aIF / lNl - s.f2 / (s.f1 - s.f2) * nWl[k] - (s.nlBiasCyc ?? 0);
+    n1Float.push(n1);
+  }
+  let ref = 0;
+  for (let k = 1; k < kept.length; k++) {
+    const e = sats[kept[k]].elevDeg ?? -1;
+    if (e > (sats[kept[ref]].elevDeg ?? -1)) ref = k;
+  }
+  const m = kept.length - 1;
+  const sd = new Float64Array(m);
+  const Qsd = new Float64Array(m * m);
+  const others = [];
+  for (let k = 0; k < kept.length; k++) if (k !== ref) others.push(k);
+  const n = sats.length;
+  const qN1 = (ka, kb) => Qaif[kept[ka] * n + kept[kb]] / (lamNl[ka] * lamNl[kb]);
+  for (let a = 0; a < m; a++) {
+    const ka = others[a];
+    sd[a] = n1Float[ka] - n1Float[ref];
+    for (let b = 0; b < m; b++) {
+      const kb = others[b];
+      Qsd[a * m + b] = qN1(ka, kb) - qN1(ka, ref) - qN1(ref, kb) + qN1(ref, ref);
+    }
+  }
+  const res = lambdaSearch(sd, Qsd, m, 2);
+  const ratio = res ? Math.min(res.ratio, 999.9) : 0;
+  const fixed = !!res && ratio >= ratioThreshold;
+  const refN1 = Math.round(n1Float[ref]);
+  const sdFix = res?.candidates[0];
+  const outSats = [];
+  for (let k = 0; k < kept.length; k++) {
+    const s = sats[kept[k]];
+    let n1;
+    if (k === ref) {
+      n1 = refN1;
+    } else {
+      const a = others.indexOf(k);
+      n1 = fixed && sdFix ? refN1 + Math.round(sdFix[a]) : Math.round(n1Float[k]);
+    }
+    const lNl = lamNl[k];
+    const aIF = lNl * n1 + lNl * (s.f2 / (s.f1 - s.f2)) * nWl[k];
+    outSats.push({ prn: s.prn, nWl: nWl[k], n1, aIF });
+  }
+  return {
+    fixed,
+    ratio,
+    refPrn: sats[kept[ref]].prn,
+    sats: outSats,
+    rejectedWl
+  };
+}
+
 // src/positioning/index.ts
 var GM_GPS2 = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -2672,11 +2765,14 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   lambdaSearch,
   moonEcef,
   niellMapping,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   solidEarthTide,
   solveDgnss,
   solvePpp,
   solveSpp,
   sunEcef,
-  toRtkEpoch
+  toRtkEpoch,
+  wlWavelength
 });

package/dist/positioning.d.cts

@@ -714,6 +714,85 @@ declare function moonEcef(timeMs: number): [number, number, number];
  */
 declare function solidEarthTide(rcv: [number, number, number], sun: [number, number, number], moon: [number, number, number]): [number, number, number];
 
+/**
+ * PPP-AR — integer ambiguity resolution for Precise Point Positioning.
+ *
+ * Float PPP estimates one *ionosphere-free* ambiguity per satellite arc, a
+ * real number in metres. Fixing those to integers is what takes PPP from the
+ * decimetre to the centimetre level. The classic two-step decomposition is
+ * used here:
+ *
+ *   1. **Wide-lane** N_WL = N1 − N2 from the Melbourne–Wübbena combination.
+ *      The WL wavelength is ~86 cm (GPS), so the arc-averaged MW rounds to an
+ *      integer reliably once the satellite WL fractional-cycle bias is removed.
+ *   2. **Narrow-lane** N1, recovered from the float IF ambiguity with N_WL
+ *      fixed, then resolved with LAMBDA on *between-satellite single
+ *      differences* (which cancel the receiver phase bias) after removing the
+ *      satellite NL fractional-cycle bias.
+ *
+ * The satellite WL/NL fractional-cycle biases come from a phase-bias product
+ * (Bias-SINEX OSB from an AR-capable analysis centre such as CNES or CODE).
+ * Without them, a single receiver cannot separate the satellite bias from the
+ * integer, so those biases are required inputs — pass 0 only for a network
+ * datum where they have already been removed.
+ *
+ * The relation between the IF ambiguity (metres) and the integers is
+ *   A_IF = λ_NL·N1 + λ_NL·(f2/(f1−f2))·N_WL,   λ_NL = c/(f1+f2).
+ */
+/** Wide-lane wavelength c/(f1−f2) (m). */
+declare function wlWavelength(f1: number, f2: number): number;
+/** Narrow-lane wavelength c/(f1+f2) (m). */
+declare function nlWavelength(f1: number, f2: number): number;
+/** One satellite's float ambiguity inputs for AR. */
+interface ArSat {
+    prn: string;
+    /** Float ionosphere-free ambiguity (m) from the PPP EKF. */
+    aIF: number;
+    /** Arc-averaged Melbourne–Wübbena wide-lane estimate (cycles). */
+    mwCyc: number;
+    f1: number;
+    f2: number;
+    /** Satellite wide-lane fractional-cycle bias (cycles). Default 0. */
+    wlBiasCyc?: number;
+    /** Satellite narrow-lane fractional-cycle bias (cycles). Default 0. */
+    nlBiasCyc?: number;
+    /** Elevation (deg), used only to pick the reference satellite. */
+    elevDeg?: number;
+}
+interface ArOptions {
+    /** Max |MW − round| to accept a wide-lane (cycles). Default 0.25. */
+    wlThreshold?: number;
+    /** Min LAMBDA ratio to accept the narrow-lane fix. Default 3. */
+    ratioThreshold?: number;
+}
+interface ArFixedSat {
+    prn: string;
+    nWl: number;
+    n1: number;
+    /** Fixed IF ambiguity (m). */
+    aIF: number;
+}
+interface ArResult {
+    /** Whether the narrow-lane fix passed the ratio test. */
+    fixed: boolean;
+    /** LAMBDA ratio statistic (second-best / best). */
+    ratio: number;
+    /** Reference PRN whose absolute N1 stays float (SD datum). */
+    refPrn: string | null;
+    /** Fixed satellites (includes the reference with its float-rounded N1). */
+    sats: ArFixedSat[];
+    /** PRNs whose wide-lane did not round confidently (excluded). */
+    rejectedWl: string[];
+}
+/**
+ * Resolve PPP integer ambiguities from the float IF ambiguities and their
+ * covariance. `Qaif` is the n×n float covariance (m²) in the same order as
+ * `sats`. Returns the fixed narrow-lane / wide-lane integers and the ratio
+ * test result. A `fixed: false` result means the data did not support a
+ * confident fix (keep the float solution).
+ */
+declare function resolvePppAmbiguities(sats: ArSat[], Qaif: Float64Array, opts?: ArOptions): ArResult;
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -803,4 +882,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 CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -714,6 +714,85 @@ declare function moonEcef(timeMs: number): [number, number, number];
  */
 declare function solidEarthTide(rcv: [number, number, number], sun: [number, number, number], moon: [number, number, number]): [number, number, number];
 
+/**
+ * PPP-AR — integer ambiguity resolution for Precise Point Positioning.
+ *
+ * Float PPP estimates one *ionosphere-free* ambiguity per satellite arc, a
+ * real number in metres. Fixing those to integers is what takes PPP from the
+ * decimetre to the centimetre level. The classic two-step decomposition is
+ * used here:
+ *
+ *   1. **Wide-lane** N_WL = N1 − N2 from the Melbourne–Wübbena combination.
+ *      The WL wavelength is ~86 cm (GPS), so the arc-averaged MW rounds to an
+ *      integer reliably once the satellite WL fractional-cycle bias is removed.
+ *   2. **Narrow-lane** N1, recovered from the float IF ambiguity with N_WL
+ *      fixed, then resolved with LAMBDA on *between-satellite single
+ *      differences* (which cancel the receiver phase bias) after removing the
+ *      satellite NL fractional-cycle bias.
+ *
+ * The satellite WL/NL fractional-cycle biases come from a phase-bias product
+ * (Bias-SINEX OSB from an AR-capable analysis centre such as CNES or CODE).
+ * Without them, a single receiver cannot separate the satellite bias from the
+ * integer, so those biases are required inputs — pass 0 only for a network
+ * datum where they have already been removed.
+ *
+ * The relation between the IF ambiguity (metres) and the integers is
+ *   A_IF = λ_NL·N1 + λ_NL·(f2/(f1−f2))·N_WL,   λ_NL = c/(f1+f2).
+ */
+/** Wide-lane wavelength c/(f1−f2) (m). */
+declare function wlWavelength(f1: number, f2: number): number;
+/** Narrow-lane wavelength c/(f1+f2) (m). */
+declare function nlWavelength(f1: number, f2: number): number;
+/** One satellite's float ambiguity inputs for AR. */
+interface ArSat {
+    prn: string;
+    /** Float ionosphere-free ambiguity (m) from the PPP EKF. */
+    aIF: number;
+    /** Arc-averaged Melbourne–Wübbena wide-lane estimate (cycles). */
+    mwCyc: number;
+    f1: number;
+    f2: number;
+    /** Satellite wide-lane fractional-cycle bias (cycles). Default 0. */
+    wlBiasCyc?: number;
+    /** Satellite narrow-lane fractional-cycle bias (cycles). Default 0. */
+    nlBiasCyc?: number;
+    /** Elevation (deg), used only to pick the reference satellite. */
+    elevDeg?: number;
+}
+interface ArOptions {
+    /** Max |MW − round| to accept a wide-lane (cycles). Default 0.25. */
+    wlThreshold?: number;
+    /** Min LAMBDA ratio to accept the narrow-lane fix. Default 3. */
+    ratioThreshold?: number;
+}
+interface ArFixedSat {
+    prn: string;
+    nWl: number;
+    n1: number;
+    /** Fixed IF ambiguity (m). */
+    aIF: number;
+}
+interface ArResult {
+    /** Whether the narrow-lane fix passed the ratio test. */
+    fixed: boolean;
+    /** LAMBDA ratio statistic (second-best / best). */
+    ratio: number;
+    /** Reference PRN whose absolute N1 stays float (SD datum). */
+    refPrn: string | null;
+    /** Fixed satellites (includes the reference with its float-rounded N1). */
+    sats: ArFixedSat[];
+    /** PRNs whose wide-lane did not round confidently (excluded). */
+    rejectedWl: string[];
+}
+/**
+ * Resolve PPP integer ambiguities from the float IF ambiguities and their
+ * covariance. `Qaif` is the n×n float covariance (m²) in the same order as
+ * `sats`. Returns the fixed narrow-lane / wide-lane integers and the ratio
+ * test result. A `fixed: false` result means the data did not support a
+ * confident fix (keep the float solution).
+ */
+declare function resolvePppAmbiguities(sats: ArSat[], Qaif: Float64Array, opts?: ArOptions): ArResult;
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -803,4 +882,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 CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type PppAntennaModel, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SatGeom, type SppOptions, type SppSolution, type TropoMapping, buildPppAntenna, createPppCorrections, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -12,14 +12,17 @@ import {
   lambdaSearch,
   moonEcef,
   niellMapping,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   solidEarthTide,
   solveDgnss,
   solvePpp,
   solveSpp,
   sunEcef,
-  toRtkEpoch
-} from "./chunk-CRFFUU5Y.js";
+  toRtkEpoch,
+  wlWavelength
+} from "./chunk-EVFD4UNZ.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
@@ -41,11 +44,14 @@ export {
   lambdaSearch,
   moonEcef,
   niellMapping,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   solidEarthTide,
   solveDgnss,
   solvePpp,
   solveSpp,
   sunEcef,
-  toRtkEpoch
+  toRtkEpoch,
+  wlWavelength
 };

package/dist/rinex.cjs

@@ -23,15 +23,18 @@ __export(rinex_exports, {
   EMPTY_WARNINGS: () => EMPTY_WARNINGS,
   SYSTEM_ORDER: () => SYSTEM_ORDER,
   WarningAccumulator: () => WarningAccumulator,
+  biasMetres: () => biasMetres,
   clkBias: () => clkBias,
   createGzipLineSink: () => createGzipLineSink,
   crxDecompress: () => crxDecompress,
   crxRepair: () => crxRepair,
+  findSatBias: () => findSatBias,
   fmtD: () => fmtD,
   fmtF: () => fmtF,
   hdrLine: () => hdrLine,
   padL: () => padL,
   padR: () => padR,
+  parseBiasSinex: () => parseBiasSinex,
   parseClk: () => parseClk,
   parseCrxDataLine: () => parseCrxDataLine,
   parseIonex: () => parseIonex,
@@ -597,6 +600,7 @@ var SYSTEM_NAMES = {
   I: "NavIC",
   S: "SBAS"
 };
+var C_LIGHT = 299792458;
 
 // src/rinex/parser.ts
 var CHUNK_SIZE = 2 * 1024 * 1024;
@@ -2450,6 +2454,81 @@ function clkBias(clk, prn, tMs) {
   return c[lo] + f * (c[lo + 1] - c[lo]);
 }
 
+// src/rinex/bias-sinex.ts
+function parseEpoch(s, openTo) {
+  const m = /^(\d{4}):(\d{3}):(\d{5})$/.exec(s.trim());
+  if (!m) return openTo;
+  const y = +m[1];
+  if (y === 0) return openTo;
+  const doy = +m[2];
+  const sod = +m[3];
+  return Date.UTC(y, 0, 1) + (doy - 1) * 864e5 + sod * 1e3;
+}
+function parseBiasSinex(text) {
+  const records = [];
+  let inBlock = false;
+  for (const line of text.split("\n")) {
+    if (line.startsWith("+BIAS/SOLUTION")) {
+      inBlock = true;
+      continue;
+    }
+    if (line.startsWith("-BIAS/SOLUTION")) {
+      inBlock = false;
+      continue;
+    }
+    if (!inBlock || line.length < 65 || line.startsWith("*")) continue;
+    const type = line.slice(1, 4).trim();
+    if (type !== "DSB" && type !== "OSB" && type !== "ISB") continue;
+    const prn = line.slice(11, 14).trim();
+    const station = line.slice(15, 24).trim();
+    const obs1 = line.slice(25, 29).trim();
+    const obs2 = line.slice(30, 34).trim();
+    const startMs = parseEpoch(line.slice(35, 49), -Infinity);
+    const endMs = parseEpoch(line.slice(50, 64), Infinity);
+    const tail = line.slice(65).trim().split(/\s+/);
+    if (tail.length < 2) continue;
+    const unit = tail[0];
+    const value = Number(tail[1]);
+    if (!Number.isFinite(value)) continue;
+    records.push({
+      type,
+      prn,
+      station,
+      obs1,
+      obs2,
+      startMs,
+      endMs,
+      value,
+      unit
+    });
+  }
+  const bySat = /* @__PURE__ */ new Map();
+  for (const r of records) {
+    if (!r.prn || r.station) continue;
+    let a = bySat.get(r.prn);
+    if (!a) {
+      a = [];
+      bySat.set(r.prn, a);
+    }
+    a.push(r);
+  }
+  return { records, bySat };
+}
+function biasMetres(r) {
+  return r.unit === "ns" ? r.value * 1e-9 * C_LIGHT : r.value;
+}
+function findSatBias(bias, prn, obs1, obs2, tMs) {
+  const recs = bias.bySat.get(prn);
+  if (!recs) return null;
+  for (const r of recs) {
+    if (r.obs1 !== obs1) continue;
+    if (obs2 === void 0 ? r.obs2 !== "" : r.obs2 !== obs2) continue;
+    if (tMs != null && (tMs < r.startMs || tMs > r.endMs)) continue;
+    return r;
+  }
+  return null;
+}
+
 // src/rinex/obs-writer-common.ts
 function createGzipLineSink() {
   const encoder = new TextEncoder();
@@ -2828,15 +2907,18 @@ function writeRinex4ObsBlob(header, epochs, obsTypes) {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,
   WarningAccumulator,
+  biasMetres,
   clkBias,
   createGzipLineSink,
   crxDecompress,
   crxRepair,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
   parseClk,
   parseCrxDataLine,
   parseIonex,

package/dist/rinex.d.cts

@@ -182,6 +182,54 @@ interface Nav4Input {
  */
 declare function writeRinexNav4(nav: Nav4Input): string;
 
+/**
+ * Bias-SINEX (`.BIA`) parser — GNSS code/phase biases published by the IGS
+ * analysis centres. Two record kinds matter here:
+ *
+ *   - **DSB** differential signal bias (OBS1−OBS2), the classic inter-code
+ *     bias used for ionosphere/DCB corrections in SPP.
+ *   - **OSB** observable-specific (absolute) bias for a single OBS, the
+ *     ingredient PPP-AR needs: the *phase* OSBs are the satellite
+ *     fractional-cycle biases that make the ambiguities integer-recoverable.
+ *
+ * Values keep their file unit (`ns` for code/most phase biases, `cyc` for
+ * cycle-based phase biases); `biasMetres` converts ns → m via c.
+ */
+interface BiasRecord {
+    /** 'DSB' | 'OSB' | 'ISB' | … */
+    type: string;
+    /** Satellite PRN (e.g. 'G01', 'E11'); '' for a pure station bias. */
+    prn: string;
+    /** Station (9-char), '' for satellite biases. */
+    station: string;
+    /** First observable code (e.g. 'C1W', 'L1C'). */
+    obs1: string;
+    /** Second observable (DSB only; '' for OSB). */
+    obs2: string;
+    /** Validity start / end (GPS-scale ms); ±Infinity when open. */
+    startMs: number;
+    endMs: number;
+    /** Estimated value in the file unit. */
+    value: number;
+    /** File unit ('ns' | 'cyc' | …). */
+    unit: string;
+}
+interface BiasSinex {
+    records: BiasRecord[];
+    /** All satellite records grouped by PRN (station === ''). */
+    bySat: Map<string, BiasRecord[]>;
+}
+/** Parse a Bias-SINEX file. */
+declare function parseBiasSinex(text: string): BiasSinex;
+/** A satellite bias record's value in metres (ns → m via c; cyc left as-is). */
+declare function biasMetres(r: BiasRecord): number;
+/**
+ * Look up a satellite bias: a DSB by its (obs1,obs2) pair, or an OSB by obs1
+ * alone (pass obs2 = undefined), valid at `tMs` (default: any). Returns the
+ * record or null.
+ */
+declare function findSatBias(bias: BiasSinex, prn: string, obs1: string, obs2?: string, tMs?: number): BiasRecord | null;
+
 /**
  * RINEX 3.04 / 4.01 observation file writer.
  *
@@ -263,4 +311,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.d.ts

@@ -182,6 +182,54 @@ interface Nav4Input {
  */
 declare function writeRinexNav4(nav: Nav4Input): string;
 
+/**
+ * Bias-SINEX (`.BIA`) parser — GNSS code/phase biases published by the IGS
+ * analysis centres. Two record kinds matter here:
+ *
+ *   - **DSB** differential signal bias (OBS1−OBS2), the classic inter-code
+ *     bias used for ionosphere/DCB corrections in SPP.
+ *   - **OSB** observable-specific (absolute) bias for a single OBS, the
+ *     ingredient PPP-AR needs: the *phase* OSBs are the satellite
+ *     fractional-cycle biases that make the ambiguities integer-recoverable.
+ *
+ * Values keep their file unit (`ns` for code/most phase biases, `cyc` for
+ * cycle-based phase biases); `biasMetres` converts ns → m via c.
+ */
+interface BiasRecord {
+    /** 'DSB' | 'OSB' | 'ISB' | … */
+    type: string;
+    /** Satellite PRN (e.g. 'G01', 'E11'); '' for a pure station bias. */
+    prn: string;
+    /** Station (9-char), '' for satellite biases. */
+    station: string;
+    /** First observable code (e.g. 'C1W', 'L1C'). */
+    obs1: string;
+    /** Second observable (DSB only; '' for OSB). */
+    obs2: string;
+    /** Validity start / end (GPS-scale ms); ±Infinity when open. */
+    startMs: number;
+    endMs: number;
+    /** Estimated value in the file unit. */
+    value: number;
+    /** File unit ('ns' | 'cyc' | …). */
+    unit: string;
+}
+interface BiasSinex {
+    records: BiasRecord[];
+    /** All satellite records grouped by PRN (station === ''). */
+    bySat: Map<string, BiasRecord[]>;
+}
+/** Parse a Bias-SINEX file. */
+declare function parseBiasSinex(text: string): BiasSinex;
+/** A satellite bias record's value in metres (ns → m via c; cyc left as-is). */
+declare function biasMetres(r: BiasRecord): number;
+/**
+ * Look up a satellite bias: a DSB by its (obs1,obs2) pair, or an OSB by obs1
+ * alone (pass obs2 = undefined), valid at `tMs` (default: any). Returns the
+ * record or null.
+ */
+declare function findSatBias(bias: BiasSinex, prn: string, obs1: string, obs2?: string, tMs?: number): BiasRecord | null;
+
 /**
  * RINEX 3.04 / 4.01 observation file writer.
  *
@@ -263,4 +311,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type BiasRecord, type BiasSinex, type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, biasMetres, createGzipLineSink, crxDecompress, crxRepair, findSatBias, fmtD, fmtF, hdrLine, padL, padR, parseBiasSinex, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };

package/dist/rinex.js

@@ -1,12 +1,15 @@
 import {
   EMPTY_WARNINGS,
   WarningAccumulator,
+  biasMetres,
   createGzipLineSink,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
   parseIonex,
   parseNavFile,
   stationHeaderLines,
@@ -17,7 +20,7 @@ import {
   writeRinexNav,
   writeRinexNav4,
   writeRinexObsBlob
-} from "./chunk-KDM2NVBO.js";
+} from "./chunk-5RNDQQ6J.js";
 import "./chunk-726LQBGM.js";
 import "./chunk-IQ5OIMQD.js";
 import {
@@ -40,15 +43,18 @@ export {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,
   WarningAccumulator,
+  biasMetres,
   clkBias,
   createGzipLineSink,
   crxDecompress,
   crxRepair,
+  findSatBias,
   fmtD,
   fmtF,
   hdrLine,
   padL,
   padR,
+  parseBiasSinex,
   parseClk,
   parseCrxDataLine,
   parseIonex,

package/package.json

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