gnss-js 2.0.0 → 2.1.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, C as ClkFile, S as Sp3File } from './clk-DGVBGRzA.cjs';
+import { a as HasOrbitCorrection, b as HasClockCorrection } from './has-BaXnQarR.cjs';
 import { AntexFile } from './antex.cjs';
 import './ephemeris-BAnBLLB9.cjs';
 import './ecef-CF0uAysr.cjs';
@@ -706,6 +707,77 @@ declare class SbasProcessor {
     private searchIgp;
 }
 
+/**
+ * Apply SSR / HAS orbit + clock corrections to a broadcast ephemeris to obtain
+ * a precise satellite ECEF position and clock offset — the bridge between the
+ * decode-only SSR/HAS decoders (`rtcm3/ssr`, `rtcm3/igs-ssr`, `navbits/has`)
+ * and a PPP/SPP solve.
+ *
+ * Convention (internal, "additive"): the orbit deltas are stored in the
+ * satellite radial / along-track / cross-track (RAC) frame such that
+ *
+ *     correctedPos   = broadcastPos   + R·Δ(radial, along, cross)
+ *     correctedClock = broadcastClock + (c0 + c1·dt + c2·dt²)/c
+ *
+ * This matches Galileo HAS directly (ICD §7.2–7.3). RTCM-SSR and IGS-SSR use
+ * the **opposite** orbit sign, so their adapters negate the orbit deltas when
+ * building an {@link OrbitClockCorrection} (see `ssrToOrbitClock`, added with
+ * the RTCM/IGS wiring). Clock and bias sign conventions are handled per source
+ * in the adapters, not here.
+ */
+
+/**
+ * A normalized per-satellite orbit + clock correction, ready to apply to a
+ * broadcast ephemeris. Additive convention (see file header). Orbit deltas are
+ * in the RAC frame in metres, referenced to the broadcast issue-of-data `iod`
+ * at `refTimeMs`; optional `dot*` give their linear rate (RTCM/IGS SSR carry
+ * these, HAS does not). Clock is a metres polynomial about `refTimeMs`.
+ */
+interface OrbitClockCorrection {
+    prn: string;
+    /** Broadcast IODE (GPS/…) or IODnav (Galileo) the deltas refer to. */
+    iod: number;
+    /** Correction reference epoch, GPS-scale ms (matches `computeSatPosition`). */
+    refTimeMs: number;
+    radial: number;
+    along: number;
+    cross: number;
+    dotRadial?: number;
+    dotAlong?: number;
+    dotCross?: number;
+    /** Clock polynomial (m): c0 + c1·dt + c2·dt². */
+    clockC0: number;
+    clockC1?: number;
+    clockC2?: number;
+    /** Additive high-rate clock term (m), when a separate high-rate stream is combined. */
+    highRateClock?: number;
+    /** Validity window (s) from `refTimeMs`; null = unbounded / unknown. */
+    validitySeconds?: number | null;
+}
+/** A corrected satellite state: ECEF position (m) + the clock offset (s) to ADD
+ *  to the broadcast satellite clock. The caller forms the total clock as
+ *  `satClockCorrection(eph, t) + clockOffsetS`. */
+interface AppliedOrbitClock {
+    x: number;
+    y: number;
+    z: number;
+    clockOffsetS: number;
+}
+/**
+ * Corrected satellite ECEF position + clock offset at emission time `timeMs`,
+ * from a broadcast ephemeris and an SSR/HAS orbit-clock correction. Returns
+ * null when the correction's IOD doesn't match the ephemeris, the correction
+ * has expired, or the orbit can't be computed.
+ */
+declare function applyOrbitClock(eph: Ephemeris, corr: OrbitClockCorrection, timeMs: number): AppliedOrbitClock | null;
+/**
+ * Build a normalized {@link OrbitClockCorrection} from a Galileo-HAS orbit
+ * correction (+ optional matching clock correction). HAS already uses the
+ * additive convention, so no sign flip. Returns null when the orbit deltas are
+ * flagged "not available". `refTimeMs` is the message ToH as GPS-scale ms.
+ */
+declare function hasToOrbitClock(orbit: HasOrbitCorrection, clock: HasClockCorrection | undefined, refTimeMs: number, validitySeconds: number | null): OrbitClockCorrection | null;
+
 /**
  * ANTEX phase-centre application for PPP: satellite and receiver antenna
  * PCO (offset) + PCV (elevation-dependent variation), ionosphere-free
@@ -830,6 +902,10 @@ interface PppSatObs {
      * the antenna corrections; optional if corrections are off. */
     band1?: string;
     band2?: string;
+    /** RINEX-3 tracking codes for the two bands (e.g. '1C','2W'). Used to look
+     * up satellite code (OSB) biases from an SSR/HAS source; optional. */
+    code1?: string;
+    code2?: string;
     /** True if a cycle slip is flagged on either band (LLI or detected). */
     slip: boolean;
 }
@@ -952,7 +1028,90 @@ interface PppArc {
     startMs: number;
     endMs: number;
 }
-declare function solvePpp(epochs: PppEpoch[], sp3: Sp3File, opts: PppOptions): PppSolution;
+/** A raw satellite sample at an emission time: ECEF position + velocity
+ *  (m, m/s) and the clock offset (s) WITHOUT the periodic relativistic term —
+ *  the PPP light-time loop adds relativity + Sagnac. Null when unavailable. */
+interface SatSample {
+    x: number;
+    y: number;
+    z: number;
+    vx: number;
+    vy: number;
+    vz: number;
+    clkS: number;
+}
+/** Supplies precise satellite state to the PPP filter. Implemented by the SP3
+ *  path (default) and by an SSR/HAS source (broadcast eph + corrections), so
+ *  the same filter runs on precise products or a real-time correction stream. */
+interface PppEphemerisSource {
+    satState(prn: string, tEmitMs: number): SatSample | null;
+    /** Optional satellite code (OSB) bias (m) for a RINEX-3 tracking code, to
+     *  make the pseudorange consistent with the (SSR/HAS) clock reference. The
+     *  corrected code is `observed − codeBias`. SP3/CLK sources omit this (the
+     *  float ambiguity + broadcast clock reference absorb it). */
+    codeBias?(prn: string, signalCode: string, tEmitMs: number): number | null;
+}
+/** Default source: SP3 orbits + clocks, with optional high-rate CLK. Velocity
+ *  is a central finite difference (±0.5 s); the clock prefers the high-rate CLK
+ *  and falls back to the SP3 5-minute clock — matching the original solver. */
+declare class Sp3EphemerisSource implements PppEphemerisSource {
+    private readonly sp3;
+    private readonly clk?;
+    constructor(sp3: Sp3File, clk?: ClkFile | undefined);
+    satState(prn: string, tEmitMs: number): SatSample | null;
+}
+/**
+ * Streaming float PPP — the forward EKF as a stateful engine fed one epoch at a
+ * time (mirrors {@link RtkFloatEngine}). `process(epoch)` returns that epoch's
+ * result and carries the filter forward; `solution()` finalises the session.
+ * `solvePpp` (below) is a thin batch wrapper. Same estimator as before — a
+ * broadcast/SP3/SSR source supplies satellite state via {@link PppEphemerisSource}.
+ */
+declare class PppEngine {
+    private readonly source;
+    private readonly opts;
+    private readonly elevMask;
+    private readonly codeSigma;
+    private readonly phaseSigma;
+    private readonly ztdQ;
+    private readonly kinematic;
+    private readonly corrections?;
+    private readonly truth?;
+    private x;
+    private P;
+    private clkIdx;
+    private clockSeeded;
+    private ambIdx;
+    private lastSeen;
+    private arcs;
+    private arcSnap;
+    private mwState;
+    private series;
+    private ei;
+    private t0;
+    private convergenceSec;
+    private convergedStreak;
+    constructor(opts: PppOptions, source: PppEphemerisSource);
+    /** Reinitialise the filter (fresh session). */
+    reset(): void;
+    private flushArc;
+    /** Detect a cycle slip via the Melbourne–Wübbena wide-lane ambiguity. */
+    private detectSlip;
+    private growState;
+    private dropState;
+    /** State index of a constellation's receiver clock, creating it lazily. */
+    private clockOf;
+    /** Process one epoch, carrying the filter forward; returns its result. */
+    process(epoch: PppEpoch): PppEpochResult;
+    /** Finalise the session: capture the AR state, flush open arcs, assemble. */
+    solution(): PppSolution;
+}
+/**
+ * Static/kinematic float PPP over a whole session — a thin batch wrapper that
+ * drives {@link PppEngine} epoch by epoch. Accepts an Sp3File (wrapped as the
+ * default source) or any {@link PppEphemerisSource} (e.g. SSR/HAS-fed).
+ */
+declare function solvePpp(epochs: PppEpoch[], sp3OrSource: Sp3File | PppEphemerisSource, opts: PppOptions): PppSolution;
 
 /**
  * Troposphere for PPP: a-priori Saastamoinen zenith hydrostatic delay
@@ -1408,4 +1567,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 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 NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkTrackPoint, type SatGeom, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, 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, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, 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 RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, 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, applyOrbitClock, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

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, C as ClkFile, S as Sp3File } from './clk-DGVBGRzA.js';
+import { a as HasOrbitCorrection, b as HasClockCorrection } from './has-BaXnQarR.js';
 import { AntexFile } from './antex.js';
 import './ephemeris-BAnBLLB9.js';
 import './ecef-CF0uAysr.js';
@@ -706,6 +707,77 @@ declare class SbasProcessor {
     private searchIgp;
 }
 
+/**
+ * Apply SSR / HAS orbit + clock corrections to a broadcast ephemeris to obtain
+ * a precise satellite ECEF position and clock offset — the bridge between the
+ * decode-only SSR/HAS decoders (`rtcm3/ssr`, `rtcm3/igs-ssr`, `navbits/has`)
+ * and a PPP/SPP solve.
+ *
+ * Convention (internal, "additive"): the orbit deltas are stored in the
+ * satellite radial / along-track / cross-track (RAC) frame such that
+ *
+ *     correctedPos   = broadcastPos   + R·Δ(radial, along, cross)
+ *     correctedClock = broadcastClock + (c0 + c1·dt + c2·dt²)/c
+ *
+ * This matches Galileo HAS directly (ICD §7.2–7.3). RTCM-SSR and IGS-SSR use
+ * the **opposite** orbit sign, so their adapters negate the orbit deltas when
+ * building an {@link OrbitClockCorrection} (see `ssrToOrbitClock`, added with
+ * the RTCM/IGS wiring). Clock and bias sign conventions are handled per source
+ * in the adapters, not here.
+ */
+
+/**
+ * A normalized per-satellite orbit + clock correction, ready to apply to a
+ * broadcast ephemeris. Additive convention (see file header). Orbit deltas are
+ * in the RAC frame in metres, referenced to the broadcast issue-of-data `iod`
+ * at `refTimeMs`; optional `dot*` give their linear rate (RTCM/IGS SSR carry
+ * these, HAS does not). Clock is a metres polynomial about `refTimeMs`.
+ */
+interface OrbitClockCorrection {
+    prn: string;
+    /** Broadcast IODE (GPS/…) or IODnav (Galileo) the deltas refer to. */
+    iod: number;
+    /** Correction reference epoch, GPS-scale ms (matches `computeSatPosition`). */
+    refTimeMs: number;
+    radial: number;
+    along: number;
+    cross: number;
+    dotRadial?: number;
+    dotAlong?: number;
+    dotCross?: number;
+    /** Clock polynomial (m): c0 + c1·dt + c2·dt². */
+    clockC0: number;
+    clockC1?: number;
+    clockC2?: number;
+    /** Additive high-rate clock term (m), when a separate high-rate stream is combined. */
+    highRateClock?: number;
+    /** Validity window (s) from `refTimeMs`; null = unbounded / unknown. */
+    validitySeconds?: number | null;
+}
+/** A corrected satellite state: ECEF position (m) + the clock offset (s) to ADD
+ *  to the broadcast satellite clock. The caller forms the total clock as
+ *  `satClockCorrection(eph, t) + clockOffsetS`. */
+interface AppliedOrbitClock {
+    x: number;
+    y: number;
+    z: number;
+    clockOffsetS: number;
+}
+/**
+ * Corrected satellite ECEF position + clock offset at emission time `timeMs`,
+ * from a broadcast ephemeris and an SSR/HAS orbit-clock correction. Returns
+ * null when the correction's IOD doesn't match the ephemeris, the correction
+ * has expired, or the orbit can't be computed.
+ */
+declare function applyOrbitClock(eph: Ephemeris, corr: OrbitClockCorrection, timeMs: number): AppliedOrbitClock | null;
+/**
+ * Build a normalized {@link OrbitClockCorrection} from a Galileo-HAS orbit
+ * correction (+ optional matching clock correction). HAS already uses the
+ * additive convention, so no sign flip. Returns null when the orbit deltas are
+ * flagged "not available". `refTimeMs` is the message ToH as GPS-scale ms.
+ */
+declare function hasToOrbitClock(orbit: HasOrbitCorrection, clock: HasClockCorrection | undefined, refTimeMs: number, validitySeconds: number | null): OrbitClockCorrection | null;
+
 /**
  * ANTEX phase-centre application for PPP: satellite and receiver antenna
  * PCO (offset) + PCV (elevation-dependent variation), ionosphere-free
@@ -830,6 +902,10 @@ interface PppSatObs {
      * the antenna corrections; optional if corrections are off. */
     band1?: string;
     band2?: string;
+    /** RINEX-3 tracking codes for the two bands (e.g. '1C','2W'). Used to look
+     * up satellite code (OSB) biases from an SSR/HAS source; optional. */
+    code1?: string;
+    code2?: string;
     /** True if a cycle slip is flagged on either band (LLI or detected). */
     slip: boolean;
 }
@@ -952,7 +1028,90 @@ interface PppArc {
     startMs: number;
     endMs: number;
 }
-declare function solvePpp(epochs: PppEpoch[], sp3: Sp3File, opts: PppOptions): PppSolution;
+/** A raw satellite sample at an emission time: ECEF position + velocity
+ *  (m, m/s) and the clock offset (s) WITHOUT the periodic relativistic term —
+ *  the PPP light-time loop adds relativity + Sagnac. Null when unavailable. */
+interface SatSample {
+    x: number;
+    y: number;
+    z: number;
+    vx: number;
+    vy: number;
+    vz: number;
+    clkS: number;
+}
+/** Supplies precise satellite state to the PPP filter. Implemented by the SP3
+ *  path (default) and by an SSR/HAS source (broadcast eph + corrections), so
+ *  the same filter runs on precise products or a real-time correction stream. */
+interface PppEphemerisSource {
+    satState(prn: string, tEmitMs: number): SatSample | null;
+    /** Optional satellite code (OSB) bias (m) for a RINEX-3 tracking code, to
+     *  make the pseudorange consistent with the (SSR/HAS) clock reference. The
+     *  corrected code is `observed − codeBias`. SP3/CLK sources omit this (the
+     *  float ambiguity + broadcast clock reference absorb it). */
+    codeBias?(prn: string, signalCode: string, tEmitMs: number): number | null;
+}
+/** Default source: SP3 orbits + clocks, with optional high-rate CLK. Velocity
+ *  is a central finite difference (±0.5 s); the clock prefers the high-rate CLK
+ *  and falls back to the SP3 5-minute clock — matching the original solver. */
+declare class Sp3EphemerisSource implements PppEphemerisSource {
+    private readonly sp3;
+    private readonly clk?;
+    constructor(sp3: Sp3File, clk?: ClkFile | undefined);
+    satState(prn: string, tEmitMs: number): SatSample | null;
+}
+/**
+ * Streaming float PPP — the forward EKF as a stateful engine fed one epoch at a
+ * time (mirrors {@link RtkFloatEngine}). `process(epoch)` returns that epoch's
+ * result and carries the filter forward; `solution()` finalises the session.
+ * `solvePpp` (below) is a thin batch wrapper. Same estimator as before — a
+ * broadcast/SP3/SSR source supplies satellite state via {@link PppEphemerisSource}.
+ */
+declare class PppEngine {
+    private readonly source;
+    private readonly opts;
+    private readonly elevMask;
+    private readonly codeSigma;
+    private readonly phaseSigma;
+    private readonly ztdQ;
+    private readonly kinematic;
+    private readonly corrections?;
+    private readonly truth?;
+    private x;
+    private P;
+    private clkIdx;
+    private clockSeeded;
+    private ambIdx;
+    private lastSeen;
+    private arcs;
+    private arcSnap;
+    private mwState;
+    private series;
+    private ei;
+    private t0;
+    private convergenceSec;
+    private convergedStreak;
+    constructor(opts: PppOptions, source: PppEphemerisSource);
+    /** Reinitialise the filter (fresh session). */
+    reset(): void;
+    private flushArc;
+    /** Detect a cycle slip via the Melbourne–Wübbena wide-lane ambiguity. */
+    private detectSlip;
+    private growState;
+    private dropState;
+    /** State index of a constellation's receiver clock, creating it lazily. */
+    private clockOf;
+    /** Process one epoch, carrying the filter forward; returns its result. */
+    process(epoch: PppEpoch): PppEpochResult;
+    /** Finalise the session: capture the AR state, flush open arcs, assemble. */
+    solution(): PppSolution;
+}
+/**
+ * Static/kinematic float PPP over a whole session — a thin batch wrapper that
+ * drives {@link PppEngine} epoch by epoch. Accepts an Sp3File (wrapped as the
+ * default source) or any {@link PppEphemerisSource} (e.g. SSR/HAS-fed).
+ */
+declare function solvePpp(epochs: PppEpoch[], sp3OrSource: Sp3File | PppEphemerisSource, opts: PppOptions): PppSolution;
 
 /**
  * Troposphere for PPP: a-priori Saastamoinen zenith hydrostatic delay
@@ -1408,4 +1567,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 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 NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkTrackPoint, type SatGeom, type SbasCorrectionSource, type SbasIgp, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, 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, type CorrectionResult, type Degradation, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, 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 RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, 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, applyOrbitClock, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -1,7 +1,10 @@
 import {
   IONO_L1_M_PER_TECU,
+  PppEngine,
   RtkFloatEngine,
   SbasProcessor,
+  Sp3EphemerisSource,
+  applyOrbitClock,
   buildPppAntenna,
   createPppCorrections,
   estimateNarrowlaneFcb,
@@ -12,6 +15,7 @@ import {
   gimSlantIonoDelayL1,
   gimVerticalTec,
   gmst,
+  hasToOrbitClock,
   ionoFree,
   klobucharDelay,
   lambdaReduction,
@@ -32,7 +36,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-2CMAFEHI.js";
+} from "./chunk-DWUQRJ37.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-JYWYQ7DW.js";
 import "./chunk-HVXYFUCB.js";
@@ -46,8 +50,11 @@ import "./chunk-V2D5KUKH.js";
 import "./chunk-FIEWO4J4.js";
 export {
   IONO_L1_M_PER_TECU,
+  PppEngine,
   RtkFloatEngine,
   SbasProcessor,
+  Sp3EphemerisSource,
+  applyOrbitClock,
   buildPppAntenna,
   createPppCorrections,
   estimateNarrowlaneFcb,
@@ -58,6 +65,7 @@ export {
   gimSlantIonoDelayL1,
   gimVerticalTec,
   gmst,
+  hasToOrbitClock,
   ionoFree,
   klobucharDelay,
   lambdaReduction,

package/dist/sbf.d.cts

@@ -1,6 +1,7 @@
 export { S as SbfAlmanac, a as SbfAlmanacResult, b as SbfGlonassAlmanac, c as SbfKeplerAlmanac, d as SbfNavResult, p as parseSbfAlmanac, e as parseSbfNav } from './nav-BYmYoh9A.cjs';
 import { R as ReceiverPvt } from './receiver-pvt-CVFvi4kT.cjs';
 import { K as KeplerEphemeris, C as CnavEphemeris, G as GlonassEphemeris, E as Ephemeris } from './nav-DImXUdbp.cjs';
+import { H as HasMessage } from './has-BaXnQarR.cjs';
 
 /**
  * Septentrio SBF framing shared by the measurement and navigation
@@ -496,128 +497,6 @@ declare function sbasL5MessageName(type: number): string;
 /** True for L5 types that carry actual corrections (vs test/null/spare). */
 declare function isSbasL5Correction(type: number): boolean;
 
-/** Constellations HAS corrects: GPS (GNSS ID 0), Galileo (GNSS ID 2). */
-type HasSystem = 'G' | 'E';
-/** One system's satellite/signal mask (ICD §6.4.3). */
-interface HasSystemMask {
-    /** GNSS ID as broadcast: 0 GPS, 2 Galileo. */
-    gnssId: number;
-    system: HasSystem;
-    /** Masked satellites in mask order, e.g. ["E03", "E05", ...]. */
-    prns: string[];
-    /** Set bits of the 16-bit signal mask, ascending. */
-    signalIndices: number[];
-    /** ICD names for `signalIndices`. */
-    signals: string[];
-    /**
-     * Cell mask, [satellite][signal] in mask order — true when the
-     * sat×signal combination carries biases. Null when the cell-mask
-     * availability flag is 0 (all combinations available).
-     */
-    cellMask: boolean[][] | null;
-    /** Nav-message index the corrections refer to (0 = LNAV / I/NAV). */
-    navMessage: number;
-}
-/** The satellite/signal mask block of one mask ID. */
-interface HasMasks {
-    maskId: number;
-    systems: HasSystemMask[];
-}
-/** Orbit correction for one satellite (ICD §6.5). */
-interface HasOrbitCorrection {
-    system: HasSystem;
-    prn: string;
-    /** Broadcast issue of data the deltas refer to (IODE / IODNav). */
-    gnssIod: number;
-    /** Delta radial in m (LSB 0.0025), null = data not available. */
-    deltaRadial: number | null;
-    /** Delta in-track in m (LSB 0.008), null = data not available. */
-    deltaInTrack: number | null;
-    /** Delta cross-track in m (LSB 0.008), null = data not available. */
-    deltaCrossTrack: number | null;
-}
-interface HasOrbitBlock {
-    validityIndex: number;
-    /** Validity interval in s (null = reserved index 15). */
-    validitySeconds: number | null;
-    corrections: HasOrbitCorrection[];
-}
-/** Clock correction for one satellite (ICD §6.6-6.7). */
-interface HasClockCorrection {
-    system: HasSystem;
-    prn: string;
-    /**
-     * Delta clock C0 in m, multiplier folded in; positive means the
-     * satellite clock offset grows (add Δclock/c to the broadcast
-     * clock). Null when data not available or the satellite is flagged
-     * "shall not be used" (see `notUsable`).
-     */
-    deltaClock: number | null;
-    /** True for the reserved "insufficient accuracy, do not use" value. */
-    notUsable: boolean;
-}
-interface HasClockBlock {
-    validityIndex: number;
-    validitySeconds: number | null;
-    /** Delta-clock C0 multiplier per system (1-4), in mask order. */
-    multipliers: {
-        system: HasSystem;
-        multiplier: number;
-    }[];
-    corrections: HasClockCorrection[];
-}
-/** Code or phase bias of one satellite × signal cell (ICD §6.8-6.9). */
-interface HasSignalBias {
-    system: HasSystem;
-    prn: string;
-    signalIndex: number;
-    signal: string;
-    /** Bias in m (code, LSB 0.02) or cycles (phase, LSB 0.01); null = N/A. */
-    bias: number | null;
-    /** Phase-discontinuity indicator (phase biases only). */
-    discontinuity?: number;
-}
-interface HasBiasBlock {
-    validityIndex: number;
-    validitySeconds: number | null;
-    biases: HasSignalBias[];
-}
-/** Which sub-blocks the MT1 header flags as present (ICD §6.2). */
-interface HasContentFlags {
-    mask: boolean;
-    orbit: boolean;
-    clockFullSet: boolean;
-    clockSubset: boolean;
-    codeBias: boolean;
-    phaseBias: boolean;
-}
-/** One decoded HAS MT1 message. */
-interface HasMessage {
-    /** Receiver time of week (s) of the page completing the message. */
-    tow?: number;
-    /** HAS status of the carrying pages (0 test, 1 operational). */
-    status: number;
-    messageId: number;
-    /** Message size in pages (= pages used for the RS decode). */
-    messageSize: number;
-    /** Time of hour of the corrections, s (0-3599, ICD §6.2). */
-    toh: number;
-    maskId: number;
-    iodSetId: number;
-    flags: HasContentFlags;
-    /** The mask block in force (own or cached, see `maskFromCache`). */
-    masks?: HasMasks;
-    /** True when `masks` came from an earlier message with this mask ID. */
-    maskFromCache: boolean;
-    orbit?: HasOrbitBlock;
-    clockFullSet?: HasClockBlock;
-    clockSubset?: HasClockBlock;
-    codeBias?: HasBiasBlock;
-    phaseBias?: HasBiasBlock;
-    /** Null when fully parsed, else why the sub-blocks are missing. */
-    parseError: string | null;
-}
-
 /**
  * Septentrio SBF GALRawCNAV (4024) blocks: one 492-bit Galileo E6-B
  * C/NAV page per block (after deinterleaving and Viterbi decoding),
@@ -731,4 +610,4 @@ interface SbfParseResult {
  */
 declare function parseSbfMeas(data: Uint8Array): SbfParseResult;
 
-export { CnavEphemeris, type DfmcCensus, type DfmcMessageCb, type GeoL5Result, type HasMessage, ReceiverPvt, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, type SbasMessageCb, type SbfBdsNavResult, type SbfCnavEphemeris, type SbfCnavResult, type SbfGalEphemeris, type SbfGalNavResult, type SbfGeoL5Result, type SbfGeoNavResult, type SbfGloNavResult, type SbfGpsNavResult, type SbfHasMessage, type SbfHasResult, type SbfIonoUtcResult, type SbfMeasEpoch, type SbfMeasurement, type SbfNavCounts, type SbfNavigation, type SbfParseResult, type SbfPvt, type SbfPvtResult, type SbfReceiverSetup, decodeSbfNavigation, feedGeoBlock, feedGeoL5Block, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, crc16 as sbfCrc16, scanSbfFrames };
+export { CnavEphemeris, type DfmcCensus, type DfmcMessageCb, type GeoL5Result, HasMessage, ReceiverPvt, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, type SbasMessageCb, type SbfBdsNavResult, type SbfCnavEphemeris, type SbfCnavResult, type SbfGalEphemeris, type SbfGalNavResult, type SbfGeoL5Result, type SbfGeoNavResult, type SbfGloNavResult, type SbfGpsNavResult, type SbfHasMessage, type SbfHasResult, type SbfIonoUtcResult, type SbfMeasEpoch, type SbfMeasurement, type SbfNavCounts, type SbfNavigation, type SbfParseResult, type SbfPvt, type SbfPvtResult, type SbfReceiverSetup, decodeSbfNavigation, feedGeoBlock, feedGeoL5Block, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, crc16 as sbfCrc16, scanSbfFrames };