gnss-js 1.10.0 → 1.12.0

package/dist/positioning.d.cts

@@ -3,6 +3,36 @@ import { DopValues } from './orbit.cjs';
 import './ephemeris-Ohl6NAfw.cjs';
 import './ecef-CF0uAysr.cjs';
 
+/**
+ * Klobuchar broadcast ionospheric model (IS-GPS-200, 20.3.3.5.2.5).
+ *
+ * Evaluates the single-frequency vertical→slant ionospheric delay from
+ * the eight coefficients broadcast in the GPS navigation message (and
+ * carried in RINEX nav headers as GPSA/GPSB, or ION ALPHA/BETA in
+ * RINEX 2). The model captures ~50% RMS of the actual delay — crude,
+ * but exactly what a single-frequency receiver applies, which makes it
+ * both a useful SPP correction and an honest "what broadcast users
+ * get" reference against measured or GIM TEC.
+ */
+/** The eight Klobuchar coefficients from the navigation message. */
+interface KlobucharCoeffs {
+    /** GPSA: alpha0..alpha3 (s, s/sc, s/sc², s/sc³). */
+    alpha: number[];
+    /** GPSB: beta0..beta3 (s, s/sc, s/sc², s/sc³). */
+    beta: number[];
+}
+/**
+ * Slant ionospheric group delay on GPS L1, in seconds.
+ *
+ * @param coeffs Broadcast alpha/beta coefficients
+ * @param latRad Receiver geodetic latitude (radians)
+ * @param lonRad Receiver geodetic longitude (radians)
+ * @param azRad Satellite azimuth (radians)
+ * @param elRad Satellite elevation (radians)
+ * @param gpsTowSec GPS time of week (seconds)
+ */
+declare function klobucharDelay(coeffs: KlobucharCoeffs, latRad: number, lonRad: number, azRad: number, elRad: number, gpsTowSec: number): number;
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -11,9 +41,11 @@ import './ecef-CF0uAysr.cjs';
  * constellation (absorbing inter-system time offsets), satellite clock
  * polynomial + relativistic correction, broadcast group delay
  * (TGD/BGD), Sagnac (Earth-rotation) correction, elevation
- * masking/weighting, and a simple tropospheric model. Ionospheric delay is NOT modelled — use the iono-free
- * combination (ionoFree helper) with dual-frequency pseudoranges for
- * metre-level results, or expect ~2–10 m of iono bias on L1-only.
+ * masking/weighting, and a simple tropospheric model. Ionospheric
+ * delay can be modelled from the broadcast Klobuchar coefficients
+ * (`iono` option) for single-frequency measurements; for metre-level
+ * results prefer the iono-free combination (ionoFree helper) with
+ * dual-frequency pseudoranges.
  */
 
 interface SppOptions {
@@ -25,6 +57,14 @@ interface SppOptions {
     maxIterations?: number;
     /** Convergence threshold on the position update (m). Default 1e-4. */
     convergenceM?: number;
+    /**
+     * Broadcast Klobuchar coefficients (RINEX nav header GPSA/GPSB —
+     * `NavResult.header.ionoCorrections`). When given, the modelled
+     * slant delay is removed from every single-frequency pseudorange,
+     * scaled to each system's primary frequency. Omit when feeding
+     * iono-free combinations.
+     */
+    iono?: KlobucharCoeffs;
     /**
      * Apply the broadcast group-delay correction (GPS TGD, Galileo
      * BGD E5a/E1, BeiDou TGD1) to the satellite clock. Correct for
@@ -73,4 +113,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 SppOptions, type SppSolution, ionoFree, satClockCorrection, solveSpp };
+export { type KlobucharCoeffs, type SppOptions, type SppSolution, ionoFree, klobucharDelay, satClockCorrection, solveSpp };

package/dist/positioning.d.ts

@@ -3,6 +3,36 @@ import { DopValues } from './orbit.js';
 import './ephemeris-Ohl6NAfw.js';
 import './ecef-CF0uAysr.js';
 
+/**
+ * Klobuchar broadcast ionospheric model (IS-GPS-200, 20.3.3.5.2.5).
+ *
+ * Evaluates the single-frequency vertical→slant ionospheric delay from
+ * the eight coefficients broadcast in the GPS navigation message (and
+ * carried in RINEX nav headers as GPSA/GPSB, or ION ALPHA/BETA in
+ * RINEX 2). The model captures ~50% RMS of the actual delay — crude,
+ * but exactly what a single-frequency receiver applies, which makes it
+ * both a useful SPP correction and an honest "what broadcast users
+ * get" reference against measured or GIM TEC.
+ */
+/** The eight Klobuchar coefficients from the navigation message. */
+interface KlobucharCoeffs {
+    /** GPSA: alpha0..alpha3 (s, s/sc, s/sc², s/sc³). */
+    alpha: number[];
+    /** GPSB: beta0..beta3 (s, s/sc, s/sc², s/sc³). */
+    beta: number[];
+}
+/**
+ * Slant ionospheric group delay on GPS L1, in seconds.
+ *
+ * @param coeffs Broadcast alpha/beta coefficients
+ * @param latRad Receiver geodetic latitude (radians)
+ * @param lonRad Receiver geodetic longitude (radians)
+ * @param azRad Satellite azimuth (radians)
+ * @param elRad Satellite elevation (radians)
+ * @param gpsTowSec GPS time of week (seconds)
+ */
+declare function klobucharDelay(coeffs: KlobucharCoeffs, latRad: number, lonRad: number, azRad: number, elRad: number, gpsTowSec: number): number;
+
 /**
  * Single-point positioning (SPP) from pseudoranges and broadcast
  * ephemerides.
@@ -11,9 +41,11 @@ import './ecef-CF0uAysr.js';
  * constellation (absorbing inter-system time offsets), satellite clock
  * polynomial + relativistic correction, broadcast group delay
  * (TGD/BGD), Sagnac (Earth-rotation) correction, elevation
- * masking/weighting, and a simple tropospheric model. Ionospheric delay is NOT modelled — use the iono-free
- * combination (ionoFree helper) with dual-frequency pseudoranges for
- * metre-level results, or expect ~2–10 m of iono bias on L1-only.
+ * masking/weighting, and a simple tropospheric model. Ionospheric
+ * delay can be modelled from the broadcast Klobuchar coefficients
+ * (`iono` option) for single-frequency measurements; for metre-level
+ * results prefer the iono-free combination (ionoFree helper) with
+ * dual-frequency pseudoranges.
  */
 
 interface SppOptions {
@@ -25,6 +57,14 @@ interface SppOptions {
     maxIterations?: number;
     /** Convergence threshold on the position update (m). Default 1e-4. */
     convergenceM?: number;
+    /**
+     * Broadcast Klobuchar coefficients (RINEX nav header GPSA/GPSB —
+     * `NavResult.header.ionoCorrections`). When given, the modelled
+     * slant delay is removed from every single-frequency pseudorange,
+     * scaled to each system's primary frequency. Omit when feeding
+     * iono-free combinations.
+     */
+    iono?: KlobucharCoeffs;
     /**
      * Apply the broadcast group-delay correction (GPS TGD, Galileo
      * BGD E5a/E1, BeiDou TGD1) to the satellite clock. Correct for
@@ -73,4 +113,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 SppOptions, type SppSolution, ionoFree, satClockCorrection, solveSpp };
+export { type KlobucharCoeffs, type SppOptions, type SppSolution, ionoFree, klobucharDelay, satClockCorrection, solveSpp };

package/dist/positioning.js

@@ -1,9 +1,10 @@
 import {
   ionoFree,
+  klobucharDelay,
   satClockCorrection,
   solveSpp
-} from "./chunk-ZTFKRZMM.js";
-import "./chunk-AK6LNWZX.js";
+} from "./chunk-LKFSXVIZ.js";
+import "./chunk-DMNFB7Q6.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
@@ -11,6 +12,7 @@ import "./chunk-LEEU5OIO.js";
 import "./chunk-FIEWO4J4.js";
 export {
   ionoFree,
+  klobucharDelay,
   satClockCorrection,
   solveSpp
 };