brahe 1.4.0

Brahe is a modern satellite dynamics library for research and engineering applications designed to be easy-to-learn, high-performance, and quick-to-deploy. The north-star of the development is enabling users to solve meaningful problems and answer questions quickly, easily, and correctly.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63

# /// script
# dependencies = ["brahe", "plotly", "kaleido"]
# ///
from pathlib import Path
import brahe as bh


bh.initialize_eop()

# Download GP data for all GPS satellites from CelesTrak
client = bh.celestrak.CelestrakClient()
gp_records = client.get_gp(group="gps-ops")
propagators = [rec.to_sgp_propagator(step_size=60.0) for rec in gp_records]

# Propagate each satellite one orbit
for prop in propagators:
    orbital_period = bh.orbital_period(prop.semi_major_axis)
    prop.propagate_to(prop.epoch + orbital_period)

# Create interactive 3D plot with Earth texture
fig = bh.plot_trajectory_3d(
    [
        {
            "trajectory": prop.trajectory,
            "mode": "markers",
            "size": 2,
            "label": prop.get_name(),
        }
        for prop in propagators
    ],
    units="km",
    show_earth=True,
    earth_texture="natural_earth_50m",
    backend="plotly",
    view_azimuth=45.0,
    view_elevation=30.0,
    view_distance=2.4,
)

# Enable grid and axis lines for better visualization
fig.update_layout(
    scene=dict(
        xaxis=dict(
            showbackground=False, gridcolor="lightgray", zerolinecolor="lightgray"
        ),
        yaxis=dict(
            showbackground=False, gridcolor="lightgray", zerolinecolor="lightgray"
        ),
        zaxis=dict(
            showbackground=False, gridcolor="lightgray", zerolinecolor="lightgray"
        ),
    )
)

# Disable the legend for clarity
fig.update_layout(showlegend=False)

# Remove title
fig.update_layout(title_text="")

# Save off plotly figure to svg and pdf files
fig.write_image(Path(__file__).parent / "gps_satellites.svg")
fig.write_image(Path(__file__).parent / "gps_satellites.pdf")