gnss-js 2.2.0 → 2.4.0

package/dist/positioning.d.cts

@@ -1,8 +1,8 @@
 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 { a as HasOrbitCorrection, b as HasClockCorrection } from './has-BaXnQarR.cjs';
 import './ephemeris-BAnBLLB9.cjs';
 import './ecef-CF0uAysr.cjs';
 import './nav-BYmYoh9A.cjs';
@@ -83,6 +83,87 @@ declare function gimVerticalTec(grid: IonexGrid, timeMs: number, latDeg: number,
  */
 declare function gimSlantIonoDelayL1(grid: IonexGrid, latRad: number, lonRad: number, azRad: number, elRad: number, timeMs: number): number | 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;
+
+/**
+ * Per-frequency receiver antenna corrections for RTK.
+ *
+ * RTK forms double differences per signal (L1, L2, … separately), so — unlike
+ * PPP, which works on the ionosphere-free combination — it needs the antenna
+ * phase-centre offset (PCO) and no-azimuth phase-centre variation (PCV) for
+ * each frequency in isolation, not IF-combined. This module is therefore a
+ * deliberately separate model from {@link buildPppAntenna}.
+ *
+ * The satellite antenna cancels in the base↔rover double difference, so only
+ * the *receiver* antennas matter — and only when base and rover differ (same
+ * antenna type + orientation cancels too). The correction is the same range
+ * quantity the PPP receiver term applies (PCO projected on the line of sight +
+ * elevation PCV + marker→ARP delta), added to the modelled geometric range.
+ */
+
+/** A single-frequency receiver antenna offset (PCO in ENU metres + no-azimuth
+ *  PCV over zenith). Shares {@link AntennaOffset}'s shape so {@link interpPcv}
+ *  applies directly. */
+type RtkAntennaOffset = AntennaOffset;
+interface RtkAntennaModel {
+    /** Receiver PCO + PCV for one ANTEX frequency code (e.g. 'G01', 'E05'), or
+     *  null when the antenna type or that frequency isn't calibrated. */
+    rcvOffset(antType: string, freq: string): RtkAntennaOffset | null;
+}
+/**
+ * Build a per-frequency receiver antenna model from an ANTEX file (e.g.
+ * igs20.atx). Receiver antennas are keyed by their 20-char type (RINEX header
+ * "ANT # / TYPE"); satellite entries are ignored (they cancel in the DD).
+ */
+declare function buildRtkAntenna(antex: AntexFile): RtkAntennaModel;
+/**
+ * Additive range correction (m) for one receiver observing one satellite,
+ * matching the PPP receiver-antenna convention: it is *added to the modelled
+ * geometric range* (marker→satellite). `losRcvToSat` is the unit line of sight
+ * from the receiver to the satellite; `antDeltaEnu` is the marker→ARP offset
+ * in local ENU (m).
+ */
+declare function rcvAntennaRangeM(off: RtkAntennaOffset, losRcvToSat: [number, number, number], latRad: number, lonRad: number, elRad: number, antDeltaEnu?: readonly [number, number, number]): number;
+/** Per-receiver antenna configuration for {@link RtkFloatOptions.antenna}. */
+interface RtkReceiverAntenna {
+    /** RINEX antenna type (20-char "ANT # / TYPE"), keyed into the model. */
+    type: string;
+    /** Marker→ARP offset in local ENU (m). RINEX "ANTENNA: DELTA H/E/N" is
+     *  [Up, East, North] → pass as [E, N, U]. Default [0,0,0]. */
+    deltaEnu?: [number, number, number];
+}
+/** RTK antenna configuration: a per-frequency model + the base/rover antennas.
+ *  Corrections are applied only for a receiver whose antenna is supplied and
+ *  found in the model; the satellite antenna is never needed (cancels in the
+ *  double difference). */
+interface RtkAntennaConfig {
+    model: RtkAntennaModel;
+    base?: RtkReceiverAntenna;
+    rover?: RtkReceiverAntenna;
+}
+/** ANTEX frequency code for an RTK signal group ("G1C" → "G01"): system letter
+ *  + '0' + band digit. */
+declare function antexFreqOfGroup(group: string): string;
+
 /**
  * RTK positioning, stage 1: epoch-wise double-differenced (DD)
  * processing of a base/rover receiver pair.
@@ -306,6 +387,14 @@ interface RtkFloatOptions {
     partialFixing?: boolean;
     /** Minimum DD ambiguities to attempt/accept a fix. Default 4. */
     minFixAmbiguities?: number;
+    /**
+     * Receiver antenna phase-centre corrections (per-frequency PCO + PCV) for
+     * the base and rover. Only worthwhile on mixed-antenna baselines — with the
+     * same antenna type + orientation at both ends the correction cancels in the
+     * double difference. The satellite antenna is never needed (it cancels too).
+     * Build the model with {@link buildRtkAntenna}. Default: none.
+     */
+    antenna?: RtkAntennaConfig | null;
 }
 interface RtkFloatSolution {
     /** Epoch, GPS-scale milliseconds. */
@@ -579,6 +668,8 @@ interface Degradation {
     cionoRamp: number;
     rssUdre: boolean;
     rssIono: boolean;
+    /** Cᶜᵒᵛᵃʳⁱᵃⁿᶜᵉ — the MT28 δUDRE quantization term (m), 0 when no MT28. */
+    ccovariance: number;
 }
 /** A per-satellite SBAS correction, ready to apply in an SPP solve. */
 interface SbasSatCorrection {
@@ -626,6 +717,10 @@ declare class SbasProcessor {
     private tlat;
     private degr;
     private mt27;
+    /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+     *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+     *  scale exponent. Keyed by PRN; cleared on a mask change. */
+    private mt28;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -650,6 +745,17 @@ declare class SbasProcessor {
      *  set: the highest-priority containing region wins (ties → lower δUDRE),
      *  else the outside value. 1 when no MT27 has been received. */
     deltaUdre(latDeg: number, lonDeg: number): number;
+    /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+     *  two satellites per message, each an upper-triangular Cholesky factor E and a
+     *  scale exponent that together define the relative clock-ephemeris covariance
+     *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+     *  the current MT1 mask, so a mismatched IODP is dropped. */
+    private decodeMt28;
+    /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+     *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+     *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+     *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+    private mt28DeltaUdre;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -693,7 +799,11 @@ declare class SbasProcessor {
      * no valid (unexpired) long-term + fast correction is available. Add `dPos`
      * to the satellite ECEF position and `dClkS` to its clock bias.
      */
-    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number): SbasSatCorrection | null;
+    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number, 
+    /** Unit line-of-sight vector user→satellite in ECEF, for the MT28
+     *  geometry-dependent δUDRE (final protection-level pass). When given and
+     *  MT28 data exists for the PRN, it supersedes the MT27 region δUDRE. */
+    losUnitEcef?: readonly [number, number, number]): SbasSatCorrection | null;
     /** The last-decoded MT10 degradation factors, or null if none seen. */
     get degradation(): Readonly<Degradation> | null;
     /**
@@ -778,27 +888,6 @@ declare function applyOrbitClock(eph: Ephemeris, corr: OrbitClockCorrection, tim
  */
 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
- * 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
@@ -1489,7 +1578,10 @@ interface SbasCorrectionSource {
     satCorrection(prn: string, week: number, tow: number, 
     /** User geodetic lat/lon (rad) for the MT27 service-region δUDRE, when
      *  known (the final protection-level pass). Omitted during iteration. */
-    userLatRad?: number, userLonRad?: number): {
+    userLatRad?: number, userLonRad?: number, 
+    /** Unit LOS user→satellite (ECEF) for the MT28 geometry-dependent δUDRE
+     *  (final protection-level pass). */
+    losUnitEcef?: readonly [number, number, number]): {
         dPos: [number, number, number];
         dClkS: number;
         varM2: number;
@@ -1567,4 +1659,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 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 };
+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 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, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -1,8 +1,8 @@
 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 { a as HasOrbitCorrection, b as HasClockCorrection } from './has-BaXnQarR.js';
 import './ephemeris-BAnBLLB9.js';
 import './ecef-CF0uAysr.js';
 import './nav-CN8rxL4h.js';
@@ -83,6 +83,87 @@ declare function gimVerticalTec(grid: IonexGrid, timeMs: number, latDeg: number,
  */
 declare function gimSlantIonoDelayL1(grid: IonexGrid, latRad: number, lonRad: number, azRad: number, elRad: number, timeMs: number): number | 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;
+
+/**
+ * Per-frequency receiver antenna corrections for RTK.
+ *
+ * RTK forms double differences per signal (L1, L2, … separately), so — unlike
+ * PPP, which works on the ionosphere-free combination — it needs the antenna
+ * phase-centre offset (PCO) and no-azimuth phase-centre variation (PCV) for
+ * each frequency in isolation, not IF-combined. This module is therefore a
+ * deliberately separate model from {@link buildPppAntenna}.
+ *
+ * The satellite antenna cancels in the base↔rover double difference, so only
+ * the *receiver* antennas matter — and only when base and rover differ (same
+ * antenna type + orientation cancels too). The correction is the same range
+ * quantity the PPP receiver term applies (PCO projected on the line of sight +
+ * elevation PCV + marker→ARP delta), added to the modelled geometric range.
+ */
+
+/** A single-frequency receiver antenna offset (PCO in ENU metres + no-azimuth
+ *  PCV over zenith). Shares {@link AntennaOffset}'s shape so {@link interpPcv}
+ *  applies directly. */
+type RtkAntennaOffset = AntennaOffset;
+interface RtkAntennaModel {
+    /** Receiver PCO + PCV for one ANTEX frequency code (e.g. 'G01', 'E05'), or
+     *  null when the antenna type or that frequency isn't calibrated. */
+    rcvOffset(antType: string, freq: string): RtkAntennaOffset | null;
+}
+/**
+ * Build a per-frequency receiver antenna model from an ANTEX file (e.g.
+ * igs20.atx). Receiver antennas are keyed by their 20-char type (RINEX header
+ * "ANT # / TYPE"); satellite entries are ignored (they cancel in the DD).
+ */
+declare function buildRtkAntenna(antex: AntexFile): RtkAntennaModel;
+/**
+ * Additive range correction (m) for one receiver observing one satellite,
+ * matching the PPP receiver-antenna convention: it is *added to the modelled
+ * geometric range* (marker→satellite). `losRcvToSat` is the unit line of sight
+ * from the receiver to the satellite; `antDeltaEnu` is the marker→ARP offset
+ * in local ENU (m).
+ */
+declare function rcvAntennaRangeM(off: RtkAntennaOffset, losRcvToSat: [number, number, number], latRad: number, lonRad: number, elRad: number, antDeltaEnu?: readonly [number, number, number]): number;
+/** Per-receiver antenna configuration for {@link RtkFloatOptions.antenna}. */
+interface RtkReceiverAntenna {
+    /** RINEX antenna type (20-char "ANT # / TYPE"), keyed into the model. */
+    type: string;
+    /** Marker→ARP offset in local ENU (m). RINEX "ANTENNA: DELTA H/E/N" is
+     *  [Up, East, North] → pass as [E, N, U]. Default [0,0,0]. */
+    deltaEnu?: [number, number, number];
+}
+/** RTK antenna configuration: a per-frequency model + the base/rover antennas.
+ *  Corrections are applied only for a receiver whose antenna is supplied and
+ *  found in the model; the satellite antenna is never needed (cancels in the
+ *  double difference). */
+interface RtkAntennaConfig {
+    model: RtkAntennaModel;
+    base?: RtkReceiverAntenna;
+    rover?: RtkReceiverAntenna;
+}
+/** ANTEX frequency code for an RTK signal group ("G1C" → "G01"): system letter
+ *  + '0' + band digit. */
+declare function antexFreqOfGroup(group: string): string;
+
 /**
  * RTK positioning, stage 1: epoch-wise double-differenced (DD)
  * processing of a base/rover receiver pair.
@@ -306,6 +387,14 @@ interface RtkFloatOptions {
     partialFixing?: boolean;
     /** Minimum DD ambiguities to attempt/accept a fix. Default 4. */
     minFixAmbiguities?: number;
+    /**
+     * Receiver antenna phase-centre corrections (per-frequency PCO + PCV) for
+     * the base and rover. Only worthwhile on mixed-antenna baselines — with the
+     * same antenna type + orientation at both ends the correction cancels in the
+     * double difference. The satellite antenna is never needed (it cancels too).
+     * Build the model with {@link buildRtkAntenna}. Default: none.
+     */
+    antenna?: RtkAntennaConfig | null;
 }
 interface RtkFloatSolution {
     /** Epoch, GPS-scale milliseconds. */
@@ -579,6 +668,8 @@ interface Degradation {
     cionoRamp: number;
     rssUdre: boolean;
     rssIono: boolean;
+    /** Cᶜᵒᵛᵃʳⁱᵃⁿᶜᵉ — the MT28 δUDRE quantization term (m), 0 when no MT28. */
+    ccovariance: number;
 }
 /** A per-satellite SBAS correction, ready to apply in an SPP solve. */
 interface SbasSatCorrection {
@@ -626,6 +717,10 @@ declare class SbasProcessor {
     private tlat;
     private degr;
     private mt27;
+    /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+     *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+     *  scale exponent. Keyed by PRN; cleared on a mask change. */
+    private mt28;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -650,6 +745,17 @@ declare class SbasProcessor {
      *  set: the highest-priority containing region wins (ties → lower δUDRE),
      *  else the outside value. 1 when no MT27 has been received. */
     deltaUdre(latDeg: number, lonDeg: number): number;
+    /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+     *  two satellites per message, each an upper-triangular Cholesky factor E and a
+     *  scale exponent that together define the relative clock-ephemeris covariance
+     *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+     *  the current MT1 mask, so a mismatched IODP is dropped. */
+    private decodeMt28;
+    /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+     *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+     *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+     *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+    private mt28DeltaUdre;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -693,7 +799,11 @@ declare class SbasProcessor {
      * no valid (unexpired) long-term + fast correction is available. Add `dPos`
      * to the satellite ECEF position and `dClkS` to its clock bias.
      */
-    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number): SbasSatCorrection | null;
+    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number, 
+    /** Unit line-of-sight vector user→satellite in ECEF, for the MT28
+     *  geometry-dependent δUDRE (final protection-level pass). When given and
+     *  MT28 data exists for the PRN, it supersedes the MT27 region δUDRE. */
+    losUnitEcef?: readonly [number, number, number]): SbasSatCorrection | null;
     /** The last-decoded MT10 degradation factors, or null if none seen. */
     get degradation(): Readonly<Degradation> | null;
     /**
@@ -778,27 +888,6 @@ declare function applyOrbitClock(eph: Ephemeris, corr: OrbitClockCorrection, tim
  */
 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
- * 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
@@ -1489,7 +1578,10 @@ interface SbasCorrectionSource {
     satCorrection(prn: string, week: number, tow: number, 
     /** User geodetic lat/lon (rad) for the MT27 service-region δUDRE, when
      *  known (the final protection-level pass). Omitted during iteration. */
-    userLatRad?: number, userLonRad?: number): {
+    userLatRad?: number, userLonRad?: number, 
+    /** Unit LOS user→satellite (ECEF) for the MT28 geometry-dependent δUDRE
+     *  (final protection-level pass). */
+    losUnitEcef?: readonly [number, number, number]): {
         dPos: [number, number, number];
         dClkS: number;
         varM2: number;
@@ -1567,4 +1659,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 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 };
+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 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, antexFreqOfGroup, applyOrbitClock, buildPppAntenna, buildRtkAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, hasToOrbitClock, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, rcvAntennaRangeM, resolvePppAmbiguities, satClockCorrection, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -4,8 +4,10 @@ import {
   RtkFloatEngine,
   SbasProcessor,
   Sp3EphemerisSource,
+  antexFreqOfGroup,
   applyOrbitClock,
   buildPppAntenna,
+  buildRtkAntenna,
   createPppCorrections,
   estimateNarrowlaneFcb,
   estimateNetworkFcbs,
@@ -24,6 +26,7 @@ import {
   niellMapping,
   nlWavelength,
   postProcessRtk,
+  rcvAntennaRangeM,
   resolvePppAmbiguities,
   satClockCorrection,
   sbasIonoDeg,
@@ -36,7 +39,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-DWUQRJ37.js";
+} from "./chunk-B2XWO4RP.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-JYWYQ7DW.js";
 import "./chunk-HVXYFUCB.js";
@@ -54,8 +57,10 @@ export {
   RtkFloatEngine,
   SbasProcessor,
   Sp3EphemerisSource,
+  antexFreqOfGroup,
   applyOrbitClock,
   buildPppAntenna,
+  buildRtkAntenna,
   createPppCorrections,
   estimateNarrowlaneFcb,
   estimateNetworkFcbs,
@@ -74,6 +79,7 @@ export {
   niellMapping,
   nlWavelength,
   postProcessRtk,
+  rcvAntennaRangeM,
   resolvePppAmbiguities,
   satClockCorrection,
   sbasIonoDeg,

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "2.2.0",
+  "version": "2.4.0",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,