gnss-js 1.27.0 → 1.28.0

package/dist/positioning.d.cts

@@ -1,6 +1,7 @@
 import { E as Ephemeris } from './nav-DImXUdbp.cjs';
 import { DopValues } from './orbit.cjs';
-import { I as IonexGrid } from './ionex-DaW5x9nq.cjs';
+import { I as IonexGrid, S as Sp3File } from './sp3-BAq8l39b.cjs';
+import { AntexFile } from './antex.cjs';
 import './ephemeris-Ohl6NAfw.cjs';
 import './ecef-CF0uAysr.cjs';
 import './nav-BYmYoh9A.cjs';
@@ -480,6 +481,234 @@ declare function lambdaSearch(a: Float64Array, Q: Float64Array, n: number, m?: n
  */
 declare function lambdaReduction(Q: Float64Array, n: number): Float64Array | null;
 
+/**
+ * ANTEX phase-centre application for PPP: satellite and receiver antenna
+ * PCO (offset) + PCV (elevation-dependent variation), ionosphere-free
+ * combined, with the nominal yaw-steering satellite attitude used to
+ * rotate the satellite PCO from the body frame into ECEF.
+ */
+
+interface AntennaOffset {
+    /** IF-combined PCO (m): satellite → body-frame [X,Y,Z]; receiver → ENU. */
+    pco: [number, number, number];
+    /** IF-combined no-azimuth PCV samples (m) over zenith/nadir angle. */
+    pcvZen1Deg: number;
+    pcvDzenDeg: number;
+    pcvNoazi: number[];
+}
+interface PppAntennaModel {
+    satOffset(prn: string, dateMs: number, band1: string, band2: string, f1: number, f2: number): AntennaOffset | null;
+    rcvOffset(antType: string, band1: string, band2: string, f1: number, f2: number): AntennaOffset | null;
+}
+declare function buildPppAntenna(antex: AntexFile): PppAntennaModel;
+
+/**
+ * Optional PPP range corrections: satellite & receiver antenna phase-centre
+ * offsets, phase wind-up, and solid-earth tides. These take a converged
+ * float PPP solution from ~decimetre toward centimetre.
+ *
+ * Everything here is opt-in via `PppOptions.corrections`; the core solver
+ * runs without it. Corrections needing Sun/Moon geometry are computed once
+ * per epoch in `epochContext` and reused across satellites.
+ */
+/** Per-satellite geometry handed to the correction models. */
+interface SatGeom {
+    prn: string;
+    satEcef: [number, number, number];
+    rcvEcef: [number, number, number];
+    /** Unit line of sight, receiver → satellite. */
+    los: [number, number, number];
+    elRad: number;
+    azRad: number;
+    f1: number;
+    f2: number;
+    /** ANTEX frequency codes for the two bands (e.g. 'G01','G02'). */
+    band1?: string;
+    band2?: string;
+    /** Iono-free coefficient g = f1²/(f1²−f2²). */
+    g: number;
+}
+/** Epoch-level context (Sun/Moon positions, tide displacement). */
+interface PppEpochContext {
+    timeMs: number;
+    /** Sun position, ECEF (m). */
+    sun: [number, number, number];
+    /** Moon position, ECEF (m). */
+    moon: [number, number, number];
+    /** Solid-earth-tide receiver displacement, ECEF (m). */
+    tideDisp: [number, number, number];
+}
+interface CorrectionResult {
+    /** Additive range correction applied to BOTH code and phase (m). */
+    rangeM: number;
+    /** Additional phase-only correction (wind-up), metres. */
+    phaseWindupM: number;
+}
+/** Correction configuration + hooks. Supplied by the caller. */
+interface PppCorrections {
+    /** Build the per-epoch Sun/Moon/tide context. */
+    epochContext?: (timeMs: number, rcvEcef: [number, number, number]) => PppEpochContext;
+    /** Per-satellite correction evaluator. */
+    evaluate?: (geom: SatGeom, ctx: PppEpochContext | null) => CorrectionResult;
+}
+interface PppCorrectionConfig {
+    antenna: PppAntennaModel;
+    /** Receiver antenna type (RINEX header "ANT # / TYPE"). */
+    rcvAntType: string;
+    /** Marker→ARP offset in local ENU (m). RINEX "ANTENNA: DELTA H/E/N"
+     * is [Up, East, North] → pass as [E, N, U]. */
+    antDeltaEnu?: [number, number, number];
+    satPco?: boolean;
+    rcvPco?: boolean;
+    tides?: boolean;
+    windup?: boolean;
+}
+/** Wire satellite/receiver antenna + tides + wind-up into a PppCorrections. */
+declare function createPppCorrections(cfg: PppCorrectionConfig): PppCorrections;
+
+/**
+ * Static float Precise Point Positioning (PPP).
+ *
+ * Dual-frequency ionosphere-free code + carrier phase, precise satellite
+ * orbits and clocks (SP3), a forward Extended Kalman Filter estimating
+ * [position, receiver clock, zenith wet troposphere, one float ambiguity
+ * per satellite arc]. No base station — absolute cm–dm coordinates from a
+ * single observation file.
+ *
+ * Corrections applied (see PppOptions to toggle the optional ones):
+ *  - precise satellite orbit + clock (SP3, Lagrange / linear)
+ *  - periodic relativistic satellite-clock term (−2·r·v/c²)
+ *  - Earth-rotation (Sagnac) during signal transit
+ *  - dry troposphere (Saastamoinen ZHD, a priori) + estimated wet ZWD,
+ *    mapped by the Niell mapping function
+ *  - elevation-dependent weighting, cycle-slip re-initialisation
+ *  - (optional, added by the caller supplying products) satellite &
+ *    receiver antenna phase-centre offsets, phase wind-up, solid-earth tides
+ *
+ * The ionosphere-free combination removes the first-order ionosphere; the
+ * float ambiguity absorbs the phase biases and the integer cycles.
+ */
+
+/** One satellite's dual-frequency observation at an epoch. */
+interface PppSatObs {
+    prn: string;
+    /** Band-1 / band-2 centre frequencies (Hz). */
+    f1: number;
+    f2: number;
+    /** Pseudorange on band 1 / band 2 (metres). */
+    c1: number;
+    c2: number;
+    /** Carrier phase on band 1 / band 2 (cycles). */
+    l1: number;
+    l2: number;
+    /** ANTEX frequency codes for the two bands (e.g. 'G01','G02'). Used by
+     * the antenna corrections; optional if corrections are off. */
+    band1?: string;
+    band2?: string;
+    /** True if a cycle slip is flagged on either band (LLI or detected). */
+    slip: boolean;
+}
+interface PppEpoch {
+    /** GPS-scale epoch time (ms). */
+    timeMs: number;
+    obs: PppSatObs[];
+}
+interface PppOptions {
+    /** A priori receiver ECEF position (m). Required — from the RINEX header
+     * approx position or an SPP solution. */
+    aprioriPos: [number, number, number];
+    /** Ground-truth ECEF (m) for reporting ENU error in the series. */
+    groundTruth?: [number, number, number];
+    /** Elevation cutoff (deg). Default 10. */
+    elevationMaskDeg?: number;
+    /** Optional corrections (antenna PCO/PCV, wind-up, tides). */
+    corrections?: PppCorrections;
+    /** IF code σ at zenith (m). Default 3.0 (already IF-inflated). */
+    codeSigma?: number;
+    /** IF phase σ at zenith (m). Default 0.01. */
+    phaseSigma?: number;
+    /** Wet-troposphere random-walk process noise (m²/epoch). Default 1e-8. */
+    ztdProcessNoise?: number;
+    /** Model + estimate the troposphere. Default true. */
+    troposphere?: boolean;
+}
+interface PppEpochResult {
+    timeMs: number;
+    /** ENU error vs ground truth (m), if truth given. */
+    enu: [number, number, number] | null;
+    /** 3D error vs ground truth (m), if truth given. */
+    error3d: number | null;
+    position: [number, number, number];
+    nSats: number;
+    /** Estimated zenith wet delay (m). */
+    ztdWet: number;
+    /** RMS of post-fit phase residuals this epoch (m). */
+    phaseResRms: number;
+}
+interface PppSolution {
+    /** Final estimated ECEF position (m). */
+    position: [number, number, number];
+    /** Final geodetic [latDeg, lonDeg, heightM]. */
+    llh: [number, number, number];
+    /** Per-epoch series (convergence). */
+    series: PppEpochResult[];
+    /** Estimated zenith wet delay at the end (m). */
+    ztdWet: number;
+    /** Epochs processed. */
+    epochsUsed: number;
+    /** Seconds from first epoch until 3D error first stays < 0.1 m
+     * (null if never, or no ground truth). */
+    convergenceSec: number | null;
+    /** Final 3D error vs ground truth (m), if given. */
+    finalError3d: number | null;
+}
+declare function solvePpp(epochs: PppEpoch[], sp3: Sp3File, opts: PppOptions): PppSolution;
+
+/**
+ * Troposphere for PPP: a-priori Saastamoinen zenith hydrostatic delay
+ * (ZHD) and the Niell (NMF) hydrostatic + wet mapping functions.
+ *
+ * The estimator carries the zenith WET delay as a state; the dry (ZHD)
+ * part is fixed from the standard atmosphere. Niell 1996 ("Global mapping
+ * functions for the atmosphere delay at radio wavelengths", JGR 101).
+ */
+interface TropoMapping {
+    /** A-priori zenith hydrostatic delay (m). */
+    zhd: number;
+    /** Hydrostatic mapping factor at the elevation. */
+    mh: number;
+    /** Wet mapping factor at the elevation. */
+    mw: number;
+}
+/**
+ * Saastamoinen ZHD + Niell hydrostatic/wet mapping at an elevation.
+ * @param elRad elevation (rad), @param latRad geodetic latitude (rad),
+ * @param heightM ellipsoidal height (m), @param timeMs epoch (GPS ms).
+ */
+declare function niellMapping(elRad: number, latRad: number, heightM: number, timeMs: number): TropoMapping;
+
+/**
+ * Low-precision Sun/Moon ephemerides, sidereal time, and the IERS degree-2
+ * solid-earth tide — enough for the PPP corrections (satellite attitude,
+ * solid-earth tides, phase wind-up) at the millimetre–centimetre level.
+ *
+ * Sun ~0.01°, Moon ~0.1–0.3° accuracy; solid-earth-tide step-1 (in-phase)
+ * displacement good to ~mm. Precession/nutation/polar-motion are neglected
+ * (their effect on the Sun/Moon *direction* used here is well below the
+ * correction accuracy needed).
+ */
+/** Greenwich mean sidereal time (rad). */
+declare function gmst(timeMs: number): number;
+/** Sun position in ECEF (m). Low-precision almanac (Meeus). */
+declare function sunEcef(timeMs: number): [number, number, number];
+/** Moon position in ECEF (m). Low-precision series (Montenbruck & Gill). */
+declare function moonEcef(timeMs: number): [number, number, number];
+/**
+ * Solid-earth-tide displacement of a station (ECEF, m). IERS conventions
+ * step-1 (degree 2, in-phase) — the dominant few-cm/dm periodic term.
+ */
+declare function solidEarthTide(rcv: [number, number, number], sun: [number, number, number], moon: [number, number, number]): [number, number, number];
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -569,4 +798,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 DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SppOptions, type SppSolution, gimSlantIonoDelayL1, gimVerticalTec, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, satClockCorrection, solveDgnss, solveSpp, toRtkEpoch };
+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 };

package/dist/positioning.d.ts

@@ -1,6 +1,7 @@
 import { E as Ephemeris } from './nav-DImXUdbp.js';
 import { DopValues } from './orbit.js';
-import { I as IonexGrid } from './ionex-DaW5x9nq.js';
+import { I as IonexGrid, S as Sp3File } from './sp3-BAq8l39b.js';
+import { AntexFile } from './antex.js';
 import './ephemeris-Ohl6NAfw.js';
 import './ecef-CF0uAysr.js';
 import './nav-CN8rxL4h.js';
@@ -480,6 +481,234 @@ declare function lambdaSearch(a: Float64Array, Q: Float64Array, n: number, m?: n
  */
 declare function lambdaReduction(Q: Float64Array, n: number): Float64Array | null;
 
+/**
+ * ANTEX phase-centre application for PPP: satellite and receiver antenna
+ * PCO (offset) + PCV (elevation-dependent variation), ionosphere-free
+ * combined, with the nominal yaw-steering satellite attitude used to
+ * rotate the satellite PCO from the body frame into ECEF.
+ */
+
+interface AntennaOffset {
+    /** IF-combined PCO (m): satellite → body-frame [X,Y,Z]; receiver → ENU. */
+    pco: [number, number, number];
+    /** IF-combined no-azimuth PCV samples (m) over zenith/nadir angle. */
+    pcvZen1Deg: number;
+    pcvDzenDeg: number;
+    pcvNoazi: number[];
+}
+interface PppAntennaModel {
+    satOffset(prn: string, dateMs: number, band1: string, band2: string, f1: number, f2: number): AntennaOffset | null;
+    rcvOffset(antType: string, band1: string, band2: string, f1: number, f2: number): AntennaOffset | null;
+}
+declare function buildPppAntenna(antex: AntexFile): PppAntennaModel;
+
+/**
+ * Optional PPP range corrections: satellite & receiver antenna phase-centre
+ * offsets, phase wind-up, and solid-earth tides. These take a converged
+ * float PPP solution from ~decimetre toward centimetre.
+ *
+ * Everything here is opt-in via `PppOptions.corrections`; the core solver
+ * runs without it. Corrections needing Sun/Moon geometry are computed once
+ * per epoch in `epochContext` and reused across satellites.
+ */
+/** Per-satellite geometry handed to the correction models. */
+interface SatGeom {
+    prn: string;
+    satEcef: [number, number, number];
+    rcvEcef: [number, number, number];
+    /** Unit line of sight, receiver → satellite. */
+    los: [number, number, number];
+    elRad: number;
+    azRad: number;
+    f1: number;
+    f2: number;
+    /** ANTEX frequency codes for the two bands (e.g. 'G01','G02'). */
+    band1?: string;
+    band2?: string;
+    /** Iono-free coefficient g = f1²/(f1²−f2²). */
+    g: number;
+}
+/** Epoch-level context (Sun/Moon positions, tide displacement). */
+interface PppEpochContext {
+    timeMs: number;
+    /** Sun position, ECEF (m). */
+    sun: [number, number, number];
+    /** Moon position, ECEF (m). */
+    moon: [number, number, number];
+    /** Solid-earth-tide receiver displacement, ECEF (m). */
+    tideDisp: [number, number, number];
+}
+interface CorrectionResult {
+    /** Additive range correction applied to BOTH code and phase (m). */
+    rangeM: number;
+    /** Additional phase-only correction (wind-up), metres. */
+    phaseWindupM: number;
+}
+/** Correction configuration + hooks. Supplied by the caller. */
+interface PppCorrections {
+    /** Build the per-epoch Sun/Moon/tide context. */
+    epochContext?: (timeMs: number, rcvEcef: [number, number, number]) => PppEpochContext;
+    /** Per-satellite correction evaluator. */
+    evaluate?: (geom: SatGeom, ctx: PppEpochContext | null) => CorrectionResult;
+}
+interface PppCorrectionConfig {
+    antenna: PppAntennaModel;
+    /** Receiver antenna type (RINEX header "ANT # / TYPE"). */
+    rcvAntType: string;
+    /** Marker→ARP offset in local ENU (m). RINEX "ANTENNA: DELTA H/E/N"
+     * is [Up, East, North] → pass as [E, N, U]. */
+    antDeltaEnu?: [number, number, number];
+    satPco?: boolean;
+    rcvPco?: boolean;
+    tides?: boolean;
+    windup?: boolean;
+}
+/** Wire satellite/receiver antenna + tides + wind-up into a PppCorrections. */
+declare function createPppCorrections(cfg: PppCorrectionConfig): PppCorrections;
+
+/**
+ * Static float Precise Point Positioning (PPP).
+ *
+ * Dual-frequency ionosphere-free code + carrier phase, precise satellite
+ * orbits and clocks (SP3), a forward Extended Kalman Filter estimating
+ * [position, receiver clock, zenith wet troposphere, one float ambiguity
+ * per satellite arc]. No base station — absolute cm–dm coordinates from a
+ * single observation file.
+ *
+ * Corrections applied (see PppOptions to toggle the optional ones):
+ *  - precise satellite orbit + clock (SP3, Lagrange / linear)
+ *  - periodic relativistic satellite-clock term (−2·r·v/c²)
+ *  - Earth-rotation (Sagnac) during signal transit
+ *  - dry troposphere (Saastamoinen ZHD, a priori) + estimated wet ZWD,
+ *    mapped by the Niell mapping function
+ *  - elevation-dependent weighting, cycle-slip re-initialisation
+ *  - (optional, added by the caller supplying products) satellite &
+ *    receiver antenna phase-centre offsets, phase wind-up, solid-earth tides
+ *
+ * The ionosphere-free combination removes the first-order ionosphere; the
+ * float ambiguity absorbs the phase biases and the integer cycles.
+ */
+
+/** One satellite's dual-frequency observation at an epoch. */
+interface PppSatObs {
+    prn: string;
+    /** Band-1 / band-2 centre frequencies (Hz). */
+    f1: number;
+    f2: number;
+    /** Pseudorange on band 1 / band 2 (metres). */
+    c1: number;
+    c2: number;
+    /** Carrier phase on band 1 / band 2 (cycles). */
+    l1: number;
+    l2: number;
+    /** ANTEX frequency codes for the two bands (e.g. 'G01','G02'). Used by
+     * the antenna corrections; optional if corrections are off. */
+    band1?: string;
+    band2?: string;
+    /** True if a cycle slip is flagged on either band (LLI or detected). */
+    slip: boolean;
+}
+interface PppEpoch {
+    /** GPS-scale epoch time (ms). */
+    timeMs: number;
+    obs: PppSatObs[];
+}
+interface PppOptions {
+    /** A priori receiver ECEF position (m). Required — from the RINEX header
+     * approx position or an SPP solution. */
+    aprioriPos: [number, number, number];
+    /** Ground-truth ECEF (m) for reporting ENU error in the series. */
+    groundTruth?: [number, number, number];
+    /** Elevation cutoff (deg). Default 10. */
+    elevationMaskDeg?: number;
+    /** Optional corrections (antenna PCO/PCV, wind-up, tides). */
+    corrections?: PppCorrections;
+    /** IF code σ at zenith (m). Default 3.0 (already IF-inflated). */
+    codeSigma?: number;
+    /** IF phase σ at zenith (m). Default 0.01. */
+    phaseSigma?: number;
+    /** Wet-troposphere random-walk process noise (m²/epoch). Default 1e-8. */
+    ztdProcessNoise?: number;
+    /** Model + estimate the troposphere. Default true. */
+    troposphere?: boolean;
+}
+interface PppEpochResult {
+    timeMs: number;
+    /** ENU error vs ground truth (m), if truth given. */
+    enu: [number, number, number] | null;
+    /** 3D error vs ground truth (m), if truth given. */
+    error3d: number | null;
+    position: [number, number, number];
+    nSats: number;
+    /** Estimated zenith wet delay (m). */
+    ztdWet: number;
+    /** RMS of post-fit phase residuals this epoch (m). */
+    phaseResRms: number;
+}
+interface PppSolution {
+    /** Final estimated ECEF position (m). */
+    position: [number, number, number];
+    /** Final geodetic [latDeg, lonDeg, heightM]. */
+    llh: [number, number, number];
+    /** Per-epoch series (convergence). */
+    series: PppEpochResult[];
+    /** Estimated zenith wet delay at the end (m). */
+    ztdWet: number;
+    /** Epochs processed. */
+    epochsUsed: number;
+    /** Seconds from first epoch until 3D error first stays < 0.1 m
+     * (null if never, or no ground truth). */
+    convergenceSec: number | null;
+    /** Final 3D error vs ground truth (m), if given. */
+    finalError3d: number | null;
+}
+declare function solvePpp(epochs: PppEpoch[], sp3: Sp3File, opts: PppOptions): PppSolution;
+
+/**
+ * Troposphere for PPP: a-priori Saastamoinen zenith hydrostatic delay
+ * (ZHD) and the Niell (NMF) hydrostatic + wet mapping functions.
+ *
+ * The estimator carries the zenith WET delay as a state; the dry (ZHD)
+ * part is fixed from the standard atmosphere. Niell 1996 ("Global mapping
+ * functions for the atmosphere delay at radio wavelengths", JGR 101).
+ */
+interface TropoMapping {
+    /** A-priori zenith hydrostatic delay (m). */
+    zhd: number;
+    /** Hydrostatic mapping factor at the elevation. */
+    mh: number;
+    /** Wet mapping factor at the elevation. */
+    mw: number;
+}
+/**
+ * Saastamoinen ZHD + Niell hydrostatic/wet mapping at an elevation.
+ * @param elRad elevation (rad), @param latRad geodetic latitude (rad),
+ * @param heightM ellipsoidal height (m), @param timeMs epoch (GPS ms).
+ */
+declare function niellMapping(elRad: number, latRad: number, heightM: number, timeMs: number): TropoMapping;
+
+/**
+ * Low-precision Sun/Moon ephemerides, sidereal time, and the IERS degree-2
+ * solid-earth tide — enough for the PPP corrections (satellite attitude,
+ * solid-earth tides, phase wind-up) at the millimetre–centimetre level.
+ *
+ * Sun ~0.01°, Moon ~0.1–0.3° accuracy; solid-earth-tide step-1 (in-phase)
+ * displacement good to ~mm. Precession/nutation/polar-motion are neglected
+ * (their effect on the Sun/Moon *direction* used here is well below the
+ * correction accuracy needed).
+ */
+/** Greenwich mean sidereal time (rad). */
+declare function gmst(timeMs: number): number;
+/** Sun position in ECEF (m). Low-precision almanac (Meeus). */
+declare function sunEcef(timeMs: number): [number, number, number];
+/** Moon position in ECEF (m). Low-precision series (Montenbruck & Gill). */
+declare function moonEcef(timeMs: number): [number, number, number];
+/**
+ * Solid-earth-tide displacement of a station (ECEF, m). IERS conventions
+ * step-1 (degree 2, in-phase) — the dominant few-cm/dm periodic term.
+ */
+declare function solidEarthTide(rcv: [number, number, number], sun: [number, number, number], moon: [number, number, number]): [number, number, number];
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -569,4 +798,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 DgnssOptions, type DgnssSolution, type EphemerisSource, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type SppOptions, type SppSolution, gimSlantIonoDelayL1, gimVerticalTec, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, satClockCorrection, solveDgnss, solveSpp, toRtkEpoch };
+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 };

package/dist/positioning.js

@@ -1,34 +1,51 @@
 import {
   IONO_L1_M_PER_TECU,
   RtkFloatEngine,
+  buildPppAntenna,
+  createPppCorrections,
   gimSlantIonoDelayL1,
   gimVerticalTec,
+  gmst,
   ionoFree,
   klobucharDelay,
   lambdaReduction,
   lambdaSearch,
+  moonEcef,
+  niellMapping,
   satClockCorrection,
+  solidEarthTide,
   solveDgnss,
+  solvePpp,
   solveSpp,
+  sunEcef,
   toRtkEpoch
-} from "./chunk-4XGKCVE3.js";
+} from "./chunk-O65WA6CM.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
+import "./chunk-4KMWWCUR.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";
 import "./chunk-FIEWO4J4.js";
 export {
   IONO_L1_M_PER_TECU,
   RtkFloatEngine,
+  buildPppAntenna,
+  createPppCorrections,
   gimSlantIonoDelayL1,
   gimVerticalTec,
+  gmst,
   ionoFree,
   klobucharDelay,
   lambdaReduction,
   lambdaSearch,
+  moonEcef,
+  niellMapping,
   satClockCorrection,
+  solidEarthTide,
   solveDgnss,
+  solvePpp,
   solveSpp,
+  sunEcef,
   toRtkEpoch
 };

package/dist/rinex.d.cts

@@ -2,7 +2,7 @@ import { R as RinexHeader } from './parser-JPjjFgeP.cjs';
 export { E as EpochSummary, P as ParseOptions, 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';
 import { a as NavResult, N as NavHeader, E as Ephemeris, R as RinexCnavEphemeris } from './nav-DImXUdbp.cjs';
 export { G as GlonassEphemeris, K as KeplerEphemeris, p as parseNavFile } from './nav-DImXUdbp.cjs';
-export { I as IonexGrid, p as parseIonex } from './ionex-DaW5x9nq.cjs';
+export { I as IonexGrid, S as Sp3File, a as Sp3Sample, p as parseIonex, b as parseSp3, s as sp3Position } from './sp3-BAq8l39b.cjs';
 
 /**
  * Compact RINEX (Hatanaka / CRX) decompression utilities.
@@ -182,53 +182,6 @@ interface Nav4Input {
  */
 declare function writeRinexNav4(nav: Nav4Input): string;
 
-/**
- * SP3-c/-d precise orbit and clock parser, with Lagrange interpolation
- * for evaluating positions between the tabulated epochs.
- *
- * Positions are converted to meters (SP3 tabulates km), clocks to
- * seconds (SP3 tabulates microseconds). Missing positions (all-zero
- * coordinates) and the 999999.999999 bad-clock sentinel become null.
- * Epoch tags are read in the file's own time scale (GPS for IGS/MGEX
- * products) as epoch milliseconds, matching the convention of the
- * RINEX parsers in this library.
- */
-interface Sp3Sample {
-    /** ECEF meters. */
-    x: number;
-    y: number;
-    z: number;
-    /** Satellite clock offset in seconds, or null when flagged bad. */
-    clk: number | null;
-}
-interface Sp3File {
-    version: string;
-    /** Epoch timestamps (ms). */
-    epochs: number[];
-    /** Tabulated interval in seconds (from the ## header line). */
-    intervalSec: number;
-    /** Time system from the %c line, e.g. 'GPS'. */
-    timeSystem: string;
-    /** prn → per-epoch samples (null where the satellite is missing). */
-    satellites: Record<string, (Sp3Sample | null)[]>;
-}
-declare function parseSp3(text: string): Sp3File;
-/**
- * Satellite position at an arbitrary time by Lagrange interpolation
- * over the `order` nearest tabulated samples (default 9, the standard
- * choice for 5–15 min SP3 tables; centimetre-level between nodes).
- * Clock is interpolated linearly between the bracketing samples.
- *
- * Returns null outside the covered span, near data gaps, or when the
- * satellite is absent.
- */
-declare function sp3Position(sp3: Sp3File, prn: string, tMs: number, order?: number): {
-    x: number;
-    y: number;
-    z: number;
-    clk: number | null;
-} | null;
-
 /**
  * RINEX 3.04 / 4.01 observation file writer.
  *
@@ -310,4 +263,4 @@ declare function createGzipLineSink(): {
  */
 declare function stationHeaderLines(header: RinexHeader): string[];
 
-export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, type Sp3File, type Sp3Sample, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseSp3, sp3Position, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };
+export { type CompactEpoch, type CrxField, type DiffState, EMPTY_WARNINGS, Ephemeris, type Nav4Input, NavHeader, NavResult, RinexCnavEphemeris, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, type WriteCrxOptions, createGzipLineSink, crxDecompress, crxRepair, fmtD, fmtF, hdrLine, padL, padR, parseCrxDataLine, stationHeaderLines, writeCrx, writeModernObsBlob, writeRinex2ObsBlob, writeRinex4ObsBlob, writeRinexNav, writeRinexNav4, writeRinexObsBlob };