keplemon 1.3.1__cp310-cp310-macosx_11_0_arm64.whl → 1.5.0__cp310-cp310-macosx_11_0_arm64.whl

keplemon/__init__.py

@@ -4,6 +4,7 @@ from keplemon._keplemon import (  # type: ignore
     set_thread_count,
     set_license_file_path,
     get_license_file_path,
+    set_jpl_ephemeris_file_path,
 )
 from pathlib import Path
 
@@ -18,10 +19,15 @@ set_license_file_path(ASSETS_DIRECTORY.as_posix())
 TIME_CONSTANTS_PATH = ASSETS_DIRECTORY / "time_constants.dat"
 load_time_constants(TIME_CONSTANTS_PATH.as_posix())
 
+# Load the JPL path
+JPL_EPHEMERIS_PATH = ASSETS_DIRECTORY / "JPLcon_1950_2050.405"
+set_jpl_ephemeris_file_path(JPL_EPHEMERIS_PATH.as_posix())
+
 __all__ = [
     "get_thread_count",
     "set_thread_count",
     "TIME_CONSTANTS_PATH",
+    "JPL_EPHEMERIS_PATH",
     "set_license_file_path",
     "PACKAGE_DIRECTORY",
     "ASSETS_DIRECTORY",

keplemon/bodies.pyi

@@ -7,6 +7,8 @@ from keplemon.elements import (
     GeodeticPosition,
     OrbitPlotData,
     TopocentricElements,
+    RelativeState,
+    BoreToBodyAngles,
 )
 from keplemon.catalogs import TLECatalog
 from keplemon.time import Epoch, TimeSpan
@@ -80,6 +82,26 @@ class Satellite:
         """
         ...
 
+    def get_relative_state_at_epoch(self, other: Satellite, epoch: Epoch) -> RelativeState | None:
+        """
+        Calculate the relative state between this satellite and another satellite at a given epoch.
+
+        Args:
+            other: Secondary satellite to calculate the relative state against
+            epoch: UTC epoch at which the relative state will be calculated
+        """
+        ...
+
+    def get_body_angles_at_epoch(self, other: Satellite, epoch: Epoch) -> BoreToBodyAngles | None:
+        """
+        Calculate the bore-to-body angles between this satellite and another satellite at a given epoch.
+
+        Args:
+            other: Secondary satellite to calculate the bore-to-body angles against
+            epoch: UTC epoch at which the bore-to-body angles will be calculated
+        """
+        ...
+
     def get_plot_data(self, start: Epoch, end: Epoch, step: TimeSpan) -> OrbitPlotData | None: ...
     def get_observatory_access_report(
         self,

keplemon/elements.py

@@ -12,6 +12,8 @@ from keplemon._keplemon.elements import (  # type: ignore
     GeodeticPosition,
     OrbitPlotData,
     OrbitPlotState,
+    RelativeState,
+    BoreToBodyAngles,
 )
 
 __all__ = [
@@ -28,4 +30,6 @@ __all__ = [
     "GeodeticPosition",
     "OrbitPlotData",
     "OrbitPlotState",
+    "RelativeState",
+    "BoreToBodyAngles",
 ]

keplemon/elements.pyi

@@ -4,6 +4,18 @@ from keplemon.time import Epoch
 from keplemon.enums import Classification, KeplerianType, ReferenceFrame
 from keplemon.events import CloseApproach
 
+class RelativeState:
+    epoch: Epoch
+    position: CartesianVector
+    velocity: CartesianVector
+    origin_satellite_id: str
+    secondary_satellite_id: str
+
+class BoreToBodyAngles:
+    earth_angle: float
+    sun_angle: float
+    moon_angle: float
+
 class OrbitPlotData:
     satellite_id: str
     epochs: list[str]

keplemon/saal/astro_func_interface.py

@@ -3,6 +3,9 @@ from keplemon._keplemon.saal.astro_func_interface import (  # type: ignore
     ra_dec_to_az_el_time,
     ra_dec_to_az_el,
     mean_motion_to_sma,
+    sma_to_mean_motion,
+    kozai_to_brouwer,
+    brouwer_to_kozai,
     topo_date_to_equinox,
     topo_equinox_to_date,
     theta_teme_to_lla,
@@ -20,6 +23,7 @@ from keplemon._keplemon.saal.astro_func_interface import (  # type: ignore
     XA_TOPO_SIZE,
     YROFEQNX_2000,
     YROFEQNX_CURR,
+    get_jpl_sun_and_moon_position,
 )
 
 __all__ = [
@@ -44,4 +48,8 @@ __all__ = [
     "YROFEQNX_2000",
     "YROFEQNX_CURR",
     "topo_equinox_to_date",
+    "get_jpl_sun_and_moon_position",
+    "kozai_to_brouwer",
+    "brouwer_to_kozai",
+    "sma_to_mean_motion",
 ]

keplemon/saal/astro_func_interface.pyi

@@ -2,13 +2,41 @@
 
 def mean_motion_to_sma(mean_motion: float) -> float:
     """
-    Convert mean motion to semi-major axis.
+    Convert mean motion to semi-major axis in kilometers.
 
     Args:
         mean_motion: Mean motion in revolutions/day.
+    """
+    ...
 
-    Returns:
-        Semi-major axis in kilometers.
+def sma_to_mean_motion(sma: float) -> float:
+    """
+    Convert semi-major axis to mean motion in revolutions/day.
+
+    Args:
+        sma: Semi-major axis in kilometers.
+    """
+    ...
+
+def kozai_to_brouwer(e_kozai: float, i_kozai: float, n_kozai: float) -> float:
+    """
+    Convert Kozai orbital elements to Brouwer orbital elements.
+
+    Args:
+        e_kozai: Eccentricity (unitless).
+        i_kozai: Inclination in degrees.
+        n_kozai: Mean motion in revolutions/day.
+    """
+    ...
+
+def brouwer_to_kozai(e_brouwer: float, i_brouwer: float, n_brouwer: float) -> float:
+    """
+    Convert Brouwer orbital elements to Kozai orbital elements.
+
+    Args:
+        e_brouwer: Eccentricity (unitless).
+        i_brouwer: Inclination in degrees.
+        n_brouwer: Mean motion in revolutions/day.
     """
     ...
 
@@ -121,6 +149,20 @@ def theta_teme_to_lla(theta: float, lat: float, long: float, ra: float, dec: flo
     """
     ...
 
+def get_jpl_sun_and_moon_position(ds50_utc: float) -> tuple[list[float], list[float]]:
+    """
+    Get the JPL ephemeris positions of the Sun and Moon in TEME coordinates.
+
+    Args:
+        ds50_utc: Epoch in DS50 UTC format.
+
+    Returns:
+        A tuple containing two lists:
+            - Sun position in TEME coordinates [x, y, z] in kilometers.
+            - Moon position in TEME coordinates [x, y, z] in kilometers.
+    """
+    ...
+
 def time_teme_to_lla(ds50_utc: float, lat: float, long: float, ra: float, dec: float) -> tuple[float, float]:
     """
     Convert TEME coordinates to latitude, longitude, and altitude.

keplemon/saal/sat_state_interface.py

@@ -0,0 +1,5 @@
+from keplemon._keplemon.saal.sat_state_interface import (  # type: ignore
+    get_relative_state,
+)
+
+__all__ = ["get_relative_state"]

keplemon/saal/sat_state_interface.pyi

@@ -0,0 +1,9 @@
+def get_relative_state(state_1: list[float], state_2: list[float], utc_ds50: float) -> list[float]:
+    """Calculate the relative state between two satellites
+
+    Args:
+        state_1: primary satellite TEME state vector [x, y, z, vx, vy, vz] in km and km/s
+        state_2: secondary satellite TEME state vector [x, y, z, vx, vy, vz] in km and km/s
+        utc_ds50: UTC time in days since 1950
+    """
+    ...

keplemon/time.py

@@ -27,8 +27,14 @@ def _now() -> Epoch:
     return _from_datetime(datetime.now(timezone.utc))
 
 
+def _to_datetime(epoch: Epoch) -> datetime:
+    iso_str = epoch.to_iso()
+    return datetime.fromisoformat(iso_str.replace("Z", "+00:00"))
+
+
 Epoch.now = staticmethod(_now)
 Epoch.from_datetime = staticmethod(_from_datetime)
+Epoch.to_datetime = _to_datetime
 
 
 def request_time_constants_update(output_path: str) -> None:

keplemon/time.pyi

@@ -146,6 +146,13 @@ class Epoch:
     day_of_year: float
     """Decimal day of the year (1-365.999...)"""
 
+    def to_datetime(self) -> datetime:
+        """
+        Returns:
+            Aware datetime object in UTC time system
+        """
+        ...
+
     @classmethod
     def from_datetime(cls, dt: datetime) -> Epoch:
         """

{keplemon-1.3.1.dist-info → keplemon-1.5.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: keplemon
-Version: 1.3.1
+Version: 1.5.0
 Requires-Dist: requests
 Requires-Dist: click
 Requires-Dist: maturin>=1.0,<2.0 ; extra == 'dev'

{keplemon-1.3.1.dist-info → keplemon-1.5.0.dist-info}/RECORD

@@ -1,10 +1,10 @@
-keplemon-1.3.1.dist-info/METADATA,sha256=fl7n1Iefcs9qhxB00Y617U8cbU94YAwmuswu_482oJc,863
-keplemon-1.3.1.dist-info/WHEEL,sha256=PmVieto1wuHPE0V9Yj-HDpDcXvcgT7RQ6xfDnzOpcS8,104
-keplemon-1.3.1.dist-info/entry_points.txt,sha256=eYbCkvQvWfRDQ0LzaCELov1xeLAxQEHlfdgNq-LXyb0,49
-keplemon/__init__.py,sha256=M9q5lNYh_BE6l4xCGJ5IH5PQH9aNm4q_r67ljsNkKvM,832
+keplemon-1.5.0.dist-info/METADATA,sha256=5zcpgXY79lEuvKCf89MaAfoXWmf0xeEsqGEu-zLVHzo,863
+keplemon-1.5.0.dist-info/WHEEL,sha256=PmVieto1wuHPE0V9Yj-HDpDcXvcgT7RQ6xfDnzOpcS8,104
+keplemon-1.5.0.dist-info/entry_points.txt,sha256=eYbCkvQvWfRDQ0LzaCELov1xeLAxQEHlfdgNq-LXyb0,49
+keplemon/__init__.py,sha256=G4RZhxNf7pu0I_2qASrVAe4aZ6ezY9YvbmnwZP7KJx8,1034
 keplemon/__init__.pyi,sha256=uE60ln_KJgcfvKburVmbcKT0h_wLPgjBWuyNLgI8ETI,1295
 keplemon/__main__.py,sha256=-3GVkDOA0lV0MIqU9gPb4zbVimg2lA8HMkvdPDw1O28,669
-keplemon/_keplemon.cpython-310-darwin.so,sha256=YRRttNqWFmMhZsARLbwbdQdOFhgAPLxwJ1uYQxXbZ7c,2076352
+keplemon/_keplemon.cpython-310-darwin.so,sha256=DIoXR8PmmsJj3pmK4207aONEeJAGykBpaVASZEyna3I,2103024
 keplemon/assets/EGM-2008.GEO,sha256=K2nG8HGLATIHZYMfw3GSClYOTCuZ7rq4RdCeUNgCw5A,148770
 keplemon/assets/EGM-96.GEO,sha256=VBkILuvEMwAPuWmUHy2PeyEfULOwJ4PEJLNf5hr84mU,148770
 keplemon/assets/GEM_5-22.GEO,sha256=stemYLn1ChXa-VdLGHYfa15AXZa_xxGZQ65p4c3gffI,6852
@@ -17,11 +17,11 @@ keplemon/assets/SGP4_Open_License.txt,sha256=0WofOXQb5YJqnYhXWXnBdCajiTJQAT60UAk
 keplemon/assets/WGS84-70.GEO,sha256=ARjEC_5s2SVd0Kh9udbTy1ztBwTeuBYPOhUVJgIqit8,148510
 keplemon/assets/time_constants.dat,sha256=qDpJ2UrQvIDfxsBc4P2AdLS-b2lyR7RCzjqmeG4Ypl8,1226736
 keplemon/bodies.py,sha256=XnaY6XTuj8CHM3XOwOSY3E8nSo0RWwCcAY0FGxAVWa8,208
-keplemon/bodies.pyi,sha256=xpciSiscdGgXJ_RTLCd08U6vJIdNX0dYAnlBzPW8240,10091
+keplemon/bodies.pyi,sha256=zLfQQPBP-6ormpm6dAObZaTgnUPJnoBW4wChgoJVm-w,10959
 keplemon/catalogs.py,sha256=lw71NiXlVtb-z3pQR03afxtkLca4HJcnpZ6kDCcR-Lk,102
 keplemon/catalogs.pyi,sha256=NU_6Mc_JY8xTYAKXOqubtPnt91YA9ZHd1hMqMFf6onc,589
-keplemon/elements.py,sha256=P4tlpQpCUG5jkgAT_yIeL2V14AhXtekMYxRBOsPlOgc,609
-keplemon/elements.pyi,sha256=L5s_2AJz6HYYIncKmpkbYcejMVht09WWVrXRqo-ZV-Y,10240
+keplemon/elements.py,sha256=aVWooAU5RB2RTPirjOrhIyaEohfTc1hLa0Fgc83fYy8,695
+keplemon/elements.pyi,sha256=pgFwj8AzKSWr0QI9WQLWjIAbTv5QtGA9tJ9EQCR-k5o,10491
 keplemon/enums.py,sha256=Jh0tFHg_rZXnOyLUXmHjSm3MSZFbvQKTBcP0BqHXeMY,308
 keplemon/enums.pyi,sha256=s9uqkm3Zrx6HLV4dQBScRiUy0CT4QoQwFjaCLOEMW1c,2330
 keplemon/estimation.py,sha256=2K87pFhOfaFohbAtBXv2SA6m_oIN-56toJOftu350fY,228
@@ -51,14 +51,16 @@ keplemon/propagation.pyi,sha256=qJOCbtgYma4MOg3-aw9_r1IvBDnknwsJlMwlgaJiEZk,688
 keplemon/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 keplemon/saal/__init__.py,sha256=aoTB13q7IvXuQ80jgGCnwXA29yjpOqLiTdrtyjyRqLE,123
 keplemon/saal/__init__.pyi,sha256=GPB5f0dcK7srvSDq2i5wvHMyi-OYZakMvlrstKdDwkk,143
-keplemon/saal/astro_func_interface.py,sha256=Z3UEbviqkGCU4OcpvgShvMoT5_WeenYZkec6A5hyDgs,984
-keplemon/saal/astro_func_interface.pyi,sha256=it0olsVygZtIY4ekh_Hbu1_q_o87IvE7Ezslzr0_hZo,4782
+keplemon/saal/astro_func_interface.py,sha256=zf56Sgw8WnRZFt-jQtiYUSxbzgPZlGPW6gA496bRSeM,1198
+keplemon/saal/astro_func_interface.pyi,sha256=6nr0xlb5r3PL4GFzRHH_-lFEMHvJQzh5x5WY3R3-DdM,6006
 keplemon/saal/obs_interface.py,sha256=EtsaPKixjMWfCUpMd79SvhCs3f9sRcBaal6-ickj3Vs,248
 keplemon/saal/obs_interface.pyi,sha256=eIXbFnZSF3cX3MyXaCYLUAp0wUAbJiQ4rosqpXdf2I0,228
+keplemon/saal/sat_state_interface.py,sha256=XSV6HEoVGytePa71MdlTbXuAUCgDFLDR04r92yO67So,134
+keplemon/saal/sat_state_interface.pyi,sha256=UbZvnEt3sEc5Fydwm9QKvIKXnMzOZL6pQVYQEv2L5UY,414
 keplemon/saal/sgp4_prop_interface.py,sha256=BBKlIvXdUypcZspC8GJ2rPqfFoDKy0RM7w5qsqvobrE,127
 keplemon/saal/sgp4_prop_interface.pyi,sha256=Fv41UrqSPuITtYy_CcbyMxhysqtN3oeY38Ov4t81hvw,235
 keplemon/saal/time_func_interface.py,sha256=cshqJ15p_gcenMdmVuXTIoLeij1gsgVi0tujRQ4O6kA,421
 keplemon/saal/time_func_interface.pyi,sha256=GCj_EOmOceJorYQLGQQj1fE2cHxPvNrYml1DLvsaMy4,1508
-keplemon/time.py,sha256=OuJYhv9WFi0SF3P_6YKLSoLWDnbhQ-W8CDlNUvQJsxQ,3236
-keplemon/time.pyi,sha256=tFbv5qdVMGexvSHr_a6RxkXe0aDKcNWOxW3GmKf3eRc,6628
-keplemon-1.3.1.dist-info/RECORD,,
+keplemon/time.py,sha256=1VTb9RbMvgEkJkRFeuPrFXnvO_ZA9HuSGbuYzh76d_4,3410
+keplemon/time.pyi,sha256=r37R56aJCPQJONTYn1stT6_s_XBI_BAPjSb17l_Gox8,6774
+keplemon-1.5.0.dist-info/RECORD,,