gnss-js 1.29.0 → 1.31.0

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,
@@ -119,6 +120,7 @@ __export(src_exports, {
   buildObsIndices: () => buildObsIndices,
   buildPppAntenna: () => buildPppAntenna,
   clampUnit: () => clampUnit,
+  clkBias: () => clkBias,
   computeAllPositions: () => computeAllPositions,
   computeDop: () => computeDop,
   computeIonoRate: () => computeIonoRate,
@@ -149,6 +151,7 @@ __export(src_exports, {
   ephInfoToEphemeris: () => ephInfoToEphemeris,
   euclidean3D: () => euclidean3D,
   fetchSourcetable: () => fetchSourcetable,
+  findSatBias: () => findSatBias,
   fmtD: () => fmtD,
   fmtF: () => fmtF,
   formatUTCTime: () => formatUTCTime,
@@ -231,12 +234,15 @@ __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,
   parseIonex: () => parseIonex,
   parseNavFile: () => parseNavFile,
@@ -272,6 +278,7 @@ __export(src_exports, {
   readString: () => readString,
   reportDecodeError: () => reportDecodeError,
   resetGloFreqCache: () => resetGloFreqCache,
+  resolvePppAmbiguities: () => resolvePppAmbiguities,
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   satClockCorrection: () => satClockCorrection,
@@ -297,6 +304,7 @@ __export(src_exports, {
   updateStreamStats: () => updateStreamStats,
   verifyChecksum: () => verifyChecksum,
   vincenty: () => vincenty,
+  wlWavelength: () => wlWavelength,
   writeCrx: () => writeCrx,
   writeModernObsBlob: () => writeModernObsBlob,
   writeRinex2ObsBlob: () => writeRinex2ObsBlob,
@@ -3896,6 +3904,141 @@ function sp3Position(sp3, prn, tMs, order = 9) {
   return { x, y, z, clk };
 }
 
+// src/rinex/clk.ts
+var GPS_EPOCH_MS4 = Date.UTC(1980, 0, 6);
+function parseClk(text) {
+  const sats = {};
+  let inData = false;
+  for (const line of text.split("\n")) {
+    if (!inData) {
+      if (line.includes("END OF HEADER")) inData = true;
+      continue;
+    }
+    if (line.length < 3 || line[0] !== "A" || line[1] !== "S") continue;
+    const p = line.slice(2).trim().split(/\s+/);
+    if (p.length < 9) continue;
+    const prn = p[0];
+    const y = +p[1];
+    const mo = +p[2];
+    const d = +p[3];
+    const h = +p[4];
+    const mi = +p[5];
+    const s = +p[6];
+    const bias = +p[8];
+    if (!Number.isFinite(bias)) continue;
+    const t = Date.UTC(y, mo - 1, d, h, mi, Math.floor(s)) + s % 1 * 1e3;
+    let e = sats[prn];
+    if (!e) {
+      e = { t: [], c: [] };
+      sats[prn] = e;
+    }
+    e.t.push(t);
+    e.c.push(bias);
+  }
+  let intervalSec = 30;
+  for (const e of Object.values(sats)) {
+    if (e.t.length >= 2) {
+      intervalSec = Math.round((e.t[1] - e.t[0]) / 1e3);
+      break;
+    }
+  }
+  return { sats, intervalSec };
+}
+function clkBias(clk, prn, tMs) {
+  const e = clk.sats[prn];
+  if (!e || e.t.length === 0) return null;
+  const { t, c } = e;
+  let lo = 0;
+  let hi = t.length - 1;
+  if (tMs < t[0] || tMs > t[hi]) return null;
+  while (lo < hi) {
+    const mid = lo + hi + 1 >> 1;
+    if (t[mid] <= tMs) lo = mid;
+    else hi = mid - 1;
+  }
+  if (t[lo] === tMs || lo === t.length - 1) return c[lo];
+  const t0 = t[lo];
+  const t1 = t[lo + 1];
+  if (t1 - t0 > clk.intervalSec * 1e3 * 4) return null;
+  const f = (tMs - t0) / (t1 - t0);
+  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();
@@ -6147,7 +6290,7 @@ var BdsAssembler = class {
 };
 
 // src/navbits/gal.ts
-var GPS_EPOCH_MS4 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS5 = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK5 = 7 * 86400;
 var HALF_WEEK3 = 302400;
 var GST_GPS_WEEK_OFFSET = 1024;
@@ -6160,7 +6303,7 @@ function galFnavPageCrcOk(page) {
   return crc24q(page, 214) === getBitU(page, 214, 24);
 }
 function gpsDate2(week, sec) {
-  return new Date(GPS_EPOCH_MS4 + (week * SEC_PER_WEEK5 + sec) * 1e3);
+  return new Date(GPS_EPOCH_MS5 + (week * SEC_PER_WEEK5 + sec) * 1e3);
 }
 function weekOfToe(week, tow, toe) {
   if (toe - tow > HALF_WEEK3) return week - 1;
@@ -6455,7 +6598,7 @@ var GalFnavAssembler = class {
 };
 
 // src/navbits/glo.ts
-var GPS_EPOCH_MS5 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS6 = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK6 = 7 * 86400;
 var SEC_PER_DAY = 86400;
 var GLO_STRING_BYTES = 10;
@@ -6549,7 +6692,7 @@ function decodeGloStrings(strings, refDate, opts = {}) {
   const slot = getBitU(b, i, 5);
   if (frn1 !== 1 || frn2 !== 2 || frn3 !== 3 || frn4 !== 4) return null;
   if (slot < 1 || slot > 27) return null;
-  const utcSec = (getUtcDate(refDate).getTime() - GPS_EPOCH_MS5) / 1e3;
+  const utcSec = (getUtcDate(refDate).getTime() - GPS_EPOCH_MS6) / 1e3;
   const week = Math.floor(utcSec / SEC_PER_WEEK6);
   let tow = utcSec - week * SEC_PER_WEEK6;
   const tod = tow % SEC_PER_DAY;
@@ -6564,7 +6707,7 @@ function decodeGloStrings(strings, refDate, opts = {}) {
     system: "R",
     prn: `R${String(slot).padStart(2, "0")}`,
     // RINEX GLONASS epochs are UTC: build the UTC toe Date directly.
-    tocDate: new Date(GPS_EPOCH_MS5 + (week * SEC_PER_WEEK6 + tow + toe) * 1e3),
+    tocDate: new Date(GPS_EPOCH_MS6 + (week * SEC_PER_WEEK6 + tow + toe) * 1e3),
     tauN: -tauN,
     // RINEX stores −τn; the SIS carries τn (ICD sign)
     gammaN,
@@ -6618,7 +6761,7 @@ var GloStringAssembler = class {
 };
 
 // src/ubx/rawnav.ts
-var GPS_EPOCH_MS6 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS7 = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK7 = 7 * 86400;
 var two2 = (n) => String(n).padStart(2, "0");
 function parseUbxRawNav(data, opts = {}) {
@@ -6638,7 +6781,7 @@ function parseUbxRawNav(data, opts = {}) {
       if (week === 0 && opts.refWeek !== void 0) week = opts.refWeek;
       if (week > 0) {
         refDate = new Date(
-          GPS_EPOCH_MS6 + (week * SEC_PER_WEEK7 + rcvTow) * 1e3
+          GPS_EPOCH_MS7 + (week * SEC_PER_WEEK7 + rcvTow) * 1e3
         );
       }
       continue;
@@ -6718,7 +6861,7 @@ function parseUbxRawNav(data, opts = {}) {
 }
 
 // src/ubx/index.ts
-var GPS_EPOCH_MS7 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS8 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK3 = 7 * 864e5;
 var SIGNALS = {
   0: { 0: ["G", "1C"], 3: ["G", "2X"], 4: ["G", "2X"] },
@@ -6801,7 +6944,7 @@ function parseUbxRawx(data) {
         if (!codes.includes(sig[1])) codes.push(sig[1]);
       }
       epochs.push({
-        timeMs: GPS_EPOCH_MS7 + week * MS_PER_WEEK3 + Math.round(rcvTow * 1e3),
+        timeMs: GPS_EPOCH_MS8 + week * MS_PER_WEEK3 + Math.round(rcvTow * 1e3),
         leapS: leapValid ? leapS : null,
         meas
       });
@@ -6872,7 +7015,7 @@ function svidToPrn(svid) {
 
 // src/sbf/nav.ts
 var PI2 = Math.PI;
-var GPS_EPOCH_MS8 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS9 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS2 = Date.UTC(2006, 0, 1);
 var SEC_PER_WEEK8 = 7 * 86400;
 var MS_PER_WEEK4 = SEC_PER_WEEK8 * 1e3;
@@ -6884,7 +7027,7 @@ function adjustWeek(refWeek, wn, mod) {
   if (offset < -(mod / 2 - 1)) offset += mod;
   return refWeek - offset;
 }
-var gpsMs = (week, sec) => GPS_EPOCH_MS8 + week * MS_PER_WEEK4 + sec * 1e3;
+var gpsMs = (week, sec) => GPS_EPOCH_MS9 + week * MS_PER_WEEK4 + sec * 1e3;
 var sowOf2 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK8;
 function decodeGpsQzsNav(view, b, sys) {
   const svid = view.getUint8(b + 14);
@@ -7034,7 +7177,7 @@ function decodeGloNav(view, b) {
   const tocDate = new Date(toeGpsMs - leapMs);
   const tofGpsMs = gpsMs(wnc, 0) + view.getUint32(b + 8, true);
   const tofLeapMs = getGpsLeap(new Date(tofGpsMs)) * 1e3;
-  const messageFrameTime = (tofGpsMs - tofLeapMs - GPS_EPOCH_MS8) / 1e3 % SEC_PER_WEEK8;
+  const messageFrameTime = (tofGpsMs - tofLeapMs - GPS_EPOCH_MS9) / 1e3 % SEC_PER_WEEK8;
   return {
     system: "R",
     prn,
@@ -7306,7 +7449,7 @@ function isFnavDummy(page) {
 }
 
 // src/sbf/rawnav-bds.ts
-var GPS_EPOCH_MS9 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS10 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK5 = 7 * 86400 * 1e3;
 var TOW_DNU = 4294967295;
 var WNC_DNU = 65535;
@@ -7376,7 +7519,7 @@ function parseSbfGloNav(data) {
     const eph = assembler.push(
       prn,
       str,
-      new Date(GPS_EPOCH_MS9 + wnc * MS_PER_WEEK5 + towMs),
+      new Date(GPS_EPOCH_MS10 + wnc * MS_PER_WEEK5 + towMs),
       view.getUint8(b + 18) - 8
       // FreqNr, offset by 8
     );
@@ -7910,7 +8053,7 @@ function parseSbfHas(data) {
 
 // src/sbf/index.ts
 var CLIGHT = 299792458;
-var GPS_EPOCH_MS10 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS11 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK6 = 7 * 864e5;
 var M3_SYS = ["G", "R", "E", "C", "S", "J", "I"];
 var M3_PR_BASE = [19e6, 19e6, 22e6, 2e7, 34e6, 34e6, 34e6];
@@ -8564,7 +8707,7 @@ function parseSbfMeas(data) {
       const tow = view.getUint32(i + 8, true);
       const wnc = view.getUint16(i + 12, true);
       if (tow !== 4294967295 && wnc !== 65535) {
-        ensureEpoch(GPS_EPOCH_MS10 + wnc * MS_PER_WEEK6 + tow);
+        ensureEpoch(GPS_EPOCH_MS11 + wnc * MS_PER_WEEK6 + tow);
         if (id === 4109) decodeMeas3Ranges(i, tow);
         else if (id === 4110) decodeMeas3CN0(i, len);
         else if (id === 4111) decodeMeas3Doppler(i, len);
@@ -8629,15 +8772,15 @@ var ID_GALEPHEMERIS = 1122;
 var ID_QZSSEPHEMERIS = 1336;
 var ID_GPSEPHEM = 7;
 var ID_BDSEPHEMERIS = 1696;
-var GPS_EPOCH_MS11 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS12 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS3 = Date.UTC(2006, 0, 1);
 var SEC_PER_WEEK9 = 7 * 86400;
 var SEC_PER_DAY2 = 86400;
 var MS_PER_WEEK7 = SEC_PER_WEEK9 * 1e3;
-var gpsMs2 = (week, sec) => GPS_EPOCH_MS11 + week * MS_PER_WEEK7 + sec * 1e3;
+var gpsMs2 = (week, sec) => GPS_EPOCH_MS12 + week * MS_PER_WEEK7 + sec * 1e3;
 var sowOf3 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK9;
 function nearWeekMs(refMs, sow) {
-  const refWeek = Math.floor((refMs - GPS_EPOCH_MS11) / MS_PER_WEEK7);
+  const refWeek = Math.floor((refMs - GPS_EPOCH_MS12) / MS_PER_WEEK7);
   let ms = gpsMs2(refWeek, sow);
   if (ms < refMs - MS_PER_WEEK7 / 2) ms += MS_PER_WEEK7;
   else if (ms > refMs + MS_PER_WEEK7 / 2) ms -= MS_PER_WEEK7;
@@ -8664,9 +8807,9 @@ function decodeGloEphemeris(view, p) {
   tof += Math.floor(tow / SEC_PER_DAY2) * SEC_PER_DAY2;
   if (tof < tow - 43200) tof += SEC_PER_DAY2;
   else if (tof > tow + 43200) tof -= SEC_PER_DAY2;
-  const toeGpsMs = GPS_EPOCH_MS11 + (week * SEC_PER_WEEK9 + tow) * 1e3;
-  const tofGpsMs = GPS_EPOCH_MS11 + (week * SEC_PER_WEEK9 + tof) * 1e3;
-  const tofUtcSec = (getUtcDate(new Date(tofGpsMs)).getTime() - GPS_EPOCH_MS11) / 1e3;
+  const toeGpsMs = GPS_EPOCH_MS12 + (week * SEC_PER_WEEK9 + tow) * 1e3;
+  const tofGpsMs = GPS_EPOCH_MS12 + (week * SEC_PER_WEEK9 + tof) * 1e3;
+  const tofUtcSec = (getUtcDate(new Date(tofGpsMs)).getTime() - GPS_EPOCH_MS12) / 1e3;
   return {
     system: "R",
     prn: `R${String(slot).padStart(2, "0")}`,
@@ -8731,7 +8874,7 @@ function decodeGalEphemeris(view, p, headerMs) {
       omega: view.getFloat64(p + 60, true),
       omegaDot: view.getFloat64(p + 140, true),
       idot: view.getFloat64(p + 124, true),
-      week: Math.floor((toeMs - GPS_EPOCH_MS11) / MS_PER_WEEK7),
+      week: Math.floor((toeMs - GPS_EPOCH_MS12) / MS_PER_WEEK7),
       // RINEX Galileo SVH bit layout (E1B DVS/HS in bits 0-2, E5a in
       // 3-5, E5b in 6-8) — same packing as RTKLIB.
       svHealth: svhE5b << 7 | dvsE5b << 6 | svhE5a << 4 | dvsE5a << 3 | svhE1b << 1 | dvsE1b,
@@ -8903,7 +9046,7 @@ function parseNovatelNav(data) {
 }
 
 // src/novatel/index.ts
-var GPS_EPOCH_MS12 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS13 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK8 = 7 * 864e5;
 var HDR = 28;
 var ID_RANGE = 43;
@@ -9073,7 +9216,7 @@ function parseNovatelRange(data) {
           if (!codes.includes(code)) codes.push(code);
         }
         epochs.push({
-          timeMs: GPS_EPOCH_MS12 + week * MS_PER_WEEK8 + towMs,
+          timeMs: GPS_EPOCH_MS13 + week * MS_PER_WEEK8 + towMs,
           meas
         });
       }
@@ -9125,7 +9268,7 @@ function parseNovatelRange(data) {
           if (!codes.includes(code)) codes.push(code);
         }
         epochs.push({
-          timeMs: GPS_EPOCH_MS12 + week * MS_PER_WEEK8 + towMs,
+          timeMs: GPS_EPOCH_MS13 + week * MS_PER_WEEK8 + towMs,
           meas
         });
       }
@@ -9176,11 +9319,11 @@ function* trimbleFrames(data, stats) {
 }
 
 // src/trimble/nav.ts
-var GPS_EPOCH_MS13 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS14 = Date.UTC(1980, 0, 6);
 var SEC_PER_WEEK10 = 7 * 86400;
 var MS_PER_WEEK9 = SEC_PER_WEEK10 * 1e3;
 var GPS_PI2 = 3.1415926535898;
-var gpsMs3 = (week, sec) => GPS_EPOCH_MS13 + week * MS_PER_WEEK9 + sec * 1e3;
+var gpsMs3 = (week, sec) => GPS_EPOCH_MS14 + week * MS_PER_WEEK9 + sec * 1e3;
 var sowOf4 = (dateMs) => dateMs / 1e3 % SEC_PER_WEEK10;
 var SUB_GPS_EPHEMERIS = 1;
 var SUB_ION_UTC = 3;
@@ -9282,7 +9425,7 @@ function parseTrimbleNav(data) {
 }
 
 // src/trimble/index.ts
-var GPS_EPOCH_MS14 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS15 = Date.UTC(1980, 0, 6);
 var MS_PER_WEEK10 = 7 * 864e5;
 var REC_TYPE_27 = 6;
 var SYS = {
@@ -9480,7 +9623,7 @@ function decodeEpoch(epoch, epochs, obsCodes) {
     o += c.length;
   }
   const dv = new DataView(M.buffer, M.byteOffset, M.byteLength);
-  const timeMs = GPS_EPOCH_MS14 + epoch.week * MS_PER_WEEK10 + epoch.towMs;
+  const timeMs = GPS_EPOCH_MS15 + epoch.week * MS_PER_WEEK10 + epoch.towMs;
   const meas = [];
   let q = 0;
   for (let sv = 0; sv < epoch.nsv && q + 8 <= M.length; sv++) {
@@ -9661,13 +9804,13 @@ function* binexRecords(data, view, stats) {
 }
 
 // src/binex/nav.ts
-var GPS_EPOCH_MS15 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS16 = Date.UTC(1980, 0, 6);
 var GLO_DAY0_MS = Date.UTC(1980, 0, 1);
 var BDT_EPOCH_MS4 = Date.UTC(2006, 0, 1);
 var SEC_PER_WEEK11 = 7 * 86400;
 var MS_PER_WEEK11 = SEC_PER_WEEK11 * 1e3;
 var SC2RAD = GPS_PI;
-var gpsMs4 = (week, sec) => GPS_EPOCH_MS15 + week * MS_PER_WEEK11 + sec * 1e3;
+var gpsMs4 = (week, sec) => GPS_EPOCH_MS16 + week * MS_PER_WEEK11 + sec * 1e3;
 var sowOf5 = (dateMs) => (dateMs / 1e3 % SEC_PER_WEEK11 + SEC_PER_WEEK11) % SEC_PER_WEEK11;
 var two5 = (n) => String(n).padStart(2, "0");
 var Reader = class {
@@ -9955,7 +10098,7 @@ function decodeBinexEph(view, body, len, le) {
 }
 
 // src/binex/index.ts
-var GPS_EPOCH_MS16 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS17 = Date.UTC(1980, 0, 6);
 var CLIGHT2 = 299792458;
 var OBS_CODES = {
   0: [
@@ -10213,7 +10356,7 @@ function decodeObsRecord(data, view, body, len, obsCodes) {
   p += 4;
   const msec = view.getUint16(p);
   p += 2;
-  const timeMs = GPS_EPOCH_MS16 + min * 6e4 + msec;
+  const timeMs = GPS_EPOCH_MS17 + min * 6e4 + msec;
   const flag = data[p++];
   const nsat = (flag & 63) + 1;
   if (flag & 128) p += 3;
@@ -10299,7 +10442,7 @@ function parseBinex(data) {
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
-var GPS_EPOCH_MS17 = Date.UTC(1980, 0, 6);
+var GPS_EPOCH_MS18 = Date.UTC(1980, 0, 6);
 var BDT_EPOCH_MS5 = Date.UTC(2006, 0, 1);
 var WEEK_MS = 6048e5;
 var BDT_MINUS_GPST_MS = 14e3;
@@ -10308,7 +10451,7 @@ function almanacEpochMs(alm) {
   if (alm.system === "C") {
     return BDT_EPOCH_MS5 + alm.weekAlm * WEEK_MS + alm.toaSec * 1e3 + BDT_MINUS_GPST_MS;
   }
-  return GPS_EPOCH_MS17 + alm.weekAlm * WEEK_MS + alm.toaSec * 1e3;
+  return GPS_EPOCH_MS18 + alm.weekAlm * WEEK_MS + alm.toaSec * 1e3;
 }
 function almanacSatPosition(alm, timeMs) {
   const GM = alm.system === "E" ? GM_GAL : alm.system === "C" ? GM_BDS : GM_GPS;
@@ -10817,10 +10960,10 @@ function selectBest(ephs, timeMs) {
   }
   return best;
 }
-var GPS_EPOCH_MS18 = START_GPS_TIME.getTime();
+var GPS_EPOCH_MS19 = START_GPS_TIME.getTime();
 var MS_PER_WEEK12 = 7 * 864e5;
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
-var GAL_EPOCH_MS = GPS_EPOCH_MS18 + 1024 * MS_PER_WEEK12;
+var GAL_EPOCH_MS = GPS_EPOCH_MS19 + 1024 * MS_PER_WEEK12;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
   if (sys === "R") {
@@ -10871,8 +11014,8 @@ function ephInfoToEphemeris(info) {
   } else if (sys === "E") {
     epochMs = GAL_EPOCH_MS;
   } else {
-    epochMs = GPS_EPOCH_MS18;
-    const weeksNow = (info.lastReceived - GPS_EPOCH_MS18) / MS_PER_WEEK12;
+    epochMs = GPS_EPOCH_MS19;
+    const weeksNow = (info.lastReceived - GPS_EPOCH_MS19) / MS_PER_WEEK12;
     week += 1024 * Math.round((weeksNow - week) / 1024);
   }
   const tocDate = new Date(epochMs + week * MS_PER_WEEK12 + tocSec * 1e3);
@@ -12834,7 +12977,7 @@ function sagnacRotate(x, y, z, travelS) {
   const s = Math.sin(theta);
   return [x * c + y * s, -x * s + y * c, z];
 }
-function satStateAtEmission(sp3, prn, recvTimeMs, rcv) {
+function satStateAtEmission(sp3, prn, recvTimeMs, rcv, clk) {
   let travelS = 0.075;
   let sx = 0;
   let sy = 0;
@@ -12853,7 +12996,7 @@ function satStateAtEmission(sp3, prn, recvTimeMs, rcv) {
       const vz = (pPlus.z - pMinus.z) / 1;
       velDotPos = p.x * vx + p.y * vy + p.z * vz;
     }
-    clkS = p.clk;
+    clkS = (clk && clkBias(clk, prn, tEmit)) ?? p.clk;
     const [rx, ry, rz] = sagnacRotate(p.x, p.y, p.z, travelS);
     sx = rx;
     sy = ry;
@@ -12967,7 +13110,7 @@ function solvePpp(epochs, sp3, opts) {
     for (const o of epoch.obs) {
       if (o.c1 === 0 || o.c2 === 0 || o.l1 === 0 || o.l2 === 0) continue;
       const slip = detectSlip(o, ei);
-      const sat = satStateAtEmission(sp3, o.prn, epoch.timeMs, rcv);
+      const sat = satStateAtEmission(sp3, o.prn, epoch.timeMs, rcv, opts.clk);
       if (!sat) continue;
       const [elRad, azRad] = getAer(
         sat.state.x,
@@ -13162,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;
@@ -15741,6 +15974,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   analyzeQuality,
   applyHelmert,
   applyIonoDcb,
+  biasMetres,
   binexCrc16,
   binexCrc32,
   binexCsum8,
@@ -15749,6 +15983,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   buildObsIndices,
   buildPppAntenna,
   clampUnit,
+  clkBias,
   computeAllPositions,
   computeDop,
   computeIonoRate,
@@ -15779,6 +16014,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   ephInfoToEphemeris,
   euclidean3D,
   fetchSourcetable,
+  findSatBias,
   fmtD,
   fmtF,
   formatUTCTime,
@@ -15861,12 +16097,15 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   msmEpochToDate,
   navTimesFromEph,
   niellMapping,
+  nlWavelength,
   nmeaCoordToDecimal,
   oem4Frames,
   padL,
   padR,
   parseAntex,
+  parseBiasSinex,
   parseBinex,
+  parseClk,
   parseCrxDataLine,
   parseIonex,
   parseNavFile,
@@ -15902,6 +16141,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolvePppAmbiguities,
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
@@ -15927,6 +16167,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   updateStreamStats,
   verifyChecksum,
   vincenty,
+  wlWavelength,
   writeCrx,
   writeModernObsBlob,
   writeRinex2ObsBlob,