siderust 0.6.0

High-precision astronomy and satellite mechanics in Rust.
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77

// SPDX-License-Identifier: AGPL-3.0-or-later
// Copyright (C) 2026 Vallés Puig, Ramon

//! Star observability in altitude + azimuth ranges.
//!
//! Run with: `cargo run --example 09_star_observability`

use siderust::bodies::catalog::SIRIUS;
use siderust::calculus::altitude::{AltitudePeriodsProvider, AltitudeQuery};
use siderust::calculus::azimuth::{AzimuthProvider, AzimuthQuery};
use siderust::observatories::ROQUE_DE_LOS_MUCHACHOS;
use siderust::time::intersect_periods;
use siderust::time::{ModifiedJulianDate, Period};

use qtty::*;

fn main() {
    println!("Star observability: altitude + azimuth constraints\n");

    let observer = ROQUE_DE_LOS_MUCHACHOS;
    let target = SIRIUS;

    // One-night search window (MJD TT).
    let t_0 = ModifiedJulianDate::new(60000.0);
    let window = Period::new(t_0, t_0 + Days::new(1.0));

    // Constraint 1: altitude between 25° and 65°.
    let altitude_query = AltitudeQuery {
        observer,
        window,
        min_altitude: Degrees::new(25.0),
        max_altitude: Degrees::new(65.0),
    };
    let altitude_periods = target.altitude_periods(&altitude_query);

    // Constraint 2: azimuth between 110° and 220° (ESE -> SW sector).
    let azimuth_query = AzimuthQuery {
        observer,
        window,
        min_azimuth: Degrees::new(110.0),
        max_azimuth: Degrees::new(220.0),
    };
    let azimuth_periods = target.azimuth_periods(&azimuth_query);

    // Final observability: periods satisfying both constraints simultaneously.
    let observable_periods = intersect_periods(&altitude_periods, &azimuth_periods);

    println!("Observer: {}", observer);
    println!("Target: Sirius");
    println!("Window: {}\n", window);

    println!(
        "Altitude range: {}..{}",
        altitude_query.min_altitude, altitude_query.max_altitude
    );
    println!(
        "Azimuth range:  {}..{}\n",
        azimuth_query.min_azimuth, azimuth_query.max_azimuth
    );

    println!("Matched periods: {}", observable_periods.len());
    for (idx, period) in observable_periods.iter().enumerate() {
        let hours = period.duration_days().to::<Hour>();
        println!(
            "  {}. {} -> {}  ({})",
            idx + 1,
            period.start,
            period.end,
            hours
        );
    }

    let total_hours = observable_periods.iter().fold(Hours::new(0.0), |acc, p| {
        acc + p.duration_days().to::<Hour>()
    });
    println!("\nTotal observable time in both ranges: {}", total_hours);
}