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
49
50
51
52
53
54
55
56
57
58

use itertools::multizip;
use nalgebra::Rotation3;
use serde_json::json;
use tool::{Vector, Vectors};

fn main() -> Result<(), spice::KernelError> {
    let mut kernel = spice::Kernel::new("rsc/krn/hera_study_PO_EA_2026.tm")?;
    let frame_inertial = "J2000";
    let primary = "DIDYMOS";
    let secondary = "DIMORPHOS";
    let primary_frame_name = "DIDYMOS_FIXED";
    let secondary_frame_name = "DIMORPHOS_FIXED";
    let date = "2027-MAR-23 16:00:00";
    let abcorr = "NONE";
    let time_step = 1.0 * tool::MINUTE;
    let duration = 2.5 * tool::HOUR;

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

    let mut secondary_positions = Vectors::<f64>::zeros(size);
    let mut primary_frames = vec![Rotation3::new(Vector::zeros()); size];
    let mut secondary_frames = vec![Rotation3::new(Vector::zeros()); size];

    for (et, mut secondary_pos, primary_frame, secondary_frame) in multizip((
        ephemeris_times.iter().cloned(),
        secondary_positions.column_iter_mut(),
        primary_frames.iter_mut(),
        secondary_frames.iter_mut(),
    )) {
        secondary_pos.copy_from(&spice::spkpos(secondary, et, frame_inertial, abcorr, primary).0);
        *primary_frame = spice::pxform(frame_inertial, primary_frame_name, et);
        *secondary_frame = spice::pxform(frame_inertial, secondary_frame_name, et);
    }

    // Export.
    let json = json!({
        "context": {
            "frame_inertial": frame_inertial,
            "primary": primary,
            "secondary": secondary,
            "primary_frame_name": primary_frame_name,
            "secondary_frame_name": secondary_frame_name,
            "abcorr": abcorr,
            "date": date,
            "duration": duration,
            "time_step": time_step,
        },
        "secondary_positions": secondary_positions,
        "primary_frames": primary_frames,
        "secondary_frames": secondary_frames,
    });
    tool::writejs!("rsc/data/tmp/binary_inertial.json", json);

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