gnss-js 1.29.0 → 1.31.0

package/dist/index.d.cts

@@ -4,10 +4,10 @@ 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 { I as IonexGrid, S as Sp3File, a as Sp3Sample, p as parseIonex, b as parseSp3, s as sp3Position } from './sp3-BAq8l39b.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';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.cjs';
 export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.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,10 +4,10 @@ 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 { I as IonexGrid, S as Sp3File, a as Sp3Sample, p as parseIonex, b as parseSp3, s as sp3Position } from './sp3-BAq8l39b.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';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.js';
 export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.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-5WHKMKIB.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-6GTS5MVI.js";
+} from "./chunk-5RNDQQ6J.js";
 import "./chunk-726LQBGM.js";
 import "./chunk-IQ5OIMQD.js";
 import {
@@ -230,9 +236,11 @@ import {
   systemName
 } from "./chunk-3U5AX7PY.js";
 import {
+  clkBias,
+  parseClk,
   parseSp3,
   sp3Position
-} from "./chunk-4KMWWCUR.js";
+} from "./chunk-GPEP67TG.js";
 import {
   BitReader,
   RTCM3_MESSAGE_NAMES,
@@ -426,6 +434,7 @@ export {
   analyzeQuality,
   applyHelmert,
   applyIonoDcb,
+  biasMetres,
   binexCrc16,
   binexCrc32,
   binexCsum8,
@@ -434,6 +443,7 @@ export {
   buildObsIndices,
   buildPppAntenna,
   clampUnit,
+  clkBias,
   computeAllPositions,
   computeDop,
   computeIonoRate,
@@ -464,6 +474,7 @@ export {
   ephInfoToEphemeris,
   euclidean3D,
   fetchSourcetable,
+  findSatBias,
   fmtD,
   fmtF,
   formatUTCTime,
@@ -546,12 +557,15 @@ export {
   msmEpochToDate,
   navTimesFromEph,
   niellMapping,
+  nlWavelength,
   nmeaCoordToDecimal,
   oem4Frames,
   padL,
   padR,
   parseAntex,
+  parseBiasSinex,
   parseBinex,
+  parseClk,
   parseCrxDataLine,
   parseIonex,
   parseNavFile,
@@ -587,6 +601,7 @@ export {
   readString,
   reportDecodeError,
   resetGloFreqCache,
+  resolvePppAmbiguities,
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
@@ -612,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);
 
@@ -1597,6 +1600,28 @@ function sp3Position(sp3, prn, tMs, order = 9) {
   return { x, y, z, clk };
 }
 
+// src/rinex/clk.ts
+var GPS_EPOCH_MS2 = Date.UTC(1980, 0, 6);
+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/positioning/ppp-tropo.ts
 var HYD_AVG_A = [
   0.0012769934,
@@ -2080,7 +2105,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;
@@ -2099,7 +2124,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;
@@ -2213,7 +2238,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,
@@ -2408,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;
@@ -2650,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

@@ -1,6 +1,6 @@
 import { E as Ephemeris } from './nav-DImXUdbp.cjs';
 import { DopValues } from './orbit.cjs';
-import { I as IonexGrid, S as Sp3File } from './sp3-BAq8l39b.cjs';
+import { I as IonexGrid, C as ClkFile, S as Sp3File } from './clk-DGVBGRzA.cjs';
 import { AntexFile } from './antex.cjs';
 import './ephemeris-Ohl6NAfw.cjs';
 import './ecef-CF0uAysr.cjs';
@@ -631,6 +631,11 @@ interface PppOptions {
     ztdProcessNoise?: number;
     /** Model + estimate the troposphere. Default true. */
     troposphere?: boolean;
+    /** High-rate precise satellite clocks (RINEX CLOCK, typically 30 s). When
+     * given, satellite clock offsets come from this by linear interpolation
+     * instead of the SP3 5-minute clocks — the dominant PPP accuracy term.
+     * Satellites/epochs absent from the CLK fall back to the SP3 clock. */
+    clk?: ClkFile;
 }
 interface PppEpochResult {
     timeMs: number;
@@ -709,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.
@@ -798,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

@@ -1,6 +1,6 @@
 import { E as Ephemeris } from './nav-DImXUdbp.js';
 import { DopValues } from './orbit.js';
-import { I as IonexGrid, S as Sp3File } from './sp3-BAq8l39b.js';
+import { I as IonexGrid, C as ClkFile, S as Sp3File } from './clk-DGVBGRzA.js';
 import { AntexFile } from './antex.js';
 import './ephemeris-Ohl6NAfw.js';
 import './ecef-CF0uAysr.js';
@@ -631,6 +631,11 @@ interface PppOptions {
     ztdProcessNoise?: number;
     /** Model + estimate the troposphere. Default true. */
     troposphere?: boolean;
+    /** High-rate precise satellite clocks (RINEX CLOCK, typically 30 s). When
+     * given, satellite clock offsets come from this by linear interpolation
+     * instead of the SP3 5-minute clocks — the dominant PPP accuracy term.
+     * Satellites/epochs absent from the CLK fall back to the SP3 clock. */
+    clk?: ClkFile;
 }
 interface PppEpochResult {
     timeMs: number;
@@ -709,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.
@@ -798,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,18 +12,21 @@ import {
   lambdaSearch,
   moonEcef,
   niellMapping,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   solidEarthTide,
   solveDgnss,
   solvePpp,
   solveSpp,
   sunEcef,
-  toRtkEpoch
-} from "./chunk-5WHKMKIB.js";
+  toRtkEpoch,
+  wlWavelength
+} from "./chunk-EVFD4UNZ.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
-import "./chunk-4KMWWCUR.js";
+import "./chunk-GPEP67TG.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";
 import "./chunk-FIEWO4J4.js";
@@ -41,11 +44,14 @@ export {
   lambdaSearch,
   moonEcef,
   niellMapping,
+  nlWavelength,
+  resolvePppAmbiguities,
   satClockCorrection,
   solidEarthTide,
   solveDgnss,
   solvePpp,
   solveSpp,
   sunEcef,
-  toRtkEpoch
+  toRtkEpoch,
+  wlWavelength
 };