ootk 4.0.9 → 5.0.0

package/dist/observation/ObservationUtils.js.map

@@ -1 +1 @@
-{"version":3,"file":"ObservationUtils.js","sourceRoot":"","sources":["../../src/observation/ObservationUtils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AAEnC,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE;IAC7D,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAEhB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9C,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,qBAAqB,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,IAAY,EAAE,OAAO,GAAG,KAAK,EAAU,EAAE,CACrG,CAAC,OAAO,CAAC,CAAC,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;AAEtD,MAAM,CAAC,MAAM,0BAA0B,GAAG,CAAC,MAAgB,EAAU,EAAE;IACrE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAEnE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAEpB,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7B,CAAC;IAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;AAC5B,CAAC,CAAC"}
+{"version":3,"file":"ObservationUtils.js","sourceRoot":"","sources":["../../src/observation/ObservationUtils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAuB,QAAQ,EAAyC,MAAM,EAAE,MAAM,YAAY,CAAC;AAE1G,MAAM,CAAC,MAAM,eAAe,GAAG,CAAC,EAAW,EAAE,GAAY,EAAE,CAAa,EAAwB,EAAE;IAChG,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAEzB,OAAO,IAAI,QAAQ,CAAC,CAAC,GAAG,EAAE,GAAG,EAAgB,EAAE,CAAC,GAAG,EAAE,GAAG,EAAgB,EAAE,CAAC,GAAG,EAAgB,CAAC,CAAC;AAClG,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,eAAe,GAAG,CAC7B,EAAW,EACX,GAAY,EACZ,CAAa,EACb,KAAuB,EACvB,MAAwB,EACxB,IAAyB,EACM,EAAE;IACjC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAEzB,OAAO,IAAI,QAAQ,CACjB,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,MAAM,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,KAA4B,EAClF,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,MAAM,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,KAA4B,EAClF,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,MAA6B,CACnD,CAAC;AACJ,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE;IAC7D,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAEhB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9C,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,qBAAqB,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,IAAY,EAAE,OAAO,GAAG,KAAK,EAAU,EAAE,CACrG,CAAC,OAAO,CAAC,CAAC,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;AAEtD,MAAM,CAAC,MAAM,0BAA0B,GAAG,CAAC,MAAgB,EAAU,EAAE;IACrE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAEnE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAEpB,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7B,CAAC;IAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;AAC5B,CAAC,CAAC"}

package/dist/observation/RAE.d.ts

@@ -0,0 +1,111 @@
+/**
+ * @author Theodore Kruczek.
+ * @license MIT
+ * @copyright (c) 2022-2025 Theodore Kruczek Permission is
+ * hereby granted, free of charge, to any person obtaining a copy of this
+ * software and associated documentation files (the "Software"), to deal in the
+ * Software without restriction, including without limitation the rights to use,
+ * copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
+ * the Software, and to permit persons to whom the Software is furnished to do
+ * so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ */
+import { J2000 } from '../coordinate/J2000.js';
+import { AngularDistanceMethod } from '../enums/AngularDistanceMethod.js';
+import { Degrees, Kilometers, Radians } from '../main.js';
+import { Vector3D } from '../operations/Vector3D.js';
+import { EpochUTC } from '../time/EpochUTC.js';
+export declare class RAE {
+    epoch: EpochUTC;
+    rng: Kilometers;
+    azRad: Radians;
+    elRad: Radians;
+    /** The range rate of the satellite relative to the observer in kilometers per second. */
+    rngRate?: number | undefined;
+    /** The azimuth rate of the satellite relative to the observer in radians per second. */
+    azRateRad?: number | undefined;
+    /** The elevation rate of the satellite relative to the observer in radians per second. */
+    elRateRad?: number | undefined;
+    constructor(epoch: EpochUTC, rng: Kilometers, azRad: Radians, elRad: Radians, 
+    /** The range rate of the satellite relative to the observer in kilometers per second. */
+    rngRate?: number | undefined, 
+    /** The azimuth rate of the satellite relative to the observer in radians per second. */
+    azRateRad?: number | undefined, 
+    /** The elevation rate of the satellite relative to the observer in radians per second. */
+    elRateRad?: number | undefined);
+    static fromDegrees(epoch: EpochUTC, range: Kilometers, azimuth: Degrees, elevation: Degrees, rangeRate?: number, azimuthRate?: number, elevationRate?: number): RAE;
+    /**
+     * Create a [Razel] object from an inertial [state] and [site] vector.
+     * @param state The inertial [state] vector.
+     * @param site The observer [site] vector.
+     * @returns A new [Razel] object.
+     */
+    static fromStateVector(state: J2000, site: J2000): RAE;
+    /**
+     * Gets the azimuth in degrees.
+     * @returns The azimuth in degrees.
+     */
+    get az(): Degrees;
+    /**
+     * Gets the elevation angle in degrees.
+     * @returns The elevation angle in degrees.
+     */
+    get el(): Degrees;
+    /**
+     * Gets the azimuth rate in degrees per second.
+     * @returns The azimuth rate in degrees per second, or undefined if it is not available.
+     */
+    get azRate(): number | undefined;
+    /**
+     * Gets the elevation rate in degrees per second.
+     * @returns The elevation rate in degrees per second, or undefined if the elevation rate is not set.
+     */
+    get elRate(): number | undefined;
+    toString(): string;
+    /**
+     * Return the position relative to the observer [site].
+     *
+     * An optional azimuth [az] _(rad)_ and elevation [el] _(rad)_ value can be
+     * passed to override the values contained in this observation.
+     * @param site The observer [site].
+     * @param azRad Azimuth _(rad)_.
+     * @param elRad Elevation _(rad)_.
+     * @returns A [Vector3D] object.
+     */
+    position(site: J2000, azRad?: Radians, elRad?: Radians): Vector3D<Kilometers>;
+    /**
+     * Convert this observation into a [J2000] state vector.
+     *
+     * This will throw an error if the [rangeRate], [elevationRate], or
+     * [azimuthRate] are not defined.
+     * @param site The observer [site].
+     * @returns A [J2000] state vector.
+     */
+    toStateVector(site: J2000): J2000;
+    /**
+     * Calculate the angular distance _(rad)_ between this and another [Razel]
+     * object.
+     * @param razel The other [Razel] object.
+     * @param method The angular distance method to use.
+     * @returns The angular distance _(rad)_.
+     */
+    angle(razel: RAE, method?: AngularDistanceMethod): number;
+    /**
+     * Calculate the angular distance _(°)_ between this and another [Razel]
+     * object.
+     * @param razel The other [Razel] object.
+     * @param method The angular distance method to use.
+     * @returns The angular distance _(°)_.
+     */
+    angleDegrees(razel: RAE, method?: AngularDistanceMethod): number;
+}

package/dist/observation/RAE.js

@@ -0,0 +1,215 @@
+/**
+ * @author Theodore Kruczek.
+ * @license MIT
+ * @copyright (c) 2022-2025 Theodore Kruczek Permission is
+ * hereby granted, free of charge, to any person obtaining a copy of this
+ * software and associated documentation files (the "Software"), to deal in the
+ * Software without restriction, including without limitation the rights to use,
+ * copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
+ * the Software, and to permit persons to whom the Software is furnished to do
+ * so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ */
+/* eslint-disable no-undefined */
+import { ITRF } from '../coordinate/ITRF.js';
+import { AngularDistanceMethod } from '../enums/AngularDistanceMethod.js';
+import { Vector3D } from '../operations/Vector3D.js';
+import { DEG2RAD, halfPi, RAD2DEG, TAU } from '../utils/constants.js';
+import { angularDistance } from '../utils/functions.js';
+// / Range, azimuth, and elevation.
+export class RAE {
+    epoch;
+    rng;
+    azRad;
+    elRad;
+    rngRate;
+    azRateRad;
+    elRateRad;
+    constructor(epoch, rng, azRad, elRad, 
+    /** The range rate of the satellite relative to the observer in kilometers per second. */
+    rngRate, 
+    /** The azimuth rate of the satellite relative to the observer in radians per second. */
+    azRateRad, 
+    /** The elevation rate of the satellite relative to the observer in radians per second. */
+    elRateRad) {
+        this.epoch = epoch;
+        this.rng = rng;
+        this.azRad = azRad;
+        this.elRad = elRad;
+        this.rngRate = rngRate;
+        this.azRateRad = azRateRad;
+        this.elRateRad = elRateRad;
+        // Do nothing
+    }
+    // / Create a new [Razel] object, using degrees for the angular values.
+    static fromDegrees(epoch, range, azimuth, elevation, rangeRate, azimuthRate, elevationRate) {
+        const azimuthRateRad = azimuthRate ? azimuthRate * DEG2RAD : undefined;
+        const elevationRateRad = elevationRate ? elevationRate * DEG2RAD : undefined;
+        return new RAE(epoch, range, (azimuth * DEG2RAD), (elevation * DEG2RAD), rangeRate, azimuthRateRad, elevationRateRad);
+    }
+    /**
+     * Create a [Razel] object from an inertial [state] and [site] vector.
+     * @param state The inertial [state] vector.
+     * @param site The observer [site] vector.
+     * @returns A new [Razel] object.
+     */
+    static fromStateVector(state, site) {
+        const stateEcef = state.toITRF();
+        const siteEcef = site.toITRF();
+        const po2 = halfPi;
+        const r = stateEcef.position.subtract(siteEcef.position);
+        const rDot = stateEcef.velocity;
+        const geo = siteEcef.toGeodetic();
+        const p = r.rotZ(geo.lon).rotY((po2 - geo.lat));
+        const pDot = rDot.rotZ(geo.lon).rotY((po2 - geo.lat));
+        const pS = p.x;
+        const pE = p.y;
+        const pZ = p.z;
+        const pSDot = pDot.x;
+        const pEDot = pDot.y;
+        const pZDot = pDot.z;
+        const pMag = p.magnitude();
+        const pSEMag = Math.sqrt(pS * pS + pE * pE);
+        const elevation = Math.asin(pZ / pMag);
+        let azimuth;
+        if (elevation !== po2) {
+            azimuth = Math.atan2(-pE, pS) + Math.PI;
+        }
+        else {
+            azimuth = Math.atan2(-pEDot, pSDot) + Math.PI;
+        }
+        const rangeRate = p.dot(pDot) / pMag;
+        const azimuthRate = (pSDot * pE - pEDot * pS) / (pS * pS + pE * pE);
+        const elevationRate = (pZDot - rangeRate * Math.sin(elevation)) / pSEMag;
+        return new RAE(state.epoch, pMag, (azimuth % TAU), elevation, rangeRate, azimuthRate, elevationRate);
+    }
+    /**
+     * Gets the azimuth in degrees.
+     * @returns The azimuth in degrees.
+     */
+    get az() {
+        return this.azRad * RAD2DEG;
+    }
+    /**
+     * Gets the elevation angle in degrees.
+     * @returns The elevation angle in degrees.
+     */
+    get el() {
+        return this.elRad * RAD2DEG;
+    }
+    /**
+     * Gets the azimuth rate in degrees per second.
+     * @returns The azimuth rate in degrees per second, or undefined if it is not available.
+     */
+    get azRate() {
+        return this.azRateRad ? this.azRateRad * RAD2DEG : undefined;
+    }
+    /**
+     * Gets the elevation rate in degrees per second.
+     * @returns The elevation rate in degrees per second, or undefined if the elevation rate is not set.
+     */
+    get elRate() {
+        return this.elRateRad ? this.elRateRad * RAD2DEG : undefined;
+    }
+    toString() {
+        return [
+            '[RazEl]',
+            `  Epoch:     ${this.epoch}`,
+            `  Azimuth:   ${this.az.toFixed(4)}°`,
+            `  Elevation: ${this.el.toFixed(4)}°`,
+            `  Range:     ${this.rng.toFixed(3)} km`,
+        ].join('\n');
+    }
+    /**
+     * Return the position relative to the observer [site].
+     *
+     * An optional azimuth [az] _(rad)_ and elevation [el] _(rad)_ value can be
+     * passed to override the values contained in this observation.
+     * @param site The observer [site].
+     * @param azRad Azimuth _(rad)_.
+     * @param elRad Elevation _(rad)_.
+     * @returns A [Vector3D] object.
+     */
+    position(site, azRad, elRad) {
+        const ecef = site.toITRF();
+        const geo = ecef.toGeodetic();
+        const po2 = halfPi;
+        const newAz = azRad ?? this.azRad;
+        const newEl = elRad ?? this.elRad;
+        const sAz = Math.sin(newAz);
+        const cAz = Math.cos(newAz);
+        const sEl = Math.sin(newEl);
+        const cEl = Math.cos(newEl);
+        const pSez = new Vector3D((-this.rng * cEl * cAz), (this.rng * cEl * sAz), (this.rng * sEl));
+        const rEcef = pSez
+            .rotY(-(po2 - geo.lat))
+            .rotZ(-geo.lon)
+            .add(ecef.position);
+        return new ITRF(this.epoch, rEcef, Vector3D.origin).toJ2000().position;
+    }
+    /**
+     * Convert this observation into a [J2000] state vector.
+     *
+     * This will throw an error if the [rangeRate], [elevationRate], or
+     * [azimuthRate] are not defined.
+     * @param site The observer [site].
+     * @returns A [J2000] state vector.
+     */
+    toStateVector(site) {
+        // If the rates are not defined then assume stationary
+        this.rngRate ??= 0;
+        this.elRateRad ??= 0;
+        this.azRateRad ??= 0;
+        const ecef = site.toITRF();
+        const geo = ecef.toGeodetic();
+        const po2 = halfPi;
+        const sAz = Math.sin(this.azRad);
+        const cAz = Math.cos(this.azRad);
+        const sEl = Math.sin(this.elRad);
+        const cEl = Math.cos(this.elRad);
+        const pSez = new Vector3D((-this.rng * cEl * cAz), (this.rng * cEl * sAz), (this.rng * sEl));
+        const pDotSez = new Vector3D((-this.rngRate * cEl * cAz +
+            this.rng * sEl * cAz * this.elRateRad +
+            this.rng * cEl * sAz * this.azRateRad), (this.rngRate * cEl * sAz -
+            this.rng * sEl * sAz * this.elRateRad +
+            this.rng * cEl * cAz * this.azRateRad), (this.rngRate * sEl + this.rng * cEl * this.elRateRad));
+        const pEcef = pSez.rotY(-(po2 - geo.lat)).rotZ(-geo.lon);
+        const pDotEcef = pDotSez
+            .rotY(-(po2 - geo.lat))
+            // TODO: #13 Intermediate unit type is incorrect.
+            .rotZ(-geo.lon);
+        const rEcef = pEcef.add(ecef.position);
+        return new ITRF(this.epoch, rEcef, pDotEcef).toJ2000();
+    }
+    /**
+     * Calculate the angular distance _(rad)_ between this and another [Razel]
+     * object.
+     * @param razel The other [Razel] object.
+     * @param method The angular distance method to use.
+     * @returns The angular distance _(rad)_.
+     */
+    angle(razel, method = AngularDistanceMethod.Cosine) {
+        return angularDistance(this.azRad, this.elRad, razel.azRad, razel.elRad, method);
+    }
+    /**
+     * Calculate the angular distance _(°)_ between this and another [Razel]
+     * object.
+     * @param razel The other [Razel] object.
+     * @param method The angular distance method to use.
+     * @returns The angular distance _(°)_.
+     */
+    angleDegrees(razel, method = AngularDistanceMethod.Cosine) {
+        return this.angle(razel, method) * RAD2DEG;
+    }
+}
+//# sourceMappingURL=RAE.js.map

package/dist/observation/RAE.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RAE.js","sourceRoot":"","sources":["../../src/observation/RAE.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,iCAAiC;AACjC,OAAO,EAAE,IAAI,EAAE,MAAM,uBAAuB,CAAC;AAE7C,OAAO,EAAE,qBAAqB,EAAE,MAAM,mCAAmC,CAAC;AAE1E,OAAO,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AAErD,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,uBAAuB,CAAC;AACtE,OAAO,EAAE,eAAe,EAAE,MAAM,uBAAuB,CAAC;AAExD,mCAAmC;AACnC,MAAM,OAAO,GAAG;IAEL;IACA;IACA;IACA;IAEA;IAEA;IAEA;IAVT,YACS,KAAe,EACf,GAAe,EACf,KAAc,EACd,KAAc;IACrB,yFAAyF;IAClF,OAAgB;IACvB,wFAAwF;IACjF,SAAkB;IACzB,0FAA0F;IACnF,SAAkB;QATlB,UAAK,GAAL,KAAK,CAAU;QACf,QAAG,GAAH,GAAG,CAAY;QACf,UAAK,GAAL,KAAK,CAAS;QACd,UAAK,GAAL,KAAK,CAAS;QAEd,YAAO,GAAP,OAAO,CAAS;QAEhB,cAAS,GAAT,SAAS,CAAS;QAElB,cAAS,GAAT,SAAS,CAAS;QAEzB,aAAa;IACf,CAAC;IAED,uEAAuE;IACvE,MAAM,CAAC,WAAW,CAChB,KAAe,EACf,KAAiB,EACjB,OAAgB,EAChB,SAAkB,EAClB,SAAkB,EAClB,WAAoB,EACpB,aAAsB;QAEtB,MAAM,cAAc,GAAG,WAAW,CAAC,CAAC,CAAC,WAAW,GAAG,OAAO,CAAC,CAAC,CAAC,SAAS,CAAC;QACvE,MAAM,gBAAgB,GAAG,aAAa,CAAC,CAAC,CAAC,aAAa,GAAG,OAAO,CAAC,CAAC,CAAC,SAAS,CAAC;QAE7E,OAAO,IAAI,GAAG,CACZ,KAAK,EACL,KAAK,EACL,CAAC,OAAO,GAAG,OAAO,CAAY,EAC9B,CAAC,SAAS,GAAG,OAAO,CAAY,EAChC,SAAS,EACT,cAAc,EACd,gBAAgB,CACjB,CAAC;IACJ,CAAC;IAED;;;;;OAKG;IACH,MAAM,CAAC,eAAe,CAAC,KAAY,EAAE,IAAW;QAC9C,MAAM,SAAS,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;QACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,SAAS,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;QACzD,MAAM,IAAI,GAAG,SAAS,CAAC,QAAQ,CAAC;QAChC,MAAM,GAAG,GAAG,QAAQ,CAAC,UAAU,EAAE,CAAC;QAClC,MAAM,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAY,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAY,CAAC,CAAC;QACjE,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;QACrB,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;QACrB,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5C,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC;QACvC,IAAI,OAAO,CAAC;QAEZ,IAAI,SAAS,KAAK,GAAG,EAAE,CAAC;YACtB,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC1C,CAAC;aAAM,CAAC;YACN,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAChD,CAAC;QACD,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;QACrC,MAAM,WAAW,GAAG,CAAC,KAAK,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACpE,MAAM,aAAa,GAAG,CAAC,KAAK,GAAG,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,MAAM,CAAC;QAEzE,OAAO,IAAI,GAAG,CACZ,KAAK,CAAC,KAAK,EACX,IAAI,EACJ,CAAC,OAAO,GAAG,GAAG,CAAY,EAC1B,SAAoB,EACpB,SAAS,EACT,WAAW,EACX,aAAa,CACd,CAAC;IACJ,CAAC;IAED;;;OAGG;IACH,IAAI,EAAE;QACJ,OAAO,IAAI,CAAC,KAAK,GAAG,OAAkB,CAAC;IACzC,CAAC;IAED;;;OAGG;IACH,IAAI,EAAE;QACJ,OAAO,IAAI,CAAC,KAAK,GAAG,OAAkB,CAAC;IACzC,CAAC;IAED;;;OAGG;IACH,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,GAAG,OAAO,CAAC,CAAC,CAAC,SAAS,CAAC;IAC/D,CAAC;IAED;;;OAGG;IACH,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,GAAG,OAAO,CAAC,CAAC,CAAC,SAAS,CAAC;IAC/D,CAAC;IAED,QAAQ;QACN,OAAO;YACL,SAAS;YACT,gBAAgB,IAAI,CAAC,KAAK,EAAE;YAC5B,gBAAgB,IAAI,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG;YACrC,gBAAgB,IAAI,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG;YACrC,gBAAgB,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK;SACzC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACf,CAAC;IAED;;;;;;;;;OASG;IACH,QAAQ,CAAC,IAAW,EAAE,KAAe,EAAE,KAAe;QACpD,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,MAAM,CAAC;QACnB,MAAM,KAAK,GAAG,KAAK,IAAI,IAAI,CAAC,KAAK,CAAC;QAClC,MAAM,KAAK,GAAG,KAAK,IAAI,IAAI,CAAC,KAAK,CAAC;QAClC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC5B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC5B,MAAM,IAAI,GAAG,IAAI,QAAQ,CACvB,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAe,EACrC,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAe,EACpC,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAe,CAC/B,CAAC;QACF,MAAM,KAAK,GAAG,IAAI;aACf,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAY,CAAC;aACjC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAc,CAAC;aACzB,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAEtB,OAAO,IAAI,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,QAAQ,CAAC,MAAuC,CAAC,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC;IAC1G,CAAC;IAED;;;;;;;OAOG;IACH,aAAa,CAAC,IAAW;QACvB,sDAAsD;QACtD,IAAI,CAAC,OAAO,KAAK,CAAC,CAAC;QACnB,IAAI,CAAC,SAAS,KAAK,CAAC,CAAC;QACrB,IAAI,CAAC,SAAS,KAAK,CAAC,CAAC;QAErB,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,MAAM,CAAC;QACnB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjC,MAAM,IAAI,GAAG,IAAI,QAAQ,CACvB,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAe,EACrC,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAe,EACpC,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAe,CAC/B,CAAC;QACF,MAAM,OAAO,GAAG,IAAI,QAAQ,CAC1B,CAAC,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,GAAG,GAAG;YACxB,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,SAAS;YACrC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,SAAS,CAAe,EACtD,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,GAAG,GAAG;YACvB,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,SAAS;YACrC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,SAAS,CAAe,EACtD,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,SAAS,CAAe,CACrE,CAAC;QACF,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAY,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAc,CAAC,CAAC;QAC/E,MAAM,QAAQ,GAAG,OAAO;aACrB,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAY,CAAC;YAClC,iDAAiD;aAChD,IAAI,CAAC,CAAC,GAAG,CAAC,GAAc,CAA6C,CAAC;QACzE,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAEvC,OAAO,IAAI,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC,OAAO,EAAE,CAAC;IACzD,CAAC;IAED;;;;;;OAMG;IACH,KAAK,CAAC,KAAU,EAAE,SAAgC,qBAAqB,CAAC,MAAM;QAC5E,OAAO,eAAe,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;IACnF,CAAC;IAED;;;;;;OAMG;IACH,YAAY,CAAC,KAAU,EAAE,SAAgC,qBAAqB,CAAC,MAAM;QACnF,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,OAAO,CAAC;IAC7C,CAAC;CACF"}

package/dist/observation/RadecGeocentric.d.ts

@@ -0,0 +1,111 @@
+/**
+ * @author Theodore Kruczek.
+ * @license MIT
+ * @copyright (c) 2022-2025 Theodore Kruczek Permission is
+ * hereby granted, free of charge, to any person obtaining a copy of this
+ * software and associated documentation files (the "Software"), to deal in the
+ * Software without restriction, including without limitation the rights to use,
+ * copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
+ * the Software, and to permit persons to whom the Software is furnished to do
+ * so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ */
+import { Degrees, DegreesPerSecond, Kilometers, KilometersPerSecond, Radians, RadiansPerSecond } from '../main.js';
+import { J2000 } from '../coordinate/J2000.js';
+import { AngularDistanceMethod } from '../enums/AngularDistanceMethod.js';
+import { Vector3D } from '../operations/Vector3D.js';
+import { EpochUTC } from '../time/EpochUTC.js';
+/**
+ * Represents a geocentric right ascension and declination observation.
+ *
+ * In geocentric coordinates, observations are considered from the Earth's center. This approach simplifies calculations
+ * for distant celestial objects, as it assumes a uniform observation point that ignores the observer's specific
+ * location on Earth.
+ */
+export declare class RadecGeocentric {
+    epoch: EpochUTC;
+    rightAscension: Radians;
+    declination: Radians;
+    range?: Kilometers | undefined;
+    rightAscensionRate?: (RadiansPerSecond | null) | undefined;
+    declinationRate?: (RadiansPerSecond | null) | undefined;
+    rangeRate?: (KilometersPerSecond | null) | undefined;
+    constructor(epoch: EpochUTC, rightAscension: Radians, declination: Radians, range?: Kilometers | undefined, rightAscensionRate?: (RadiansPerSecond | null) | undefined, declinationRate?: (RadiansPerSecond | null) | undefined, rangeRate?: (KilometersPerSecond | null) | undefined);
+    /**
+     * Creates a RadecGeocentric object from the given parameters in degrees.
+     * @param epoch - The epoch in UTC.
+     * @param rightAscensionDegrees - The right ascension in degrees.
+     * @param declinationDegrees - The declination in degrees.
+     * @param range - The range in kilometers (optional).
+     * @param rightAscensionRateDegrees - The right ascension rate in degrees per second (optional).
+     * @param declinationRateDegrees - The declination rate in degrees per second (optional).
+     * @param rangeRate - The range rate in kilometers per second (optional).
+     * @returns A new RadecGeocentric object.
+     */
+    static fromDegrees(epoch: EpochUTC, rightAscensionDegrees: Degrees, declinationDegrees: Degrees, range?: Kilometers, rightAscensionRateDegrees?: Degrees, declinationRateDegrees?: Degrees, rangeRate?: KilometersPerSecond): RadecGeocentric;
+    /**
+     * Creates a RadecGeocentric object from a state vector in J2000 coordinates.
+     * @param state - The J2000 state vector.
+     * @returns A new RadecGeocentric object.
+     */
+    static fromStateVector(state: J2000): RadecGeocentric;
+    /**
+     * Gets the right ascension in degrees.
+     * @returns The right ascension in degrees.
+     */
+    get rightAscensionDegrees(): Degrees;
+    /**
+     * Gets the declination in degrees.
+     * @returns The declination in degrees.
+     */
+    get declinationDegrees(): Degrees;
+    /**
+     * Gets the right ascension rate in degrees per second.
+     * @returns The right ascension rate in degrees per second, or null if it is not available.
+     */
+    get rightAscensionRateDegrees(): DegreesPerSecond | null;
+    /**
+     * Gets the rate of change of declination in degrees per second.
+     * @returns The rate of change of declination in degrees per second, or null if not available.
+     */
+    get declinationRateDegrees(): DegreesPerSecond | null;
+    /**
+     * Calculates the position vector in geocentric coordinates.
+     * @param range - The range in kilometers (optional). If not provided, it uses the default range or 1.0 kilometer.
+     * @returns The position vector in geocentric coordinates.
+     */
+    position(range?: Kilometers): Vector3D<Kilometers>;
+    /**
+     * Calculates the velocity vector of the celestial object.
+     * @param range - The range of the celestial object in kilometers. If not provided, it uses the stored range value.
+     * @param rangeRate - The range rate of the celestial object in kilometers per second.
+     * If not provided, it uses the stored range rate value.
+     * @returns The velocity vector of the celestial object in kilometers per second.
+     * @throws Error if the right ascension rate or declination rate is missing.
+     */
+    velocity(range?: Kilometers, rangeRate?: KilometersPerSecond): Vector3D<KilometersPerSecond>;
+    /**
+     * Calculates the angular distance between two celestial coordinates (RA and Dec).
+     * @param radec - The celestial coordinates to compare with.
+     * @param method - The method to use for calculating the angular distance. Default is `AngularDistanceMethod.Cosine`.
+     * @returns The angular distance between the two celestial coordinates in radians.
+     */
+    angle(radec: RadecGeocentric, method?: AngularDistanceMethod): Radians;
+    /**
+     * Calculates the angle in degrees between two RadecGeocentric objects.
+     * @param radec - The RadecGeocentric object to calculate the angle with.
+     * @param method - The method to use for calculating the angular distance. Default is AngularDistanceMethod.Cosine.
+     * @returns The angle in degrees.
+     */
+    angleDegrees(radec: RadecGeocentric, method?: AngularDistanceMethod): Degrees;
+}

package/dist/observation/RadecGeocentric.js

@@ -0,0 +1,169 @@
+/**
+ * @author Theodore Kruczek.
+ * @license MIT
+ * @copyright (c) 2022-2025 Theodore Kruczek Permission is
+ * hereby granted, free of charge, to any person obtaining a copy of this
+ * software and associated documentation files (the "Software"), to deal in the
+ * Software without restriction, including without limitation the rights to use,
+ * copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
+ * the Software, and to permit persons to whom the Software is furnished to do
+ * so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in
+ * all copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ */
+import { AngularDistanceMethod } from '../enums/AngularDistanceMethod.js';
+import { DEG2RAD, RAD2DEG, TAU } from '../utils/constants.js';
+import { angularDistance } from '../utils/functions.js';
+import { radecToPosition, radecToVelocity } from './ObservationUtils.js';
+/**
+ * Represents a geocentric right ascension and declination observation.
+ *
+ * In geocentric coordinates, observations are considered from the Earth's center. This approach simplifies calculations
+ * for distant celestial objects, as it assumes a uniform observation point that ignores the observer's specific
+ * location on Earth.
+ */
+export class RadecGeocentric {
+    epoch;
+    rightAscension;
+    declination;
+    range;
+    rightAscensionRate;
+    declinationRate;
+    rangeRate;
+    constructor(epoch, rightAscension, declination, range, rightAscensionRate, declinationRate, rangeRate) {
+        this.epoch = epoch;
+        this.rightAscension = rightAscension;
+        this.declination = declination;
+        this.range = range;
+        this.rightAscensionRate = rightAscensionRate;
+        this.declinationRate = declinationRate;
+        this.rangeRate = rangeRate;
+        // Nothing to do here.
+    }
+    /**
+     * Creates a RadecGeocentric object from the given parameters in degrees.
+     * @param epoch - The epoch in UTC.
+     * @param rightAscensionDegrees - The right ascension in degrees.
+     * @param declinationDegrees - The declination in degrees.
+     * @param range - The range in kilometers (optional).
+     * @param rightAscensionRateDegrees - The right ascension rate in degrees per second (optional).
+     * @param declinationRateDegrees - The declination rate in degrees per second (optional).
+     * @param rangeRate - The range rate in kilometers per second (optional).
+     * @returns A new RadecGeocentric object.
+     */
+    static fromDegrees(epoch, rightAscensionDegrees, declinationDegrees, range, rightAscensionRateDegrees, declinationRateDegrees, rangeRate) {
+        const rightAscensionRate = rightAscensionRateDegrees
+            ? rightAscensionRateDegrees * DEG2RAD
+            : null;
+        const declinationRate = declinationRateDegrees ? declinationRateDegrees * DEG2RAD : null;
+        return new RadecGeocentric(epoch, rightAscensionDegrees * DEG2RAD, declinationDegrees * DEG2RAD, range, rightAscensionRate, declinationRate, rangeRate);
+    }
+    /**
+     * Creates a RadecGeocentric object from a state vector in J2000 coordinates.
+     * @param state - The J2000 state vector.
+     * @returns A new RadecGeocentric object.
+     */
+    static fromStateVector(state) {
+        const rI = state.position.x;
+        const rJ = state.position.y;
+        const rK = state.position.z;
+        const vI = state.velocity.x;
+        const vJ = state.velocity.y;
+        const vK = state.velocity.z;
+        const rMag = state.position.magnitude();
+        const declination = Math.asin(rK / rMag);
+        const rIJMag = Math.sqrt(rI * rI + rJ * rJ);
+        let rightAscension;
+        if (rIJMag !== 0) {
+            rightAscension = Math.atan2(rJ, rI);
+        }
+        else {
+            rightAscension = Math.atan2(vJ, vI);
+        }
+        const rangeRate = state.position.dot(state.velocity) / rMag;
+        const rightAscensionRate = (vI * rJ - vJ * rI) / (-(rJ * rJ) - rI * rI);
+        const declinationRate = (vK - rangeRate * (rK / rMag)) / rIJMag;
+        return new RadecGeocentric(state.epoch, rightAscension % TAU, declination, rMag, rightAscensionRate, declinationRate, rangeRate);
+    }
+    /**
+     * Gets the right ascension in degrees.
+     * @returns The right ascension in degrees.
+     */
+    get rightAscensionDegrees() {
+        return this.rightAscension * RAD2DEG;
+    }
+    /**
+     * Gets the declination in degrees.
+     * @returns The declination in degrees.
+     */
+    get declinationDegrees() {
+        return this.declination * RAD2DEG;
+    }
+    /**
+     * Gets the right ascension rate in degrees per second.
+     * @returns The right ascension rate in degrees per second, or null if it is not available.
+     */
+    get rightAscensionRateDegrees() {
+        return this.rightAscensionRate ? this.rightAscensionRate * RAD2DEG : null;
+    }
+    /**
+     * Gets the rate of change of declination in degrees per second.
+     * @returns The rate of change of declination in degrees per second, or null if not available.
+     */
+    get declinationRateDegrees() {
+        return this.declinationRate ? this.declinationRate * RAD2DEG : null;
+    }
+    /**
+     * Calculates the position vector in geocentric coordinates.
+     * @param range - The range in kilometers (optional). If not provided, it uses the default range or 1.0 kilometer.
+     * @returns The position vector in geocentric coordinates.
+     */
+    position(range) {
+        const r = range ?? this.range ?? 1.0;
+        return radecToPosition(this.rightAscension, this.declination, r);
+    }
+    /**
+     * Calculates the velocity vector of the celestial object.
+     * @param range - The range of the celestial object in kilometers. If not provided, it uses the stored range value.
+     * @param rangeRate - The range rate of the celestial object in kilometers per second.
+     * If not provided, it uses the stored range rate value.
+     * @returns The velocity vector of the celestial object in kilometers per second.
+     * @throws Error if the right ascension rate or declination rate is missing.
+     */
+    velocity(range, rangeRate) {
+        if (!this.rightAscensionRate || !this.declinationRate) {
+            throw new Error('Velocity unsolvable, missing ra/dec rates.');
+        }
+        const r = range ?? this.range ?? 1.0;
+        const rd = rangeRate ?? this.rangeRate ?? 0.0;
+        return radecToVelocity(this.rightAscension, this.declination, r, this.rightAscensionRate, this.declinationRate, rd);
+    }
+    /**
+     * Calculates the angular distance between two celestial coordinates (RA and Dec).
+     * @param radec - The celestial coordinates to compare with.
+     * @param method - The method to use for calculating the angular distance. Default is `AngularDistanceMethod.Cosine`.
+     * @returns The angular distance between the two celestial coordinates in radians.
+     */
+    angle(radec, method = AngularDistanceMethod.Cosine) {
+        return angularDistance(this.rightAscension, this.declination, radec.rightAscension, radec.declination, method);
+    }
+    /**
+     * Calculates the angle in degrees between two RadecGeocentric objects.
+     * @param radec - The RadecGeocentric object to calculate the angle with.
+     * @param method - The method to use for calculating the angular distance. Default is AngularDistanceMethod.Cosine.
+     * @returns The angle in degrees.
+     */
+    angleDegrees(radec, method = AngularDistanceMethod.Cosine) {
+        return this.angle(radec, method) * RAD2DEG;
+    }
+}
+//# sourceMappingURL=RadecGeocentric.js.map

package/dist/observation/RadecGeocentric.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RadecGeocentric.js","sourceRoot":"","sources":["../../src/observation/RadecGeocentric.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAIH,OAAO,EAAE,qBAAqB,EAAE,MAAM,mCAAmC,CAAC;AAG1E,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,uBAAuB,CAAC;AAC9D,OAAO,EAAE,eAAe,EAAE,MAAM,uBAAuB,CAAC;AACxD,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,uBAAuB,CAAC;AAEzE;;;;;;GAMG;AACH,MAAM,OAAO,eAAe;IAEjB;IACA;IACA;IACA;IACA;IACA;IACA;IAPT,YACS,KAAe,EACf,cAAuB,EACvB,WAAoB,EACpB,KAAkB,EAClB,kBAA4C,EAC5C,eAAyC,EACzC,SAAsC;QANtC,UAAK,GAAL,KAAK,CAAU;QACf,mBAAc,GAAd,cAAc,CAAS;QACvB,gBAAW,GAAX,WAAW,CAAS;QACpB,UAAK,GAAL,KAAK,CAAa;QAClB,uBAAkB,GAAlB,kBAAkB,CAA0B;QAC5C,oBAAe,GAAf,eAAe,CAA0B;QACzC,cAAS,GAAT,SAAS,CAA6B;QAE7C,sBAAsB;IACxB,CAAC;IAED;;;;;;;;;;OAUG;IACH,MAAM,CAAC,WAAW,CAChB,KAAe,EACf,qBAA8B,EAC9B,kBAA2B,EAC3B,KAAkB,EAClB,yBAAmC,EACnC,sBAAgC,EAChC,SAA+B;QAE/B,MAAM,kBAAkB,GAAG,yBAAyB;YAClD,CAAC,CAAC,yBAAyB,GAAG,OAA2B;YACzD,CAAC,CAAC,IAAI,CAAC;QACT,MAAM,eAAe,GAAG,sBAAsB,CAAC,CAAC,CAAC,sBAAsB,GAAG,OAA2B,CAAC,CAAC,CAAC,IAAI,CAAC;QAE7G,OAAO,IAAI,eAAe,CACxB,KAAK,EACL,qBAAqB,GAAG,OAAkB,EAC1C,kBAAkB,GAAG,OAAkB,EACvC,KAAK,EACL,kBAAkB,EAClB,eAAe,EACf,SAAS,CACV,CAAC;IACJ,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,eAAe,CAAC,KAAY;QACjC,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC5B,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QACxC,MAAM,WAAW,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC;QACzC,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5C,IAAI,cAAc,CAAC;QAEnB,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;YACjB,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,CAAC;aAAM,CAAC;YACN,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,CAAC;QACD,MAAM,SAAS,GAAG,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC,QAAQ,CAAC,GAAG,IAA2B,CAAC;QACnF,MAAM,kBAAkB,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACxE,MAAM,eAAe,GAAG,CAAC,EAAE,GAAG,SAAS,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC;QAEhE,OAAO,IAAI,eAAe,CACxB,KAAK,CAAC,KAAK,EACX,cAAc,GAAG,GAAc,EAC/B,WAAsB,EACtB,IAAI,EACJ,kBAAsC,EACtC,eAAmC,EACnC,SAAS,CACV,CAAC;IACJ,CAAC;IAED;;;OAGG;IACH,IAAI,qBAAqB;QACvB,OAAO,IAAI,CAAC,cAAc,GAAG,OAAkB,CAAC;IAClD,CAAC;IAED;;;OAGG;IACH,IAAI,kBAAkB;QACpB,OAAO,IAAI,CAAC,WAAW,GAAG,OAAkB,CAAC;IAC/C,CAAC;IAED;;;OAGG;IACH,IAAI,yBAAyB;QAC3B,OAAO,IAAI,CAAC,kBAAkB,CAAC,CAAC,CAAC,IAAI,CAAC,kBAAkB,GAAG,OAA2B,CAAC,CAAC,CAAC,IAAI,CAAC;IAChG,CAAC;IAED;;;OAGG;IACH,IAAI,sBAAsB;QACxB,OAAO,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,eAAe,GAAG,OAA2B,CAAC,CAAC,CAAC,IAAI,CAAC;IAC1F,CAAC;IAED;;;;OAIG;IACH,QAAQ,CAAC,KAAkB;QACzB,MAAM,CAAC,GAAG,KAAK,IAAI,IAAI,CAAC,KAAK,IAAI,GAAiB,CAAC;QAEnD,OAAO,eAAe,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;IACnE,CAAC;IAED;;;;;;;OAOG;IACH,QAAQ,CAAC,KAAkB,EAAE,SAA+B;QAC1D,IAAI,CAAC,IAAI,CAAC,kBAAkB,IAAI,CAAC,IAAI,CAAC,eAAe,EAAE,CAAC;YACtD,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;QAChE,CAAC;QACD,MAAM,CAAC,GAAG,KAAK,IAAI,IAAI,CAAC,KAAK,IAAI,GAAiB,CAAC;QACnD,MAAM,EAAE,GAAG,SAAS,IAAI,IAAI,CAAC,SAAS,IAAI,GAA0B,CAAC;QAErE,OAAO,eAAe,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,WAAW,EAAE,CAAC,EAAE,IAAI,CAAC,kBAAkB,EAAE,IAAI,CAAC,eAAe,EAAE,EAAE,CAAC,CAAC;IACtH,CAAC;IAED;;;;;OAKG;IACH,KAAK,CAAC,KAAsB,EAAE,SAAgC,qBAAqB,CAAC,MAAM;QACxF,OAAO,eAAe,CAAC,IAAI,CAAC,cAAc,EAAE,IAAI,CAAC,WAAW,EAAE,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,WAAW,EAAE,MAAM,CAAC,CAAC;IACjH,CAAC;IAED;;;;;OAKG;IACH,YAAY,CAAC,KAAsB,EAAE,SAAgC,qBAAqB,CAAC,MAAM;QAC/F,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,OAAkB,CAAC;IACxD,CAAC;CACF"}