open-space-toolkit-astrodynamics 2.4.1__py311-none-any.whl → 5.0.0__py311-none-any.whl

ostk/astrodynamics/test/trajectory/orbit/messages/__init__.py

@@ -1 +0,0 @@
-# Apache License 2.0

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/__init__.py

@@ -1 +0,0 @@
-# Apache License 2.0

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/conftest.py

@@ -1,18 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-import pathlib
-
-from ostk.core.filesystem import Path
-from ostk.core.filesystem import File
-
-
-@pytest.fixture
-def data_directory_path() -> str:
-    return f"{pathlib.Path(__file__).parent.absolute()}/data"
-
-
-@pytest.fixture
-def opm_file(data_directory_path: str) -> File:
-    return File.path(Path.parse(f"{data_directory_path}/opm_1.yaml"))

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/data/opm_1.yaml

@@ -1,44 +0,0 @@
-# Dummy SpaceX OPM output
-
-# Notes:
-# - ECEF velocity is Earth relative
-# - Apogee/Perigee altitude assumes a spherical Earth, 6378.137 km radius
-# - Orbital elements are computed in an inertial frame realized by inertially
-#   freezing the WGS84 ECEF frame at time of current state
-# - State is post-deployment, so includes separation delta-velocity
-
-
-header:
-  generation_date: 2020-01-01T12:34:56.789Z
-  launch_date: 2020-01-02T12:34:56.789Z
-
-
-deployments:
-
-- name: satellite_a
-  sequence_number: 1
-  mission_time_s: 3600.0
-  date: 2020-01-02T13:34:56.789Z
-  r_ecef_m: [693289.644, 6876578.628, -133035.288]
-  v_ecef_m_per_s: [1305.783, 39.783, 7525.920]
-  mean_perigee_altitude_km: 526.768
-  mean_apogee_altitude_km: 568.430
-  mean_inclination_deg: 97.123
-  mean_argument_of_perigee_deg: -179.513
-  mean_longitude_ascending_node_deg: 85.057
-  mean_mean_anomaly_deg: 179.263
-  ballistic_coef_kg_per_m2: 47.55
-
-- name: satellite_b
-  sequence_number: 2
-  mission_time_s: 7200.0
-  date: 2020-01-02T14:34:56.789Z
-  r_ecef_m: [699863.059, 6875647.517, -123777.595]
-  v_ecef_m_per_s: [1504.658, 6.705, 7538.669]
-  mean_perigee_altitude_km: 536.779
-  mean_apogee_altitude_km: 529.851
-  mean_inclination_deg: 97.124
-  mean_argument_of_perigee_deg: 136.875
-  mean_longitude_ascending_node_deg: 85.032
-  mean_mean_anomaly_deg: -127.164
-  ballistic_coef_kg_per_m2: 44.26

ostk/astrodynamics/test/trajectory/orbit/messages/spacex/test_opm.py

@@ -1,108 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-from datetime import datetime
-
-from ostk.core.containers import Dictionary
-from ostk.core.filesystem import Path
-from ostk.core.filesystem import File
-
-from ostk.physics.units import Length
-from ostk.physics.units import Angle
-from ostk.physics.time import Instant
-from ostk.physics.time import Scale
-from ostk.physics.time import Duration
-from ostk.physics.coordinate import Position
-from ostk.physics.coordinate import Velocity
-from ostk.physics.coordinate import Frame
-
-from ostk.astrodynamics.trajectory.orbit.messages.spacex import OPM
-
-
-@pytest.fixture
-def opm() -> OPM:
-    return OPM(
-        header=OPM.Header(
-            generation_date=Instant.date_time(datetime(2020, 1, 2, 3, 4, 5), Scale.UTC),
-            launch_date=Instant.date_time(datetime(2020, 1, 3, 3, 4, 5), Scale.UTC),
-        ),
-        deployments=[
-            OPM.Deployment(
-                name="A",
-                sequence_number=1,
-                mission_time=Duration.hours(3600.0),
-                date=Instant.date_time(datetime(2020, 1, 4, 3, 4, 5), Scale.UTC),
-                position=Position.meters((1.0, 2.0, 3.0), Frame.ITRF()),
-                velocity=Velocity.meters_per_second((4.0, 5.0, 6.0), Frame.ITRF()),
-                mean_perigee_altitude=Length.kilometers(500.0),
-                mean_apogee_altitude=Length.kilometers(500.0),
-                mean_inclination=Angle.degrees(1.0),
-                mean_argument_of_perigee=Angle.degrees(2.0),
-                mean_longitude_ascending_node=Angle.degrees(3.0),
-                mean_mean_anomaly=Angle.degrees(4.0),
-                ballistic_coefficient=123.456,
-            ),
-        ],
-    )
-
-
-class TestOPM:
-    def test_constructor(self, opm: OPM):
-        assert opm is not None
-
-    def test_get_header(self, opm: OPM):
-        assert opm.get_header() is not None
-
-    def test_get_deployments(self, opm: OPM):
-        assert len(opm.get_deployments()) == 1
-
-    def test_get_deployment_at(self, opm: OPM):
-        assert opm.get_deployment_at(index=0).name == "A"
-
-    def test_get_deployment_with_name(self, opm: OPM):
-        assert opm.get_deployment_with_name(name="A").name == "A"
-
-    def test_undefined(self):
-        assert OPM.undefined().is_defined() is False
-
-    def test_dictionary(self):
-        dictionary = Dictionary(
-            {
-                "header": {
-                    "generation_date": "2020-01-01T12:34:56.789Z",
-                    "launch_date": "2020-01-02T12:34:56.789Z",
-                },
-                "deployments": [
-                    {
-                        "name": "satellite_a",
-                        "sequence_number": 1,
-                        "mission_time_s": 3600.0,
-                        "date": "2020-01-02T13:34:56.789Z",
-                        "r_ecef_m": [693289.644, 6876578.628, -133035.288],
-                        "v_ecef_m_per_s": [1305.783, 39.783, 7525.920],
-                        "mean_perigee_altitude_km": 526.768,
-                        "mean_apogee_altitude_km": 568.430,
-                        "mean_inclination_deg": 97.123,
-                        "mean_argument_of_perigee_deg": -179.513,
-                        "mean_longitude_ascending_node_deg": 85.057,
-                        "mean_mean_anomaly_deg": 179.263,
-                        "ballistic_coef_kg_per_m2": 47.55,
-                    }
-                ],
-            }
-        )
-
-        assert OPM.dictionary(dictionary) is not None
-
-    def test_parse(self, opm_file: File):
-        with open(str(opm_file.get_path().to_string()), "r") as stream:
-            assert OPM.parse(string=stream.read()) is not None
-
-    def test_load(self, opm_file: File):
-        assert OPM.load(file=opm_file) is not None
-
-
-class TestOPMDeployment:
-    def test_to_state(self, opm: OPM):
-        assert opm.get_deployment_with_name(name="A").to_state() is not None

ostk/astrodynamics/test/trajectory/orbit/models/__init__.py

@@ -1 +0,0 @@
-# Apache License 2.0

ostk/astrodynamics/test/trajectory/orbit/models/kepler/__init__.py

@@ -1 +0,0 @@
-# Apache License 2.0

ostk/astrodynamics/test/trajectory/orbit/models/kepler/test_coe.py

@@ -1,112 +0,0 @@
-# Apache License 2.0
-
-from ostk.physics.units import Length
-from ostk.physics.units import Angle
-from ostk.physics.environment.gravitational import Earth
-
-from ostk.astrodynamics.trajectory.orbit.models.kepler import COE
-
-
-def construct_coe():
-    a = Length.kilometers(7000.0)
-    e = 0.1
-    i = Angle.degrees(35.0)
-    raan = Angle.degrees(40.0)
-    aop = Angle.degrees(50.0)
-    nu = Angle.degrees(60.0)
-
-    coe = COE(a, e, i, raan, aop, nu)
-
-    return coe
-
-
-def test_trajectory_orbit_models_kepler_coe_constructors():
-    a = Length.kilometers(7000.0)
-    e = 0.1
-    i = Angle.degrees(35.0)
-    raan = Angle.degrees(40.0)
-    aop = Angle.degrees(50.0)
-    nu = Angle.degrees(60.0)
-
-    coe: COE = COE(a, e, i, raan, aop, nu)
-
-    assert coe is not None
-    assert isinstance(coe, COE)
-    assert coe.is_defined()
-
-    coe: COE = COE.undefined()
-
-    assert coe is not None
-    assert isinstance(coe, COE)
-    assert coe.is_defined() is False
-
-
-def test_trajectory_orbit_models_kepler_coe_comparators():
-    coe: COE = construct_coe()
-
-    assert coe == coe
-    assert (coe != coe) is False
-
-
-def test_trajectory_orbit_models_kepler_coe_getters():
-    coe: COE = construct_coe()
-
-    # get_semi_major_axis()
-
-    assert coe.get_semi_major_axis() is not None
-
-    # get_eccentricity()
-
-    assert coe.get_eccentricity() is not None
-
-    # get_inclination()
-
-    assert coe.get_inclination() is not None
-
-    # get_raan()
-
-    assert coe.get_raan() is not None
-
-    # get_aop()
-
-    assert coe.get_aop() is not None
-
-    # get_true_anomaly()
-
-    assert coe.get_true_anomaly() is not None
-
-    # get_mean_anomaly()
-
-    assert coe.get_mean_anomaly() is not None
-
-    # get_eccentric_anomaly()
-
-    assert coe.get_eccentric_anomaly() is not None
-
-
-def test_trajectory_orbit_models_kepler_coe_get_mean_motion():
-    coe: COE = construct_coe()
-
-    assert coe.get_mean_motion(Earth.EGM2008.gravitational_parameter) is not None
-
-
-def test_trajectory_orbit_models_kepler_coe_get_orbital_period():
-    coe: COE = construct_coe()
-
-    assert coe.get_orbital_period(Earth.EGM2008.gravitational_parameter) is not None
-
-
-def test_trajectory_orbit_models_kepler_coe_static_methods():
-    coe: COE = construct_coe()
-
-    assert COE.eccentric_anomaly_from_true_anomaly(Angle.degrees(0.0), 0.0) is not None
-    assert COE.true_anomaly_from_eccentric_anomaly(Angle.degrees(0.0), 0.0) is not None
-    assert COE.mean_anomaly_from_eccentric_anomaly(Angle.degrees(0.0), 0.0) is not None
-    assert (
-        COE.eccentric_anomaly_from_mean_anomaly(Angle.degrees(0.0), 0.0, 0.0) is not None
-    )
-
-
-def test_string_from_element():
-    element_str = COE.string_from_element(COE.Element.SemiMajorAxis)
-    assert element_str == "SemiMajorAxis"

ostk/astrodynamics/test/trajectory/orbit/models/sgp4/__init__.py

@@ -1 +0,0 @@
-# Apache License 2.0

ostk/astrodynamics/test/trajectory/orbit/models/sgp4/test_tle.py

@@ -1,331 +0,0 @@
-# Apache License 2.0
-
-import pytest
-
-import tempfile
-
-from ostk.core.filesystem import Path
-from ostk.core.filesystem import File
-
-from ostk.physics.units import Time
-from ostk.physics.units import Angle
-from ostk.physics.units import Derived
-from ostk.physics.time import Instant
-from ostk.physics.time import Scale
-from ostk.physics.time import DateTime
-
-from ostk.astrodynamics.trajectory.orbit.models.sgp4 import TLE
-
-
-@pytest.fixture
-def tle() -> TLE:
-    return TLE(
-        satellite_name="Satellite",
-        first_line="1 25544U 98067A   18231.17878740  .00000187  00000-0  10196-4 0  9994",
-        second_line="2 25544  51.6447  64.7824 0005971  73.1467  36.4366 15.53848234128316",
-    )
-
-
-class TestTLE:
-    def test_constructor(self):
-        tle = TLE(
-            first_line="1 25544U 98067A   18231.17878740  .00000187  00000-0  10196-4 0  9994",
-            second_line="2 25544  51.6447  64.7824 0005971  73.1467  36.4366 15.53848234128316",
-        )
-
-        assert isinstance(tle, TLE)
-
-    def test_constructor_with_satellite_name(self):
-        tle = TLE(
-            satellite_name="Satellite",
-            first_line="1 25544U 98067A   18231.17878740  .00000187  00000-0  10196-4 0  9994",
-            second_line="2 25544  51.6447  64.7824 0005971  73.1467  36.4366 15.53848234128316",
-        )
-
-        assert isinstance(tle, TLE)
-
-    def test_comparators(self, tle: TLE):
-        assert (tle == tle) is True
-        assert (tle != tle) is False
-
-    def test_is_defined(self, tle: TLE):
-        assert tle.is_defined()
-
-    def test_get_satellite_name(self, tle: TLE):
-        assert tle.get_satellite_name() == "Satellite"
-
-    def test_get_first_line(self, tle: TLE):
-        assert (
-            tle.get_first_line()
-            == "1 25544U 98067A   18231.17878740  .00000187  00000-0  10196-4 0  9994"
-        )
-
-    def test_get_second_line(self, tle: TLE):
-        assert (
-            tle.get_second_line()
-            == "2 25544  51.6447  64.7824 0005971  73.1467  36.4366 15.53848234128316"
-        )
-
-    def test_get_satellite_number(self, tle: TLE):
-        assert tle.get_satellite_number() == 25544
-
-    def test_get_classification(self, tle: TLE):
-        assert tle.get_classification() == "U"
-
-    def test_get_international_designator(self, tle: TLE):
-        assert tle.get_international_designator() == "98067A"
-
-    def test_get_epoch(self, tle: TLE):
-        assert tle.get_epoch() == Instant.date_time(
-            DateTime(2018, 8, 19, 4, 17, 27, 231, 360, 0), Scale.UTC
-        )
-
-    def test_get_mean_motion_first_time_derivative_divided_by_two(self, tle: TLE):
-        assert (
-            tle.get_mean_motion_first_time_derivative_divided_by_two()
-            == 1.8700000000000001e-06
-        )
-
-    def test_get_mean_motion_second_time_derivative_divided_by_six(self, tle: TLE):
-        assert tle.get_mean_motion_second_time_derivative_divided_by_six() == 0.0
-
-    def test_get_b_star_drag_term(self, tle: TLE):
-        assert tle.get_b_star_drag_term() == 1.0196e-05
-
-    def test_get_ephemeris_type(self, tle: TLE):
-        assert tle.get_ephemeris_type() == 0
-
-    def test_get_element_set_number(self, tle: TLE):
-        assert tle.get_element_set_number() == 999
-
-    def test_get_first_line_checksum(self, tle: TLE):
-        assert tle.get_first_line_checksum() == 4
-
-    def test_get_inclination(self, tle: TLE):
-        assert tle.get_inclination().in_degrees() == 51.6447
-
-    def test_get_raan(self, tle: TLE):
-        assert tle.get_raan().in_degrees() == 64.782399999999996
-
-    def test_get_eccentricity(self, tle: TLE):
-        assert tle.get_eccentricity() == 0.00059710000000000004
-
-    def test_get_aop(self, tle: TLE):
-        assert tle.get_aop().in_degrees() == 73.146699999999996
-
-    def test_get_mean_anomaly(self, tle: TLE):
-        assert tle.get_mean_anomaly().in_degrees() == 36.436599999999999
-
-    def test_get_mean_motion(self, tle: TLE):
-        assert (
-            tle.get_mean_motion().in_unit(
-                Derived.Unit.angular_velocity(Angle.Unit.Revolution, Time.Unit.Day)
-            )
-            == 15.53848234
-        )
-
-    def test_get_revolution_number_at_epoch(self, tle: TLE):
-        assert tle.get_revolution_number_at_epoch() == 12831
-
-    def test_get_second_line_checksum(self, tle: TLE):
-        assert tle.get_second_line_checksum() == 6
-
-    def test_set_satellite_number(self, tle: TLE):
-        tle.set_satellite_number(99959)
-
-        assert tle.get_satellite_number() == 99959
-
-        tle.set_satellite_number(25544)
-
-        assert tle.get_satellite_number() == 25544
-
-    def test_set_epoch(self, tle: TLE):
-        tle.set_epoch(
-            Instant.date_time(DateTime(2019, 9, 20, 5, 18, 28, 232, 361, 0), Scale.UTC)
-        )
-
-        assert tle.get_epoch() == Instant.date_time(
-            DateTime(2019, 9, 20, 5, 18, 28, 231, 776, 0), Scale.UTC
-        )
-
-        tle.set_epoch(
-            Instant.date_time(DateTime(2018, 8, 19, 4, 17, 27, 231, 360, 0), Scale.UTC)
-        )
-
-        assert tle.get_epoch() == Instant.date_time(
-            DateTime(2018, 8, 19, 4, 17, 27, 231, 360, 0), Scale.UTC
-        )
-
-    def test_set_revolution_number_at_epoch(self, tle: TLE):
-        tle.set_revolution_number_at_epoch(2345)
-
-        assert tle.get_revolution_number_at_epoch() == 2345
-        assert (
-            tle.get_first_line()
-            == "1 25544U 98067A   18231.17878740  .00000187  00000-0  10196-4 0  9994"
-        )
-        assert (
-            tle.get_second_line()
-            == "2 25544  51.6447  64.7824 0005971  73.1467  36.4366 15.53848234 23455"
-        )
-
-        tle.set_revolution_number_at_epoch(6)
-
-        assert tle.get_revolution_number_at_epoch() == 6
-        assert (
-            tle.get_first_line()
-            == "1 25544U 98067A   18231.17878740  .00000187  00000-0  10196-4 0  9994"
-        )
-        assert (
-            tle.get_second_line()
-            == "2 25544  51.6447  64.7824 0005971  73.1467  36.4366 15.53848234    67"
-        )
-
-        other_tle = TLE(
-            "1 99993U 21990ZZZ 21182.62513889  .00000763  00000-0  42347-4 0 99995",
-            "2 99993 097.5132 311.4037 0016005 231.4378 006.3908 15.13696975009992",
-        )
-
-        assert other_tle.get_revolution_number_at_epoch() == 999
-
-        other_tle.set_revolution_number_at_epoch(999)
-
-        assert other_tle.get_revolution_number_at_epoch() == 999
-        assert (
-            other_tle.get_first_line()
-            == "1 99993U 21990ZZZ 21182.62513889  .00000763  00000-0  42347-4 0 99995"
-        )
-        assert (
-            other_tle.get_second_line()
-            == "2 99993 097.5132 311.4037 0016005 231.4378 006.3908 15.13696975  9992"
-        )
-
-        other_tle.set_revolution_number_at_epoch(9909)
-
-        assert other_tle.get_revolution_number_at_epoch() == 9909
-        assert (
-            other_tle.get_first_line()
-            == "1 99993U 21990ZZZ 21182.62513889  .00000763  00000-0  42347-4 0 99995"
-        )
-        assert (
-            other_tle.get_second_line()
-            == "2 99993 097.5132 311.4037 0016005 231.4378 006.3908 15.13696975 99092"
-        )
-
-        other_tle.set_revolution_number_at_epoch(99)
-
-        assert other_tle.get_revolution_number_at_epoch() == 99
-        assert (
-            other_tle.get_first_line()
-            == "1 99993U 21990ZZZ 21182.62513889  .00000763  00000-0  42347-4 0 99995"
-        )
-        assert (
-            other_tle.get_second_line()
-            == "2 99993 097.5132 311.4037 0016005 231.4378 006.3908 15.13696975   993"
-        )
-
-    def test_can_parse(self, tle: TLE):
-        assert TLE.can_parse(tle.get_first_line(), tle.get_second_line())
-
-    def test_parse(self, tle: TLE):
-        tle_lines: str = "\n".join(
-            [
-                str(tle.get_satellite_name()),
-                str(tle.get_first_line()),
-                str(tle.get_second_line()),
-            ]
-        )
-
-        assert TLE.parse(tle_lines) == tle
-
-    def test_load(self, tle: TLE):
-        tle_lines: str = "\n".join(
-            [
-                str(tle.get_satellite_name()),
-                str(tle.get_first_line()),
-                str(tle.get_second_line()),
-            ]
-        )
-
-        tle_file = tempfile.NamedTemporaryFile(suffix=".tle", delete=False)
-        tle_file.write(tle_lines.encode())
-        tle_file.close()
-
-        try:
-            assert TLE.load(File.path(Path.parse(tle_file.name))) == tle
-
-        finally:
-            File.path(Path.parse(tle_file.name)).remove()
-
-    def test_construct_with_satellite_name(self, tle: TLE):
-        constructed_tle = TLE.construct(
-            satellite_name="Satellite",
-            satellite_number=25544,
-            classification="U",
-            international_designator="98067A",
-            epoch=Instant.date_time(
-                DateTime(2018, 8, 19, 4, 17, 27, 231, 360, 0), Scale.UTC
-            ),
-            mean_motion_first_time_derivative_divided_by_two=1.8700000000000001e-06,
-            mean_motion_second_time_derivative_divided_by_six=0.0,
-            b_star_drag_term=1.0196e-05,
-            ephemeris_type=0,
-            element_set_number=999,
-            inclination=Angle.degrees(51.6447),
-            raan=Angle.degrees(64.782399999999996),
-            eccentricity=0.00059710000000000004,
-            aop=Angle.degrees(73.146699999999996),
-            mean_anomaly=Angle.degrees(36.436599999999999),
-            mean_motion=Derived(
-                15.53848234,
-                Derived.Unit.angular_velocity(Angle.Unit.Revolution, Time.Unit.Day),
-            ),
-            revolution_number_at_epoch=12831,
-        )
-
-        assert constructed_tle == tle
-
-    def test_construct_without_satellite_name(self, tle: TLE):
-        constructed_tle = TLE.construct(
-            satellite_number=25544,
-            classification="U",
-            international_designator="98067A",
-            epoch=Instant.date_time(
-                DateTime(2018, 8, 19, 4, 17, 27, 231, 360, 0), Scale.UTC
-            ),
-            mean_motion_first_time_derivative_divided_by_two=1.8700000000000001e-06,
-            mean_motion_second_time_derivative_divided_by_six=0.0,
-            b_star_drag_term=1.0196e-05,
-            ephemeris_type=0,
-            element_set_number=999,
-            inclination=Angle.degrees(51.6447),
-            raan=Angle.degrees(64.782399999999996),
-            eccentricity=0.00059710000000000004,
-            aop=Angle.degrees(73.146699999999996),
-            mean_anomaly=Angle.degrees(36.436599999999999),
-            mean_motion=Derived(
-                15.53848234,
-                Derived.Unit.angular_velocity(Angle.Unit.Revolution, Time.Unit.Day),
-            ),
-            revolution_number_at_epoch=12831,
-        )
-
-        assert constructed_tle.get_first_line() == tle.get_first_line()
-        assert constructed_tle.get_second_line() == tle.get_second_line()
-
-    def test_generate_checksum(self, tle: TLE):
-        assert (
-            TLE.generate_checksum(tle.get_first_line()) == tle.get_first_line_checksum()
-        )
-        assert (
-            TLE.generate_checksum(tle.get_second_line()) == tle.get_second_line_checksum()
-        )
-
-        assert (
-            TLE.generate_checksum(f"{str(tle.get_first_line())[:-1]}0")
-            == tle.get_first_line_checksum()
-        )
-        assert (
-            TLE.generate_checksum(f"{str(tle.get_second_line())[:-1]}9")
-            == tle.get_second_line_checksum()
-        )