nrl-tracker 0.19.1__tar.gz → 1.6.0__tar.gz

{nrl_tracker-0.19.1/nrl_tracker.egg-info → nrl_tracker-1.6.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 0.19.1
+Version: 1.6.0
 Summary: Python port of the U.S. Naval Research Laboratory's Tracker Component Library for target tracking algorithms
 Author: Original: David F. Crouse, Naval Research Laboratory
 Maintainer: Python Port Contributors
@@ -45,6 +45,7 @@ Provides-Extra: dev
 Requires-Dist: pytest>=7.0.0; extra == "dev"
 Requires-Dist: pytest-cov>=4.0.0; extra == "dev"
 Requires-Dist: pytest-xdist>=3.0.0; extra == "dev"
+Requires-Dist: pytest-benchmark>=4.0.0; extra == "dev"
 Requires-Dist: hypothesis>=6.0.0; extra == "dev"
 Requires-Dist: black>=23.0.0; extra == "dev"
 Requires-Dist: isort>=5.12.0; extra == "dev"
@@ -55,27 +56,34 @@ Requires-Dist: sphinx>=6.0.0; extra == "dev"
 Requires-Dist: sphinx-rtd-theme>=1.2.0; extra == "dev"
 Requires-Dist: myst-parser>=1.0.0; extra == "dev"
 Requires-Dist: nbsphinx>=0.9.0; extra == "dev"
+Provides-Extra: benchmark
+Requires-Dist: pytest-benchmark>=4.0.0; extra == "benchmark"
 Provides-Extra: all
 Requires-Dist: nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualization]; extra == "all"
 
 # Tracker Component Library (Python)
 
+[![PyPI version](https://img.shields.io/badge/pypi-v1.6.0-blue.svg)](https://pypi.org/project/nrl-tracker/)
 [![Python 3.10+](https://img.shields.io/badge/python-3.10+-blue.svg)](https://www.python.org/downloads/)
 [![License: Public Domain](https://img.shields.io/badge/License-Public%20Domain-brightgreen.svg)](https://en.wikipedia.org/wiki/Public_domain)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)
+[![Tests](https://img.shields.io/badge/tests-1922%20passing-success.svg)](https://github.com/nedonatelli/TCL)
+[![MATLAB Parity](https://img.shields.io/badge/MATLAB%20Parity-99%25-brightgreen.svg)](docs/gap_analysis.rst)
 
 A Python port of the [U.S. Naval Research Laboratory's Tracker Component Library](https://github.com/USNavalResearchLaboratory/TrackerComponentLibrary), a comprehensive collection of algorithms for target tracking, estimation, coordinate systems, and related mathematical functions.
 
+**1,070+ functions** | **150+ modules** | **1,922 tests** | **99% MATLAB parity**
+
 ## Overview
 
 The Tracker Component Library provides building blocks for developing target tracking algorithms, including:
 
 - **Coordinate Systems**: Conversions between Cartesian, spherical, geodetic, and other coordinate systems
 - **Dynamic Models**: State transition matrices for constant velocity, coordinated turn, and other motion models
-- **Estimation Algorithms**: Kalman filters (EKF, UKF, etc.), particle filters, and batch estimation
+- **Estimation Algorithms**: Kalman filters (EKF, UKF, CKF, H-infinity), particle filters, and batch estimation
 - **Assignment Algorithms**: Hungarian algorithm, auction algorithms, and multi-dimensional assignment
 - **Mathematical Functions**: Special functions, statistics, numerical integration, and more
-- **Astronomical Code**: Ephemeris calculations, time systems, celestial mechanics
+- **Astronomical Code**: SGP4/SDP4 propagation, TLE parsing, reference frames (GCRF, ITRF, TEME, TOD, MOD), ephemerides
 - **Navigation**: Geodetic calculations, INS algorithms, GNSS utilities
 - **Geophysical Models**: Gravity, magnetism, atmosphere, and terrain models
 

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/README.md

@@ -1,21 +1,26 @@
 # Tracker Component Library (Python)
 
+[![PyPI version](https://img.shields.io/badge/pypi-v1.6.0-blue.svg)](https://pypi.org/project/nrl-tracker/)
 [![Python 3.10+](https://img.shields.io/badge/python-3.10+-blue.svg)](https://www.python.org/downloads/)
 [![License: Public Domain](https://img.shields.io/badge/License-Public%20Domain-brightgreen.svg)](https://en.wikipedia.org/wiki/Public_domain)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)
+[![Tests](https://img.shields.io/badge/tests-1922%20passing-success.svg)](https://github.com/nedonatelli/TCL)
+[![MATLAB Parity](https://img.shields.io/badge/MATLAB%20Parity-99%25-brightgreen.svg)](docs/gap_analysis.rst)
 
 A Python port of the [U.S. Naval Research Laboratory's Tracker Component Library](https://github.com/USNavalResearchLaboratory/TrackerComponentLibrary), a comprehensive collection of algorithms for target tracking, estimation, coordinate systems, and related mathematical functions.
 
+**1,070+ functions** | **150+ modules** | **1,922 tests** | **99% MATLAB parity**
+
 ## Overview
 
 The Tracker Component Library provides building blocks for developing target tracking algorithms, including:
 
 - **Coordinate Systems**: Conversions between Cartesian, spherical, geodetic, and other coordinate systems
 - **Dynamic Models**: State transition matrices for constant velocity, coordinated turn, and other motion models
-- **Estimation Algorithms**: Kalman filters (EKF, UKF, etc.), particle filters, and batch estimation
+- **Estimation Algorithms**: Kalman filters (EKF, UKF, CKF, H-infinity), particle filters, and batch estimation
 - **Assignment Algorithms**: Hungarian algorithm, auction algorithms, and multi-dimensional assignment
 - **Mathematical Functions**: Special functions, statistics, numerical integration, and more
-- **Astronomical Code**: Ephemeris calculations, time systems, celestial mechanics
+- **Astronomical Code**: SGP4/SDP4 propagation, TLE parsing, reference frames (GCRF, ITRF, TEME, TOD, MOD), ephemerides
 - **Navigation**: Geodetic calculations, INS algorithms, GNSS utilities
 - **Geophysical Models**: Gravity, magnetism, atmosphere, and terrain models
 

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0/nrl_tracker.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: nrl-tracker
-Version: 0.19.1
+Version: 1.6.0
 Summary: Python port of the U.S. Naval Research Laboratory's Tracker Component Library for target tracking algorithms
 Author: Original: David F. Crouse, Naval Research Laboratory
 Maintainer: Python Port Contributors
@@ -45,6 +45,7 @@ Provides-Extra: dev
 Requires-Dist: pytest>=7.0.0; extra == "dev"
 Requires-Dist: pytest-cov>=4.0.0; extra == "dev"
 Requires-Dist: pytest-xdist>=3.0.0; extra == "dev"
+Requires-Dist: pytest-benchmark>=4.0.0; extra == "dev"
 Requires-Dist: hypothesis>=6.0.0; extra == "dev"
 Requires-Dist: black>=23.0.0; extra == "dev"
 Requires-Dist: isort>=5.12.0; extra == "dev"
@@ -55,27 +56,34 @@ Requires-Dist: sphinx>=6.0.0; extra == "dev"
 Requires-Dist: sphinx-rtd-theme>=1.2.0; extra == "dev"
 Requires-Dist: myst-parser>=1.0.0; extra == "dev"
 Requires-Dist: nbsphinx>=0.9.0; extra == "dev"
+Provides-Extra: benchmark
+Requires-Dist: pytest-benchmark>=4.0.0; extra == "benchmark"
 Provides-Extra: all
 Requires-Dist: nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualization]; extra == "all"
 
 # Tracker Component Library (Python)
 
+[![PyPI version](https://img.shields.io/badge/pypi-v1.6.0-blue.svg)](https://pypi.org/project/nrl-tracker/)
 [![Python 3.10+](https://img.shields.io/badge/python-3.10+-blue.svg)](https://www.python.org/downloads/)
 [![License: Public Domain](https://img.shields.io/badge/License-Public%20Domain-brightgreen.svg)](https://en.wikipedia.org/wiki/Public_domain)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)
+[![Tests](https://img.shields.io/badge/tests-1922%20passing-success.svg)](https://github.com/nedonatelli/TCL)
+[![MATLAB Parity](https://img.shields.io/badge/MATLAB%20Parity-99%25-brightgreen.svg)](docs/gap_analysis.rst)
 
 A Python port of the [U.S. Naval Research Laboratory's Tracker Component Library](https://github.com/USNavalResearchLaboratory/TrackerComponentLibrary), a comprehensive collection of algorithms for target tracking, estimation, coordinate systems, and related mathematical functions.
 
+**1,070+ functions** | **150+ modules** | **1,922 tests** | **99% MATLAB parity**
+
 ## Overview
 
 The Tracker Component Library provides building blocks for developing target tracking algorithms, including:
 
 - **Coordinate Systems**: Conversions between Cartesian, spherical, geodetic, and other coordinate systems
 - **Dynamic Models**: State transition matrices for constant velocity, coordinated turn, and other motion models
-- **Estimation Algorithms**: Kalman filters (EKF, UKF, etc.), particle filters, and batch estimation
+- **Estimation Algorithms**: Kalman filters (EKF, UKF, CKF, H-infinity), particle filters, and batch estimation
 - **Assignment Algorithms**: Hungarian algorithm, auction algorithms, and multi-dimensional assignment
 - **Mathematical Functions**: Special functions, statistics, numerical integration, and more
-- **Astronomical Code**: Ephemeris calculations, time systems, celestial mechanics
+- **Astronomical Code**: SGP4/SDP4 propagation, TLE parsing, reference frames (GCRF, ITRF, TEME, TOD, MOD), ephemerides
 - **Navigation**: Geodetic calculations, INS algorithms, GNSS utilities
 - **Geophysical Models**: Gravity, magnetism, atmosphere, and terrain models
 

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/nrl_tracker.egg-info/SOURCES.txt

@@ -9,6 +9,7 @@ nrl_tracker.egg-info/dependency_links.txt
 nrl_tracker.egg-info/requires.txt
 nrl_tracker.egg-info/top_level.txt
 pytcl/__init__.py
+pytcl/logging_config.py
 pytcl/assignment_algorithms/__init__.py
 pytcl/assignment_algorithms/data_association.py
 pytcl/assignment_algorithms/gating.py
@@ -19,11 +20,16 @@ pytcl/assignment_algorithms/two_dimensional/__init__.py
 pytcl/assignment_algorithms/two_dimensional/assignment.py
 pytcl/assignment_algorithms/two_dimensional/kbest.py
 pytcl/astronomical/__init__.py
+pytcl/astronomical/ephemerides.py
 pytcl/astronomical/lambert.py
 pytcl/astronomical/orbital_mechanics.py
 pytcl/astronomical/reference_frames.py
+pytcl/astronomical/relativity.py
+pytcl/astronomical/sgp4.py
 pytcl/astronomical/time_systems.py
+pytcl/astronomical/tle.py
 pytcl/atmosphere/__init__.py
+pytcl/atmosphere/ionosphere.py
 pytcl/atmosphere/models.py
 pytcl/clustering/__init__.py
 pytcl/clustering/dbscan.py
@@ -31,6 +37,7 @@ pytcl/clustering/gaussian_mixture.py
 pytcl/clustering/hierarchical.py
 pytcl/clustering/kmeans.py
 pytcl/containers/__init__.py
+pytcl/containers/base.py
 pytcl/containers/cluster_set.py
 pytcl/containers/covertree.py
 pytcl/containers/kd_tree.py
@@ -59,8 +66,11 @@ pytcl/dynamic_estimation/smoothers.py
 pytcl/dynamic_estimation/batch_estimation/__init__.py
 pytcl/dynamic_estimation/kalman/__init__.py
 pytcl/dynamic_estimation/kalman/extended.py
+pytcl/dynamic_estimation/kalman/h_infinity.py
 pytcl/dynamic_estimation/kalman/linear.py
 pytcl/dynamic_estimation/kalman/square_root.py
+pytcl/dynamic_estimation/kalman/sr_ukf.py
+pytcl/dynamic_estimation/kalman/ud_filter.py
 pytcl/dynamic_estimation/kalman/unscented.py
 pytcl/dynamic_estimation/measurement_update/__init__.py
 pytcl/dynamic_estimation/particle_filters/__init__.py
@@ -106,9 +116,13 @@ pytcl/mathematical_functions/signal_processing/filters.py
 pytcl/mathematical_functions/signal_processing/matched_filter.py
 pytcl/mathematical_functions/special_functions/__init__.py
 pytcl/mathematical_functions/special_functions/bessel.py
+pytcl/mathematical_functions/special_functions/debye.py
 pytcl/mathematical_functions/special_functions/elliptic.py
 pytcl/mathematical_functions/special_functions/error_functions.py
 pytcl/mathematical_functions/special_functions/gamma_functions.py
+pytcl/mathematical_functions/special_functions/hypergeometric.py
+pytcl/mathematical_functions/special_functions/lambert_w.py
+pytcl/mathematical_functions/special_functions/marcum_q.py
 pytcl/mathematical_functions/statistics/__init__.py
 pytcl/mathematical_functions/statistics/distributions.py
 pytcl/mathematical_functions/statistics/estimators.py
@@ -119,8 +133,10 @@ pytcl/mathematical_functions/transforms/wavelets.py
 pytcl/misc/__init__.py
 pytcl/navigation/__init__.py
 pytcl/navigation/geodesy.py
+pytcl/navigation/great_circle.py
 pytcl/navigation/ins.py
 pytcl/navigation/ins_gnss.py
+pytcl/navigation/rhumb.py
 pytcl/performance_evaluation/__init__.py
 pytcl/performance_evaluation/estimation_metrics.py
 pytcl/performance_evaluation/track_metrics.py
@@ -153,11 +169,15 @@ tests/test_astronomical.py
 tests/test_clustering.py
 tests/test_coordinate_systems.py
 tests/test_coverage_boost.py
+tests/test_coverage_boost_2.py
 tests/test_dynamic_models.py
 tests/test_egm.py
 tests/test_emm.py
+tests/test_ephemerides.py
 tests/test_gaussian_mixtures.py
 tests/test_geophysical.py
+tests/test_great_circle.py
+tests/test_h_infinity.py
 tests/test_ins.py
 tests/test_ins_gnss.py
 tests/test_kalman_filters.py
@@ -168,16 +188,23 @@ tests/test_performance_evaluation.py
 tests/test_phase6_specialized.py
 tests/test_plotting.py
 tests/test_projections.py
+tests/test_relativity.py
+tests/test_rhumb.py
+tests/test_sgp4.py
 tests/test_signal_processing.py
 tests/test_smoothers.py
+tests/test_spatial_containers_parametrized.py
 tests/test_spatial_structures.py
+tests/test_special_functions_phase12.py
 tests/test_static_estimation.py
 tests/test_terrain.py
 tests/test_terrain_loaders.py
 tests/test_tides.py
+tests/test_tod_mod.py
 tests/test_trackers.py
 tests/test_tracking_containers.py
 tests/test_transforms.py
 tests/test_v030_comprehensive.py
 tests/test_v030_features.py
+tests/test_validation.py
 tests/unit/test_core.py

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/nrl_tracker.egg-info/requires.txt

@@ -9,10 +9,14 @@ nrl-tracker[astronomy,dev,geodesy,optimization,signal,visualization]
 astropy>=5.0
 jplephem>=2.18
 
+[benchmark]
+pytest-benchmark>=4.0.0
+
 [dev]
 pytest>=7.0.0
 pytest-cov>=4.0.0
 pytest-xdist>=3.0.0
+pytest-benchmark>=4.0.0
 hypothesis>=6.0.0
 black>=23.0.0
 isort>=5.12.0

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "nrl-tracker"
-version = "0.19.1"
+version = "1.6.0"
 description = "Python port of the U.S. Naval Research Laboratory's Tracker Component Library for target tracking algorithms"
 readme = "README.md"
 authors = [
@@ -67,6 +67,7 @@ dev = [
     "pytest>=7.0.0",
     "pytest-cov>=4.0.0",
     "pytest-xdist>=3.0.0",
+    "pytest-benchmark>=4.0.0",
     "hypothesis>=6.0.0",
     "black>=23.0.0",
     "isort>=5.12.0",
@@ -78,6 +79,9 @@ dev = [
     "myst-parser>=1.0.0",
     "nbsphinx>=0.9.0",
 ]
+benchmark = [
+    "pytest-benchmark>=4.0.0",
+]
 all = [
     "nrl-tracker[astronomy,geodesy,visualization,optimization,signal,dev]",
 ]
@@ -100,6 +104,9 @@ testpaths = ["tests"]
 markers = [
     "slow: marks tests as slow (deselect with '-m \"not slow\"')",
     "matlab_validated: tests validated against MATLAB reference output",
+    "benchmark: marks tests as benchmarks",
+    "light: marks benchmark as part of light (PR) suite",
+    "full: marks benchmark as part of full (main) suite",
 ]
 filterwarnings = [
     "ignore::DeprecationWarning",

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/pytcl/__init__.py

@@ -20,7 +20,7 @@ References
        no. 5, pp. 18-27, May 2017.
 """
 
-__version__ = "0.18.0"
+__version__ = "1.6.0"
 __author__ = "Python Port Contributors"
 __original_author__ = "David F. Crouse, Naval Research Laboratory"
 

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/pytcl/assignment_algorithms/data_association.py

@@ -10,13 +10,8 @@ from typing import List, NamedTuple, Optional
 import numpy as np
 from numpy.typing import ArrayLike, NDArray
 
-from pytcl.assignment_algorithms.gating import (
-    mahalanobis_batch,
-    mahalanobis_distance,
-)
-from pytcl.assignment_algorithms.two_dimensional import (
-    assign2d,
-)
+from pytcl.assignment_algorithms.gating import mahalanobis_batch, mahalanobis_distance
+from pytcl.assignment_algorithms.two_dimensional import assign2d
 
 
 class AssociationResult(NamedTuple):

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/pytcl/assignment_algorithms/jpda.py

@@ -12,6 +12,7 @@ as JPDA can handle measurement origin uncertainty in cluttered environments.
 from typing import List, NamedTuple, Optional, Tuple
 
 import numpy as np
+from numba import njit
 from numpy.typing import ArrayLike, NDArray
 from scipy.stats import chi2
 
@@ -314,45 +315,34 @@ def _jpda_exact(
     return beta
 
 
-def _jpda_approximate(
-    likelihood_matrix: NDArray,
-    gated: NDArray,
+@njit(cache=True)
+def _jpda_approximate_core(
+    likelihood_matrix: np.ndarray,
+    gated: np.ndarray,
     detection_prob: float,
     clutter_density: float,
-) -> NDArray:
-    """
-    Approximate JPDA using parametric approach.
-
-    Uses the approach from [1] which is O(n_tracks * n_meas^2).
-
-    References
-    ----------
-    .. [1] Fitzgerald, R.J., "Development of Practical PDA Logic for
-           Multitarget Tracking by Microprocessor", American Control
-           Conference, 1986.
-    """
-    n_tracks, n_meas = likelihood_matrix.shape
-    beta = np.zeros((n_tracks, n_meas + 1))
-
-    # For each track, compute association probabilities independently
-    # then apply correction for shared measurements
+) -> np.ndarray:
+    """JIT-compiled core of approximate JPDA computation."""
+    n_tracks = likelihood_matrix.shape[0]
+    n_meas = likelihood_matrix.shape[1]
+    beta = np.zeros((n_tracks, n_meas + 1), dtype=np.float64)
 
     # Likelihood ratio for each measurement given each track
-    # L[i,j] = likelihood(z_j | track i) / clutter_density
-    L = np.zeros((n_tracks, n_meas))
+    L = np.zeros((n_tracks, n_meas), dtype=np.float64)
     for i in range(n_tracks):
         for j in range(n_meas):
             if gated[i, j] and clutter_density > 0:
                 L[i, j] = likelihood_matrix[i, j] / clutter_density
             elif gated[i, j]:
-                L[i, j] = likelihood_matrix[i, j] * 1e10  # Large value
+                L[i, j] = likelihood_matrix[i, j] * 1e10
 
     # Compute delta factors (accounts for other tracks)
-    delta = np.ones((n_tracks, n_meas))
+    delta = np.ones((n_tracks, n_meas), dtype=np.float64)
 
     for j in range(n_meas):
-        # Sum of likelihood ratios for measurement j
-        sum_L = np.sum(L[:, j])
+        sum_L = 0.0
+        for i in range(n_tracks):
+            sum_L += L[i, j]
         for i in range(n_tracks):
             if sum_L > 0:
                 delta[i, j] = 1.0 / (1.0 + sum_L - L[i, j])
@@ -361,7 +351,6 @@ def _jpda_approximate(
 
     # Compute association probabilities
     for i in range(n_tracks):
-        # Denominator for normalization
         denom = 1.0 - detection_prob
 
         for j in range(n_meas):
@@ -369,9 +358,9 @@ def _jpda_approximate(
                 beta[i, j] = detection_prob * L[i, j] * delta[i, j]
                 denom += beta[i, j]
 
-        # Normalize
         if denom > 0:
-            beta[i, :n_meas] /= denom
+            for j in range(n_meas):
+                beta[i, j] /= denom
             beta[i, n_meas] = (1.0 - detection_prob) / denom
         else:
             beta[i, n_meas] = 1.0
@@ -379,6 +368,31 @@ def _jpda_approximate(
     return beta
 
 
+def _jpda_approximate(
+    likelihood_matrix: NDArray,
+    gated: NDArray,
+    detection_prob: float,
+    clutter_density: float,
+) -> NDArray:
+    """
+    Approximate JPDA using parametric approach.
+
+    Uses the approach from [1] which is O(n_tracks * n_meas^2).
+
+    References
+    ----------
+    .. [1] Fitzgerald, R.J., "Development of Practical PDA Logic for
+           Multitarget Tracking by Microprocessor", American Control
+           Conference, 1986.
+    """
+    return _jpda_approximate_core(
+        likelihood_matrix.astype(np.float64),
+        gated.astype(np.bool_),
+        detection_prob,
+        clutter_density,
+    )
+
+
 def jpda_update(
     track_states: List[ArrayLike],
     track_covariances: List[ArrayLike],

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/pytcl/assignment_algorithms/two_dimensional/assignment.py

@@ -211,14 +211,21 @@ def auction(
             # Compute values: value[j] = -cost[i,j] - prices[j]
             values = -cost[i, :] - prices
 
-            # Find best and second best
-            sorted_idx = np.argsort(values)[::-1]
-            best_j = sorted_idx[0]
-            best_value = values[best_j]
-
-            if len(sorted_idx) > 1:
-                second_value = values[sorted_idx[1]]
+            # Find best and second best using argpartition (O(n) vs O(n log n))
+            if len(values) >= 2:
+                # Get indices of top 2 values
+                top2_idx = np.argpartition(values, -2)[-2:]
+                # Determine which is best and second best
+                if values[top2_idx[0]] > values[top2_idx[1]]:
+                    best_j = top2_idx[0]
+                    second_value = values[top2_idx[1]]
+                else:
+                    best_j = top2_idx[1]
+                    second_value = values[top2_idx[0]]
+                best_value = values[best_j]
             else:
+                best_j = np.argmax(values)
+                best_value = values[best_j]
                 second_value = -np.inf
 
             # Bid increment

{nrl_tracker-0.19.1 → nrl_tracker-1.6.0}/pytcl/astronomical/__init__.py

@@ -2,7 +2,8 @@
 Astronomical calculations for target tracking.
 
 This module provides time system conversions, orbital mechanics,
-Lambert problem solvers, and reference frame transformations.
+Lambert problem solvers, reference frame transformations, and high-precision
+ephemerides for celestial bodies.
 
 Examples
 --------
@@ -18,8 +19,19 @@ Examples
 >>> r1 = np.array([5000, 10000, 2100])
 >>> r2 = np.array([-14600, 2500, 7000])
 >>> sol = lambert_universal(r1, r2, 3600)
+
+>>> # Query Sun position with high precision
+>>> from pytcl.astronomical import sun_position
+>>> r_sun, v_sun = sun_position(2451545.0)  # J2000.0
 """
 
+from pytcl.astronomical.ephemerides import (
+    DEEphemeris,
+    barycenter_position,
+    moon_position,
+    planet_position,
+    sun_position,
+)
 from pytcl.astronomical.lambert import (
     LambertSolution,
     bi_elliptic_transfer,
@@ -28,8 +40,10 @@ from pytcl.astronomical.lambert import (
     lambert_universal,
     minimum_energy_transfer,
 )
-from pytcl.astronomical.orbital_mechanics import (  # Constants; Types; Anomaly conversions; Element conversions; Propagation; Orbital quantities
-    GM_EARTH,
+from pytcl.astronomical.orbital_mechanics import (
+    GM_EARTH,  # Constants; Types; Anomaly conversions; Element conversions; Propagation; Orbital quantities
+)
+from pytcl.astronomical.orbital_mechanics import (
     GM_JUPITER,
     GM_MARS,
     GM_MOON,
@@ -62,8 +76,10 @@ from pytcl.astronomical.orbital_mechanics import (  # Constants; Types; Anomaly
     true_to_mean_anomaly,
     vis_viva,
 )
-from pytcl.astronomical.reference_frames import (  # Time utilities; Precession; Nutation; Earth rotation; Polar motion; Full transformations; Ecliptic/equatorial
-    earth_rotation_angle,
+from pytcl.astronomical.reference_frames import (
+    earth_rotation_angle,  # Time utilities; Precession; Nutation; Earth rotation; Polar motion; Full transformations; Ecliptic/equatorial
+)
+from pytcl.astronomical.reference_frames import (
     ecef_to_eci,
     eci_to_ecef,
     ecliptic_to_equatorial,
@@ -71,20 +87,59 @@ from pytcl.astronomical.reference_frames import (  # Time utilities; Precession;
     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.time_systems import (  # Julian dates; Time scales; Unix time; GPS week; Sidereal time; Leap seconds; Constants
-    JD_GPS_EPOCH,
+from pytcl.astronomical.relativity import (
+    C_LIGHT,  # Physical constants; Schwarzschild metric; Time dilation; Shapiro delay; Precession; PN effects; Range corrections
+)
+from pytcl.astronomical.relativity import (
+    G_GRAV,
+    geodetic_precession,
+    gravitational_time_dilation,
+    lense_thirring_precession,
+    post_newtonian_acceleration,
+    proper_time_rate,
+    relativistic_range_correction,
+    schwarzschild_precession_per_orbit,
+    schwarzschild_radius,
+    shapiro_delay,
+)
+from pytcl.astronomical.sgp4 import (
+    SGP4Satellite,
+    SGP4State,
+    sgp4_propagate,
+    sgp4_propagate_batch,
+)
+from pytcl.astronomical.time_systems import (
+    JD_GPS_EPOCH,  # Julian dates; Time scales; Unix time; GPS week; Sidereal time; Leap seconds; Constants
+)
+from pytcl.astronomical.time_systems import (
     JD_J2000,
     JD_UNIX_EPOCH,
     MJD_OFFSET,
@@ -112,6 +167,17 @@ from pytcl.astronomical.time_systems import (  # Julian dates; Time scales; Unix
     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
@@ -218,4 +284,57 @@ __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
+    "sun_position",
+    "moon_position",
+    "planet_position",
+    "barycenter_position",
+    # Relativity - Constants
+    "C_LIGHT",
+    "G_GRAV",
+    # Relativity - Functions
+    "schwarzschild_radius",
+    "gravitational_time_dilation",
+    "proper_time_rate",
+    "shapiro_delay",
+    "schwarzschild_precession_per_orbit",
+    "post_newtonian_acceleration",
+    "geodetic_precession",
+    "lense_thirring_precession",
+    "relativistic_range_correction",
 ]