rust-spice 0.4.7

WOW! The complete NASA/NAIF Spice toolkit is actually usable on 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

use itertools::multizip;
use serde_json::json;
use spice;
use tool::Vectors;

fn main() -> Result<(), spice::KernelError> {
    let mut kernel = spice::Kernel::new("rsc/krn/hera_study_PO_EMA_2024.tm")?;
    let frame = "J2000";
    let target = "DIDYMOS_BARYCENTER";
    let observer = "SUN";
    let date = "2027-MAR-23 16:00:00";
    let abcorr = "NONE";
    let time_step = 1.0 * tool::DAY;
    let duration = 192.0 * tool::DAY;

    let et_start = spice::str2et(date);
    let ephemeris_times = tool::linspace(et_start, et_start + duration, time_step);
    let size = ephemeris_times.len();

    println!("{}", spice::timout(et_start + duration, spice::TIME_FORMAT));

    let mut positions = Vectors::<f64>::zeros(size);

    for (et, mut position) in
        multizip((ephemeris_times.iter().cloned(), positions.column_iter_mut()))
    {
        let (pos, _) = spice::spkpos(target, et, frame, abcorr, observer);
        position.copy_from(&pos);
    }

    // Export
    let json = json!({
        "context": {
            "frame": frame,
            "observer": observer,
            "target": target,
            "abcorr": abcorr,
            "date": date,
            "duration": duration,
            "time_step": time_step,
        },
        "positions": positions,
    });
    tool::writejs!("rsc/data/tmp/position.json", json);

    kernel.unload()?;
    Ok(())
}