gnss-js 1.3.0 → 1.3.2

package/dist/index.d.ts

@@ -1,6 +1,6 @@
 export { D as DAYS_MJD2000_MJD, a as DAYS_MJD_JULIAN, H as HourCode, M as MILLISECONDS_GPS_TAI, b as MILLISECONDS_IN_DAY, c as MILLISECONDS_IN_HOUR, d as MILLISECONDS_IN_MINUTE, e as MILLISECONDS_IN_SECOND, f as MILLISECONDS_IN_WEEK, g as MILLISECONDS_TT_TAI, R as RINEX_CODES, S as SECONDS_IN_DAY, h as SECONDS_IN_HOUR, i as SECONDS_IN_MINUTE, j as SECONDS_IN_WEEK, k as START_BDS_TIME, l as START_GAL_TIME, m as START_GLO_LEAP, n as START_GPS_TIME, o as START_JULIAN_TAI, p as START_MJD_UNIX_SECONDS, q as START_NTP_TIME, r as START_TAI_TIME, s as START_UNIX_TIME, T as TimeDifference } from './time-DnI1VpE8.js';
 export { WGS84_ECCENTRICITY_SQUARED, WGS84_FLATTENING, WGS84_SEMI_MAJOR_AXIS, WGS84_SEMI_MINOR_AXIS } from './constants.js';
-export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OnSlipDetected, S as SYSTEM_NAMES, i as SYS_SHORT, j as buildGloChannelMap, k as buildObsIndices, l as formatUTCTime, m as getFreq, n as gloFreq } from './gnss-C-tgoYNa.js';
+export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSatellite, C as C_LIGHT, D as DEFAULT_ELEV_MASK_DEG, c as DUAL_FREQ_PAIRS, F as FREQ, G as GLO_CHANNEL_FALLBACK, d as GLO_F1_BASE, e as GLO_F1_STEP, f as GLO_F2_BASE, g as GLO_F2_STEP, h as GLO_F3, O as OMEGA_E, i as OnSlipDetected, S as SYSTEM_NAMES, j as SYS_SHORT, k as buildGloChannelMap, l as buildObsIndices, m as formatUTCTime, n as getFreq, o as gloFreq } from './gnss-DhnEr0Dy.js';
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.js';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.js';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.js';
@@ -9,7 +9,9 @@ export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResul
 export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, p as parseNavFile } from './nav-BAI1a9vK.js';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.js';
 export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
-export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.js';
+export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.js';
+export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
+export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } from './positioning.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';

package/dist/index.js

@@ -8,7 +8,7 @@ import {
   CycleSlipAccumulator,
   MultipathAccumulator,
   analyzeQuality
-} from "./chunk-G3N4S3DM.js";
+} from "./chunk-AW6BORTZ.js";
 import {
   frequencyLabel,
   parseAntex
@@ -53,7 +53,7 @@ import {
   phiBOCc,
   phiBOCs,
   phiBPSK
-} from "./chunk-5SPJH4MG.js";
+} from "./chunk-2K3FCJBX.js";
 import {
   Scale,
   getDateFromDayOfWeek,
@@ -87,6 +87,17 @@ import {
   rhumbLine,
   vincenty
 } from "./chunk-MIIM4LDY.js";
+import {
+  REFERENCE_FRAMES,
+  applyHelmert,
+  dateToEpoch,
+  transformFrame
+} from "./chunk-4KOXMHUW.js";
+import {
+  ionoFree,
+  satClockCorrection,
+  solveSpp
+} from "./chunk-PX7TPPTQ.js";
 import {
   computeAllPositions,
   computeDop,
@@ -94,11 +105,12 @@ import {
   computeSatPosition,
   computeVisibility,
   ecefToAzEl,
+  ephInfoToEphemeris,
   glonassPosition,
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-ZCNXERQ2.js";
+} from "./chunk-PQYBTE42.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -166,7 +178,7 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-OZCYOM5D.js";
+} from "./chunk-3U5AX7PY.js";
 import {
   BitReader,
   RTCM3_MESSAGE_NAMES,
@@ -184,7 +196,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-7EEWQ5DU.js";
+} from "./chunk-65CQINSB.js";
 import {
   DAYS_MJD2000_MJD,
   DAYS_MJD_JULIAN,
@@ -224,6 +236,7 @@ import {
   GLO_F2_BASE,
   GLO_F2_STEP,
   GLO_F3,
+  OMEGA_E,
   SYSTEM_NAMES,
   SYS_SHORT,
   buildGloChannelMap,
@@ -231,7 +244,7 @@ import {
   formatUTCTime,
   getFreq,
   gloFreq
-} from "./chunk-W5WKEV7U.js";
+} from "./chunk-FIEWO4J4.js";
 export {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -288,6 +301,8 @@ export {
   MILLISECONDS_IN_WEEK,
   MILLISECONDS_TT_TAI,
   MultipathAccumulator,
+  OMEGA_E,
+  REFERENCE_FRAMES,
   RINEX_CODES,
   RTCM3_MESSAGE_NAMES,
   Rtcm3Decoder,
@@ -314,6 +329,7 @@ export {
   WGS84_SEMI_MINOR_AXIS,
   WarningAccumulator,
   analyzeQuality,
+  applyHelmert,
   buildGloChannelMap,
   buildObsIndices,
   clampUnit,
@@ -329,12 +345,14 @@ export {
   createStreamStats,
   crxDecompress,
   crxRepair,
+  dateToEpoch,
   decodeEphemeris,
   decodeMsmFull,
   deg2dms,
   deg2rad,
   ecefToAzEl,
   ecefToGeodetic,
+  ephInfoToEphemeris,
   euclidean3D,
   fetchSourcetable,
   fmtF,
@@ -402,6 +420,7 @@ export {
   greatCircleMidpoint,
   hdrLine,
   horizonDistance,
+  ionoFree,
   keplerPosition,
   msmEpochToDate,
   navTimesFromEph,
@@ -424,11 +443,14 @@ export {
   resetGloFreqCache,
   rhumbLine,
   rtcm3Constellation,
+  satClockCorrection,
   selectEphemeris,
   setGloFreqNumber,
   setRtcm3DebugHandler,
+  solveSpp,
   systemCmp,
   systemName,
+  transformFrame,
   updateStationMeta,
   updateStreamStats,
   verifyChecksum,

package/dist/orbit.cjs

@@ -27,6 +27,7 @@ __export(orbit_exports, {
   computeVisibility: () => computeVisibility,
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
+  ephInfoToEphemeris: () => ephInfoToEphemeris,
   geodeticToEcef: () => geodeticToEcef,
   glonassPosition: () => glonassPosition,
   keplerPosition: () => keplerPosition,
@@ -79,6 +80,9 @@ var START_GPS_TIME = /* @__PURE__ */ new Date("1980-01-06T00:00:00Z");
 var START_BDS_TIME = /* @__PURE__ */ new Date("2006-01-01T00:00:14Z");
 var START_NTP_TIME = /* @__PURE__ */ new Date("1900-01-01T00:00:00Z");
 
+// src/constants/gnss.ts
+var OMEGA_E = 72921151467e-15;
+
 // src/time/gnss.ts
 function getTaiDate(date) {
   return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
@@ -166,7 +170,6 @@ var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
-var OMEGA_E = 72921151467e-15;
 var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
 var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
@@ -375,7 +378,8 @@ function selectEphemeris(ephemerides, prn, timeMs) {
     if (eph.prn !== prn) continue;
     const ephTime = eph.tocDate.getTime();
     const dt = Math.abs(timeMs - ephTime);
-    if (dt > 4 * 3600 * 1e3) continue;
+    const maxAge = eph.system === "R" || eph.system === "S" ? 18e5 : 4 * 36e5;
+    if (dt > maxAge) continue;
     if (dt < bestDt) {
       bestDt = dt;
       best = eph;
@@ -483,13 +487,13 @@ function ephInfoToEphemeris(info) {
       return null;
     if (info.vx === void 0 || info.vy === void 0 || info.vz === void 0)
       return null;
-    const now = new Date(info.lastReceived);
+    const moscow = new Date(info.lastReceived + 3 * 3600 * 1e3);
     const utcMidnight = Date.UTC(
-      now.getUTCFullYear(),
-      now.getUTCMonth(),
-      now.getUTCDate()
+      moscow.getUTCFullYear(),
+      moscow.getUTCMonth(),
+      moscow.getUTCDate()
     );
-    const tbSec = (info.tb ?? 0) * 900;
+    const tbSec = (info.tb ?? 0) * 60;
     const moscowOffset = 3 * 3600;
     const tocMs = utcMidnight + (tbSec - moscowOffset) * 1e3;
     return {
@@ -660,6 +664,7 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
   computeVisibility,
   ecefToAzEl,
   ecefToGeodetic,
+  ephInfoToEphemeris,
   geodeticToEcef,
   glonassPosition,
   keplerPosition,

package/dist/orbit.d.cts

@@ -108,6 +108,8 @@ interface AllPositionsData {
  * When rxPos is provided, az/el are computed; otherwise they are 0.
  */
 declare function computeAllPositions(ephemerides: Ephemeris[], times: number[], rxPos?: [number, number, number]): AllPositionsData;
+/** Convert an RTCM3 EphemerisInfo to the Ephemeris type used by orbit computation. */
+declare function ephInfoToEphemeris(info: EphemerisInfo): Ephemeris | null;
 /**
  * Compute current satellite az/el from live RTCM3 ephemeris data.
  * Returns array of SatAzEl for all satellites with valid ephemeris.
@@ -151,4 +153,4 @@ interface VisibilityResult {
  */
 declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number): VisibilityResult;
 
-export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };
+export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };

package/dist/orbit.d.ts

@@ -108,6 +108,8 @@ interface AllPositionsData {
  * When rxPos is provided, az/el are computed; otherwise they are 0.
  */
 declare function computeAllPositions(ephemerides: Ephemeris[], times: number[], rxPos?: [number, number, number]): AllPositionsData;
+/** Convert an RTCM3 EphemerisInfo to the Ephemeris type used by orbit computation. */
+declare function ephInfoToEphemeris(info: EphemerisInfo): Ephemeris | null;
 /**
  * Compute current satellite az/el from live RTCM3 ephemeris data.
  * Returns array of SatAzEl for all satellites with valid ephemeris.
@@ -151,4 +153,4 @@ interface VisibilityResult {
  */
 declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number): VisibilityResult;
 
-export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };
+export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };

package/dist/orbit.js

@@ -5,11 +5,12 @@ import {
   computeSatPosition,
   computeVisibility,
   ecefToAzEl,
+  ephInfoToEphemeris,
   glonassPosition,
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-ZCNXERQ2.js";
+} from "./chunk-PQYBTE42.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
@@ -17,6 +18,7 @@ import {
 } from "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-LEEU5OIO.js";
+import "./chunk-FIEWO4J4.js";
 export {
   computeAllPositions,
   computeDop,
@@ -25,6 +27,7 @@ export {
   computeVisibility,
   ecefToAzEl,
   ecefToGeodetic,
+  ephInfoToEphemeris,
   geodeticToEcef,
   glonassPosition,
   keplerPosition,

package/dist/positioning.cjs

@@ -63,6 +63,10 @@ var START_GPS_TIME = /* @__PURE__ */ new Date("1980-01-06T00:00:00Z");
 var START_BDS_TIME = /* @__PURE__ */ new Date("2006-01-01T00:00:14Z");
 var START_NTP_TIME = /* @__PURE__ */ new Date("1900-01-01T00:00:00Z");
 
+// src/constants/gnss.ts
+var C_LIGHT = 299792458;
+var OMEGA_E = 72921151467e-15;
+
 // src/time/gnss.ts
 function getTaiDate(date) {
   return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
@@ -150,7 +154,6 @@ var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
-var OMEGA_E = 72921151467e-15;
 var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
 var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
@@ -366,9 +369,6 @@ function computeSatPosition(eph, timeMs) {
 var GPS_EPOCH_MS = START_GPS_TIME.getTime();
 var BDS_EPOCH_MS = START_BDS_TIME.getTime();
 
-// src/constants/gnss.ts
-var C_LIGHT = 299792458;
-
 // src/positioning/index.ts
 var GM_GPS2 = 3986005e8;
 var F_REL = -4442807633e-19;
@@ -400,8 +400,7 @@ function tropoDelay(elevationRad) {
   return 2.47 / (sinEl + 0.0121);
 }
 function sagnac(pos, travelTimeS) {
-  const OMEGA_E2 = 72921151467e-15;
-  const theta = OMEGA_E2 * travelTimeS;
+  const theta = OMEGA_E * travelTimeS;
   const c = Math.cos(theta);
   const s = Math.sin(theta);
   return [pos.x * c + pos.y * s, -pos.x * s + pos.y * c, pos.z];
@@ -429,7 +428,8 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
     elevationMaskDeg = 10,
     troposphere = true,
     maxIterations = 15,
-    convergenceM = 1e-4
+    convergenceM = 1e-4,
+    tgd = true
   } = opts;
   const all = [...pseudoranges.keys()].filter((p) => ephemerides.has(p));
   const minSats = (list) => 3 + new Set(list.map((p) => p[0])).size;
@@ -460,9 +460,11 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         const sys = prn[0];
         const clk = clock2.get(sys);
         const tTx = timeMs - psr / C_LIGHT * 1e3;
-        const sat = computeSatPosition(eph, tTx);
+        const dtsClock = satClockCorrection(eph, tTx);
+        const sat = computeSatPosition(eph, tTx - dtsClock * 1e3);
         if (!Number.isFinite(sat.x)) continue;
-        const dts = satClockCorrection(eph, tTx);
+        const isKepler = eph.system !== "R" && eph.system !== "S";
+        const dts = dtsClock - (tgd && isKepler ? eph.tgd : 0);
         const travel = Math.hypot(sat.x - x2, sat.y - y2, sat.z - z2) / C_LIGHT;
         const [sx, sy, sz] = sagnac(sat, travel);
         const rho = Math.hypot(sx - x2, sy - y2, sz - z2);

package/dist/positioning.d.cts

@@ -9,9 +9,9 @@ import './ecef-CF0uAysr.cjs';
  *
  * Weighted iterative least squares with one receiver-clock unknown per
  * constellation (absorbing inter-system time offsets), satellite clock
- * polynomial + relativistic correction, Sagnac (Earth-rotation)
- * correction, elevation masking/weighting, and a simple tropospheric
- * model. Ionospheric delay is NOT modelled — use the iono-free
+ * 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.
  */
@@ -25,6 +25,14 @@ interface SppOptions {
     maxIterations?: number;
     /** Convergence threshold on the position update (m). Default 1e-4. */
     convergenceM?: number;
+    /**
+     * Apply the broadcast group-delay correction (GPS TGD, Galileo
+     * BGD E5a/E1, BeiDou TGD1) to the satellite clock. Correct for
+     * single-frequency measurements on the primary frequency (C1C/E1/B1I)
+     * — the default. Disable when feeding iono-free combinations, whose
+     * reference the broadcast clock already matches. Default true.
+     */
+    tgd?: boolean;
 }
 interface SppSolution {
     /** Receiver position, ECEF meters. */

package/dist/positioning.d.ts

@@ -9,9 +9,9 @@ import './ecef-CF0uAysr.js';
  *
  * Weighted iterative least squares with one receiver-clock unknown per
  * constellation (absorbing inter-system time offsets), satellite clock
- * polynomial + relativistic correction, Sagnac (Earth-rotation)
- * correction, elevation masking/weighting, and a simple tropospheric
- * model. Ionospheric delay is NOT modelled — use the iono-free
+ * 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.
  */
@@ -25,6 +25,14 @@ interface SppOptions {
     maxIterations?: number;
     /** Convergence threshold on the position update (m). Default 1e-4. */
     convergenceM?: number;
+    /**
+     * Apply the broadcast group-delay correction (GPS TGD, Galileo
+     * BGD E5a/E1, BeiDou TGD1) to the satellite clock. Correct for
+     * single-frequency measurements on the primary frequency (C1C/E1/B1I)
+     * — the default. Disable when feeding iono-free combinations, whose
+     * reference the broadcast clock already matches. Default true.
+     */
+    tgd?: boolean;
 }
 interface SppSolution {
     /** Receiver position, ECEF meters. */

package/dist/positioning.js

@@ -1,204 +1,14 @@
 import {
-  computeDop,
-  computeSatPosition,
-  ecefToAzEl
-} from "./chunk-ZCNXERQ2.js";
+  ionoFree,
+  satClockCorrection,
+  solveSpp
+} from "./chunk-PX7TPPTQ.js";
+import "./chunk-PQYBTE42.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-LEEU5OIO.js";
-import {
-  C_LIGHT
-} from "./chunk-W5WKEV7U.js";
-
-// src/positioning/index.ts
-var GM_GPS = 3986005e8;
-var F_REL = -4442807633e-19;
-function satClockCorrection(eph, tMs) {
-  if (eph.system === "R" || eph.system === "S") {
-    const dt2 = (tMs - eph.tocDate.getTime()) / 1e3;
-    return eph.tauN + eph.gammaN * dt2;
-  }
-  const k = eph;
-  const dt = (tMs - k.tocDate.getTime()) / 1e3;
-  const poly = k.af0 + k.af1 * dt + k.af2 * dt * dt;
-  const a = k.sqrtA * k.sqrtA;
-  const n = Math.sqrt(GM_GPS / (a * a * a)) + k.deltaN;
-  const tk = dt;
-  const Mk = k.m0 + n * tk;
-  let Ek = Mk;
-  for (let i = 0; i < 8; i++) {
-    Ek = Mk + k.e * Math.sin(Ek);
-  }
-  const rel = F_REL * k.e * k.sqrtA * Math.sin(Ek);
-  return poly + rel;
-}
-function ionoFree(p1, p2, f1, f2) {
-  const g = f1 * f1 / (f1 * f1 - f2 * f2);
-  return g * p1 - (g - 1) * p2;
-}
-function tropoDelay(elevationRad) {
-  const sinEl = Math.sin(elevationRad);
-  return 2.47 / (sinEl + 0.0121);
-}
-function sagnac(pos, travelTimeS) {
-  const OMEGA_E = 72921151467e-15;
-  const theta = OMEGA_E * travelTimeS;
-  const c = Math.cos(theta);
-  const s = Math.sin(theta);
-  return [pos.x * c + pos.y * s, -pos.x * s + pos.y * c, pos.z];
-}
-function solveLinear(A, b) {
-  const n = b.length;
-  const M = A.map((row, i) => [...row, b[i]]);
-  for (let col = 0; col < n; col++) {
-    let pivot = col;
-    for (let r = col + 1; r < n; r++) {
-      if (Math.abs(M[r][col]) > Math.abs(M[pivot][col])) pivot = r;
-    }
-    if (Math.abs(M[pivot][col]) < 1e-12) return null;
-    [M[col], M[pivot]] = [M[pivot], M[col]];
-    for (let r = 0; r < n; r++) {
-      if (r === col) continue;
-      const f = M[r][col] / M[col][col];
-      for (let c = col; c <= n; c++) M[r][c] -= f * M[col][c];
-    }
-  }
-  return M.map((row, i) => row[n] / M[i][i]);
-}
-function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
-  const {
-    elevationMaskDeg = 10,
-    troposphere = true,
-    maxIterations = 15,
-    convergenceM = 1e-4
-  } = opts;
-  const all = [...pseudoranges.keys()].filter((p) => ephemerides.has(p));
-  const minSats = (list) => 3 + new Set(list.map((p) => p[0])).size;
-  if (all.length < minSats(all)) return null;
-  const gaussNewton = (prns) => {
-    const systems = [...new Set(prns.map((p) => p[0]))].sort();
-    const dim = 3 + systems.length;
-    const sysIndex = new Map(systems.map((s, i) => [s, 3 + i]));
-    let x2 = 0;
-    let y2 = 0;
-    let z2 = 0;
-    const clock2 = new Map(systems.map((s) => [s, 0]));
-    const residuals2 = {};
-    let iterations2 = 0;
-    let converged2 = false;
-    for (let iter = 0; iter < maxIterations; iter++) {
-      iterations2 = iter + 1;
-      const HtH = Array.from(
-        { length: dim },
-        () => new Array(dim).fill(0)
-      );
-      const Htv = new Array(dim).fill(0);
-      let rows = 0;
-      for (const k of Object.keys(residuals2)) delete residuals2[k];
-      for (const prn of prns) {
-        const psr = pseudoranges.get(prn);
-        const eph = ephemerides.get(prn);
-        const sys = prn[0];
-        const clk = clock2.get(sys);
-        const tTx = timeMs - psr / C_LIGHT * 1e3;
-        const sat = computeSatPosition(eph, tTx);
-        if (!Number.isFinite(sat.x)) continue;
-        const dts = satClockCorrection(eph, tTx);
-        const travel = Math.hypot(sat.x - x2, sat.y - y2, sat.z - z2) / C_LIGHT;
-        const [sx, sy, sz] = sagnac(sat, travel);
-        const rho = Math.hypot(sx - x2, sy - y2, sz - z2);
-        const ux = (x2 - sx) / rho;
-        const uy = (y2 - sy) / rho;
-        const uz = (z2 - sz) / rho;
-        let elev = Math.PI / 2;
-        let weight = 1;
-        const positionSane = x2 * x2 + y2 * y2 + z2 * z2 > 1e12;
-        if (positionSane) {
-          const azel = ecefToAzEl(x2, y2, z2, sx, sy, sz);
-          elev = azel.el;
-          if (iter >= 2 && elev < elevationMaskDeg * Math.PI / 180) continue;
-          const sinEl = Math.max(Math.sin(elev), 0.1);
-          weight = sinEl * sinEl;
-        }
-        const tropo = troposphere && positionSane ? tropoDelay(elev) : 0;
-        const predicted = rho + clk - C_LIGHT * dts + tropo;
-        const v = psr - predicted;
-        const h = new Array(dim).fill(0);
-        h[0] = ux;
-        h[1] = uy;
-        h[2] = uz;
-        h[sysIndex.get(sys)] = 1;
-        for (let i = 0; i < dim; i++) {
-          for (let j = 0; j < dim; j++) HtH[i][j] += weight * h[i] * h[j];
-          Htv[i] += weight * h[i] * v;
-        }
-        residuals2[prn] = v;
-        rows++;
-      }
-      if (rows < dim) return null;
-      const dx = solveLinear(HtH, Htv);
-      if (!dx) return null;
-      x2 += dx[0];
-      y2 += dx[1];
-      z2 += dx[2];
-      for (const s of systems)
-        clock2.set(s, clock2.get(s) + dx[sysIndex.get(s)]);
-      if (Math.hypot(dx[0], dx[1], dx[2]) < convergenceM) {
-        converged2 = true;
-        break;
-      }
-    }
-    return { x: x2, y: y2, z: z2, clock: clock2, residuals: residuals2, iterations: iterations2, converged: converged2 };
-  };
-  const REJECT_THRESHOLD_M = 50;
-  let candidates = all;
-  const rejected = [];
-  let inner = null;
-  for (let round = 0; round <= all.length; round++) {
-    inner = gaussNewton(candidates);
-    if (!inner) return null;
-    const vals = Object.entries(inner.residuals);
-    let worst = null;
-    let worstAbs = REJECT_THRESHOLD_M;
-    for (const [prn, v] of vals) {
-      if (Math.abs(v) > worstAbs) {
-        worstAbs = Math.abs(v);
-        worst = prn;
-      }
-    }
-    if (!worst) break;
-    const remaining = candidates.filter((p) => p !== worst);
-    if (remaining.length < minSats(remaining)) break;
-    rejected.push(worst);
-    candidates = remaining;
-  }
-  if (!inner) return null;
-  const { x, y, z, clock, residuals, iterations, converged } = inner;
-  const used = [];
-  const azels = [];
-  for (const prn of candidates) {
-    if (!(prn in residuals)) continue;
-    const eph = ephemerides.get(prn);
-    const sat = computeSatPosition(eph, timeMs);
-    if (!Number.isFinite(sat.x)) continue;
-    const azel = ecefToAzEl(x, y, z, sat.x, sat.y, sat.z);
-    if (azel.el >= elevationMaskDeg * Math.PI / 180) {
-      used.push(prn);
-      azels.push(azel);
-    }
-  }
-  return {
-    position: [x, y, z],
-    clockBias: Object.fromEntries(clock),
-    usedSatellites: used,
-    rejectedSatellites: rejected,
-    residuals,
-    dop: computeDop(azels),
-    iterations,
-    converged
-  };
-}
+import "./chunk-FIEWO4J4.js";
 export {
   ionoFree,
   satClockCorrection,

package/dist/rinex.js

@@ -16,8 +16,8 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-OZCYOM5D.js";
-import "./chunk-W5WKEV7U.js";
+} from "./chunk-3U5AX7PY.js";
+import "./chunk-FIEWO4J4.js";
 export {
   EMPTY_WARNINGS,
   SYSTEM_ORDER,

package/dist/rtcm3.cjs

@@ -1319,6 +1319,7 @@ var RTCM3_MESSAGE_NAMES = {
   1031: "GLONASS Network RTK Residual",
   1032: "Physical Reference Station",
   1033: "Receiver+Antenna Descriptor",
+  1041: "NavIC Ephemeris",
   1042: "BeiDou Ephemeris",
   1043: "SBAS Ephemeris",
   1044: "QZSS Ephemeris",
@@ -1383,6 +1384,8 @@ var RTCM3_MESSAGE_NAMES = {
   1064: "GLONASS SSR Clock",
   1065: "GLONASS SSR Code Bias",
   1066: "GLONASS SSR Orbit+Clock",
+  // Biases
+  1230: "GLONASS Code-Phase Biases",
   // IGS SSR
   4076: "IGS SSR"
 };
@@ -1391,6 +1394,7 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1009 && msgType <= 1012) return "GLONASS";
   if (msgType === 1019) return "GPS";
   if (msgType === 1020) return "GLONASS";
+  if (msgType === 1041) return "NavIC";
   if (msgType === 1042) return "BeiDou";
   if (msgType === 1043) return "SBAS";
   if (msgType === 1044) return "QZSS";
@@ -1402,6 +1406,7 @@ function rtcm3Constellation(msgType) {
   if (msgType >= 1111 && msgType <= 1117) return "QZSS";
   if (msgType >= 1121 && msgType <= 1127) return "BeiDou";
   if (msgType >= 1131 && msgType <= 1137) return "NavIC";
+  if (msgType === 1230) return "GLONASS";
   return null;
 }
 function createStreamStats() {

package/dist/rtcm3.js

@@ -15,9 +15,9 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-7EEWQ5DU.js";
+} from "./chunk-65CQINSB.js";
 import "./chunk-LEEU5OIO.js";
-import "./chunk-W5WKEV7U.js";
+import "./chunk-FIEWO4J4.js";
 export {
   BitReader,
   RTCM3_MESSAGE_NAMES,