gnss-js 2.13.1 → 2.14.0

package/dist/index.d.cts

@@ -20,9 +20,9 @@ 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, EPHEMERIS_VALIDITY_MS, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, isEphemerisValid, keplerPosition, maskRadForAzimuth, navTimesFromEph, nearestValidEphemeris, rangeRate, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
-export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, MsmResolveOptions, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.cjs';
+export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, BLQ_CONSTITUENTS, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, MsmResolveOptions, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OceanLoadingCoefficients, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, PppTroposphereModel, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, oceanLoadingEnu, parseBlq, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, 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 { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CoordinateSample, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, HelmertFit, HelmertModel, HelmertOptions, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, SeriesComponent, SeriesFit, SeriesFitOptions, SeriesPoint, analyzeQuality, applyHelmertFit, applyIonoDcb, computeIonoRate, detrendIonoArcs, estimateHelmert, fitCoordinateSeries, parseSinexBiasDcb } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';
 export { NmeaFix, NmeaStats, NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum } from './nmea.cjs';
 export { BAND_PRESETS, BandPreset, CONSTELLATIONS, ConstellationRow, DELTA_GLO_L1, DELTA_GLO_L2, FREQ_BDS_B1A, FREQ_BDS_B1C, FREQ_BDS_B1I, FREQ_BDS_B2I, FREQ_BDS_B3A, FREQ_BDS_B3I, FREQ_GAL_E1, FREQ_GAL_E5, FREQ_GAL_E5a, FREQ_GAL_E5b, FREQ_GAL_E6, FREQ_GLO_L1, FREQ_GLO_L1OC, FREQ_GLO_L2, FREQ_GLO_L2OC, FREQ_GLO_L3OC, FREQ_GPS_L1, FREQ_GPS_L2, FREQ_GPS_L5, FREQ_NAVIC_L5, FREQ_QZS_L1, FREQ_QZS_L2, FREQ_QZS_L5, FREQ_QZS_L6, Modulation, SignalDef, computePsdDb, phiAltBOC, phiBOCc, phiBOCs, phiBPSK } from './signals.cjs';

package/dist/index.d.ts

@@ -20,9 +20,9 @@ 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, EPHEMERIS_VALIDITY_MS, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, isEphemerisValid, keplerPosition, maskRadForAzimuth, navTimesFromEph, nearestValidEphemeris, rangeRate, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
-export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, MsmResolveOptions, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.js';
+export { AntennaOffset, AppliedOrbitClock, ArFixedSat, ArOptions, ArResult, ArSat, BLQ_CONSTITUENTS, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, MsmResolveOptions, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, OceanLoadingCoefficients, OrbitClockCorrection, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEngine, PppEphemerisSource, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, PppTroposphereModel, RawObservation, ReceiverBias, ReceiverClockSample, RtkAntennaConfig, RtkAntennaModel, RtkAntennaOffset, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkReceiverAntenna, RtkTrackPoint, SatGeom, SatSample, SbasCorrectionSource, SbasIgp, SbasIonoDelay, SbasProcessor, SbasSatCorrection, Sp3EphemerisSource, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, ZeroBaselineEngine, ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, oceanLoadingEnu, parseBlq, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, 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 { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CoordinateSample, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, HelmertFit, HelmertModel, HelmertOptions, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, SeriesComponent, SeriesFit, SeriesFitOptions, SeriesPoint, analyzeQuality, applyHelmertFit, applyIonoDcb, computeIonoRate, detrendIonoArcs, estimateHelmert, fitCoordinateSeries, parseSinexBiasDcb } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';
 export { NmeaFix, NmeaStats, NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum } from './nmea.js';
 export { BAND_PRESETS, BandPreset, CONSTELLATIONS, ConstellationRow, DELTA_GLO_L1, DELTA_GLO_L2, FREQ_BDS_B1A, FREQ_BDS_B1C, FREQ_BDS_B1I, FREQ_BDS_B2I, FREQ_BDS_B3A, FREQ_BDS_B3I, FREQ_GAL_E1, FREQ_GAL_E5, FREQ_GAL_E5a, FREQ_GAL_E5b, FREQ_GAL_E6, FREQ_GLO_L1, FREQ_GLO_L1OC, FREQ_GLO_L2, FREQ_GLO_L2OC, FREQ_GLO_L3OC, FREQ_GPS_L1, FREQ_GPS_L2, FREQ_GPS_L5, FREQ_NAVIC_L5, FREQ_QZS_L1, FREQ_QZS_L2, FREQ_QZS_L5, FREQ_QZS_L6, Modulation, SignalDef, computePsdDb, phiAltBOC, phiBOCc, phiBOCs, phiBPSK } from './signals.js';

package/dist/index.js

@@ -4,11 +4,14 @@ import {
   IonoAccumulator,
   MultipathAccumulator,
   analyzeQuality,
+  applyHelmertFit,
   applyIonoDcb,
   computeIonoRate,
   detrendIonoArcs,
+  estimateHelmert,
+  fitCoordinateSeries,
   parseSinexBiasDcb
-} from "./chunk-MFCPYUEI.js";
+} from "./chunk-GGWK7XLZ.js";
 import {
   frequencyLabel,
   parseAntex
@@ -176,6 +179,7 @@ import {
   transformFrame
 } from "./chunk-4KOXMHUW.js";
 import {
+  BLQ_CONSTITUENTS,
   IONO_L1_M_PER_TECU,
   PppEngine,
   RtkFloatEngine,
@@ -203,6 +207,8 @@ import {
   moonEcef,
   niellMapping,
   nlWavelength,
+  oceanLoadingEnu,
+  parseBlq,
   postProcessRtk,
   rcvAntennaRangeM,
   resolveMsmRanges,
@@ -218,7 +224,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PQLPPFIK.js";
+} from "./chunk-TETPRPZ2.js";
 import "./chunk-H3AKFD25.js";
 import {
   EPHEMERIS_VALIDITY_MS,
@@ -414,6 +420,7 @@ export {
   BAND_LABELS,
   BAND_PRESETS,
   BDS_SATELLITES,
+  BLQ_CONSTITUENTS,
   BitReader,
   CONSTELLATIONS,
   C_LIGHT,
@@ -511,6 +518,7 @@ export {
   analyzeQuality,
   antexFreqOfGroup,
   applyHelmert,
+  applyHelmertFit,
   applyIonoDcb,
   applyOrbitClock,
   biasMetres,
@@ -558,6 +566,7 @@ export {
   ecefToAzEl,
   ecefToGeodetic,
   ephInfoToEphemeris,
+  estimateHelmert,
   estimateNarrowlaneFcb,
   estimateNetworkFcbs,
   estimateWidelaneFcb,
@@ -569,6 +578,7 @@ export {
   feedNavicBlock,
   fetchSourcetable,
   findSatBias,
+  fitCoordinateSeries,
   fixPppPosition,
   fmtD,
   fmtF,
@@ -663,12 +673,14 @@ export {
   nlWavelength,
   nmeaCoordToDecimal,
   obsEpochToDate,
+  oceanLoadingEnu,
   oem4Frames,
   padL,
   padR,
   parseAntex,
   parseBiasSinex,
   parseBinex,
+  parseBlq,
   parseClk,
   parseCrxDataLine,
   parseIonex,

package/dist/positioning.cjs

@@ -20,6 +20,7 @@ var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: tru
 // src/positioning/index.ts
 var positioning_exports = {};
 __export(positioning_exports, {
+  BLQ_CONSTITUENTS: () => BLQ_CONSTITUENTS,
   IONO_L1_M_PER_TECU: () => IONO_L1_M_PER_TECU,
   PppEngine: () => PppEngine,
   RtkFloatEngine: () => RtkFloatEngine,
@@ -47,6 +48,8 @@ __export(positioning_exports, {
   moonEcef: () => moonEcef,
   niellMapping: () => niellMapping,
   nlWavelength: () => nlWavelength,
+  oceanLoadingEnu: () => oceanLoadingEnu,
+  parseBlq: () => parseBlq,
   postProcessRtk: () => postProcessRtk,
   rcvAntennaRangeM: () => rcvAntennaRangeM,
   resolveMsmRanges: () => resolveMsmRanges,
@@ -3544,6 +3547,103 @@ function solidEarthTide(rcv, sun, moon) {
   return disp;
 }
 
+// src/positioning/ocean-loading.ts
+var BLQ_CONSTITUENTS = [
+  "M2",
+  "S2",
+  "N2",
+  "K2",
+  "K1",
+  "O1",
+  "P1",
+  "Q1",
+  "Mf",
+  "Mm",
+  "Ssa"
+];
+var ARGS = [
+  [140519e-9, 2, -2, 0, 0],
+  // M2
+  [145444e-9, 0, 0, 0, 0],
+  // S2
+  [13788e-8, 2, -3, 1, 0],
+  // N2
+  [145842e-9, 2, 0, 0, 0],
+  // K2
+  [72921e-9, 1, 0, 0, 0.25],
+  // K1
+  [67598e-9, 1, -2, 0, -0.25],
+  // O1
+  [72523e-9, -1, 0, 0, -0.25],
+  // P1
+  [64959e-9, 1, -3, 1, -0.25],
+  // Q1
+  [53234e-10, 0, 2, 0, 0],
+  // Mf
+  [26392e-10, 0, 1, -1, 0],
+  // Mm
+  [3982e-10, 2, 0, 0, 0]
+  // Ssa
+];
+var D2R2 = Math.PI / 180;
+function oceanLoadingEnu(c, utcMs) {
+  const day = Math.floor(utcMs / 864e5);
+  const fday = (utcMs - day * 864e5) / 1e3;
+  const days = day - Date.UTC(1975, 0, 1) / 864e5 + 1;
+  const t = (27392.500528 + 1.000000035 * days) / 36525;
+  const t2 = t * t;
+  const t3 = t2 * t;
+  const a = [
+    fday,
+    (279.69668 + 36000.768930485 * t + 303e-6 * t2) * D2R2,
+    // H0
+    (270.434358 + 481267.88314137 * t - 1133e-6 * t2 + 19e-7 * t3) * D2R2,
+    // S0
+    (334.329653 + 4069.0340329577 * t - 0.010325 * t2 - 12e-6 * t3) * D2R2,
+    // P0
+    2 * Math.PI
+  ];
+  const dp = [0, 0, 0];
+  for (let i = 0; i < ARGS.length; i++) {
+    let ang = 0;
+    for (let j = 0; j < 5; j++) ang += a[j] * ARGS[i][j];
+    for (let k = 0; k < 3; k++)
+      dp[k] += c.amplitude[i][k] * Math.cos(ang - c.phase[i][k] * D2R2);
+  }
+  return [-dp[1], -dp[2], dp[0]];
+}
+function parseBlq(text) {
+  const out = /* @__PURE__ */ new Map();
+  const lines = text.split(/\r?\n/);
+  for (let i = 0; i < lines.length; i++) {
+    const l = lines[i];
+    if (!l.trim() || l.startsWith("$$")) continue;
+    const name = l.trim().split(/\s+/)[0].toUpperCase();
+    const rows = [];
+    let j = i + 1;
+    while (j < lines.length && rows.length < 6) {
+      const r = lines[j];
+      j++;
+      if (r.startsWith("$$") || !r.trim()) continue;
+      const v = r.trim().split(/\s+/).map(Number);
+      if (v.length !== 11 || v.some((x) => !Number.isFinite(x))) break;
+      rows.push(v);
+    }
+    if (rows.length === 6) {
+      out.set(name, {
+        amplitude: BLQ_CONSTITUENTS.map(
+          (_, k) => [rows[0][k], rows[1][k], rows[2][k]]
+        ),
+        phase: BLQ_CONSTITUENTS.map(
+          (_, k) => [rows[3][k], rows[4][k], rows[5][k]]
+        )
+      });
+      i = j - 1;
+    }
+  }
+  return out;
+}
+
 // src/positioning/ppp-corrections.ts
 function applyCorrections(corrections, geom, ctx) {
   if (corrections.evaluate) return corrections.evaluate(geom, ctx);
@@ -3561,6 +3661,14 @@ function createPppCorrections(cfg) {
       const sun = sunEcef(timeMs);
       const moon = moonEcef(timeMs);
       const tideDisp = useTides ? solidEarthTide(rcvEcef, sun, moon) : [0, 0, 0];
+      if (cfg.oceanLoading) {
+        const [e, n, u] = oceanLoadingEnu(cfg.oceanLoading, timeMs);
+        const [lat, lon] = ecefToGeodetic(rcvEcef[0], rcvEcef[1], rcvEcef[2]);
+        const d = enuToEcef([e, n, u], lat, lon);
+        tideDisp[0] += d[0];
+        tideDisp[1] += d[1];
+        tideDisp[2] += d[2];
+      }
       return { timeMs, sun, moon, tideDisp };
     },
     evaluate(geom, ctx) {
@@ -3615,7 +3723,7 @@ function createPppCorrections(cfg) {
           rangeM += interpPcv(off, zenithDeg);
         }
       }
-      if (useTides && ctx) {
+      if (ctx) {
         rangeM -= dot(e, ctx.tideDisp);
       }
       if (useWindup && frame) {
@@ -4001,7 +4109,7 @@ var PppEngine = class {
     for (const ci of clkIdx.values()) P[ci][ci] = CLK_VAR;
     P[IDX_ZWD][IDX_ZWD] += ztdQ * (ei > 0 ? 1 : 0);
     const rcv = [x[IDX_X], x[IDX_Y], x[IDX_Z]];
-    const [latRad, , heightM] = ecefToGeodetic(rcv[0], rcv[1], rcv[2]);
+    const [latRad, lonRad, heightM] = ecefToGeodetic(rcv[0], rcv[1], rcv[2]);
     const epochCtx = corrections?.epochContext?.(epoch.timeMs, rcv) ?? null;
     const vis = [];
     let epochZhd = 0;
@@ -4034,7 +4142,17 @@ var PppEngine = class {
       const pIf = g * o.c1 - (g - 1) * o.c2;
       const lIf = g * (o.l1 * lam1) - (g - 1) * (o.l2 * lam2);
       if (!Number.isFinite(pIf) || !Number.isFinite(lIf)) continue;
-      const trop = opts.troposphere ?? true ? niellMapping(elRad, latRad, heightM, epoch.timeMs) : { zhd: 0, mh: 0, mw: 0 };
+      const trop = opts.troposphere ?? true ? opts.troposphereModel ? opts.troposphereModel.delays(
+        epoch.timeMs,
+        latRad,
+        lonRad,
+        heightM,
+        elRad,
+        azRad
+      ) : {
+        ...niellMapping(elRad, latRad, heightM, epoch.timeMs),
+        gradient: 0
+      } : { zhd: 0, mh: 0, mw: 0, gradient: 0 };
       const { zhd, mh, mw } = trop;
       epochZhd = zhd;
       const geom = {
@@ -4059,7 +4177,7 @@ var PppEngine = class {
         lIf,
         satClkM: sat.state.clkM,
         range,
-        tropoFixed: zhd * mh + corr.rangeM,
+        tropoFixed: zhd * mh + trop.gradient + corr.rangeM,
         mw,
         windupM: corr.phaseWindupM,
         slip,
@@ -4984,6 +5102,7 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
 }
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
+  BLQ_CONSTITUENTS,
   IONO_L1_M_PER_TECU,
   PppEngine,
   RtkFloatEngine,
@@ -5011,6 +5130,8 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   moonEcef,
   niellMapping,
   nlWavelength,
+  oceanLoadingEnu,
+  parseBlq,
   postProcessRtk,
   rcvAntennaRangeM,
   resolveMsmRanges,

package/dist/positioning.d.cts

@@ -969,6 +969,36 @@ declare function applyOrbitClock(eph: Ephemeris, corr: OrbitClockCorrection, tim
  */
 declare function hasToOrbitClock(orbit: HasOrbitCorrection, clock: HasClockCorrection | undefined, refTimeMs: number, validitySeconds: number | null): OrbitClockCorrection | null;
 
+/**
+ * Ocean tide loading: the crust deforms under the weight of the ocean
+ * tides, up to several centimetres vertically near tidal coasts. Given a
+ * station's BLQ coefficients (amplitude and phase of 11 tidal constituents
+ * for the radial, west and south displacement — the format of Onsala's
+ * ocean loading service), the displacement at any instant is the sum of the
+ * constituents with their astronomical arguments (IERS ARG.f; the method of
+ * RTKLIB's tide_oload, without nodal corrections).
+ */
+/** Constituent order of BLQ files. */
+declare const BLQ_CONSTITUENTS: readonly ["M2", "S2", "N2", "K2", "K1", "O1", "P1", "Q1", "Mf", "Mm", "Ssa"];
+/** One station's BLQ coefficients. */
+interface OceanLoadingCoefficients {
+    /** Amplitudes (m) [radial, west, south] per constituent (BLQ order). */
+    amplitude: [number, number, number][];
+    /** Phases (degrees, Greenwich lag) [radial, west, south] per constituent. */
+    phase: [number, number, number][];
+}
+/**
+ * Ocean loading displacement (m) at `utcMs`, local [East, North, Up].
+ * UTC stands in for UT1 (≤ 0.9 s, negligible for tides).
+ */
+declare function oceanLoadingEnu(c: OceanLoadingCoefficients, utcMs: number): [number, number, number];
+/**
+ * Read the coefficients of every station in a BLQ file (Onsala format:
+ * `$$` comments, a station line, then 3 amplitude rows and 3 phase rows of
+ * 11 values). Station names are upper-cased, keyed by their first word.
+ */
+declare function parseBlq(text: string): Map<string, OceanLoadingCoefficients>;
+
 /**
  * Optional PPP range corrections: satellite & receiver antenna phase-centre
  * offsets, phase wind-up, and solid-earth tides. These take a converged
@@ -1029,6 +1059,8 @@ interface PppCorrectionConfig {
     rcvPco?: boolean;
     tides?: boolean;
     windup?: boolean;
+    /** The station's ocean tide loading coefficients (BLQ), when known. */
+    oceanLoading?: OceanLoadingCoefficients;
 }
 /** Wire satellite/receiver antenna + tides + wind-up into a PppCorrections. */
 declare function createPppCorrections(cfg: PppCorrectionConfig): PppCorrections;
@@ -1115,6 +1147,10 @@ interface PppOptions {
     ztdProcessNoise?: number;
     /** Model + estimate the troposphere. Default true. */
     troposphere?: boolean;
+    /** A-priori troposphere model (zenith hydrostatic delay, mapping
+     *  functions, gradients) — e.g. GPT3 + VMF3 from `gnss-js/troposphere`.
+     *  Default: Saastamoinen with the standard atmosphere + Niell. */
+    troposphereModel?: PppTroposphereModel;
     /** Expose the final EKF state (position + ambiguities + covariance) in the
      * solution for ambiguity-resolved positioning. Default false. */
     exposeState?: boolean;
@@ -1155,6 +1191,16 @@ interface PppOptions {
         maxChi2?: number;
     };
 }
+/** A troposphere model the PPP engine can use (structurally the
+ *  `TroposphereModel` of `gnss-js/troposphere`). */
+interface PppTroposphereModel {
+    delays(utcMs: number, latRad: number, lonRad: number, hEll: number, elRad: number, azRad: number): {
+        zhd: number;
+        mh: number;
+        mw: number;
+        gradient: number;
+    };
+}
 interface PppEpochResult {
     timeMs: number;
     /** ENU error vs ground truth (m), if truth given. */
@@ -1842,4 +1888,4 @@ declare function sbasProtectionLevels(sats: readonly {
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type MsmResolveOptions, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, BLQ_CONSTITUENTS, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type MsmResolveOptions, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OceanLoadingCoefficients, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type PppTroposphereModel, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, oceanLoadingEnu, parseBlq, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -969,6 +969,36 @@ declare function applyOrbitClock(eph: Ephemeris, corr: OrbitClockCorrection, tim
  */
 declare function hasToOrbitClock(orbit: HasOrbitCorrection, clock: HasClockCorrection | undefined, refTimeMs: number, validitySeconds: number | null): OrbitClockCorrection | null;
 
+/**
+ * Ocean tide loading: the crust deforms under the weight of the ocean
+ * tides, up to several centimetres vertically near tidal coasts. Given a
+ * station's BLQ coefficients (amplitude and phase of 11 tidal constituents
+ * for the radial, west and south displacement — the format of Onsala's
+ * ocean loading service), the displacement at any instant is the sum of the
+ * constituents with their astronomical arguments (IERS ARG.f; the method of
+ * RTKLIB's tide_oload, without nodal corrections).
+ */
+/** Constituent order of BLQ files. */
+declare const BLQ_CONSTITUENTS: readonly ["M2", "S2", "N2", "K2", "K1", "O1", "P1", "Q1", "Mf", "Mm", "Ssa"];
+/** One station's BLQ coefficients. */
+interface OceanLoadingCoefficients {
+    /** Amplitudes (m) [radial, west, south] per constituent (BLQ order). */
+    amplitude: [number, number, number][];
+    /** Phases (degrees, Greenwich lag) [radial, west, south] per constituent. */
+    phase: [number, number, number][];
+}
+/**
+ * Ocean loading displacement (m) at `utcMs`, local [East, North, Up].
+ * UTC stands in for UT1 (≤ 0.9 s, negligible for tides).
+ */
+declare function oceanLoadingEnu(c: OceanLoadingCoefficients, utcMs: number): [number, number, number];
+/**
+ * Read the coefficients of every station in a BLQ file (Onsala format:
+ * `$$` comments, a station line, then 3 amplitude rows and 3 phase rows of
+ * 11 values). Station names are upper-cased, keyed by their first word.
+ */
+declare function parseBlq(text: string): Map<string, OceanLoadingCoefficients>;
+
 /**
  * Optional PPP range corrections: satellite & receiver antenna phase-centre
  * offsets, phase wind-up, and solid-earth tides. These take a converged
@@ -1029,6 +1059,8 @@ interface PppCorrectionConfig {
     rcvPco?: boolean;
     tides?: boolean;
     windup?: boolean;
+    /** The station's ocean tide loading coefficients (BLQ), when known. */
+    oceanLoading?: OceanLoadingCoefficients;
 }
 /** Wire satellite/receiver antenna + tides + wind-up into a PppCorrections. */
 declare function createPppCorrections(cfg: PppCorrectionConfig): PppCorrections;
@@ -1115,6 +1147,10 @@ interface PppOptions {
     ztdProcessNoise?: number;
     /** Model + estimate the troposphere. Default true. */
     troposphere?: boolean;
+    /** A-priori troposphere model (zenith hydrostatic delay, mapping
+     *  functions, gradients) — e.g. GPT3 + VMF3 from `gnss-js/troposphere`.
+     *  Default: Saastamoinen with the standard atmosphere + Niell. */
+    troposphereModel?: PppTroposphereModel;
     /** Expose the final EKF state (position + ambiguities + covariance) in the
      * solution for ambiguity-resolved positioning. Default false. */
     exposeState?: boolean;
@@ -1155,6 +1191,16 @@ interface PppOptions {
         maxChi2?: number;
     };
 }
+/** A troposphere model the PPP engine can use (structurally the
+ *  `TroposphereModel` of `gnss-js/troposphere`). */
+interface PppTroposphereModel {
+    delays(utcMs: number, latRad: number, lonRad: number, hEll: number, elRad: number, azRad: number): {
+        zhd: number;
+        mh: number;
+        mw: number;
+        gradient: number;
+    };
+}
 interface PppEpochResult {
     timeMs: number;
     /** ENU error vs ground truth (m), if truth given. */
@@ -1842,4 +1888,4 @@ declare function sbasProtectionLevels(sats: readonly {
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type MsmResolveOptions, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type AppliedOrbitClock, type ArFixedSat, type ArOptions, type ArResult, type ArSat, BLQ_CONSTITUENTS, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type MsmResolveOptions, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type OceanLoadingCoefficients, type OrbitClockCorrection, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, PppEngine, type PppEphemerisSource, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type PppTroposphereModel, type RawObservation, type ReceiverBias, type ReceiverClockSample, type RtkAntennaConfig, type RtkAntennaModel, type RtkAntennaOffset, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkReceiverAntenna, type RtkTrackPoint, type SatGeom, type SatSample, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, Sp3EphemerisSource, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, ZeroBaselineEngine, type ZeroBaselineOptions, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, oceanLoadingEnu, parseBlq, postProcessRtk, rcvAntennaRangeM, resolveMsmRanges, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -1,4 +1,5 @@
 import {
+  BLQ_CONSTITUENTS,
   IONO_L1_M_PER_TECU,
   PppEngine,
   RtkFloatEngine,
@@ -26,6 +27,8 @@ import {
   moonEcef,
   niellMapping,
   nlWavelength,
+  oceanLoadingEnu,
+  parseBlq,
   postProcessRtk,
   rcvAntennaRangeM,
   resolveMsmRanges,
@@ -41,7 +44,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-PQLPPFIK.js";
+} from "./chunk-TETPRPZ2.js";
 import "./chunk-H3AKFD25.js";
 import "./chunk-7T5DLVZL.js";
 import "./chunk-HVXYFUCB.js";
@@ -54,6 +57,7 @@ import "./chunk-HUXQABJI.js";
 import "./chunk-V2D5KUKH.js";
 import "./chunk-FIEWO4J4.js";
 export {
+  BLQ_CONSTITUENTS,
   IONO_L1_M_PER_TECU,
   PppEngine,
   RtkFloatEngine,
@@ -81,6 +85,8 @@ export {
   moonEcef,
   niellMapping,
   nlWavelength,
+  oceanLoadingEnu,
+  parseBlq,
   postProcessRtk,
   rcvAntennaRangeM,
   resolveMsmRanges,