nrl-tracker 0.22.5__py3-none-any.whl → 1.8.0__py3-none-any.whl

pytcl/astronomical/__init__.py

@@ -77,30 +77,46 @@ from pytcl.astronomical.orbital_mechanics import (
     vis_viva,
 )
 from pytcl.astronomical.reference_frames import (
-    earth_rotation_angle,  # Time utilities; Precession; Nutation; Earth rotation; Polar motion; Full transformations; Ecliptic/equatorial
+    earth_rotation_angle,  # Time utilities; Precession; Nutation
 )
+from pytcl.astronomical.reference_frames import ecef_to_eci  # Time utilities
 from pytcl.astronomical.reference_frames import (
-    ecef_to_eci,
     eci_to_ecef,
     ecliptic_to_equatorial,
     equation_of_equinoxes,
     equatorial_to_ecliptic,
     gast_iau82,
     gcrf_to_itrf,
+    gcrf_to_mod,
+    gcrf_to_teme,
+    gcrf_to_tod,
     gmst_iau82,
     itrf_to_gcrf,
+    itrf_to_teme,
+    itrf_to_teme_with_velocity,
+    itrf_to_tod,
     julian_centuries_j2000,
     mean_obliquity_iau80,
+    mod_to_gcrf,
+    mod_to_tod,
     nutation_angles_iau80,
     nutation_matrix,
+    pef_to_teme,
     polar_motion_matrix,
     precession_angles_iau76,
     precession_matrix_iau76,
     sidereal_rotation_matrix,
+    teme_to_gcrf,
+    teme_to_itrf,
+    teme_to_itrf_with_velocity,
+    teme_to_pef,
+    tod_to_gcrf,
+    tod_to_itrf,
+    tod_to_mod,
     true_obliquity,
 )
 from pytcl.astronomical.relativity import (
-    C_LIGHT,  # Physical constants; Schwarzschild metric; Time dilation; Shapiro delay; Precession; PN effects; Range corrections
+    C_LIGHT,  # Physical constants; Schwarzschild metric; Time dilation
 )
 from pytcl.astronomical.relativity import (
     G_GRAV,
@@ -114,6 +130,30 @@ from pytcl.astronomical.relativity import (
     schwarzschild_radius,
     shapiro_delay,
 )
+from pytcl.astronomical.sgp4 import (
+    SGP4Satellite,
+    SGP4State,
+    sgp4_propagate,
+    sgp4_propagate_batch,
+)
+from pytcl.astronomical.special_orbits import (
+    OrbitType,
+    classify_orbit,
+    eccentricity_vector,
+    escape_velocity_at_radius,
+    hyperbolic_anomaly_to_true_anomaly,
+    hyperbolic_asymptote_angle,
+    hyperbolic_deflection_angle,
+    hyperbolic_excess_velocity,
+    mean_to_parabolic_anomaly,
+    mean_to_true_anomaly_parabolic,
+    parabolic_anomaly_to_true_anomaly,
+    radius_parabolic,
+    semi_major_axis_from_energy,
+    true_anomaly_to_hyperbolic_anomaly,
+    true_anomaly_to_parabolic_anomaly,
+    velocity_parabolic,
+)
 from pytcl.astronomical.time_systems import (
     JD_GPS_EPOCH,  # Julian dates; Time scales; Unix time; GPS week; Sidereal time; Leap seconds; Constants
 )
@@ -145,6 +185,17 @@ from pytcl.astronomical.time_systems import (
     utc_to_tai,
     utc_to_tt,
 )
+from pytcl.astronomical.tle import (
+    TLE,
+    format_tle,
+    is_deep_space,
+    orbital_period_from_tle,
+    parse_tle,
+    parse_tle_3line,
+    semi_major_axis_from_mean_motion,
+    tle_epoch_to_datetime,
+    tle_epoch_to_jd,
+)
 
 __all__ = [
     # Time systems - Julian dates
@@ -251,6 +302,39 @@ __all__ = [
     # Reference frames - Ecliptic/equatorial
     "ecliptic_to_equatorial",
     "equatorial_to_ecliptic",
+    # Reference frames - TEME (for SGP4/SDP4)
+    "teme_to_pef",
+    "pef_to_teme",
+    "teme_to_itrf",
+    "itrf_to_teme",
+    "teme_to_gcrf",
+    "gcrf_to_teme",
+    "teme_to_itrf_with_velocity",
+    "itrf_to_teme_with_velocity",
+    # Reference frames - TOD/MOD (legacy conventions)
+    "gcrf_to_mod",
+    "mod_to_gcrf",
+    "gcrf_to_tod",
+    "tod_to_gcrf",
+    "mod_to_tod",
+    "tod_to_mod",
+    "tod_to_itrf",
+    "itrf_to_tod",
+    # TLE parsing
+    "TLE",
+    "parse_tle",
+    "parse_tle_3line",
+    "tle_epoch_to_jd",
+    "tle_epoch_to_datetime",
+    "format_tle",
+    "is_deep_space",
+    "semi_major_axis_from_mean_motion",
+    "orbital_period_from_tle",
+    # SGP4/SDP4 propagation
+    "SGP4State",
+    "SGP4Satellite",
+    "sgp4_propagate",
+    "sgp4_propagate_batch",
     # Ephemerides - Classes
     "DEEphemeris",
     # Ephemerides - Functions
@@ -271,4 +355,21 @@ __all__ = [
     "geodetic_precession",
     "lense_thirring_precession",
     "relativistic_range_correction",
+    # Special orbits - Parabolic and hyperbolic
+    "OrbitType",
+    "classify_orbit",
+    "mean_to_parabolic_anomaly",
+    "parabolic_anomaly_to_true_anomaly",
+    "true_anomaly_to_parabolic_anomaly",
+    "mean_to_true_anomaly_parabolic",
+    "radius_parabolic",
+    "velocity_parabolic",
+    "hyperbolic_anomaly_to_true_anomaly",
+    "true_anomaly_to_hyperbolic_anomaly",
+    "escape_velocity_at_radius",
+    "hyperbolic_excess_velocity",
+    "semi_major_axis_from_energy",
+    "hyperbolic_asymptote_angle",
+    "hyperbolic_deflection_angle",
+    "eccentricity_vector",
 ]

pytcl/astronomical/ephemerides.py

@@ -49,9 +49,10 @@ References
 
 """
 
-from typing import Literal, Optional, Tuple
+from typing import Any, Literal, Optional, Tuple
 
 import numpy as np
+from numpy.typing import NDArray
 
 # Constants for unit conversion
 AU_PER_KM = 1.0 / 149597870.7  # 1 AU in km
@@ -152,10 +153,10 @@ class DEEphemeris:
         self.version = version
         self._jplephem = jplephem
         self._kernel: Optional[object] = None
-        self._cache: dict = {}
+        self._cache: dict[str, Any] = {}
 
     @property
-    def kernel(self):
+    def kernel(self) -> Optional[object]:
         """Lazy-load ephemeris kernel on first access.
 
         Note: This requires jplephem to be installed and the kernel file
@@ -204,7 +205,7 @@ class DEEphemeris:
 
     def sun_position(
         self, jd: float, frame: Literal["icrf", "ecliptic"] = "icrf"
-    ) -> Tuple[np.ndarray, np.ndarray]:
+    ) -> Tuple[NDArray[np.floating], NDArray[np.floating]]:
         """Compute Sun position and velocity.
 
         Parameters
@@ -253,7 +254,7 @@ class DEEphemeris:
 
     def moon_position(
         self, jd: float, frame: Literal["icrf", "ecliptic", "earth_centered"] = "icrf"
-    ) -> Tuple[np.ndarray, np.ndarray]:
+    ) -> Tuple[NDArray[np.floating], NDArray[np.floating]]:
         """Compute Moon position and velocity.
 
         Parameters
@@ -323,7 +324,7 @@ class DEEphemeris:
         ],
         jd: float,
         frame: Literal["icrf", "ecliptic"] = "icrf",
-    ) -> Tuple[np.ndarray, np.ndarray]:
+    ) -> Tuple[np.ndarray[Any, Any], np.ndarray[Any, Any]]:
         """Compute planet position and velocity.
 
         Parameters
@@ -379,7 +380,7 @@ class DEEphemeris:
 
     def barycenter_position(
         self, body: str, jd: float
-    ) -> Tuple[np.ndarray, np.ndarray]:
+    ) -> Tuple[NDArray[np.floating], NDArray[np.floating]]:
         """Compute position of any body relative to Solar System Barycenter.
 
         Parameters
@@ -424,7 +425,7 @@ def _get_default_ephemeris() -> DEEphemeris:
 
 def sun_position(
     jd: float, frame: Literal["icrf", "ecliptic"] = "icrf"
-) -> Tuple[np.ndarray, np.ndarray]:
+) -> Tuple[NDArray[np.floating], NDArray[np.floating]]:
     """Convenience function: Compute Sun position and velocity.
 
     Parameters
@@ -451,7 +452,7 @@ def sun_position(
 
 def moon_position(
     jd: float, frame: Literal["icrf", "ecliptic", "earth_centered"] = "icrf"
-) -> Tuple[np.ndarray, np.ndarray]:
+) -> Tuple[NDArray[np.floating], NDArray[np.floating]]:
     """Convenience function: Compute Moon position and velocity.
 
     Parameters
@@ -482,7 +483,7 @@ def planet_position(
     ],
     jd: float,
     frame: Literal["icrf", "ecliptic"] = "icrf",
-) -> Tuple[np.ndarray, np.ndarray]:
+) -> Tuple[np.ndarray[Any, Any], np.ndarray[Any, Any]]:
     """Convenience function: Compute planet position and velocity.
 
     Parameters
@@ -509,7 +510,9 @@ def planet_position(
     return _get_default_ephemeris().planet_position(planet, jd, frame=frame)
 
 
-def barycenter_position(body: str, jd: float) -> Tuple[np.ndarray, np.ndarray]:
+def barycenter_position(
+    body: str, jd: float
+) -> Tuple[NDArray[np.floating], NDArray[np.floating]]:
     """Convenience function: Position relative to Solar System Barycenter.
 
     Parameters